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.Comment: 14 pages, 8 figures, accepted for publication in PAS

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.Comment: 15 pages, 13 figures, 3 tables, accepted for publication in MNRA

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.Comment: accepted to PAS

Effects of a Supermassive Black Hole Binary on a Nuclear Gas Disk

We study influence of a galactic central supermassive black hole (SMBH) binary on gas dynamics and star formation activity in a nuclear gas disk by making three-dimensional Tree+SPH simulations. Due to orbital motions of SMBHs, there are various resonances between gas motion and the SMBH binary motion. We have shown that these resonances create some characteristic structures of gas in the nuclear gas disk, for examples, gas elongated or filament structures, formation of gaseous spiral arms, and small gas disks around SMBHs. In these gaseous dense regions, active star formations are induced. As the result, many star burst regions are formed in the nuclear region.Comment: 19 pages, 11 figures, accepted for publication in Ap