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 propose and demonstrate a new stochastic star formation model for star-by-star simulations. In this model, an assumed stellar initial mass function (IMF) is randomly assigned to newly formed stars. We then 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 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. Herein, we demonstrate that our models can be applied to simulations varying from star clusters to galaxies for a wide range of resolutions.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 $\alpha$-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.

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.

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.

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

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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)

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.

Tidal disruption of dark matter halos around proto-globular clusters in a halo of a small galaxy is studied in the context of the hierarchical clustering scenario by using semi-cosmological N-body/SPH simulations assuming the standard cold dark matter model ($\Omega_0 = 1$). Our analysis on formation and evolution of the galaxy and its substructures archives until $z = 2.0$. In such a high-redshift universe, the Einstein-de Sitter universe is still a good approximation for a recently favored $\Lambda$-dominated universe, and then our results does not depend on the choice of cosmology. In order to resolve small gravitationally-bound clumps around galaxies and consider radiative cooling below $T = 10^4 K$, we adopt a fine mass resolution ($m_{\rm SPH} = 1.12 \times 10^3 \Msun$). Because of the cooling, each clump immediately forms a `core-halo' structure which consists of a baryonic core and a dark matter halo. The tidal force from the host galaxy mainly strips the dark matter halo from clumps and, as a result, theses clumps get dominated by baryons. Once a clump is captured by the host halo, its mass drastically decreases each pericenter passage. At $z = 2$, more than half of the clumps become baryon dominated systems (baryon mass/total mass $> 0.5$). Our results support the tidal evolution scenario of the formation of globular clusters and baryon dominated dwarf galaxies in the context of the cold dark matter universe.

Combining a theoretical model of mass accretion onto a galactic center with a high-resolution $N$-body/SPH simulation, we investigate the formation of an intermediate massive black hole (IMBH) during the hierarchical formation of a small spiral galaxy (with a total mass of $10^{10}M_{\odot}$) in the high-$z$ universe. We found that the rate of average mass accretion to the nucleus due to the radiation drag exerted by newly formed stars in the forming galaxy is $\approx 10^{-5}M_{\odot}$yr$^{-1}$. As a result of this accretion, an IMBH with $\approx 10^{4}M_{\odot}$ can be formed in the center of the spiral galaxy at $z\sim 4$. We found that a central BH coevolves with the dark matter halo from $z\sim 15$ to $z\sim 2$. The mass ratio of the BH to the dark matter halo is nearly constant $\approx (1-3) \times 10^{-6}$ from $z\sim 10$ to $z\sim 2$. This is because that change in the dark matter potential enhances star formation in the central part of the galaxy, and as a result the BH evolves due to mass accretion via the radiation drag. Therefore, our model naturally predicts a correlation between massive BHs and dark matter halos. Moreover, it is found that the final BH-to-bulge mass ratio ($\approx 5\times 10^{-5}$) in a small spiral galaxy at high-$z$ is much smaller than that in the large galaxies ($\approx 10^{-3}$). Our results also suggest that the scatter in the observed scaling relations between the bulge mass and black hole mass are caused by a time lag between BH growth and growth of bulge. We also predict that the X-ray luminosity of AGN is positively correlated with the CO luminosity in the central region. By comparing our results with the properties of Lyman break galaxies (LBGs), it is predicted that some LBGs have massive BHs of $\approx 10^{6}-10^{7}M_{\odot}$.

Using high-resolution N-body/SPH simulations with $2\times 10^6$ particles, we investigate the evolution of stellar and gaseous galactic cores during the hierarchical formation of a spiral galaxy. We find that the galactic core ($r < 300 $ pc) coevolves with the host galaxy. The average mass ratio between the baryonic core and the halo is nearly constant, $\sim $ 0.04 from $z \sim 10$ to $z \sim 2$. However, there are several `rapid-growing phases' during the evolution, in which the rate of mass accretion to the central sub-kpc region is ten times higher ($\sim 1 M_\odot$ yr pc$^{-1}$) than the average accretion rate. The rapid growth of the inner core is associated with the major merger events with a time-delay. We also find that the spin-axis of the gas core frequently changes. As a result, the angular momentum vector of the central part of the galaxy is independent of the rotation of the outer part. Our results suggest that if a constant fraction of the baryonic mass in the central several 100 pc of a galaxy is converted into a massive black hole, the black hole mass should correlate with the total mass of the galaxies.

We describe a new method to accelerate neighbor searches on GRAPE, i.e. a special purpose hardware that efficiently calculates gravitational forces and potentials in $N$-body simulations. In addition to the gravitational calculations, GRAPE simultaneously constructs the lists of neighbor particles that are necessary for Smoothed Particle Hydrodynamics (SPH). However, data transfer of the neighbor lists from GRAPE to the host computer is time consuming, and can be a bottleneck. In fact, the data transfer can take about the same time as the calculations of forces themselves. Making use of GRAPE's special treatment of neighbor lists, we can reduce the amount of data transfer if we search neighbors in the order that the neighbor lists, constructed in a single GRAPE run, overlap each other. We find that the Morton-ordering requires very low additional calculation and programming costs, and results in successful speed-up on data transfer. We show some benchmark results in the case of GRAPE-5. Typical reduction in transferred data becomes as much as 90%. This method is suitable not only for GRAPE-5, but also GRAPE-3 and the other versions of GRAPE.