Showing 1–50 of 57 results for author: Saitoh, T R

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

Chemo-dynamical Evolution of Simulated Satellites for a Milky Way-like Galaxy

Authors: Yutaka Hirai, Evan N. Kirby, Masashi Chiba, Kohei Hayashi, Borja Anguiano, Takayuki R. Saitoh, Miho N. Ishigaki, Timothy C. Beers

Abstract: The chemical abundances of Milky Way's satellites reflect their star formation histories (SFHs), yet, due to the difficulty of determining the ages of old stars, the SFHs of most satellites are poorly measured. Ongoing and upcoming surveys will obtain around ten times more medium-resolution spectra for stars in satellites than are currently available. To correctly extract SFHs from large samples o… ▽ More The chemical abundances of Milky Way's satellites reflect their star formation histories (SFHs), yet, due to the difficulty of determining the ages of old stars, the SFHs of most satellites are poorly measured. Ongoing and upcoming surveys will obtain around ten times more medium-resolution spectra for stars in satellites than are currently available. To correctly extract SFHs from large samples of chemical abundances, the relationship between chemical abundances and SFHs needs to be clarified. Here, we perform a high-resolution cosmological zoom-in simulation of a Milky Way-like galaxy with detailed models of star formation, supernova feedback, and metal diffusion. We quantify SFHs, metallicity distribution functions, and the $α$-element (Mg, Ca, and Si) abundances in satellites of the host galaxy. We find that star formation in most simulated satellites is quenched before infalling to their host. Star formation episodes in simulated satellites are separated by a few hundred Myr owing to supernova feedback; each star formation event produces groups of stars with similar [$α$/Fe] and [Fe/H]. We then perform a mock observation of the upcoming Subaru Prime Focus Spectrograph (PFS) observations. We find that Subaru PFS will be able to detect distinct groups of stars in [$α$/Fe] vs. [Fe/H] space, produced by episodic star formation. This result means that episodic SFHs can be estimated from the chemical abundances of $\gtrsim$ 1,000 stars determined with medium-resolution spectroscopy. △ Less

Submitted 22 May, 2024; v1 submitted 8 May, 2024; originally announced May 2024.

Comments: 16 pages, 9 figures, accepted for publication in The Astrophysical Journal

Impulsive Gas Fueling to Galactic Center in a Barred Galaxy Due to Falls of Gas Clouds

Authors: Hidenori Matsui, Toshiyasu Masakawa, Asao Habe, Takayuki R. Saitoh

Abstract: We have studied the evolution of the central hundred pc region of barred galaxies by performing numerical simulations realizing multi-phase nature of gas. Our simulations have shown that a stellar bar produces an oval gas ring namely the $x$-2 ring within $1~{\rm kpc}$ as the bar grows. The ring is self-gravitationally unstable enough to trigger formations of gas clouds. Although the gas clouds in… ▽ More We have studied the evolution of the central hundred pc region of barred galaxies by performing numerical simulations realizing multi-phase nature of gas. Our simulations have shown that a stellar bar produces an oval gas ring namely the $x$-2 ring within $1~{\rm kpc}$ as the bar grows. The ring is self-gravitationally unstable enough to trigger formations of gas clouds. Although the gas clouds initially rotate in the $x$-2 ring, cloud-cloud collisions and/or energy injections into the gas ring by Type II supernovae deviate some of the clouds from the ring orbit. After the deviation, the deviated clouds repeat collisions with the other clouds, which rotate in the $x$-2 ring, during several rotations. These processes effectively reduce the angular momentum of the deviated gas cloud. As a result, the gas cloud finally falls into the galactic center, and episodic gas supply to the galactic center takes place. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 10 pages, 9 figures, accepted for publication in PASJ

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

Exploring the Sun's birth radius and the distribution of planet building blocks in the Milky Way galaxy: A multi-zone Galactic chemical evolution approach

Authors: Junichi Baba, Takayuki R. Saitoh, Takuji Tsujimoto

Abstract: We explore the influence of the Milky Way galaxy's chemical evolution on the formation, structure, and habitability of the Solar system. Using a multi-zone Galactic Chemical Evolution (GCE) model, we successfully reproduce key observational constraints, including the age-metallicity ([Fe/H]) relation, metallicity distribution functions, abundance gradients, and [X/Fe] ratio trends for critical ele… ▽ More We explore the influence of the Milky Way galaxy's chemical evolution on the formation, structure, and habitability of the Solar system. Using a multi-zone Galactic Chemical Evolution (GCE) model, we successfully reproduce key observational constraints, including the age-metallicity ([Fe/H]) relation, metallicity distribution functions, abundance gradients, and [X/Fe] ratio trends for critical elements involved in planetary mineralogy, including C, O, Mg, and Si. Our GCE model suggests that the Sun formed in the inner Galactic disc, $R_{\rm birth,\odot}\approx 5$ kpc. We also combined a stoichiometric model with the GCE model to examine the temporal evolution and spatial distribution of planet building blocks (PBBs) within the Milky Way galaxy, revealing trends in the condensed mass fraction ($f_{\rm cond}$), iron-to-silicon mass fraction ($f_{\rm iron}$), and water mass fraction ($f_{\rm water}$) over time and towards the inner Galactic disc regions. Specifically, our model predicts a higher $f_{\rm cond}$ in the protoplanetary disc within the inner regions of the Milky Way galaxy, as well as an increased $f_{\rm iron}$ and a decreased $f_{\rm water}$ in the inner regions. Based on these findings, we discuss the potential impact of the Sun's birth location on the overall structure and habitability of the Solar System. △ Less

Submitted 16 October, 2023; originally announced October 2023.

Comments: 15 pages, 10 figures. Accepted for publication in MNRAS

$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

SIRIUS Project. V. Formation of off-center ionized bubbles associated with Orion Nebula Cluster

Authors: Michiko S. Fujii, Kohei Hattori, Long Wang, Yutaka Hirai, Jun Kumamoto, Yoshito Shimajiri, Takayuki R. Saitoh

Abstract: Massive stars born in star clusters terminate star cluster formation by ionizing the surrounding gas. This process is considered to be prevalent in young star clusters containing massive stars. The Orion Nebula is an excellent example associated with a forming star cluster including several massive stars (the Orion Nebula Cluster; ONC) and a 2-pc size H{\sc ii} region (ionized bubble) opening towa… ▽ More Massive stars born in star clusters terminate star cluster formation by ionizing the surrounding gas. This process is considered to be prevalent in young star clusters containing massive stars. The Orion Nebula is an excellent example associated with a forming star cluster including several massive stars (the Orion Nebula Cluster; ONC) and a 2-pc size H{\sc ii} region (ionized bubble) opening toward the observer; however, the other side is still covered with dense molecular gas. Recent astrometric data acquired by the Gaia satellite revealed the stellar kinematics in this region. By comparing this data with star cluster formation simulation results, we demonstrate that massive stars born in the ONC center were ejected via three-body encounters. Further, orbit analysis indicates that $θ^2$ Ori A, the second massive star in this region, was ejected from the ONC center toward the observer and is now returning to the cluster center. Such ejected massive stars can form a hole in the dense molecular cloud and contribute to the formation of the 2-pc bubble. Our results demonstrate that the dynamics of massive stars are essential for the formation of star clusters and H{\sc ii} regions that are not always centered by massive stars. △ Less

Submitted 9 June, 2022; originally announced June 2022.

Comments: 13 pages, 12 figures, accepted for MNRAS

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

Testing the effect of resolution on gravitational fragmentation with Lagrangian hydrodynamic schemes

Authors: Yasuyoshi Yamamoto, Takashi Okamoto, Takayuki R. Saitoh

Abstract: To study the resolution required for simulating gravitational fragmentation with newly developed Lagrangian hydrodynamic schemes, Meshless Finite Volume method (MFV) and Meshless Finite Mass method (MFM), we have performed a number of simulations of the Jeans test and compared the results with both the expected analytic solution and results from the more standard Lagrangian approach: Smoothed Part… ▽ More To study the resolution required for simulating gravitational fragmentation with newly developed Lagrangian hydrodynamic schemes, Meshless Finite Volume method (MFV) and Meshless Finite Mass method (MFM), we have performed a number of simulations of the Jeans test and compared the results with both the expected analytic solution and results from the more standard Lagrangian approach: Smoothed Particle Hydrodynamics (SPH). We find that the different schemes converge to the analytic solution when the diameter of a fluid element is smaller than a quarter of the Jeans wavelength, $λ_\mathrm{J}$. Among the three schemes, SPH/MFV shows the fastest/slowest convergence to the analytic solution. Unlike the well-known behaviour of Eulerian schemes, none of the Lagrangian schemes investigated displays artificial fragmentation when the perturbation wavelength, $λ$, is shorter than $λ_\mathrm{J}$, even at low numerical resolution. For larger wavelengths ($λ> λ_\mathrm{J}$) the growth of the perturbation is delayed when it is not well resolved. Furthermore, with poor resolution, the fragmentation seen with the MFV scheme proceeds very differently compared to the converged solution. All these results suggest that, when unresolved, the ratio of the magnitude of hydrodynamic force to that of self-gravity at the sub-resolution scale is the largest/smallest in MFV/SPH, the reasons for which we discussed in detail. These tests are repeated to investigate the effect of kernels of higher-order than the fiducial cubic spline. Our results indicate that the standard deviation of the kernel is a more appropriate definition of the size of a fluid element than its compact support radius. △ Less

Submitted 20 April, 2021; v1 submitted 18 April, 2021; originally announced April 2021.

Comments: 11 pages, 11 figures, accepted for publication in MNRAS. The typos in the authors' names have been corrected

First Star Survivors as Metal-Rich Halo Stars that Experienced Supernova Explosions in Binary Systems

Authors: Takuma Suda, Takayuki R. Saitoh, Yuki Moritani, Tadafumi Matsuno, Toshikazu Shigeyama

Abstract: The search for the first stars formed from metal-free gas in the universe is one of the key issues in astronomy because it relates to many fields, such as the formation of stars and galaxies, the evolution of the universe, and the origin of elements. It is not still clear if metal-free first stars can be found in the present universe. These first stars are thought to exist among extremely metal-po… ▽ More The search for the first stars formed from metal-free gas in the universe is one of the key issues in astronomy because it relates to many fields, such as the formation of stars and galaxies, the evolution of the universe, and the origin of elements. It is not still clear if metal-free first stars can be found in the present universe. These first stars are thought to exist among extremely metal-poor stars in the halo of our Galaxy. Here we propose a new scenario for the formation of low-mass first stars that have survived until today and observational counterparts in our Galaxy. The first stars in binary systems, consisting of massive- and low-mass stars, are examined using stellar evolution models, simulations of supernova ejecta colliding with low-mass companions, and comparisons with observed data. These first star survivors will be observed as metal-rich halo stars in our Galaxy. We may have identified a candidate star in the observational database where elemental abundances and kinematic data are available. Our models also account for the existence of several solar-metallicity stars in the literature having space velocities equivalent to the halo population. The proposed scenario demands a new channel of star formation in the early universe and is a supplementary scenario for the origin of the known metal-poor stars. △ Less

Submitted 25 March, 2021; originally announced March 2021.

Comments: 41 pages, 14 figures, accepted for publication in PASJ (to be published in open access)

SIRIUS Project. III. Star-by-star simulations of star cluster formation using a direct N-body integrator with stellar feedback

Authors: Michiko S. Fujii, Takayuki R. Saitoh, Yutaka Hirai, Long Wang

Abstract: One of the computational challenges of cluster formation simulations is resolving individual stars and simulating massive clusters with masses of more than $10^4 M_{\odot}$ without gravitational softening. Combining direct $N$-body code with smoothed-particle hydrodynamics (SPH) code, we have developed a new code, \textsc{ASURA+BRIDGE}, in which we can integrate stellar particles without softening… ▽ More One of the computational challenges of cluster formation simulations is resolving individual stars and simulating massive clusters with masses of more than $10^4 M_{\odot}$ without gravitational softening. Combining direct $N$-body code with smoothed-particle hydrodynamics (SPH) code, we have developed a new code, \textsc{ASURA+BRIDGE}, in which we can integrate stellar particles without softening. We add a feedback model for \HII regions into this code, in which thermal and momentum feedback is given within the Str{ö}mgren radius. We perform $N$-body/SPH simulations of star cluster formation. Without softening, a portion of massive stars are ejected from the forming clusters. As a result, the stellar feedback works outside the clusters. This enhances/suppresses the star formation in initially sub-virial/super-virial clouds. We find that the formed star clusters are denser than currently observed open clusters, but the mass--density relation is consistent with or even higher than that is estimated as an initial cluster density. We also find that some clusters have multiple peaks in their stellar age distribution as a consequence of their hierarchical formation. △ Less

Submitted 5 July, 2021; v1 submitted 3 March, 2021; originally announced March 2021.

Comments: 25 pages, 29 figures, accepted for publication in PASJ

SIRIUS project II: a new tree-direct hybrid code for smoothed particle hydrodynamics/N-body simulations of star clusters

Authors: M. S. Fujii, T. R. Saitoh, L. Wang, Y. Hirai

Abstract: Star clusters form via clustering star formation inside molecular clouds. In order to understand the dynamical evolution of star clusters in their early phase, in which star clusters are still embedded in their surrounding gas, we need an accurate integration of individual stellar orbits without gravitational softening in the systems including both gas and stars, as well as modeling individual sta… ▽ More Star clusters form via clustering star formation inside molecular clouds. In order to understand the dynamical evolution of star clusters in their early phase, in which star clusters are still embedded in their surrounding gas, we need an accurate integration of individual stellar orbits without gravitational softening in the systems including both gas and stars, as well as modeling individual stars with a realistic mass function. We develop a new tree-direct hybrid smoothed particle hydrodynamics/N-body code, ASURA+BRIDGE, in which stars are integrated using a direct N-body scheme or PeTar, a particle-particle particle-tree scheme code, without gravitational softening. In ASURA+BRIDGE, stars are assumed to have masses randomly drawn from a given initial mass function. With this code, we perform star-cluster formation simulations starting from molecular clouds without gravitational softening. We find that artificial dense cores in star-cluster centers due to the softening disappear when we do not use softening. We further demonstrate that star clusters are built up via mergers of smaller clumps. Star clusters formed in our simulations include some dynamically formed binaries with the minimum semi-major axes of a few au, and the binary fraction is higher for more massive stars. △ Less

Submitted 15 April, 2021; v1 submitted 14 January, 2021; originally announced January 2021.

Comments: 17 pages, 18 figures, 3 tables, accepted for publication in PASJ

SIRIUS project. I. Star formation models for star-by-star simulations of star clusters and galaxy formation

Authors: Yutaka Hirai, Michiko S. Fujii, Takayuki R. Saitoh

Abstract: Most stars are formed as star clusters in galaxies, which then disperse into galactic disks. Upcoming exascale supercomputational facilities will enable performing simulations of galaxies and their formation by resolving individual stars (star-by-star simulations). This will substantially advance our understanding of star formation in galaxies, star cluster formation, and assembly histories of gal… ▽ More Most stars are formed as star clusters in galaxies, which then disperse into galactic disks. Upcoming exascale supercomputational facilities will enable performing simulations of galaxies and their formation by resolving individual stars (star-by-star simulations). This will substantially advance our understanding of star formation in galaxies, star cluster formation, and assembly histories of galaxies. In previous galaxy simulations, a simple stellar population approximation was used. It is, however, difficult to improve the mass resolution with this approximation. Therefore, a model for forming individual stars that can be used in simulations of galaxies must be established. In this first paper of a series of the SIRIUS (SImulations Resolving IndividUal Stars) project, we demonstrate a stochastic star formation model for star-by-star simulations. An assumed stellar initial mass function (IMF) is randomly assigned to newly formed stars. We introduce a maximum search radius to assemble the mass from surrounding gas particles to form star particles. In this study, we perform a series of N-body/smoothed particle hydrodynamics simulations of star cluster formations from turbulent molecular clouds and ultra-faint dwarf galaxies as test cases. The IMF can be correctly sampled if a maximum search radius that is larger than the value estimated from the threshold density for star formation is adopted. In small clouds, the formation of massive stars is highly stochastic because of the small number of stars. We confirm that the star formation efficiency and threshold density do not strongly affect the results. We find that our model can naturally reproduce the relationship between the most massive stars and the total stellar mass of star clusters. Herein, we demonstrate that our models can be applied to simulations varying from star clusters to galaxies for a wide range of resolutions. △ Less

Submitted 17 April, 2021; v1 submitted 25 May, 2020; originally announced May 2020.

Comments: 23 pages, 22 figures, accepted for publication in PASJ

Enrichment of Strontium in Dwarf Galaxies

Authors: Yutaka Hirai, Shinya Wanajo, Takayuki R. Saitoh

Abstract: Light trans-iron elements such as Sr serve as the key to understanding the astrophysical sites of heavy elements. Spectroscopic studies of metal-poor stars have revealed large star-to-star scatters in the ratios of [Sr/Ba], which indicates that there are multiple sites for the production of Sr. Here we present the enrichment history of Sr by a series of the $N$-body/smoothed particle hydrodynamics… ▽ More Light trans-iron elements such as Sr serve as the key to understanding the astrophysical sites of heavy elements. Spectroscopic studies of metal-poor stars have revealed large star-to-star scatters in the ratios of [Sr/Ba], which indicates that there are multiple sites for the production of Sr. Here we present the enrichment history of Sr by a series of the $N$-body/smoothed particle hydrodynamics simulations of a dwarf galaxy with a stellar mass of 3 $\times$ 10$^{6}$ $M_{\odot}$. We show that binary neutron star mergers (NSMs) and asymptotic giant branch (AGB) stars contribute to the enrichment of Sr in the metallicity ranges [Fe/H] $\gtrsim$ $-$3 and [Fe/H] $\gtrsim$ $-$1, respectively. It appears insufficient, however, to explain the overall observational trends of Sr by considering only these sites. We find that the models including electron-capture supernovae (ECSNe) and rotating massive stars (RMSs), in addition to NSMs and AGBs, reasonably reproduce the enrichment histories of Sr in dwarf galaxies. The contributions of both ECSNe and NSMs make scatters of $\approx$ 0.2 dex in [Sr/Fe], [Sr/Ba], and [Sr/Zn] as can be seen for observed stars in the metallicity range [Fe/H] $<$ $-2$. We also find that the mass range of ECSN progenitors should be substantially smaller than $1\, M_\odot$ (e.g., 0.1-$0.2\, M_\odot$) to avoid over-prediction of [Sr/Ba] and [Sr/Zn] ratios. Our results demonstrate that NSMs, AGBs, ECSNe, and RMSs all play roles in the enrichment histories of Local Group dwarf galaxies, although more observational data are required to disentangle the relative contributions of these sources. △ Less

Submitted 19 September, 2019; originally announced September 2019.

Comments: 19 pages, 10 figures, accepted for publication in The Astrophysical Journal

Journal ref: ApJ, 885, 33, 2019

Radial Velocity Survey of Nearby OB Stars

Authors: Yuki Moritani, Takuma Suda, Toshikazu Shigeyama, Takayuki R. Saitoh

Abstract: We report on the current status of the radial velocity monitoring of nearby OB stars to look for binaries with small mass ratios. The combined data of radial velocities using the domestic 1-2 m-class telescopes seems to confirm the variations of radial velocities in a few weeks for four out of ten target single-lined spectroscopic binaries. More data are needed to estimate the exact periods and ma… ▽ More We report on the current status of the radial velocity monitoring of nearby OB stars to look for binaries with small mass ratios. The combined data of radial velocities using the domestic 1-2 m-class telescopes seems to confirm the variations of radial velocities in a few weeks for four out of ten target single-lined spectroscopic binaries. More data are needed to estimate the exact periods and mass distributions. △ Less

Submitted 16 December, 2018; originally announced December 2018.

Comments: 7 pages, 4 figures, accepted for publication in Stars and Galaxies

Property of Young Massive Clusters in a Galaxy-Galaxy Merger Remnant

Authors: Hidenori Matsui, Ataru Tanikawa, Takayuki R Saitoh

Abstract: We investigate the properties of young massive clusters (YMCs) in a galaxy-galaxy merger remnant by analyzing the data obtained by a gas rich major merger simulations in Matsui et al. 2012. We found that the YMCs are distributed at a few $\rm kpc$ and at $\sim 10~{\rm kpc}$ from the galactic center, in other words, there are two components of their distribution. The former are formed in filamentar… ▽ More We investigate the properties of young massive clusters (YMCs) in a galaxy-galaxy merger remnant by analyzing the data obtained by a gas rich major merger simulations in Matsui et al. 2012. We found that the YMCs are distributed at a few $\rm kpc$ and at $\sim 10~{\rm kpc}$ from the galactic center, in other words, there are two components of their distribution. The former are formed in filamentary and turbulent gas generated at a few $\rm kpc$ from the center because of galaxy encounters, and the latter are formed in tidal tails which are far from the center. The YMCs are much less concentrated than galaxy stars. The mass function of the YMCs is $dN/dM \propto M^{-2}$. Most of YMCs are formed from the second encounter to the final coalescence phase of the galactic cores, and their formation rate is especially high at final coalescence phase. Most of them consists of single stellar population in age, but YMCs with multi stellar populations in age are also formed. The multiple populations are produced by the following process: a YMC captures dense gas, and another generation stars form within the cluster. There are several YMCs formed in an isolated disk before the encounter of galaxies. These candidates contain stars with various age by capturing dense gas and forming stars. YMCs in a merger remnant, have various orbits, but large fraction of candidates have circular orbits. △ Less

Submitted 11 November, 2018; originally announced November 2018.

Comments: 10 pages, 11 figures, accepted for publication in PASJ

Discovery of a Molecular Collision Front in Interacting Galaxies NGC 4567/4568 with ALMA

Authors: Hiroyuki Kaneko, Nario Kuno, Takayuki R. Saitoh

Abstract: We present results of $^{12}$CO(J = 1-0) imaging observations of NGC 4567/4568, a galaxy pair in a close encounter, with Atacama Large Millimetre/submillimetre Array. For the first time, we find clear evidence of a molecular collision front whose velocity dispersion is 16.8$\pm$1.4 km s$^{-1}$ at the overlapping region owing to high spatial and velocity resolution. By integrating over the velocity… ▽ More We present results of $^{12}$CO(J = 1-0) imaging observations of NGC 4567/4568, a galaxy pair in a close encounter, with Atacama Large Millimetre/submillimetre Array. For the first time, we find clear evidence of a molecular collision front whose velocity dispersion is 16.8$\pm$1.4 km s$^{-1}$ at the overlapping region owing to high spatial and velocity resolution. By integrating over the velocity width that corresponds to the molecular collision front, we find a long filamentary structure with a size of 1800 pc $\times$ 350 pc at the collision front. This filamentary molecular structure spatially coincides with a dark-lane seen in R-band image. We find four molecular clouds in the filament, each with a radius of 30 pc and mass of 10$^{6} M_{\odot}$; the radii matching a typical value for giant molecular clouds and the masses corresponding to those between giant molecular clouds and giant molecular associations. All four clouds are gravitationally bound. The molecular filamentary structure and its physical conditions are similar to the structure expected via numerical simulation. The filament could be a progenitor of super star clusters. △ Less

Submitted 29 May, 2018; originally announced May 2018.

Comments: 6 pages, 4 figures, 1 table; accepted for publication in ApJL

Compact object mergers driven by gas fallback

Authors: Hiromichi Tagawa, Takayuki R. Saitoh, Bence Kocsis

Abstract: Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel… ▽ More Recently several gravitational wave detections have shown evidence for compact object mergers. However, the astrophysical origin of merging binaries is not well understood. Stellar binaries are typically at much larger separations than what is needed for the binaries to merge due to gravitational wave emission, which leads to the so-called final AU problem. In this letter we propose a new channel for mergers of compact object binaries which solves the final AU problem. We examine the binary evolution following gas expansion due to a weak failed supernova explosion, neutrino mass loss, core disturbance, or envelope instability. In such situations the binary is possibly hardened by ambient gas. We investigate the evolution of the binary system after a shock has propagated by performing smoothed particle hydrodynamics simulations. We find that significant binary hardening occurs when the gas mass bound to the binary exceeds that of the compact objects. This mechanism represents a new possibility for the pathway to mergers for gravitational wave events. △ Less

Submitted 26 June, 2018; v1 submitted 1 February, 2018; originally announced February 2018.

Comments: 6 pages, 3 figures, version published in Physical Review Letters

Journal ref: Phys. Rev. Lett. 120, 261101 (2018)

Enrichment of Zinc in galactic chemodynamical evolution models

Authors: Yutaka Hirai, Takayuki R. Saitoh, Yuhri Ishimaru, Shinya Wanajo

Abstract: The heaviest iron-peak element, Zn has been used as an important tracer of cosmic chemical evolution. Spectroscopic observations of the metal-poor stars in Local Group galaxies show that an increasing trend of [Zn/Fe] ratios toward lower metallicity. However, enrichment of Zn in galaxies is not well understood due to the poor knowledge of astrophysical sites of Zn as well as metal mixing in galaxi… ▽ More The heaviest iron-peak element, Zn has been used as an important tracer of cosmic chemical evolution. Spectroscopic observations of the metal-poor stars in Local Group galaxies show that an increasing trend of [Zn/Fe] ratios toward lower metallicity. However, enrichment of Zn in galaxies is not well understood due to the poor knowledge of astrophysical sites of Zn as well as metal mixing in galaxies. Here we show possible explanations for the observed trend by taking into account electron-capture supernovae (ECSNe) as one of the sources of Zn in our chemodynamical simulations of dwarf galaxies. We find that the ejecta from ECSNe contribute to stars with [Zn/Fe] $\gtrsim$ 0.5. We also find that scatters of [Zn/Fe] in higher metallicity originate from the ejecta of type Ia supernovae. On the other hand, it appears difficult to explain the observed trends if we do not consider ECSNe as a source of Zn. These results come from inhomogeneous spatial metallicity distribution due to the inefficiency of metal mixing. We find that the optimal value of scaling factor for metal diffusion coefficient is $\sim$ 0.01 in the shear-based metal mixing model in smoothed particle hydrodynamics simulations. These results suggest that ECSNe can be one of the contributors to the enrichment of Zn in galaxies. △ Less

Submitted 24 January, 2018; originally announced January 2018.

Comments: 16 pages, 15 figures, accepted for publication in The Astrophysical Journal

Hydrodynamic simulations of the central molecular zone with realistic Galactic potential

Authors: Jihye Shin, Sungsoo S. Kim, Junichi Baba, Takayuki R. Saitoh, Jeong-Sun Hwang, Kyungwon Chun, Shunsuke Hozumi

Abstract: We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different ba… ▽ More We present hydrodynamic simulations of gas clouds inflowing from the disk to a few hundred parsec region of the Milky Way. A gravitational potential is generated to include realistic Galactic structures by using thousands of multipole expansions that describe 6.4 million stellar particles of a self-consistent Galaxy simulation. We find that a hybrid multipole expansion model, with two different basis sets and a thick disk correction, accurately reproduces the overall structures of the Milky Way. Through non-axisymmetric Galactic structures of an elongated bar and spiral arms, gas clouds in the disk inflow to the nuclear region and form a central molecular zone (CMZ)-like nuclear ring. We find that the size of the nuclear ring evolves into ~240 pc at T~1500 Myr, regardless of the initial size. For most simulation runs, the rate of gas inflow to the nuclear region is equilibrated to ~0.02 M_sun/yr. The nuclear ring is off-centered, relative to the Galactic center, by the lopsided central mass distribution of the Galaxy model, and thus an asymmetric mass distribution of the nuclear ring arises accordingly. The vertical asymmetry of the the Galaxy model also causes the nuclear ring to be tilted along the Galactic plane. During the first ~100 Myr, the vertical frequency of the gas motion is twice that of the orbital frequency, thus the projected nuclear ring shows a twisted, infinity-like shape. △ Less

Submitted 28 April, 2017; originally announced April 2017.

Comments: Accepted for publication in ApJ

Efficiency of metal mixing in dwarf galaxies

Authors: Yutaka Hirai, Takayuki R. Saitoh

Abstract: Metal mixing plays critical roles in the enrichment of metals in galaxies. The abundance of elements such as Mg, Fe, and Ba in metal-poor stars help us understand the metal mixing in galaxies. However, the efficiency of metal mixing in galaxies is not yet understood. Here we report a series of $N$-body/smoothed particle hydrodynamics simulations of dwarf galaxies with different efficiencies of met… ▽ More Metal mixing plays critical roles in the enrichment of metals in galaxies. The abundance of elements such as Mg, Fe, and Ba in metal-poor stars help us understand the metal mixing in galaxies. However, the efficiency of metal mixing in galaxies is not yet understood. Here we report a series of $N$-body/smoothed particle hydrodynamics simulations of dwarf galaxies with different efficiencies of metal mixing using turbulence-induced mixing model. We show that metal mixing apparently occurs in dwarf galaxies from Mg and Ba abundance. We find that the scaling factor for metal diffusion larger than 0.01 is necessary to reproduce the observation of Ba abundance in dwarf galaxies. This value is consistent with the value expected from turbulence theory and experiment. We also find that timescale of metal mixing is less than 40 Myr. This timescale is shorter than that of typical dynamical times of dwarf galaxies. We demonstrate that the determination of a degree of scatters of Ba abundance by the observation will help us to constrain the efficiency of metal mixing more precisely. △ Less

Submitted 16 March, 2017; originally announced March 2017.

Comments: 7 pages, 2 figures, accepted for publication in ApJL

Imprints of Zero-Age Velocity Dispersions and Dynamical Heating on the Age-Velocity dispersion Relation

Authors: Jun Kumamoto, Junichi Baba, Takayuki R. Saitoh

Abstract: Observations of stars in the the solar vicinity show a clear tendency for old stars to have larger velocity dispersions. This relation is called the age-velocity dispersion relation (AVR) and it is believed to provide insight into the heating history of the Milky Way galaxy. Here, in order to investigate the origin of the AVR, we performed smoothed particle hydrodynamic simulations of the self-gra… ▽ More Observations of stars in the the solar vicinity show a clear tendency for old stars to have larger velocity dispersions. This relation is called the age-velocity dispersion relation (AVR) and it is believed to provide insight into the heating history of the Milky Way galaxy. Here, in order to investigate the origin of the AVR, we performed smoothed particle hydrodynamic simulations of the self-gravitating multiphase gas disks in the static disk-halo potentials. Star formation from cold and dense gas is taken into account, and we analyze the evolution of these star particles. We find that exponents of simulated AVR and the ratio of the radial to vertical velocity dispersion are close to the observed values. We also find that the simulated AVR is not a simple consequence of dynamical heating. The evolution tracks of stars with different epochs evolve gradually in the age-velocity dispersion plane as a result of: (1) the decrease in velocity dispersion in star forming regions, and (2) the decrease in the number of cold/dense/gas as scattering sources. These results suggest that the AVR involves not only the heating history of a stellar disk, but also the historical evolution of the ISM in a galaxy. △ Less

Submitted 13 January, 2017; originally announced January 2017.

Comments: 13 pages, 14 figures, accepted for publication in PASJ

Early chemo-dynamical evolution of dwarf galaxies deduced from enrichment of r-process elements

Authors: Yutaka Hirai, Yuhri Ishimaru, Takayuki R. Saitoh, Michiko S. Fujii, Jun Hidaka, Toshitaka Kajino

Abstract: The abundance of elements synthesized by the rapid neutron-capture process (r-process elements) of extremely metal-poor (EMP) stars in the Local Group galaxies gives us clues to clarify the early evolutionary history of the Milky Way halo. The Local Group dwarf galaxies would have similarly evolved with building blocks of the Milky Way halo. However, how the chemo-dynamical evolution of the buildi… ▽ More The abundance of elements synthesized by the rapid neutron-capture process (r-process elements) of extremely metal-poor (EMP) stars in the Local Group galaxies gives us clues to clarify the early evolutionary history of the Milky Way halo. The Local Group dwarf galaxies would have similarly evolved with building blocks of the Milky Way halo. However, how the chemo-dynamical evolution of the building blocks affects the abundance of r-process elements is not yet clear. In this paper, we perform a series of simulations using dwarf galaxy models with various dynamical times and total mass, which determine star-formation histories. We find that galaxies with dynamical times longer than 100 Myr have star formation rates less than $10^{-3} M_{\odot}$ yr$^{-1}$ and slowly enrich metals in their early phase. These galaxies can explain the observed large scatters of r-process abundance in EMP stars in the Milky Way halo regardless of their total mass. On the other hand, the first neutron star merger appears at a higher metallicity in galaxies with a dynamical time shorter than typical neutron star merger times. The scatters of r-process elements mainly come from inhomogeneity of the metals in the interstellar medium whereas the scatters of $α$-elements are mostly due to the difference in the yield of each supernova. Our results demonstrate that the future observations of r-process elements in EMP stars will be able to constrain the early chemo-dynamical evolution of the Local Group galaxies. △ Less

Submitted 20 December, 2016; originally announced December 2016.

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

Chemical evolution library for galaxy formation simulation

Authors: Takayuki R. Saitoh

Abstract: We have developed a software library for chemical evolution simulations of galaxy formation under the simple stellar population (SSP) approximation. In this library, all of the necessary components concerning chemical evolution, such as initial mass functions, stellar lifetimes, yields from type II and Ia supernovae, asymptotic giant branch stars, and neutron star mergers, are compiled from the li… ▽ More We have developed a software library for chemical evolution simulations of galaxy formation under the simple stellar population (SSP) approximation. In this library, all of the necessary components concerning chemical evolution, such as initial mass functions, stellar lifetimes, yields from type II and Ia supernovae, asymptotic giant branch stars, and neutron star mergers, are compiled from the literature. Various models are pre-implemented in this library so that users can choose their favorite combination of models. Subroutines of this library return released energy and masses of individual elements depending on a given event type. Since the redistribution manner of these quantities depends on the implementation of users' simulation codes, this library leaves it up to the simulation code. As demonstrations, we carry out both one-zone, closed box simulations and three-dimensional simulations of a collapsing gas and dark matter system using this library. In these simulations, we can easily compare the impact of individual models on the chemical evolution of galaxies, just by changing the control flags and parameters of the library. Since this library only deals with the part of chemical evolution under the SSP approximation, any simulation codes that use the SSP approximation -- namely particle-base and mesh codes, as well as semi-analytical models -- can use it. This library is named "CELib" after the term "Chemical Evolution Library" and is made available to the community. △ Less

Submitted 7 December, 2016; originally announced December 2016.

Comments: 36 pages, 41 figrues. Accepted for publication in AJ. The distribution site of CELib is https://bitbucket.org/tsaitoh/celib

Eventful Evolution of Giant Molecular Clouds in Dynamically Evolving Spiral Arms

Authors: Junichi Baba, Kana Morokuma-Matsui, Takayuki R. Saitoh

Abstract: The formation and evolution of giant molecular clouds (GMCs) in spiral galaxies have been investigated in the traditional framework of the combined quasi-stationary density wave and galactic shock model. However, our understanding of the dynamics of spiral arms is changing from the traditional spiral model to a dynamically evolving spiral model. In this study, we investigate the structure and evol… ▽ More The formation and evolution of giant molecular clouds (GMCs) in spiral galaxies have been investigated in the traditional framework of the combined quasi-stationary density wave and galactic shock model. However, our understanding of the dynamics of spiral arms is changing from the traditional spiral model to a dynamically evolving spiral model. In this study, we investigate the structure and evolution of GMCs in a dynamically evolving spiral arm using a three-dimensional N-body/hydrodynamic simulation of a barred spiral galaxy at parsec-scale resolution. This simulation incorporated self-gravity, molecular hydrogen formation, radiative cooling, heating due to interstellar far-ultraviolet radiation, and stellar feedback by both HII regions and Type-II supernovae. In contrast to a simple expectation based on the traditional spiral model, the GMCs exhibited no systematic evolutionary sequence across the spiral arm. Our simulation showed that the GMCs behaved as highly dynamic objects with eventful lives involving collisional build-up, collision-induced star formation, and destruction via stellar feedback. The GMC lifetimes were predicted to be short, only a few tens of millions years. We also found that, at least at the resolutions and with the feedback models used in this study, most of the GMCs without HII regions were collapsing, but half of the GMCs with HII regions were expanding owing to the HII-region feedback from stars within them. Our results support the dynamic and feedback-regulated GMC evolution scenario. Although the simulated GMCs were converging rather than virial equilibrium, they followed the observed scaling relationship well. We also analysed the effects of galactic tides and external pressure on GMC evolution and suggested that GMCs cannot be regarded as isolated systems since their evolution in disc galaxies is complicated because of these environmental effects. △ Less

Submitted 15 September, 2016; v1 submitted 31 August, 2016; originally announced September 2016.

Comments: 19 pages, 16 figures. Accepted for publication in MNRAS

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/

Enrichment of r-process elements in dwarf spheroidal galaxies in chemo-dynamical evolution model

Authors: Yutaka Hirai, Yuhri Ishimaru, Takayuki R. Saitoh, Michiko S. Fujii, Jun Hidaka, Toshitaka Kajino

Abstract: The rapid neutron-capture process (r-process) is a major process to synthesize elements heavier than iron, but the astrophysical site(s) of r-process is not identified yet. Neutron star mergers (NSMs) are suggested to be a major r-process site from nucleosynthesis studies. Previous chemical evolution studies however require unlikely short merger time of NSMs to reproduce the observed large star-to… ▽ More The rapid neutron-capture process (r-process) is a major process to synthesize elements heavier than iron, but the astrophysical site(s) of r-process is not identified yet. Neutron star mergers (NSMs) are suggested to be a major r-process site from nucleosynthesis studies. Previous chemical evolution studies however require unlikely short merger time of NSMs to reproduce the observed large star-to-star scatters in the abundance ratios of r-process elements relative to iron, [Eu/Fe], of extremely metal-poor stars in the Milky Way (MW) halo. This problem can be solved by considering chemical evolution in dwarf spheroidal galaxies (dSphs) which would be building blocks of the MW and have lower star formation efficiencies than the MW halo. We demonstrate that enrichment of r-process elements in dSphs by NSMs using an N-body/smoothed particle hydrodynamics code. Our high-resolution model reproduces the observed [Eu/Fe] by NSMs with a merger time of 100 Myr when the effect of metal mixing is taken into account. This is because metallicity is not correlated with time up to ~ 300 Myr from the start of the simulation due to low star formation efficiency in dSphs. We also confirm that this model is consistent with observed properties of dSphs such as radial profiles and metallicity distribution. The merger time and the Galactic rate of NSMs are suggested to be <~ 300 Myr and ~ $10^{-4}$ yr$^{-1}$, which are consistent with the values suggested by population synthesis and nucleosynthesis studies. This study supports that NSMs are the major astrophysical site of r-process. △ Less

Submitted 29 September, 2015; originally announced September 2015.

Comments: 16 pages, 16 figures, accepted for publication in ApJ

Journal ref: ApJ 814 (2015) 41

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

Properties of thick discs formed in clumpy galaxies

Authors: Shigeki Inoue, Takayuki R. Saitoh

Abstract: We examine a possible formation scenario of galactic thick discs with numerical simulations. Thick discs have previously been argued to form in clumpy disc phase in the high-redshift Universe, which host giant clumps of <10^9 M_sun in their highly gas-rich discs. We performed SPH simulations using isolated galaxy models for the purpose of verifying whether dynamical and chemical properties of the… ▽ More We examine a possible formation scenario of galactic thick discs with numerical simulations. Thick discs have previously been argued to form in clumpy disc phase in the high-redshift Universe, which host giant clumps of <10^9 M_sun in their highly gas-rich discs. We performed SPH simulations using isolated galaxy models for the purpose of verifying whether dynamical and chemical properties of the thick discs formed in such clumpy galaxies are compatible with observations. The results of our simulations seem nearly consistent with observations in dynamical properties such as radial and vertical density profiles, significant rotation velocity lag with height and distributions of orbital eccentricities. In addition, the thick discs in our simulations indicate nearly exponential dependence of azimuthal and vertical velocity dispersions with radius, nearly isothermal kinematics in vertical direction and negligible metallicity gradients in radial and vertical directions. However, our simulations cannot reproduce altitudinal dependence of eccentricities, metallicity relations with eccentricities or rotation velocities, which shows striking discrepancy from recent observations of the Galactic thick disc. From this result, we infer that the clumpy disc scenario for thick-disc formation would not be suitable at least for the Milky Way. Our study, however, cannot reject this scenario for external galaxies if not all galaxies form their thick discs by the same process. In addition, we found that a large fraction of thick-disc stars forms in giant clumps. △ Less

Submitted 17 March, 2014; v1 submitted 24 February, 2014; originally announced February 2014.

Comments: 15 pages, 13 figures, 3 tables, 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/

Dynamics of Non-Steady Spiral Arms in Disk Galaxies

Authors: Junichi Baba, Takayuki R. Saitoh, Keiichi Wada

Abstract: In order to understand the physical mechanisms underlying non-steady stellar spiral arms in disk galaxies, we analyzed the growing and damping phases of their spiral arms using three-dimensional $N$-body simulations. We confirmed that the spiral arms are formed due to a swing amplification mechanism that reinforces density enhancement as a seeded wake. In the damping phase, the Coriolis force exer… ▽ More In order to understand the physical mechanisms underlying non-steady stellar spiral arms in disk galaxies, we analyzed the growing and damping phases of their spiral arms using three-dimensional $N$-body simulations. We confirmed that the spiral arms are formed due to a swing amplification mechanism that reinforces density enhancement as a seeded wake. In the damping phase, the Coriolis force exerted on a portion of the arm surpasses the gravitational force that acts to shrink the portion. Consequently, the stars in the portion escape from the arm, and subsequently they form a new arm at a different location. The time-dependent nature of the spiral arms are originated in the continual repetition of this non-linear phenomenon. Since a spiral arm does not rigidly rotate, but follows the galactic differential rotation, the stars in the arm rotate at almost the same rate as the arm. In other words, every single position in the arm can be regarded as the co-rotation point. Due to interaction with their host arms, the energy and angular momentum of the stars change, thereby causing the radial migration of the stars. During this process, the kinetic energy of random motion (random energy) of the stars does not significantly increase, and the disk remains dynamically cold. Owing to this low degree of disk heating, the short-lived spiral arms can recurrently develop over many rotational periods. The resultant structure of the spiral arms in the $N$-body simulations is consistent with some observational nature of spiral galaxies. We conclude that the formation and structure of spiral arms in isolated disk galaxies can be reasonably understood by non-linear interactions between a spiral arm and its constituent stars. △ Less

Submitted 22 November, 2012; originally announced November 2012.

Comments: 16 pages (in emulateapj format), incl. 13 figures and 1 table; Accepted for publication in The Astrophysical Journal

The Schmidt-Kennicutt Law of Matched-Age Star Forming Regions; Pa-alpha Observations of the Early-Phase Interacting Galaxy Taffy I

Authors: S. Komugi, K. Tateuchi, K. Motohara, T. Takagi, D. Iono, H. Kaneko, J. Ueda, T. R. Saitoh, N. Kato, M. Konishi, S. Koshida, T. Morokuma, H. Takahashi, T. Tanabe, Y. Yoshii

Abstract: In order to test a recent hypothesis that the dispersion in the Schmidt-Kennicutt law arises from variations in the evolutionary stage of star forming molecular clouds, we compared molecular gas and recent star formation in an early-phase merger galaxy pair, Taffy I (UGC\ 12915/UGC\ 12914, VV\ 254) which went through a direct collision 20 Myr ago and whose star forming regions are expected to have… ▽ More In order to test a recent hypothesis that the dispersion in the Schmidt-Kennicutt law arises from variations in the evolutionary stage of star forming molecular clouds, we compared molecular gas and recent star formation in an early-phase merger galaxy pair, Taffy I (UGC\ 12915/UGC\ 12914, VV\ 254) which went through a direct collision 20 Myr ago and whose star forming regions are expected to have similar ages. Narrow-band Pa-alpha image is obtained using the ANIR near-infrared camera on the mini-TAO 1m telescope. The image enables us to derive accurate star formation rates within the galaxy directly. The total star formation rate, 22.2 M_sun/yr, was found to be much higher than previous estimates. Ages of individual star forming blobs estimated from equivalent widths indicate that most star forming regions are ~7 Myr old, except for a giant HII region at the bridge which is much younger. Comparison between star formation rates and molecular gas masses for the regions with the same age exhibits a surprisingly tight correlation, a slope of unity, and star formation efficiencies comparable to those of starburst galaxies. These results suggest that Taffy I has just evolved into a starburst system after the collision, and the star forming sites are at a similar stage in their evolution from natal molecular clouds except for the bridge region. The tight Schmidt-Kennicutt law supports the scenario that dispersion in the star formation law is in large part due to differences in evolutionary stage of star forming regions. △ Less

Submitted 1 August, 2012; originally announced August 2012.

Comments: Accepted for publication in Astrophysical Journal

The formation of young dense star clusters through mergers

Authors: Michiko S. Fujii, Takayuki R. Saitoh, Simon F. Portegies Zwart

Abstract: Young star clusters like R136 in the Large Magellanic Cloud and NGC 3603, Westerlund 1, and 2 in the Milky Way are dynamically more evolved than expected based on their current relaxation times. In particular, the combination of a high degree of mass segregation, a relatively low central density, and the large number of massive runaway stars in their vicinity are hard to explain with the monolithi… ▽ More Young star clusters like R136 in the Large Magellanic Cloud and NGC 3603, Westerlund 1, and 2 in the Milky Way are dynamically more evolved than expected based on their current relaxation times. In particular, the combination of a high degree of mass segregation, a relatively low central density, and the large number of massive runaway stars in their vicinity are hard to explain with the monolithic formation of these clusters. Young star clusters can achieve such a mature dynamical state if they formed through the mergers of a number of less massive clusters. The shorter relaxation times of less massive clusters cause them to dynamically evolve further by the time they merge, and the merger product preserves the memory of the dynamical evolution of its constituent clusters. With a series of $N$-body simulations, we study the dynamical evolution of single massive clusters and those that are assembled through merging smaller clusters together. We find that the formation of massive star clusters through the mergers of smaller clusters can reproduce the currently observed spatial distribution of massive stars, the density, and the characteristics (number and mass distribution) of the stars ejected as runaways from young dense clusters. We therefore conclude that these clusters and possibly other young massive star clusters formed through the mergers of smaller clusters. △ Less

Submitted 7 May, 2012; originally announced May 2012.

Comments: 13 pages, 13 figures, accepted for publication in ApJ

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

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

Natures of a clump-origin bulge: a pseudobulge-like but old metal-rich bulge

Authors: Shigeki Inoue, Takayuki R. Saitoh

Abstract: Bulges in spiral galaxies have been supposed to be classified into two types: classical bulges or pseudobulges. Classical bulges are thought to form by galactic merger with bursty star formation, whereas pseudobulges are suggested to form by secular evolution due to spiral arms and a barred structure funneling gas into the galactic centre. Noguchi (1998,1999) suggested another bulge formation scen… ▽ More Bulges in spiral galaxies have been supposed to be classified into two types: classical bulges or pseudobulges. Classical bulges are thought to form by galactic merger with bursty star formation, whereas pseudobulges are suggested to form by secular evolution due to spiral arms and a barred structure funneling gas into the galactic centre. Noguchi (1998,1999) suggested another bulge formation scenario, `clump-origin bulge'. He demonstrated using a numerical simulation that a galactic disc suffers dynamical instability to form clumpy structures in the early stage of disc formation since the premature disc is expected to be highly gas-rich, then the clumps are sucked into the galactic centre by dynamical friction and merge into a single bulge at the centre. This bulge formation scenario, which is expected to happen only at the high-redshift, is different from the galactic merger and the secular evolution. Therefore, clump-origin bulges may have their own unique properties. We perform a high-resolution N-body/smoothed particle hydrodynamics (SPH) simulation for the formation of the clump-origin bulge in an isolated galaxy model and study dynamical and chemical properties of the clump-origin bulge. We find that the clump-origin bulge resembles pseudobulges in dynamical properties, a nearly exponential surface density profile, a barred boxy shape and a significant rotation. We also find that this bulge consists of old and metal-rich stars, displaying resemblance to classical bulges. These natures, old metal-rich population but pseudobulge-like structures, mean that the clump-origin bulge can not be simply classified into classical bulges nor pseudobulges. From these results, we discuss similarities of the clump-origin bulge to the Milky Way bulge. △ Less

Submitted 6 December, 2011; v1 submitted 13 September, 2011; originally announced September 2011.

Comments: 13 pages, 10 figures, Accepted for publication in MNRAS

Shallowed cusp slope of dark matter in disc galaxy formation through clump clusters

Authors: Shigeki Inoue, Takayuki R. Saitoh

Abstract: Cusp-core problem is a controversial problem on galactic dark matter haloes. Cosmological N-body simulations has demonstrated that galactic dark matter haloes have a cuspy density profile at the centre. However, baryonic physics may affect the dark matter density profile. For example, it was suggested that adiabatic contraction of baryonic gas makes the dark matter cusp steeper. However, it is sti… ▽ More Cusp-core problem is a controversial problem on galactic dark matter haloes. Cosmological N-body simulations has demonstrated that galactic dark matter haloes have a cuspy density profile at the centre. However, baryonic physics may affect the dark matter density profile. For example, it was suggested that adiabatic contraction of baryonic gas makes the dark matter cusp steeper. However, it is still an open question if the gas falls into the galactic centre in smooth adiabatic manner. Recent numerical studies suggested that disc galaxies might experience clumpy phase in their early stage of the disc formation, which could also explain clump clusters and chain galaxies observed in high redshift Universe. In this letter, using numerical simulations with an isolated model, we study how the dark matter halo responds to these clumpy nature of baryon component in disc galaxy formation through the clump cluster phase. Our simulation demonstrates that such clumpy phase leads to a shallower density profile of the dark matter halo in the central region while clumps fall into the centre due to dynamical friction. This mechanism helps to make the central dark matter density profile shallower in the galaxies whose virial mass is as large as 5.0*10^11 solar masses. This phenomenon is caused by reaction to dynamical friction of the stellar clumps against the dark matter halo. The halo draws the clumps into the galactic centre, while kinematically heated by the clumps. We additionally run a dark matter only simulation excluding baryonic component and confirm that the resultant shallower density profile is not due to numerical artifact in the simulation, such as two-body relaxation. △ Less

Submitted 18 August, 2011; v1 submitted 3 August, 2011; originally announced August 2011.

Comments: 5 pages, 2 figures, accepted for publication in MNRAS Letters

Nuclear Star-Forming Ring of the Milky Way: Simulations

Authors: Sungsoo S. Kim, Takayuki R. Saitoh, Myoungwon Jeon, Donald F. Figer, David Merritt, Keiichi Wada

Abstract: We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with… ▽ More We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past. △ Less

Submitted 2 June, 2011; originally announced June 2011.

Comments: Accepted for publication in the Astrophysical Journal Letters, vol. 734, 2011

Interplay between Stellar Spirals and the ISM in Galactic Disks

Authors: Keiichi Wada, Junichi Baba, Takayuki R. Saitoh

Abstract: We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic simulations in a `live' stellar disk. We focus especially on spiral structures excited in a isolated galactic disk without a stellar bar. Using high-resolution, 3-dimensional N-body/SPH simulations, we found that the spiral features of the gas in galactic disks are formed by essentially different mechanis… ▽ More We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic simulations in a `live' stellar disk. We focus especially on spiral structures excited in a isolated galactic disk without a stellar bar. Using high-resolution, 3-dimensional N-body/SPH simulations, we found that the spiral features of the gas in galactic disks are formed by essentially different mechanisms from the galactic shock in stellar density waves. The stellar spiral arms and the interstellar matter on average corotate in a galactic potential at any radii. Unlike the stream motions in the galactic shock, the interstellar matter flows into the local potential minima with irregular motions. The flows converge to form dense gas clouds/filaments near the bottom of the stellar spirals, whose global structures resemble dust-lanes seen in late-type spiral galaxies. The stellar arms are non-steady; they are wound and stretched by the galactic shear, and thus local densities of the arm change on a time scale of ~ 100 Myrs, due to bifurcating or merging with other arms. This makes the gas spirals associated with the stellar arms non-steady. The association of dense gas clouds are eventually dissolved into inter-arm regions with non-cirucular motions. Star clusters are formed from the cold, dense gases, whose ages are less than ~30 Myrs, and they are roughly associated with the background stellar arms without a clear spatial offset between gas spiral arms and distribution of young stars. △ Less

Submitted 7 April, 2011; originally announced April 2011.

Comments: 13 pages, 12 figures, accepted by ApJ. Higher resolution of ms.pdf is available at http://d.pr/Nvjk A targzipped Supplementary movies is available at http://d.pr/TV61

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

On the Interpretation of the l-v Features in the Milky Way Galaxy

Authors: Junichi Baba, Takayuki R. Saitoh, Keiichi Wada

Abstract: We model the gas dynamics of barred galaxies using a three-dimensional, high-resolution, $N$-body+hydrodynamical simulation and apply it to the Milky Way in an attempt to reproduce both the large-scale structure and the clumpy morphology observed in Galactic H\emissiontype{I} and CO $l-v$ diagrams. Owing to including the multi-phase interstellar medium, self-gravity, star-formation and supernovae… ▽ More We model the gas dynamics of barred galaxies using a three-dimensional, high-resolution, $N$-body+hydrodynamical simulation and apply it to the Milky Way in an attempt to reproduce both the large-scale structure and the clumpy morphology observed in Galactic H\emissiontype{I} and CO $l-v$ diagrams. Owing to including the multi-phase interstellar medium, self-gravity, star-formation and supernovae feedback, the clumpy morphology, as well as the large-scale features, in observed $l-v$ diagrams are naturally reproduced. We identify in our $l-v$ diagrams with a number of not only large-scale peculiar features such as the '3-kpc arm', '135-km s$^{-1}$ arm' and 'Connecting arm' but also clumpy features such as `Bania clumps', and then link these features in a face-on view of our model. We give suggestions on the real structure of the Milky Way and on the fate of gas clumps in the central region. △ Less

Submitted 16 September, 2010; originally announced September 2010.

Comments: accepted to PASJ

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

A Natural Symmetrization for the Plummer Potential

Authors: Takayuki R. Saitoh, Junichiro Makino

Abstract: We propose a symmetrized form of the softened gravitational potential which is a natural extension of the Plummer potential. The gravitational potential at the position of particle i (x_i,y_i,z_i), induced by particle j at (x_j,y_j,z_j), is given by: phi_ij = -G m_j/|r_ij^2+e_i^2+e_j^2|^1/2, where G is the gravitational constant, m_j is the mass of particle j, r_ij = |(x_i-x_j)^2+(y_i-y_j)^2+(z_i-… ▽ More We propose a symmetrized form of the softened gravitational potential which is a natural extension of the Plummer potential. The gravitational potential at the position of particle i (x_i,y_i,z_i), induced by particle j at (x_j,y_j,z_j), is given by: phi_ij = -G m_j/|r_ij^2+e_i^2+e_j^2|^1/2, where G is the gravitational constant, m_j is the mass of particle j, r_ij = |(x_i-x_j)^2+(y_i-y_j)^2+(z_i-z_j)^2|^1/2 and e_i and e_j are the gravitational softening lengths of particles i and j, respectively. This form is formally an extension of the Newtonian potential to five dimensions. The derivative of this equation in the x,y, and z directions correspond to the gravitational accelerations in these directions and these are always symmetric between two particles. When one applies this potential to a group of particles with different softening lengths, as is the case with a tree code, an averaged gravitational softening length for the group can be used. We find that the most suitable averaged softening length for a group of particles is <e_j^2> = sum_j^N m_j e_j^2 / M, where M = sum_j^N m_j and N are the mass and number of all particles in the group, respectively. The leading error related to the softening length is O(sum_j r_j d(e_j^2)/r_ij^3), where r_j is the distance between particle j and the center of mass of the group and d(e_j^2) = e_j^2 - <e_j^2>. Using this averaged gravitational softening length with the tree method, one can use a single tree to evaluate the gravitational forces for a system of particles with a wide variety of gravitational softening lengths. Consequently, this will reduce the calculation cost of the gravitational force for such a system with different softenings without the need for complicated forms of softening. We present the result of simple numerical tests. We found that our modification of the Plummer potential works well. △ Less

Submitted 15 January, 2011; v1 submitted 11 May, 2010; originally announced May 2010.

Comments: 8 pages, 3 figures, submitted to NewAstronomy

FAST: A Fully Asynchronous Split Time-Integrator for Self-Gravitating Fluid

Authors: Takayuki R. Saitoh, Junichiro Makino

Abstract: We describe a new algorithm for the integration of self-gravitating fluid systems using SPH method. We split the Hamiltonian of a self-gravitating fluid system to the gravitational potential and others (kinetic and internal energies) and use different time-steps for their integrations. The time integration is done in the way similar to that used in the mixed variable or multiple stepsize symplec… ▽ More We describe a new algorithm for the integration of self-gravitating fluid systems using SPH method. We split the Hamiltonian of a self-gravitating fluid system to the gravitational potential and others (kinetic and internal energies) and use different time-steps for their integrations. The time integration is done in the way similar to that used in the mixed variable or multiple stepsize symplectic schemes. We performed three test calculations. One was the spherical collapse and the other was an explosion. We also performed a realistic test, in which the initial model was taken from a simulation of merging galaxies. In all test calculations, we found that the number of time-steps for gravitational interaction were reduced by nearly an order of magnitude when we adopted our integration method. In the case of the realistic test, in which the dark matter potential dominates the total system, the total calculation time was significantly reduced. Simulation results were almost the same with those of simulations with the ordinary individual time-step method. Our new method achieves good performance without sacrificing the accuracy of the time integration. △ Less

Submitted 18 January, 2010; v1 submitted 11 August, 2009; originally announced August 2009.

Comments: 14 pages, 8 figures, accepted for publication in PASJ

The origin of large peculiar motions of star-forming regions and spiral structures of our Galaxy

Authors: Junichi Baba, Yoshiharu Asaki, Junichiro Makino, Makoto Miyoshi, Takayuki R. Saitoh, Keiichi Wada

Abstract: Recent VLBI (Very Long Baseline Interferometer) observations determined the distances and proper motions of star-forming regions in spiral arms directly. They showed that star-forming regions and young stars have large peculiar motions, as large as 30 km/s with complex structures. Such a large peculiar motion is incompatible with the prediction of the standard theory of quasi-stationary spiral a… ▽ More Recent VLBI (Very Long Baseline Interferometer) observations determined the distances and proper motions of star-forming regions in spiral arms directly. They showed that star-forming regions and young stars have large peculiar motions, as large as 30 km/s with complex structures. Such a large peculiar motion is incompatible with the prediction of the standard theory of quasi-stationary spiral arms. We use a high-resolution, self-consistent N-body+hydrodynamical simulation to explore how the spiral arms are formed and maintained, and how star-forming regions move. We found that arms are not quasi-stationary but transient and recurrent, as suggested in alternative theories of spiral structures. Because of this transient nature of the spiral arms, star-forming regions exhibit a trend of large and complex non-circular motions, which is qualitatively consistent with the VLBI observations. Owing to this large non-circular motion, a kinematically estimated gas map of our Galaxy has a large systematic errors of ~2-3 kpc in the distance from the Sun. △ Less

Submitted 29 October, 2009; v1 submitted 28 April, 2009; originally announced April 2009.

Comments: 11 pages, 11 figures. Accepted to ApJ

Journal ref: Astrophys.J.706:471-481,2009

Mass Supply to Galactic Center due to Nested Bars in the Galaxy

Authors: Daisuke Namekata, Asao Habe, Hidenori Matsui, Takayuki R. Saitoh

Abstract: We investigate rapid mass supply process by nested bars in the Galaxy by numerical simulation. We simulate gas flow in the whole galaxy disk with nested bars, which are the outer bar and the inner bar, especially with highly spatial resolution in the galactic central region. We assume two cases of inner bar size which are a smaller one and a larger one than the radius of the 200 pc gas ring whic… ▽ More We investigate rapid mass supply process by nested bars in the Galaxy by numerical simulation. We simulate gas flow in the whole galaxy disk with nested bars, which are the outer bar and the inner bar, especially with highly spatial resolution in the galactic central region. We assume two cases of inner bar size which are a smaller one and a larger one than the radius of the 200 pc gas ring which is corresponds to the Central Molecular Zone. From our numerical results, in the large size bar cases, the inner bars with large elongation induce sufficient mass inflow and destroy the 200 pc gas ring. On the other hand, in the small size bar cases, the inner bars with large elongation induce large mass inflow and do not destroy the 200 pc gas ring. This mass inflow is caused by straight shocks excited by the inner bar. In this case, nuclear gas disks of ~ 15 pc radius are formed. The nuclear gas disks are self-gravitationally unstable and we expect formation of compact star clusters under strong tidal force in the nuclear gas disks. We discuss evolution of the nuclear gas disk. △ Less

Submitted 7 October, 2008; originally announced October 2008.

Comments: 36 pages, 11 figures. Submitted to ApJ

Journal ref: Astrophys.J.691:1525-1539,2009

A Necessary Condition for Individual Time-Steps in SPH Simulations

Authors: Takayuki R. Saitoh, Junichiro Makino

Abstract: We show that the smoothed particle hydrodynamics (SPH) method, used with individual time-steps in the way described in the literature, cannot handle strong explosion problems correctly. In the individual time-step scheme, particles determine their time-steps essentially from a local Courant condition. Thus they cannot respond to a strong shock, if the pre-shock timescale is too long compared to… ▽ More We show that the smoothed particle hydrodynamics (SPH) method, used with individual time-steps in the way described in the literature, cannot handle strong explosion problems correctly. In the individual time-step scheme, particles determine their time-steps essentially from a local Courant condition. Thus they cannot respond to a strong shock, if the pre-shock timescale is too long compared to the shock timescale. This problem is not severe in SPH simulations of galaxy formation with a temperature cutoff in the cooling function at $10^4 {\rm K}$, while it is very dangerous for simulations in which the multiphase nature of the interstellar medium under $10^4 {\rm K}$ is taken into account. A solution for this problem is to introduce a time-step limiter which reduces the time-step of a particle if it is too long compared to the time-steps of its neighbor particles. Thus this kind of time-step constraint is essential for the correct treatment of explosions in high-resolution SPH simulations with individual time-steps. △ Less

Submitted 9 April, 2009; v1 submitted 6 August, 2008; originally announced August 2008.

Comments: 5 pages, 4 figures, ApJL accepted

Journal ref: Astrophys.J.697:L99-L102,2009