Gaia Data Release 2 revealed detailed structures of nearby stars in phase space. These include the Hercules stream, whose origin is still debated. Earlier numerical studies conjectured that the observed structures originate from orbits in resonance with the bar, based on static potential models for the Milky Way. We, in contrast, approach the problem via a self-consistent, dynamic, and morphologically well-resolved model, namely a full $N$-body simulation of the Milky Way. Our simulation comprises about 5.1 billion particles in the galactic stellar bulge, bar, disk, and dark-matter halo and is evolved to 10 Gyr. Our model's disk component is composed of 200 million particles, and its simulation snapshots are stored every 10 Myr, enabling us to resolve and classify resonant orbits of representative samples of stars. After identifying the Sun's position in the simulation, we compare the distribution of stars in its neighborhood with Gaia's astrometric data, thereby establishing the role of identified resonantly trapped stars in the formation of Hercules-like structures. From our orbital spectral-analysis we identify multiple resonances and conclude that the Hercules stream is dominated by the 4:1 and 5:1 outer Lindblad and corotation resonances. In total, this yields a trimodal structure of the Hercules stream. From the relation between resonances and ridges in phase space, we conclude that the pattern speed of the Milky-Way bar is 40$-$45 km s$^{-1}$ kpc$^{-1}$.

From the decades of the theoretical studies, it is well known that the formation of the bar triggers the gas funnelling into the central sub-kpc region and leads to the formation of a kinematically cold nuclear stellar disc (NSD). We demonstrate that this mechanism can be used to identify the formation epoch of the Galactic bar, using an N-body/hydrodynamics simulation of an isolated Milky Way-like galaxy. As shown in many previous literature, our simulation shows that the bar formation triggers an intense star formation for ~1 Gyr in the central region, and forms a NSD. As a result, the oldest age limit of the NSD is relatively sharp, and the oldest population becomes similar to the age of the bar. Therefore, the age distribution of the NSD tells us the formation epoch of the bar. We discuss that a major challenge in measuring the age distribution of the NSD in the Milky Way is contamination from other non-negligible stellar components in the central region, such as a classical bulge component. We demonstrate that because the NSD is kinematically colder than the other stellar populations in the Galactic central region, the NSD population can be kinematically distinguished from the other stellar populations, if the 3D velocity of tracer stars are accurately measured. Hence, in addition to the line-of-sight velocities from spectroscopic surveys, the accurate measurements of the transverse velocities of stars are necessary, and hence the near-infrared space astrometry mission, JASMINE, would play a crutial role to identify the formation epoch of the Galactic bar. We also discuss that the accuracy of stellar age estimation is also crucial to measure the oldest limit of the NSD stellar population.

Using the cross-matched data of Gaia DR2 and 2MASS Point Source Catalog, we investigated the surface density distribution of stars aged ~1 Gyr in the thin disk in the range of 90{\deg} <= l <= 270{\deg}. We selected 4,654 stars above the turnoff corresponding to the age ~1 Gyr, that fall within a small box region in the color-magnitude diagram, (J-Ks)0 versus M(Ks), for which the distance and reddening are corrected. The selected sample shows an arm-like overdensity at 90{\deg} <= l <= 190{\deg}.This overdensity is located close to the Local arm traced by high-mass star forming regions (HMSFRs), but its pitch angle is slightly larger than that of the HMSFR-defined arm. Although the significance of the overdensity we report is marginal, its structure poses questions concerning both of the competing scenarios of spiral arms, the density-wave theory and the dynamic spiral arm model. The offset between the arms traced by stars and HMSFRs, i.e., gas, is difficult to be explained by the dynamic arm scenario. On the other hand, the pitch angle of the stellar Local arm, if confirmed, larger than that of the Perseus arm is difficult to be explained with the classical density-wave scenario. The dynamic arm scenario can explain it if the Local arm is in a growing up phase, while the Perseus arm is in a disrupting phase. Our result provides a new and complex picture of the Galactic spiral arms, and encourages further studies.

Growing interests in the chemical feature of r-process elements among nearby disk stars represented by the [Eu/Fe] vs. [Fe/H] diagram have sprouted since it can assess the origin of r-process elements through the comparison with theoretical models, including a test as to if neutron star mergers can be the major site of r-process nucleosynthesis. On the other hand, recent studies reveal that local chemistry is strongly coupled with the dynamics of Galactic disk, which predicts that stars radially move on the disk where the observed elemental feature is different at various Galactocentric distances. Here, we show that radial migration of stars across the Galactic disk plays a crucial role in shaping the r-process abundance feature in the solar vicinity. In this proposed scenario, we highlight the importance of migration from the outer disk where [r-process/Fe] of some old stars is predicted to be enhanced to the level beyond the expectation from the observed Galactic Fe and Eu radial gradient, which results in a large span of [r-process/Fe] among nearby disk stars. The variation in the [r-process/Fe] ratio seen across the Galactic disk as well as in dwarf galaxies may be an outcome of different stellar initial mass functions which change the occurrence frequency between supernovae leaving behind neutron stars and ones ending with black holes. Here we propose that enhancement in [Eu/Fe] is attributed to the initial mass function lacking high-mass stars such as > 25 solar masses in the scheme for which neutron star mergers are a major source of r-process elements.

The origin of the Hercules stream, the most prominent velocity substructure in the Solar neighbour disc stars, is still under debate. Recent accurate measurements of position, velocity, and metallicity provided by Tycho Gaia Astrometric Solution (TGAS) and RAdial Velocity Experiments (RAVE) have revealed that the Hercules stream is most clearly seen in the metal-rich region ([Fe/H] > 0), while it is not clearly seen in lower metallicity region ([Fe/H] < -0.25). By using a large number of chemo-dynamical 2D test-particle simulations with a rotating bar and/or spiral arms, we find that the observed [Fe/H] dependence of the Hercules stream is a natural consequence of the inside-out formation of the stellar disc and the existence of highly non-closed orbits in the rotating frame of the bar or spiral arms. Our successful models that reproduce the observed properties of the Hercules stream include not only fast-bar-only and fast-bar+spiral models, but also slow-bar+spiral models. This indicates that it is very difficult to estimate the pattern speed of the bar or spiral arms based only on the observations of the Hercules stream in the Solar neighbourhood. As a by-product of our simulations, we make some predictions about the locations across the Galactic plane where we can observe velocity bimodality that is not associated with the Hercules stream. These predictions can be tested by the Gaia Data Release 2, and such a test will improve our understanding of the evolution of the Milky Way stellar disc.

We construct a model for the Milky Way Galaxy composed of a stellar disc and bulge embedded in a dark-matter halo. All components are modelled as $N$-body systems with up to 8 billion equal-mass particles and integrated up to an age of 10\,Gyr. We find that net angular-momentum of the dark-matter halo with a spin parameter of $\lambda=0.06$ is required to form a relatively short bar ($\sim 4$\,kpc) with a high pattern speed (40--50\,km\,s$^{-1}$). By comparing our model with observations of the Milky Way Galaxy, we conclude that a disc mass of $\sim 3.7\times10^{10}M_{\odot}$ and an initial bulge scale length and velocity of $\sim 1$\,kpc and $\sim 300$\,km\,s$^{-1}$, respectively, fit best to the observations. The disc-to-total mass fraction ($f_{\rm d}$) appears to be an important parameter for the evolution of the Galaxy and models with $f_{\rm d}\sim 0.45$ are most similar to the Milky Way Galaxy. In addition, we compare the velocity distribution in the solar neighbourhood in our simulations with observations in the Milky Way Galaxy. In our simulations the observed gap in the velocity distribution, which is expected to be caused by the outer Lindblad resonance (the so-called Hercules stream), appears to be a time-dependent structure. The velocity distribution changes on a time scale of 20--30\,Myr and therefore it is difficult to estimate the pattern speed of the bar from the shape of the local velocity distribution alone.

We apply a simple axisymmetric disc model to 218 Galactic Cepheids whose accurate measurements of the distance and velocities are obtained by cross-matching an existing Cepheids catalogue with the Gaia DR2 data. Our model fit determines the "local centrifugal speed", $V_\mathrm{c}$ $-$ defined as the rotation speed required to balance the local radial gravitational force $-$ at the Sun's location to be $V_{c}(R_0)=236\pm 3$ km s$^{-1}$ and the Sun's azimuthal and radial peculiar motions to be $V_{\odot}=12.4\pm0.7$ km s$^{-1}$ and $U_{\odot}=7.7\pm0.9$ km s$^{-1}$, respectively. These results are obtained with strong priors on the solar radius, $R_0=8.2\pm0.1$ kpc, and Sun's angular rotation velocity, $\Omega_{\odot}=30.24\pm0.12$ km s$^{-1}$ kpc$^{-1}$. We also applied the axisymmetric model to mock data from an N-body/hydrodynamic simulation of a Milky Way-like galaxy with a bar and spiral arms. We find that our axisymmetric model fit to the young stars recovers the local centrifugal speed reasonably well, even in the face of significant non-axisymmetry. However, the local centrifugal speed determined from our Cepheid sample could suffer from systematic uncertainty as large as 6 km s$^{-1}$.

We present results of the $^{12}$CO (1--0) mosaic observations of the nearby barred-spiral galaxy M83 obtained with the Atacama Large Millimeter/submillimeter Array (ALMA). The total flux is recovered by combining the ALMA data with single-dish data obtained using the Nobeyama 45-m telescope. The combined map covers a $\sim$13 kpc$^{2}$ field that includes the galactic center, eastern bar, and spiral arm with a resolution of \timeform{2''.03} $\times$ \timeform{1''.1} ($\sim$45 pc $\times$ $\sim$25 pc). With a resolution comparable to typical sizes of giant molecular clouds (GMCs), the CO distribution in the bar and arm is resolved into many clumpy peaks that form ridge-like structures. Remarkably, in the eastern arm, the CO peaks form two arc-shaped ridges that run along the arm and exhibit a distinct difference in the activity of star formation: the one on the leading side has numerous HII regions associated with it, whereas the other one on the trailing side has only a few. To see whether GMCs form stars with uniform star formation efficiency (SFE) per free-fall time (SFEff), GMCs are identified from the data cube and then cross-matched with the catalog of HII regions to estimate the star formation rate for each of them. 179 GMCs with a median mass of 1.6 $\times$ 10$^{6}$ $M_{\odot}$ are identified. The mass-weighted average SFEff of the GMCs is $\sim$9.4 $\times$ 10$^{-3}$, which is in agreement with models of turbulence regulated star formation. Meanwhile, we find that SFEff is not universal within the mapped region. In particular, one of the arm ridges shows a high SFEff with a mass-weighted value of $\sim$2.7 $\times$ 10$^{-2}$, which is higher by more than a factor of 5 compared to the inter-arm regions. This large regional variation in SFEff favors the recent interpretation that GMCs do not form stars at a constant rate within their lifetime.

By taking advantage of the superb measurements of position and velocity for an unprecedented large number of stars provided in Gaia DR2, we have generated the first maps of the rotation velocity, $V_{\rm rot}$, and vertical velocity, $V_{\rm z}$, distributions as a function of the Galactocentric radius, $R_{\rm gal}$, across a radial range of $5<R_{\rm gal}<12$~kpc. In the $R-V_{\rm rot}$ map, we have identified many diagonal ridge features, which are compared with the location of the spiral arms and the expected outer Lindblad resonance of the Galactic bar. We have detected also radial wave-like oscillations of the peak of the vertical velocity distribution.

Recent developments in computer hardware and software enable researchers to simulate the self-gravitating evolution of galaxies at a resolution comparable to the actual number of stars. Here we present the results of a series of such simulations. We performed $N$-body simulations of disk galaxies with between 100 and 500 million particles over a wide range of initial conditions. Our calculations include a live bulge, disk, and dark matter halo, each of which is represented by self-gravitating particles in the $N$-body code. The simulations are performed using the gravitational $N$-body tree-code Bonsai running on the Piz Daint supercomputer. We find that the time scale over which the bar forms increases exponentially with decreasing disk-mass fraction and that the bar formation epoch exceeds a Hubble time when the disk-mass fraction is $\sim0.35$. These results can be explained with the swing-amplification theory. The condition for the formation of $m=2$ spirals is consistent with that for the formation of the bar, which is also an $m=2$ phenomenon. We further argue that the non-barred grand-design spiral galaxies are transitional, and that they evolve to barred galaxies on a dynamical timescale. We also confirm that the disk-mass fraction and shear rate are important parameters for the morphology of disk galaxies. The former affects the number of spiral arms and the bar formation epoch, and the latter determines the pitch angle of the spiral arms.

We carried out $^{12}$CO(J = 1 -0) observations toward bar and arm regions of strongly barred galaxies, NGC 1300 and NGC 5383, with the Nobeyama 45-m radio telescope (beamsize of 1-2 kpc in the galaxies). The aim of the observations is to qualitatively examine a new scenario for the suppression of star formation in bars based on recent high-resolution numerical simulations: higher speed collisions between molecular clouds in the bar region compared with the arm region suppress the massive star formation. CO emissions were detected from all the regions, indicating the presence of the molecular gases in the strong bars without associating clear HII regions. In both galaxies, the velocity width of the CO line profile tends to be larger in the bar region than in the arm region, which is qualitatively consistent with the new scenario.

We have discovered a clear sign of the disruption phase of the Perseus arm in the Milky Way using Cepheid variables, taking advantage of the accurately measured distances of Cepheids and the proper motions from Gaia Data Release 1. Both the Galactocentric radial and rotation velocities of 77 Cepheids within 1.5 kpc of the Perseus arm are correlated with their distances from the locus of the Perseus arm, as the trailing side is rotating faster and moving inward compared to the leading side. We also found a negative vertex deviation for the Cepheids on the trailing side, $-27.6\pm2.4$ deg, in contrast to the positive vertex deviation in the solar neighborhood. This is, to our knowledge, the first direct evidence that the vertex deviation around the Perseus arm is affected by the spiral arm. We compared these observational trends with our $N$-body/hydrodynamics simulations based on a static density-wave spiral scenario and those based on a transient dynamic spiral scenario. Although our comparisons are limited to qualitative trends, they strongly favor a conclusion that the Perseus arm is in the disruption phase of a transient arm.

We have obtained 12CO(1--0) data of the nearby barred spiral galaxy M83 from Atacama Large Millimeter/submillimeter Array and Nobeyama 45m observations. By combining these two data sets, the total CO flux has been recovered, and a high angular resolution (2" corresponding to ~40 pc at the distance of M83) has been achieved. The field of view is 3' corresponding to ~3.4 kpc and covers the galactic center, bar, and spiral arm regions. In order to investigate how these galactic structures affect gas properties, we have created a probability distribution function (PDF) of the CO integrated intensity (I_CO), peak temperature, and velocity dispersion for a region with each structure. We find that the I_CO PDF for the bar shows a bright-end tail while that for the arm does not. Since the star formation efficiency is lower in the bar, this difference in PDF shape is contrary to the trend in Milky Way studies where the bright-end tail is found for star-forming molecular clouds. While the peak temperature PDFs are similar for bar and arm regions, velocity dispersion in bar is systematically larger than in arm. This large velocity dispersion is likely a major cause of the bright-end tail and of suppressed star formation. We also investigate an effect of stellar feedback to PDF profiles and find that the different I_CO PDFs between bar and arm regions cannot be explained by the feedback effect, at least at the current spatial scale.

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.

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.

Theoretical studies on the response of interstellar gas to a gravitational potential disc with a quasi-stationary spiral arm pattern suggest that the gas experiences a sudden compression due to standing shock waves at spiral arms. This mechanism, called a galactic shock wave, predicts that gas spiral arms move from downstream to upstream of stellar arms with increasing radius inside a corotation radius. In order to investigate if this mechanism is at work in the grand-design spiral galaxy M51, we have measured azimuthal offsets between the peaks of stellar mass and gas mass distributions in its two spiral arms. The stellar mass distribution is created by the spatially resolved spectral energy distribution fitting to optical and near infrared images, while the gas mass distribution is obtained by high-resolution CO and HI data. For the inner region (r < 150"), we find that one arm is consistent with the galactic shock while the other is not. For the outer region, results are less certain due to the narrower range of offset values, the weakness of stellar arms, and the smaller number of successful offset measurements. The results suggest that the nature of two inner spiral arms are different, which is likely due to an interaction with the companion galaxy.

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.

There are two theories of stellar spiral arms in isolated disc galaxies that model stellar spiral arms with different longevities: quasi-stationary density wave theory, which characterises spirals as rigidly rotating, long-lived patterns (i.e. steady spirals), and dynamic spiral theory, which characterises spirals as differentially rotating, transient, recurrent patterns (i.e. dynamic spirals). In order to discriminate between these two spiral models observationally, we investigated the differences between the gas velocity patterns predicted by these two spiral models in hydrodynamic simulations. We found that the azimuthal phases of the velocity patterns relative to the gas density peaks (i.e. gaseous arms) differ between the two models, as do the gas flows; nevertheless, the velocity patterns themselves are similar for both models. Such similarity suggests that the mere existence of streaming motions does not conclusively confirm the steady spiral model. However, we found that the steady spiral model shows that the gaseous arms have radial streaming motions well inside the co-rotation radius, whereas the dynamic spiral model predicts that the gaseous arms tend to have tangential streaming motions. These differences suggest that the gas velocity patterns around spiral arms will enable distinction between the spiral theories.

We investigate the short-term dynamical evolution of stellar grand-design spiral arms in barred spiral galaxies using a three-dimensional (3D) $N$-body/hydrodynamic simulation. Similar to previous numerical simulations of unbarred, multiple-arm spirals, we find that grand-design spiral arms in barred galaxies are not stationary, but rather dynamic. This means that the amplitudes, pitch angles, and rotational frequencies of the spiral arms are not constant, but change within a few hundred million years (i.e. the typical rotational period of a galaxy). We also find that the clear grand-design spirals in barred galaxies appear it only when the spirals connect with the ends of the bar. Furthermore, we find that the short-term behaviour of spiral arms in the outer regions ($R>$ 1.5--2 bar radius) can be explained by the swing amplification theory and that the effects of the bar are not negligible in the inner regions ($R<$ 1.5--2 bar radius). These results suggest that, although grand-design spiral arms in barred galaxies are affected by the stellar bar, the grand-design spiral arms essentially originate not as bar-driven stationary density waves, but rather as self-excited dynamic patterns. We imply that a rigidly rotating grand-design spiral could not be a reasonable dynamical model for investigating gas flows and cloud formation even in barred spiral galaxies.

We investigate the redshift evolution of the molecular gas mass fraction (f_mol=M_mol/(M_star+M_mol), where M_mol is molecular gas mass and M_star is stellar mass) of galaxies in the redshift range of 0<z<2 as a function of the stellar mass by combining CO literature data. We observe a stellar-mass dependence of the f_mol evolution where massive galaxies have largely depleted their molecular gas at z=1, whereas the f_mol value of less massive galaxies drastically decreases from z=1. We compare the observed M_star-f_mol relation with theoretical predictions from cosmological hydrodynamic simulations and semi-analytical models for galaxy formation. Although the theoretical studies approximately reproduce the observed mass dependence of f_mol evolution, they tend to underestimate the f_mol values, particularly of less massive (<10^10 Msun) and massive galaxies (>10^11 Msun) when compared with the observational values. Our result suggests the importance of the feedback models which suppress the star formation while simultaneously preserving the molecular gas in order to reproduce the observed M_star-f_mol relation.

Theories of stellar spiral arms in disk galaxies can be grouped into two classes based on the longevity of a spiral arm. Although the quasi-stationary density wave theory supposes that spirals are rigidly-rotating, long-lived patterns, the dynamic spiral theory predicts that spirals are differentially-rotating, transient, recurrent patterns. In order to distinguish between the two spiral models from observations, we performed hydrodynamic simulations with steady and dynamic spiral models. Hydrodynamics simulations in steady spiral models demonstrated that the dust lane locations relative to the stellar spiral arms (hereafter, arm-gas offsets) depend on radius, regardless of the strength and pitch angle of the spiral and the model of the inter-stellar medium (ISM). In contrast, we found that the dynamic spiral models show no systematic radial dependence of the arm-gas offsets. The arm-gas offset radial profile method, together with the other test methods, will help us to distinguish between the two spiral models in observed spiral galaxies.

We performed 12CO(J=1-0) (hereafter, CO) observations towards 12 normal star-forming galaxies with stellar mass of Mstar=10^10.6-10^11.3 Msun at z=0.1-0.2 with the 45-m telescope at the Nobeyama Radio Observatory (NRO). The samples are selected with Dn(4000) that is a strength of the 4000 \AA break, instead of commonly used far-infrared (FIR) flux. We successfully detect the CO emissions from eight galaxies with signal-to-noise ratio (S/N) larger than three, demonstrating the effectiveness of the Dn(4000)-based sample selection. For the first time, we find a tight anti-correlation between Dn(4000) and molecular gas fraction (fmol) using literature data of nearby galaxies in which the galaxies with more fuel for star formation have younger stellar populations. We find that our CO-detected galaxies at z~0.1-0.2 also follow the same relation of nearby galaxies. This implies that the galaxies evolve along this Dn(4000)-fmol relation, and that Dn(4000) seems to be used as a proxy for fmol which requires many time-consuming observations. Based on the comparison with the model calculation with a population synthesis code, we find that star formation from metal enriched gas and its quenching in the early time are necessary to reproduce galaxies with large Dn(4000) and non-zero gas fraction.

Classical Cepheids are useful tracers of the Galactic young stellar population because their distances and ages can be determined from their period-luminosity and period-age relations. In addition, the radial velocities and chemical abundance of the Cepheids can be derived from spectroscopic observations, providing further insights into the structure and evolution of the Galaxy. Here, we report the radial velocities of classical Cepheids near the Galactic Center, three of which were reported in 2011, the other reported for the first time. The velocities of these Cepheids suggest that the stars orbit within the Nuclear Stellar Disk, a group of stars and interstellar matter occupying a region of 200 pc around the Center, although the three-dimensional velocities cannot be determined until the proper motions are known. According to our simulation, these four Cepheids formed within the Nuclear Stellar Disk like younger stars and stellar clusters therein.

The majority of astrophysics involves the study of spiral galaxies, and stars and planets within them, but how spiral arms in galaxies form and evolve is still a fundamental problem. Major progress in this field was made primarily in the 1960s, and early 1970s, but since then there has been no comprehensive update on the state of the field. In this review, we discuss the progress in theory, and in particular numerical calculations, which unlike in the 1960s and 1970s, are now commonplace, as well as recent observational developments. We set out the current status for different scenarios for spiral arm formation, the nature of the spiral arms they induce, and the consequences for gas dynamics and star formation in different types of spiral galaxies. We argue that, with possible the exception of barred galaxies, spiral arms are transient, recurrent and initiated by swing amplified instabilities in the disc. We suppose that unbarred m = 2 spiral patterns are induced by tidal interactions, and slowly wind up over time. However the mechanism for generating spiral structure does not appear to have significant consequences for star formation in galaxies.

We present the results of the wide-field $^{12}$CO (1--0) observations of the nearby barred galaxy M83 carried out with the Nobeyama Millimeter Array (NMA). The interferometric data are combined with the data obtained with the Nobeyama 45-m telescope to recover the total-flux. The target fields of the observations cover the molecular bar and part of the spiral arms, with a spatial resolution of ~110 pc x 260 pc. By exploiting the resolution and sensitivity to extended CO emission, the impact of the galactic structures on the molecular gas content is investigated in terms of the gas kinematics and the star formation. By inspecting the gas kinematics, the pattern speed of the bar is estimated to be 57.4 $\pm$ 2.8 km s$^{-1}$ kpc$^{-1}$, which places the corotation radius to be about 1.7 times the semi-major radius of the bar. Within the observed field, HII regions brighter than 10$^{37.6}$ erg s$^{-1}$ in H{\alpha} luminosity are found to be preferentially located downstream of the CO emitting regions. Azimuthal angular offsets between molecular gas and star forming (SF) calculated with the angular cross-correlation method confirm the trend. By comparing with a cloud orbit model based on the derived pattern speed, the angular offsets are found to be in accordance with a time delay of about 10 Myr. Finally, to test whether the arm/bar promote star formation efficiency (SFE $\equiv$ Star Formation Rate (SFR)/H$_2$ mass), SFR is derived with the diffuse-background-subtracted H{alpha} and 24{\mu}m images. The arm-to-interarm ratio of the SFE is found to lie in the range of 2 to 5, while it is ~1 if no background-removal is performed. The CO-SF offsets and the enhancement of the SFE in the arm/bar found in the inner region of M83 are in agreement with the predictions of the classical galactic shock model.

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.

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.

Most stars in galactic disks are believed to be born as a member of star clusters or associations. Star clusters formed in disks are disrupted due to the tidal stripping and the evolution of star clusters themselves, and as a results new stars are supplied to the galactic disks. We performed $N$-body simulations of star clusters in galactic disks, in which both star clusters and galactic disks are modeled as $N$-body ("live") systems, and as a consequence the disks form transient and recurrent spiral arms. In such non-steady spiral arms, star clusters migrate radially due to the interaction with spiral arms. We found that the migration timescale is a few hundreds Myr and that the angular momentum changes of star clusters are at most $\sim 50$% in 1 Gyr. Radial migration of star clusters to the inner region of galaxies results in a fast disruption of the star clusters because of a stronger tidal field in the inner region of the galaxy. This effect is not negligible for the disruption timescale of star clusters in galactic disks. Stars stripped from clusters form tidal tails which spread over 1--2 kpc. While the spatial distribution of tidal tails change in a complicated way due to the non-steady spiral arms, the velocity distribution conserve well even if the tidal tails are located at a few kpc from their parent clusters. Tidal tails of clusters in galactic disks might be detected using velocity plots.

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 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.

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.