Researcher profile

Takayuki R. Saitoh

Takayuki R. Saitoh contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

astro-ph.GA astro-ph.CO astro-ph.IM astro-ph.EP astro-ph astro-ph.SR

Trust snapshot

Quick read

Trust 21 - EmergingVerification L1Unclaimed author

24works

0followers

6topics

4close collaborators

Actions

Decide how to stay connected

Follow researcher0

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

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

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

Building this graph slice

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

preprint2022arXiv

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

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.

preprint2016arXiv

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

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.

preprint2016arXiv

Eventful Evolution of Giant Molecular Clouds in Dynamically Evolving Spiral Arms

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.

preprint2016arXiv

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

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.

preprint2015arXiv

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

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.

preprint2015arXiv

Smoothed Particle Hydrodynamics with Smoothed Pseudo-Density

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.

preprint2014arXiv

Properties of thick discs formed in clumpy galaxies

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.

preprint2013arXiv

A Density Independent Formulation of Smoothed Particle Hydrodynamics

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.

preprint2013arXiv

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

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.

preprint2012arXiv

Dynamics of Non-Steady Spiral Arms in Disk Galaxies

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.

preprint2012arXiv

The formation of young dense star clusters through mergers

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.

preprint2011arXiv

A Natural Symmetrization for the Plummer Potential

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.

preprint2011arXiv

Interplay between Stellar Spirals and the ISM in Galactic Disks

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.

preprint2011arXiv

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

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.

preprint2011arXiv

Nuclear Star-Forming Ring of the Milky Way: Simulations

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.

preprint2011arXiv

Origin of multiple nuclei in ultraluminous infrared galaxies

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.

preprint2011arXiv

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

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.

preprint2011arXiv

Shock-induced star cluster formation in colliding galaxies

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.

preprint2010arXiv

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

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.

preprint2010arXiv

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

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.

preprint2009arXiv

A Necessary Condition for Individual Time-Steps in SPH Simulations

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.

preprint2009arXiv

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

We investigated the evolution of interacting disk galaxies using high-resolution $N$-body/SPH simulations, taking into account the multiphase nature of the interstellar medium (ISM). In our high-resolution simulations, a large-scale starburst occurred naturally at the collision interface between two gas disks at the first encounter, resulting in the formation of star clusters. This is consistent with observations of interacting galaxies. The probability distribution function (PDF) of gas density showed clear change during the galaxy-galaxy encounter. The compression of gas at the collision interface between the gas disks first appears as an excess at $n_{\rm H} \sim 10{\rm cm^{-3}}$ in the PDF, and then the excess moves to higher densities ($n_{\rm H} \gtrsim 100{\rm cm^{-3}}$) in a few times $10^7$ years where starburst takes place. After the starburst, the PDF goes back to the quasi-steady state. These results give a simple picture of starburst phenomena in galaxy-galaxy encounters.

preprint2008arXiv

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

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.

preprint2008arXiv

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?

We performed 3-dimensional N-body/SPH simulations to study how mass resolution and other model parameters such as the star formation efficiency parameter, C* and the threshold density, nth affect structures of the galactic gaseous/stellar disk in a static galactic potential. We employ 10^6 - 10^7 particles to resolve a cold and dense (T < 100 K & n_H > 100 cm^{-3}) phase. We found that structures of the ISM and the distribution of young stars are sensitive to the assumed nth. High-nth models with nth = 100 cm^{-3} yield clumpy multi-phase features in the ISM. Young stars are distributed in a thin disk of which half-mass scale height is 10 - 30 pc. In low-nth models with nth = 0.1 cm^{-3}, the stellar disk is found to be several times thicker, and the gas disk appears smoother than the high-nth models. A high-resolution simulation with high-nth is necessary to reproduce the complex structure of the gas disk. The global properties of the model galaxies in low-nth models, such as star formation histories, are similar to those in the high-nth models when we tune the value of C* so that they reproduce the observed relation between surface gas density and surface star formation rate density. We however emphasize that high-nth models automatically reproduce the relation, regardless of the values of C*. The ISM structure, phase distribution, and distributions of young star forming region are quite similar between two runs with values of C* which differ by a factor of 15. We also found that the timescale of the flow from n_H ~1 cm^{-3} to n_H > 100 cm^{-3} is about 5 times as long as the local dynamical time and is independent of the value of C*. The use of a high-nth criterion for star formation in high-resolution simulations makes numerical models fairy insensitive to the modelling of star formation. (Abridged)

Takayuki R. Saitoh

Quick read

Decide how to stay connected

How to connect with this researcher

Open a focused conversation when the fit is right

See the researcher in context

Building this graph slice

24 published item(s)

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

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

Eventful Evolution of Giant Molecular Clouds in Dynamically Evolving Spiral Arms

The Giant Impact Simulations with Density Independent Smoothed Particle Hydrodynamics

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

Smoothed Particle Hydrodynamics with Smoothed Pseudo-Density

Properties of thick discs formed in clumpy galaxies

A Density Independent Formulation of Smoothed Particle Hydrodynamics

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

Dynamics of Non-Steady Spiral Arms in Disk Galaxies

The formation of young dense star clusters through mergers

A Natural Symmetrization for the Plummer Potential

Interplay between Stellar Spirals and the ISM in Galactic Disks

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

Nuclear Star-Forming Ring of the Milky Way: Simulations

Origin of multiple nuclei in ultraluminous infrared galaxies

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

Shock-induced star cluster formation in colliding galaxies

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

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

A Necessary Condition for Individual Time-Steps in SPH Simulations

Toward First-Principle Simulations of Galaxy Formation: II. Shock-Induced Starburst at a Collision Interface During the First Encounter of Interacting Galaxies

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

Toward first-principle simulations of galaxy formation: I. How should we choose star formation criteria in high-resolution simulations of disk galaxies?