We performed new large-scale $^{12}$CO, $^{13}$CO, and C$^{18}$O $J=$1--0 observations of the W43 giant molecular cloud complex in the tangential direction of the Scutum arm ($l\sim {30^\circ}$) as a part of the FUGIN project. The low-density gas traced by $^{12}$CO is distributed over 150 pc $\times$ 100 pc ($l \times b$), and has a large velocity dispersion (20-30 km s$^{-1}$). However, the dense gas traced by C$^{18}$O is localized in the W43 Main, G30.5, and W43 South (G29.96-0.02) high-mass star-forming regions in the W43 GMC complex, which have clumpy structures. We found at least two clouds with a velocity difference of $\sim$ 10-20 km s$^{-1}$, both of which are likely to be physically associated with these high-mass star-forming regions based on the results of high $^{13}$CO $J=$ 3-2 to $J =$ 1-0 intensity ratio and morphological correspondence with the infrared dust emission. The velocity separation of these clouds in W43 Main, G30.5, and W43 South is too large for each cloud to be gravitationally bound. We also revealed that the dense gas in the W43 GMC has a high local column density, while "the current SFE" of entire the GMC is low ($\sim 4\%$) compared with the W51 and M17 GMC. We argue that the supersonic cloud-cloud collision hypothesis can explain the origin of the local mini-starbursts and dense gas formation in the W43 GMC complex.

In many barred galaxies, star formation efficiency (SFE) in the bar is lower than those in the arm and bar-end, and its cause has still not been clear. Focusing on the strongly barred galaxy NGC 1300, we investigate the possibility that the presence of a large amount of diffuse molecular gas, which would not contribute to the SF, makes the SFE low in appearance. We examine the relation between the SFE and the diffuse molecular gas fraction ($f_{\rm dif}$), which is derived using the $^{12}$CO($1-0$) flux obtained from the interferometer of ALMA 12-m array, which has no sensitivity on diffuse (extended; FWHM $\gtrapprox 700$ pc) molecular gases due to the lack of ACA, and the total $^{12}$CO($1-0$) flux obtained from Nobeyama 45-m single-dish telescope. We find that the SFE decreases with increasing $f_{\rm dif}$. The $f_{\rm dif}$ and ${\rm SFE}$ are $0.74 - 0.91$ and $(0.06 - 0.16) ~\rm Gyr^{-1}$ in the bar regions, and $0.28 - 0.65$ and $(0.23 - 0.96) ~\rm Gyr^{-1}$ in the arm and bar-end regions. This result supports the idea that the presence of a large amount of diffuse molecular gas makes the SFE low. The suppression of the SFE in the bar has also been seen even when we exclude the diffuse molecular gas components. This suggests that the low SFE appears to be caused not only by a large amount of diffuse molecular gases but also by other mechanisms such as fast cloud-cloud collisions.

Recent galaxy observations show that star formation activity changes depending on galactic environments. In order to understand the diversity of galactic-scale star formation, it is crucial to understand the formation and evolution of giant molecular clouds in an extreme environment. We focus on observational evidence that bars in strongly barred galaxies lack massive stars even though quantities of molecular gas are sufficient to form stars. In this paper, we present a hydrodynamical simulation of a strongly barred galaxy, using a stellar potential which is taken from observational results of NGC1300, and we compare cloud properties between different galactic environments: bar, bar-end and spiral arms. We find that the mean of cloud's virial parameter is ~1 and that there is no environmental dependence, indicating that the gravitationally-bound state of a cloud is not behind the observational evidence of the lack of massive stars in strong bars. Instead, we focus on cloud-cloud collisions, which have been proposed as a triggering mechanism for massive star formation. We find that the collision speed in the bar is faster than those in the other regions. We examine the collision frequency using clouds' kinematics and conclude that the fast collisions in the bar could originate from random-like motion of clouds due to elliptical gas orbits shifted by the bar potential. These results suggest that the observed regions of lack of active star-formation in the strong bar originate from the fast cloud-cloud collisions, which are inefficient in forming massive stars, due to the galactic-scale violent gas motion.

We study effect of magnetic field on massive dense core formation in colliding unequal molecular clouds by performing magnetohydrodynamic simulations with sub-parsec resolution (0.015 pc) that can resolve the molecular cores. Initial clouds with the typical gas density of the molecular clouds are immersed in various uniform magnetic fields. The turbulent magnetic fields in the clouds consistent with the observation by Crutcher et al. (2010) are generated by the internal turbulent gas motion before the collision, if the uniform magnetic field strength is 4.0 $\mu$G. The collision speed of 10 km s$^{-1}$ is adopted, which is much larger than the sound speeds and the Alfv\'{e}n speeds of the clouds. We identify gas clumps with gas densities greater than 5$\times$10$^{-20}$ g cm$^{-3}$ as the dense cores and trace them throughout the simulations to investigate their mass evolution and gravitational boundness. We show that a greater number of massive, gravitationally bound cores are formed in the strong magnetic field (4.0 $\mu$G) models than the weak magnetic field (0.1 $\mu$G) models. This is partly because the strong magnetic field suppresses the spatial shifts of the shocked layer that should be caused by the nonlinear thin shell instability. The spatial shifts promote formation of low-mass dense cores in the weak magnetic field models. The strong magnetic fields also support low-mass dense cores against gravitational collapse. We show that the numbers of massive, gravitationally bound cores formed in the strong magnetic field models are much larger than the isolated, non-colliding cloud models, which are simulated for comparison. We discuss the implications of our numerical results on massive star formation.

Star formation activity depends on galactic-scale environments. To understand the variations in star formation activity, comparing the properties of giant molecular clouds (GMCs) among environments with different star formation efficiency (SFE) is necessary. We thus focus on a strongly barred galaxy to investigate the impact of the galactic environment on the GMCs properties, because the SFE is clearly lower in bar regions than in arm regions. In this paper, we present the $^{12}$CO($1-0$) observations toward the western bar, arm and bar-end regions of the strongly barred galaxy NGC1300 with ALMA 12-m array at a high angular resolution of $\sim$40 pc. We detected GMCs associated with the dark lanes not only in the arm and bar-end regions but also in the bar region, where massive star formation is not seen. Using the CPROPS algorithm, we identified and characterized 233 GMCs across the observed regions. Based on the Kolmogorov-Smirnov test, we find that there is virtually no significant variations in GMC properties (e.g., radius, velocity dispersion, molecular gas mass, and virial parameter) among the bar, arm and bar-end region. These results suggest that systematic differences in the physical properties of the GMCs are not the cause for SFE differences with environments, and that there should be other mechanisms which control the SFE of the GMCs such as fast cloud-cloud collisions in NGC1300.

Using the NANTEN2 Observatory, we carried out a molecular line study of high-mass star forming regions with reflection nebulae, NGC 2068 and NGC 2071, in Orion in the 13CO(J=2-1) transition. The 13CO distribution shows that there are two velocity components at 9.0 and 10.5 km/s . The blue-shifted component is in the northeast associated with NGC 2071, whereas the red-shifted component is in the southwest associated with NGC 2068. The total intensity distribution of the two clouds shows a gap of ~1 pc, suggesting that they are detached at present. A detailed spatial comparison indicates that the two show complementary distributions. The blue-shifted component lies toward an intensity depression to the northwest of the red-shifted component, where we find that a displacement of 0.8 pc makes the two clouds fit well with each other. Furthermore, a new simulation of non-frontal collisions shows that observations from 60 degrees off the collisional axis agreed well with the velocity structure in this region. On the basis of these results, we hypothesize that the two components collided with each other at a projected relative velocity 3.0 km/s estimated to be 0.3 Myr for an assumed axis of the relative motion 60 degrees off the line of sight. We assume that the two most massive early B-type stars in the cloud, illuminating stars of the two reflection nebulae, were formed by collisional triggering at the interfaces between the two clouds. Given the other young high-mass star forming regions, namely, M42, M43, and NGC 2024 (Fukui et al. 2018b; Ohama et al. 2017a), it seems possible that collisional triggering has been independently working to form O-type and early B-type stars in Orion in the last Myr over a projected distance of ~80 pc.

As part our investigation into the Galactic rotation curve, we carried out Very Long Baseline Interferometry (VLBI) observations towards the star-forming region IRAS 01123+6430 using VLBI Exploration of Radio Astrometry (VERA) to measure its annual parallax and proper motion. The annual parallax was measured to be 0.151+/-0.042 mas, which corresponds to a distance of D=6.61^{+2.55}_{-1.44} kpc, and the obtained proper motion components were ({\mu}_{\alpha}cos{\delta}, {\mu}_{\delta})=(-1.44+/-0.15, -0.27+/-0.16) mas yr^{-1} in equatorial coordinates. Assuming Galactic constants of (R_0, {\Theta}_0)=(8.05+/-0.45 kpc, 238+/-14 km s^{-1}), the Galactocentric distance and rotation velocity were measured to be (R, {\Theta})=(13.04+/-2.24 kpc, 239+/-22 km s^{-1}), which are consistent with a flat Galactic rotation curve. The newly estimated distance provides a more accurate bolometric luminosity of the central young stellar object, L_Bol=(3.11+/-2.86)\times 10^3 L_solar, which corresponds to a spectral type of B1-B2. The analysis of 12CO(J=1-0) survey data obtained with the Five College Radio Astronomical Observatory (FCRAO) 14 m telescope shows that the molecular cloud associated with IRAS 01123+6430 consists of arc-like and linear components, which well matches a structure predicted by numerical simulation of the cloud-cloud collision (CCC) phenomenon. The coexistence of arc-like and linear components implies that the relative velocity of initial two clouds was as slow as 3-5 km s^{-1}, which meets the expected criteria of massive star formation where the core mass is effectively increased in the presence of low relative velocity (~3-5 km s^{-1}), as suggested by Takahira et al.(2014).

Sh2-48 is a Galactic HII region located at 3.8 kpc with an O9.5-type star identified at its center. As a part of the FOREST Unbiased Galactic plane Imaging survey using the Nobeyama 45-m telescope (FUGIN) project, we obtained the CO J=1-0 dataset for a large area of Sh2-48 at a spatial resolution of 21"(~0.4 pc), which we used to find a molecular cloud with a total molecular mass of ~3.8x10^4 Mo associated with Sh2-48. The molecular cloud has a systematic velocity shift within a velocity range ~42-47 km/s . On the lower velocity side the CO emission spatially corresponds with the bright 8 {\mu}m filament at the western rim of Sh2-48, while the CO emission at higher velocities is separated at the eastern and western sides of the 8{\mu}m filament. This velocity change forms V-shaped, east-west-oriented feature on the position-velocity diagram. We found that these lower and higher-velocity components are, unlike the infrared and radio continuum data, physically associated with Sh2-48. To interpret the observed V-shaped velocity distribution, we assessed a cloud-cloud collision scenario and found from a comparison between the observations and simulations that the velocity distribution is an expected outcome of a collision between a cylindrical cloud and a spherical cloud, with the cylindrical cloud corresponding to the lower-velocity component, and the two separated components in the higher-velocity part interpretable as the collision-broken remnants of the spherical cloud. Based on the consistency of the ~1.3Myr estimated formation timescale of the HII region with that of the collision, we concluded that the high-mass star formation in Sh2-48 was triggered by the collision.

With the advent of modern observational efforts providing extensive giant molecular cloud catalogues, understanding the evolution of such clouds in a galactic context is of prime importance. While numerous previous numerical and theoretical works have focused on the cloud properties in isolated discs, few have looked into the cloud population in an interacting disc system. We present results of the first study investigating the evolution of the cloud population in galaxy experiencing an M51-like tidal fly-by using numerical simulations including star formation, interstellar medium cooling and stellar feedback. We see the cloud population shift to large unbound clouds in the wake of the companion passage, with the largest clouds appearing as fleeting short-lived agglomerations of smaller clouds within the tidal spiral arms, brought together by large scale streaming motions. These are then sheared apart as they leave the protection of the spiral arms. Clouds appear to lead diverse lives, even within similar environments, with some being born from gas shocked by filaments streaming into the spiral arms, and others from effectively isolated smaller colliding pairs. Overall this cloud population produces a shallower mass function than the disc in isolation, especially in the arms compared to the inter-arm regions. Direct comparisons to M51 observations show similarities between cloud populations, though models tailored to the mass and orbital models of M51 appear necessary to precisely reproduce the cloud population.

We report an observational study of the giant molecular cloud (GMC) associated with the Galactic infrared ring-like structure N35 and two nearby HII regions G024.392+00.072 (HII region A) and G024.510-00.060 (HII region B), using the new CO J=1-0 data obtained as a part of the FOREST Unbiased Galactic Plane Imaging survey with the Nobeyama 45-m telescope (FUGIN) project at a spatial resolution of 21". Our CO data revealed that the GMC, with a total molecular mass of 2.1x10^6Mo, has two velocity components over ~10-15km/s. The majority of molecular gas in the GMC is included in the lower-velocity component (LVC) at ~110-114km/s, while the higher-velocity components (HVCs) at ~118-126km/s consist of three smaller molecular clouds which are located near the three HII regions. The LVC and HVCs show spatially complementary distributions along the line-of-sight, despite large velocity separations of ~5-15km/s, and are connected in velocity by the CO emission with intermediate intensities. By comparing the observations with simulations, we discuss a scenario where collisions of the three HVCs with LVC at velocities of ~10-15km/s can provide an interpretation of these two observational signatures. The intermediate velocity features between the LVC and HVCs can be understood as broad bridge features, which indicate the turbulent motion of the gas at the collision interfaces, while the spatially complementary distributions represent the cavities created in the LVC by the HVCs through the collisions. Our model indicates that the three HII regions were formed after the onset of the collisions, and it is therefore suggested that the high-mass star formation in the GMC was triggered by the collisions.

The Orion Nebula Cluster toward the HII region M42 is the most outstanding young cluster at the smallest distance 410pc among the rich high-mass stellar clusters. By newly analyzing the archival molecular data of the 12CO(J=1-0) emission at 21" resolution, we identified at least three pairs of complementary distributions between two velocity components at 8km/s and 13km/s. We present a hypothesis that the two clouds collided with each other and triggered formation of the high-mass stars, mainly toward two regions including the nearly ten O stars, theta1 Ori and theta2 Ori, in M42 and the B star, NU Ori, in M43. The timescale of the collision is estimated to be ~0.1Myr by a ratio of the cloud size and velocity corrected for projection, which is consistent with the age of the youngest cluster members less than 0.1Myr. The majority of the low-mass cluster members were formed prior to the collision in the last one Myr. We discuss implications of the present hypothesis and the scenario of high-mass star formation by comparing with the other eight cases of triggered O star formation via cloud-cloud collision.

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 performed sub-parsec ($\sim$ 0.014 pc) scale simulations of cloud-cloud collisions of two idealized turbulent molecular clouds (MCs) with different masses in the range of $0.76 - 2.67 \times 10^4$M$_{\odot}$ and with collision speeds of 5 $-$ 30 km/s. Those parameters are larger than Takahira, Tasker and Habe (2014) (paper I) in which the colliding system showed a partial gaseous arc morphology that supports the NANTEN observations of objects indicated to be colliding MCs by numerical simulations. Gas clumps with density greater than $10^{-20}$ g cm$^{-3}$ were identified as pre-stellar cores and tracked through the simulation to investigate the effect of mass of colliding clouds and collision speeds on the resulting core population. Our results demonstrate that smaller cloud property is more important for results of cloud cloud collisions. The mass function of formed cores can be approximated by a power law relation with index $\gamma$ = -1.6 in slower cloud cloud collisions ($v \sim 5 $ km/s), in good agreement with observation of MCs. A faster relative velocity increases the number of cores formed in the early stage of collisions and shortens gas accretion phase of cores in the shocked region, leading to suppression of core growth. The bending point appears in the high mass part of the core mass function and the bending point mass decreases with increasing of the collision velocity for the same combination of colliding clouds. The high mass part of the core mass function than the bending point mass can be approximated by a power law with $\gamma$ = -2.3 that is very similar to the power index of the massive part of the observed initial stellar mass function. We discuss implication of our results for the massive star formation in our Galaxy.

We investigate star formation occurring in idealised giant molecular clouds, comparing structures that evolve in isolation versus those undergoing a collision. Two different collision speeds are investigated and the impact of photoionising radiation from the stars is determined. We find that a colliding system leads to more massive star formation both with and without the addition of feedback, raising overall star formation efficiencies (SFE) by a factor of 10 and steepening the high-mass end of the stellar mass function. This rise in SFE is due to increased turbulent compression during the cloud collision. While feedback can both promote and hinder star formation in the isolated system, it increases the SFE by approximately 1.5 times in the colliding case when the thermal speed of the resulting HII regions matches the shock propagation speed in the collision.

We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud.

High-mass star formation is one of the top-priority issues in astrophysics. Recent observational studies are revealing that cloud-cloud collisions may play a role in high-mass star formation in several places in the Milky Way and the Large Magellanic Cloud. The Trifid Nebula M20 is a well known galactic HII region ionized by a single O7.5 star. In 2011, based on the CO observations with NANTEN2 we reported that the O star was formed by the collision between two molecular clouds ~0.3,Myr ago. Those observations identified two molecular clouds towards M20, traveling at a relative velocity of 7.5 km/s. This velocity separation implies that the clouds cannot be gravitationally bound to M20, but since the clouds show signs of heating by the stars there they must be spatially coincident with it. A collision is therefore highly possible. In this paper we present the new CO J=1-0 and J=3-2 observations of the colliding clouds in M20 performed with the Mopra and ASTE telescopes. The high resolution observations revealed the two molecular clouds have peculiar spatial and velocity structures, i.e., the spatially complementary distribution between the two clouds and the bridge feature which connects the two clouds in velocity space. Based on a new comparison with numerical models, we find that this complementary distribution is an expected outcome of cloud-cloud collisions, and that the bridge feature can be interpreted as the turbulent gas excited at the interface of the collision. Our results reinforce the cloud-cloud collision scenario in M20.

We investigated the effect of photoionising feedback inside turbulent star-forming clouds, comparing the resultant star formation in both idealised profiles and more realistic cloud structures drawn from a global galaxy simulation. We performed a series of numerical simulations which compared the effect of star formation alone, photoionisation and photoionisation plus supernovae feedback. In the idealised cloud, photoionisation suppresses gas fragmentation at early times, resulting in the formation of more massive stars and an increase in the star formation efficiency. At later times, the dispersal of the dense gas causes the radiative feedback effect to switch from positive to negative as the star formation efficiency drops. In the cloud extracted from the global simulation, the initial cloud is heavily fragmented prior to the stellar feedback beginning and is largely structurally unaffected by the late injection of radiation energy. The result is a suppression of the star formation. We conclude that the efficiency of feedback is heavily dependent on the gas structure, with negative feedback dominating when the density is high.

We explore the impact of star formation and thermal stellar feedback on the giant molecular cloud (GMC) population forming in a M83-type barred spiral galaxy. We compare three high-resolution simulations (1.5 pc cell size) with different star formation/feedback models: one with no star formation, one with star formation but no feedback, and one with star formation and thermal energy injection. We analyze the resulting population of clouds, finding that we can identify the same population of massive, virialized clouds and transient, low-surface density clouds found in our previous work (that did not include star formation or feedback). Star formation and feedback can affect the mix of clouds we identify. In particular, star formation alone simply converts dense cloud gas into stars with only a small change to the cloud populations, principally resulting in a slight decrease in the transient population. Feedback, however, has a stronger impact: while it is not generally sufficient to entirely destroy the clouds, it does eject gas out of them, increasing the gas density in the inter-cloud region. This decreases the number of massive clouds, but substantially increases the transient cloud population. We also find that feedback tends to drive a net radial inflow of massive clouds, leading to an increase in the star formation rate in the bar region. We examine a number of possible reasons for this and conclude that it is possible that the drag force from the enhanced intercloud density could be responsible.

The bipolar-hyper shell (BHS) model for the North Polar Spurs (NPS-E, -W, and Loop I) and counter southern spurs (SPS-E and -W) is revisited based on numerical hydrodynamical simulations. Propagations of shock waves produced by energetic explosive events in the Galactic Center are examined. Distributions of soft X-ray brightness on the sky at 0.25, 0.7, and 1.5 keV in a +/-50 deg x +/-50 deg region around the Galactic Center are modeled by thermal emission from high-temperature plasma in the shock-compressed shell considering shadowing by the interstellar HI and H2 gases. The result is compared with the ROSAT wide field X-ray images in R2, 4 and 6 bands. The NPS and southern spurs are well reproduced by the simulation as shadowed dumbbell-shaped shock waves. We discuss the origin and energetics of the event in relation to the starburst and/or AGN activities in the Galactic Center. [ High resolution pdf is available at http://www.ioa.s.u-tokyo.ac.jp/~sofue/htdocs/2016bhs/ ]

We perform dry merger simulations to investigate the role of dry mergers in the size growth of early-type galaxies in high density environments. We replace the virialized dark matter haloes obtained by a large cosmological $N$-body simulation with $N$-body galaxy models consisting of two components, a stellar bulge and a dark matter halo, which have higher mass resolution than the cosmological simulation. We then re-simulate nine cluster forming regions, whose masses range from 1e+14 Msun to 5e+14 Msun. Masses and sizes of stellar bulges are also assumed to satisfy the stellar mass--size relation of high-z compact massive early-type galaxies. We find that dry major mergers considerably contribute to the mass and size growth of central massive galaxies. One or two dry major mergers double the average stellar mass and quadruple the average size between $z=2$ and $z=0$. These growths favorably agree with observations. Moreover, the density distributions of our simulated central massive galaxies grow from the inside-out, which is consistent with recent observations. The mass--size evolution is approximated as R propto M_{*}^{alpha}, with alpha ~ 2.24. Most of our simulated galaxies are efficiently grown by dry mergers, and their stellar mass--size relations match the ones observed in the local Universe. Our results show that the central galaxies in the cluster haloes are potential descendants of high-z (z ~ 2-3) compact massive early-type galaxies. This conclusion is consistent with previous numerical studies which investigate the formation and evolution of compact massive early-type galaxies.

We investigate the longevity of broad bridge features in position-velocity diagrams that appear as a result of cloud-cloud collisions. Broad bridges will have a finite lifetime due to the action of feedback, conversion of gas into stars and the timescale of the collision. We make a series of analytic arguments with which to estimate these lifetimes. Our simple analytic arguments suggest that for collisions between clouds larger than R~10 pc the lifetime of the broad bridge is more likely to be determined by the lifetime of the collision rather than the radiative or wind feedback disruption timescale. However for smaller clouds feedback becomes much more effective. This is because the radiative feedback timescale scales with the ionising flux Nly as R^{7/4}Nly^{-1/4} so a reduction in cloud size requires a relatively large decrease in ionising photons to maintain a given timescale. We find that our analytic arguments are consistent with new synthetic observations of numerical simulations of cloud-cloud collisions (including star formation and radiative feedback). We also argue that if the number of observable broad bridges remains ~ constant, then the disruption timescale must be roughly equivalent to the collision rate. If this is the case our analytic arguments also provide collision rate estimates, which we find are readily consistent with previous theoretical models at the scales they consider (clouds larger than about 10 pc) but are much higher for smaller clouds.

We investigate the evolution of two bars formed in fully self-consistent hydrodynamic simulations of the formation of Milky Way-mass galaxies. One galaxy shows higher central mass concentration and has a longer and stronger bar than the other at $z = 0$. The stronger bar evolves by transferring its angular momentum mainly to the dark halo. Consequently the rotation speed of the bar decreases with time, while the amplitude of the bar increases with time. These features qualitatively agree with the results obtained by idealized simulations. The pattern speed of the stronger bar largely goes up and down within a half revolution in its early evolutionary stage. These oscillations occur when the bar is misaligned with the $m = 4$ mode Fourier component. These oscillations correlate with the oscillations in the triaxilality of the dark matter halo, but differently from the way identified by idealized simulations. The amplitude of the weaker bar does not increase despite the fact that its rotation slows down with time.This result contradicts what is expected from idealized simulations and is caused by the decline of the central density associated with the mass loss and feedback from the stellar populations. The amplitude of the weaker bar is further weakens by the angular momentum injection by the interactions with stellar clumps in the disk. In the both galaxies, the bars are terminated around the 4:1 resonance.

Collisions between giant molecular clouds are a potential mechanism for triggering the formation of massive stars, or even super star clusters. The trouble is identifying this process observationally and distinguishing it from other mechanisms. We produce synthetic position-velocity diagrams from models of: cloud-cloud collisions, non-interacting clouds along the line of sight, clouds with internal radiative feedback and a more complex cloud evolving in a galactic disc, to try and identify unique signatures of collision. We find that a broad bridge feature connecting two intensity peaks, spatially correlated but separated in velocity, is a signature of a high velocity cloud-cloud collision. We show that the broad bridge feature is resilient to the effects of radiative feedback, at least to around 2.5Myr after the formation of the first massive (ionising) star. However for a head on 10km/s collision we find that this will only be observable from 20-30 per cent of viewing angles. Such broad-bridge features have been identified towards M20, a very young region of massive star formation that was concluded to be a site of cloud-cloud collision by Torii et al (2011), and also towards star formation in the outer Milky Way by Izumi et al (2014).

Cloud collision have been proposed as a way to link the small-scale star formation process with the observed global relation between the surface star formation rate and gas surface density. We suggest that this model can be improved further by allowing the productivity of such collisions to depend on the relative velocity of the two clouds. Our adjustment implements a simple step function that results in the most successful collisions being at the observed velocities for triggered star formation. By applying this to a high resolution simulation of a barred galaxy, we successfully reproduce the observational result that the star formation efficiency (SFE) in the bar is lower than that in the spiral arms. This is not possible when we use an efficiency dependent on the internal turbulence properties of the clouds. Our results suggest that high velocity collisions driven by the gravitational pull of the clouds are responsible for the low bar SFE.

We performed sub-parsec (~0.06pc) scale simulations of two idealised molecular clouds with different masses undergoing a collision. Gas clumps with density greater than 1e-20 g/cm3 (0.3e4 cm-3) were identified as pre-stellar cores and tracked through the simulation. The colliding system showed a partial gas arc morphology with core formation in the oblique shock-front at the collision interface. These characteristics support NANTEN observations of objects suspected to be colliding giant molecular clouds (GMCs). We investigated the effect of turbulence and collision speed on the resulting core population and compared the cumulative mass distribution to cores in observed GMCs. Our results suggest that a faster relative velocity increases the number of cores formed but that cores grow via accretion predominately while in the shock front, leading to a slower shock being more important for core growth. The core masses obey a power law relation with index gamma = -1.6, in good agreement with observations. This suggests that core production through collisions should follow a similar mass distribution as quiescent formation, albeit at a higher mass range. If cores can be supported against collapse during their growth, the estimated ram pressure from gas infall is of the right order to counter the radiation pressure and form a star of 100Msun.

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.

We present the results of sulfur monoxide, SO line emission observations of G0.253+0.016 with the Atacama Large Millimeter/submillimeter Array (ALMA) at an angular resolution of 1.7 arcsec. The dense and massive molecular cloud of G0.253+0.016 is highly sub-structured, yet shows no obvious signs of cluster formation. We found three outstanding features of the cloud from the SO emission, namely, shell structure of radius 1.3 pc, large velocity gradients of 20 km/s/pc with the cloud, and cores with large velocity dispersions (30-40 km/s) around the shell structure. We suggest that these large-velocity dispersion cores will form high-mass stars in the future. In attempt to explore the formation scenario of the dense cloud, we compared our results with numerical simulations, thus, we propose that G0.253+0.016 may have formed due to a cloud-cloud collision process.

We investigate the impact of galactic environment on the properties of simulated giant molecular clouds formed in a M83-type barred spiral galaxy. Our simulation uses a rotating stellar potential to create the grand design features and resolves down to 1.5 pc. From the comparison of clouds found in the bar, spiral and disc regions, we find that the typical GMC is environment independent, with a mass of 5e+5 Msun and radius 11 pc. However, the fraction of clouds in the property distribution tails varies between regions, with larger, more massive clouds with a higher velocity dispersion being found in greatest proportions in the bar, spiral and then disc. The bar clouds also show a bimodality that is not reflected in the spiral and disc clouds except in the surface density, where all three regions show two distinct peaks. We identify these features as being due to the relative proportion of three cloud types, classified via the mass-radius scaling relation, which we label A, B and C. Type A clouds have the typical values listed above and form the largest fraction in each region. Type B clouds are massive giant molecular associations while Type C clouds are unbound, transient clouds that form in dense filaments and tidal tails. The fraction of each clouds type depends on the cloud-cloud interactions, which cause mergers to build up the GMA Type Bs and tidal features in which the Type C clouds are formed. The number of cloud interactions is greatest in the bar, followed by the spiral, causing a higher fraction of both cloud types compared to the disc. While the cloud types also exist in lower resolution simulations, their identification becomes more challenging as they are not well separated populations on the mass-radius relation or distribution plots. Finally, we compare the results for three star formation models to estimate the star formation rate and efficiency in each region.

We report the discovery of a new ultra-bright submillimeter galaxy (SMG) behind the Small Magellanic Cloud (SMC). This SMG is detected as a 43.3+-8.4 mJy point source (MM J01071-7302, hereafter MMJ0107) in the 1.1 mm continuum survey of the SMC by AzTEC on the ASTE telescope. MMJ0107 is also detected in the radio (843 MHz), Herschel/SPIRE, Spitzer MIPS 24 {\mu}m, all IRAC bands, Wide-field Infrared Survey Explorer, and near-infrared (J, H, KS). We find an optical (U, B, V) source, which might be the lensing object, at a distance of 1.4 arcsec from near-infrared and IRAC sources. Photometric redshift estimates for the SMG using representative spectral energy distribution templates show the redshifts of 1.4-3.9. We estimate total far-infrared luminosity of (0.3-2.2)x10^14 {\mu}^-1 L_sun and a star formation rate of 5600-39, 000 {\mu}^-1 M_sun yr^-1, where {\mu} is the gravitational magnification factor. This apparent extreme star formation activity is likely explained by a highly magnified gravitational lens system.

We investigate the properties of massive, dense clouds formed in a barred galaxy and their possible relation to star formation, performing a two-dimensional hydrodynamical simulation with the gravitational potential obtained from the 2Mass data from the barred spiral galaxy, M83. Since the environment for cloud formation and evolution in the bar region is expected to be different from that in the spiral arm region, barred galaxies are a good target to study the environmental effects on cloud formation and the subsequent star formation. Our simulation uses for an initial 80 Myr an isothermal flow of non-self gravitating gas in the barred potential, then including radiative cooling, heating and self-gravitation of the gas for the next 40 Myr, during which dense clumps are formed. We identify many cold, dense gas clumps for which the mass is more than $10^4M_{\odot}$ (a value corresponding to the molecular clouds) and study the physical properties of these clumps. The relation of the velocity dispersion of the identified clump's internal motion with the clump size is similar to that observed in the molecular clouds of our Galaxy. We find that the virial parameters for clumps in the bar region are larger than that in the spiral arm region. From our numerical results, we estimate star formation in the bar and spiral arm regions by applying the simple model of Krumholtz and McKee (2005). The mean relation between star formation rate and gas surface density agrees well with the observed Kennicutt-Schmidt relation. The SFE in the bar region is $\sim 60 %$ of the spiral arm region. This trend is consistent with observations of barred galaxies.

Recent observations show evidence that high-z (z\sim 2 - 3) early-type galaxies (ETGs) are more compact than those with comparable mass at z\sim 0. Such a size evolution is most likely explained by the `Dry Merger Sceanario'. However, previous studies based on this scenario are not able to consistantly explain both the properties of the high-z compact massive ETGs and the local ETGs. We investigate the effect of multiple sequential dry minor mergers on the size evolution of the compact massive ETGs. From an analysis of the Millennium Simulation Database, we show that such minor (stellar mass ratio $M_{2}/M_{1} < 1/4$) mergers are extremely common during hierarchical structure formation. We perform N-body simulations of sequential minor mergers with parabolic and head-on orbits, including a dark matter component and a stellar component. Typical mass ratios of the minor mergers are $1/20 < M_{2}/M_{1} < 1/10$. We show that sequential minor mergers of compact satellite galaxies are the most efficient at promoting size growth and decreasing the velocity dispersion of the compact massive ETGs in our simulations. The change of stellar size and density of the merger remnants is consistent with recent observations. Furthermore, we construct the merger histories of candidates for the high-z compact massive ETGs using the Millennium Simulation Database, and estimate the size growth of the galaxies by the dry minor merger scenario. We can reproduce the mean size growth factor between $z=2$ and $z=0$, assuming the most efficient size growth obtained during sequential minor mergers in our simulations. However, we note that our numerical result is only valid for merger histories with typical mass ratios between 1/20 and 1/10 with parabolic and head-on orbits, and that our most efficient size growth efficiency is likely to an upper limit.

We have performed VLBI and single-dish monitoring of 22 GHz H$_{2}$O maser emission from the high mass young stellar object G353.273+0.641 with VERA (VLBI Exploration of Radio Astrometry) and Tomakamai 11-m radio telescope. Two maser flares have been detected, separated almost two years. Frequent VLBI monitoring has revealed that these flare activities have been accompanied by structural change of the prominent shock front traced by H2O maser alignments. We have detected only blue-shifted emissions and all maser features have been distributed within very small area of 200 $\times$ 200 au$^{2}$ in spite of wide velocity range (> 100 km s$^{-1}$). The light curve shows notably intermittent variation and suggests that the H$_{2}$O masers in G353.273+0.641 are excited by episodic radio jet. The time-scale of \sim2 yr and characteristic velocity of \sim500 km s$^{-1}$ also support this interpretation. Two isolated velocity components of C50 (-53 \pm 7 km s$^{-1}$) and C70 (-73 \pm 7 km s$^{-1}$) have shown synchronised linear acceleration of the flux weighted $V_{\rmn{LSR}}$ values (\sim-5 km s$^{-1}$ yr$^{-1}$) during the flare phase. This can be converted to the lower-limit momentum rate of 1.1 \times 10$^{-3}$ M$_{\sun}$ km s$^{-1}$ yr$^{-1}$. Maser properties are quite similar to that of IRAS 20126+4104 especially. This corroborates the previous suggestion that G353.273+0.641 is a candidate of high mass protostellar object. The possible pole-on geometry of disc-jet system can be suitable for direct imaging of the accretion disc in this case.

We develop one-zone galaxy formation models in the early Universe, taking into account dust formation and evolution by supernova (SN) explosions. Especially we focus on the time evolution of dust size distribution, because ${\rm H}_{2}$ formation on the dust surface plays a critical role in the star formation process in the early Universe. In the model we assume that star formation rate (SFR) is proportional to the total amount of ${\rm H}_{2}$. We consistently treat (i) the formation and size evolution of dust, (ii) the chemical reaction networks including ${\rm H}_{2}$ formation both on the surface of dust and in gas phase, and (iii) the SFR in the model. First, we find that, because of dust destruction due to both reverse and forward shocks driven by SNe, H$_{2}$ formation is more suppressed than that without dust destruction. At the galaxy age of $sim0.8 {\rm Gyr}$, for galaxy models with virial mass $M_{\rm vir}=10^{9} M_{\odot}$ and formation redshift $z_{\rm vir}=10$, the molecular fraction is 2.5 orders of magnitude less in the model with dust destruction by both shocks than that in the model without dust destruction. Second, we show that the H$_{2}$ formation rate strongly depends on the ISM density around SN progenitors. The SFR in higher ISM density is lower, since dust destruction by reverse shocks is more effective in higher ISM density. We conclude that not only the amount but also the size distribution of dust being related with the star formation activity strongly affects the evolution of galaxies in the early Universe.

We report on the astrometric observations of the 22 GHz H2O masers in the high mass star-forming region G5.89-0.39 with VERA (VLBI Exploration of Radio Astrometry). Newly derived distance of 1.28^{+0.09}_{-0.08} kpc is the most precise and significantly nearer than previous values. We revised physical parameters and reconsidered nature of G5.89-0.39 based on the new distance as follows. (1) The ionizing star of the ultra compact (UC) HII region is a late O-type (O8 - 8.5) zero age main sequence (ZAMS) star, consistent with previously established limits based on its infrared spectral line emission. (2) Crescent-like maser alignment at the position of the O type ZAMS star may trace accretion disk (or its remnant), which suggests that the star is still young and before complete evaporation of circumstellar materials. (3) Although the revised mass for the east-west outflow has been reduced, it still quite large (100 Msun) which indicates that a significant fraction of the mass is entrained material and that the dynamical age significantly underestimates the actual outflow age. Our newly-derived distance emphasizes that G5.89-0.39 is one of the nearest targets to investigate ongoing high-mass star formation and evolution in a compact cluster containing a young O-type star.

We investigate the star formation and chemical evolution in the early universe by considering the merging history of the Galaxy in the {\Lambda}CDM scenario according to the extended Press-Schechter theory. We give some possible constraints from comparisons with observation of extremely metal-poor (EMP) stars. We demonstrate that (1) The hierarchical structure formation can explain the characteristics of the observed metallicity distribution function (MDF) including a break around [Fe/H]~-4. (2) A high mass IMF of peak mass ~10Msun with the contribution of binaries, derived from the statistics of carbon enhanced EMP stars (Komiya et al. 2007), predicts the frequency of low-mass survivors consistent with the number of EMP stars observed for -4~<[Fe/H]~<-2.5. (3) The stars formed from primordial gas before the first supernova explosions in their host mini-halos are assigned to the HMP stars with [Fe/H]~-5. (4) There is no indication of significant changes in the IMF and the binary contribution at metallicity -4~<[Fe/H]~<-2.5, or even larger as long as the field stars of Galactic halo are concerned. We further study the effects of the surface pollution through the accretion of ISM along the chemical and dynamical evolution of the Galaxy for low-mass Pop.III and EMP survivors. Because of shallower potential of smaller halos, the accretion of ISM in the mini-halos in which these stars were born dominates the surface metal pollution. This can account for the surface iron abundances as observed for the HMP stars if the cooling and concentration of gas in their birth mini-halos is taken into account. We also study the feedback effect from the very massive Pop. III stars. The metal pre-pollution by PISNe is shown to be compatible with the observed lack of their nucleosynthetic signatures when some positive feedback on gas cooling works and changes IMF from being very massive to being high mass.

We present recent advances in theoretical studies of the formation and evolution of dust in primordial supernovae (SNe) that are considered to be the main sources of dust in the early universe. Being combined with the results of calculations of dust formation in the ejecta of Population III SNe, the investigations of the evolution of newly formed dust within supernova remnants (SNRs) show that smaller grains are predominantly destroyed by sputtering in the shocked gas, while larger grains are injected into the ambient medium. The mass of dust grains surviving the destruction in SNRs reaches up to 0.1--15 $M_\odot$, which is high enough to account for the content of dust observed for the host galaxies of quasars at $z > 5$. In addition, the transport of dust formed in the ejecta causes the formation of low-mass stars in the dense shell of primordial SNRs and affects the elemental composition of those stars. We also show that the flat extinction curve is expected in the high-redshift universe where SNe are the possible sources of dust.

We study the possibility that minor mergers resolve the loss cone depletion problem, which is the difficulty occured in the coalescence process of the supermassive black hole (SMBH) binary, by performing numerical simulations with a highly accurate $N$-body code. We show that the minor merger of a dwarf galaxy disturbs stellar orbits in the galactic central region of the host galaxy where the loss cone depletion is already caused by the SMBH binary. The disturbed stars are supplied into the loss cone. Stars of the dwarf galaxy are also supplied into the loss cone. The gravitational interactions between the SMBH binary and these stars become very effective. The gravitational interaction decreases the binding energy of the SMBH binary effectively. As a result, the shrink of the separation of the SMBH binary is accelerated. Our numerical results strongly suggest that the minor mergers are one of the important processes to reduce the coalescence time of the SMBH binary much less than the Hubble time.

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 investigate the evolution of dust formed in Population III supernovae (SNe) by considering its transport and processing by sputtering within the SN remnants (SNRs). We find that the fates of dust grains within SNRs heavily depend on their initial radii $a_{\rm ini}$. For Type II SNRs expanding into the ambient medium with density of $n_{\rm H,0} = 1$ cm$^{-3}$, grains of $a_{\rm ini} < 0.05$ $\mu$m are detained in the shocked hot gas and are completely destroyed, while grains of $a_{\rm ini} > 0.2$ $\mu$m are injected into the surrounding medium without being destroyed significantly. Grains with $a_{\rm ini}$ = 0.05-0.2 $\mu$m are finally trapped in the dense shell behind the forward shock. We show that the grains piled up in the dense shell enrich the gas up to 10$^{-6}-10^{-4}$ $Z_\odot$, high enough to form low-mass stars with 0.1-1 $M_\odot$. In addition, [Fe/H] in the dense shell ranges from -6 to -4.5, which is in good agreement with the ultra-metal-poor stars with [Fe/H] < -4. We suggest that newly formed dust in a Population III SN can have great impacts on the stellar mass and elemental composition of Population II.5 stars formed in the shell of the SNR.

We have made multi-epoch VLBI observations of H2O maser emission in the massive star forming region IRAS 06061+2151 with the Japanese VLBI network (JVN) from 2005 May to 2007 October. The detected maser features are distributed within an 1\arcsec$\times$1\arcsec (2000 au$\times$2000 au at the source position) around the ultra-compact H {\small\bf II} region seen in radio continuum emission. Their bipolar morphology and expanding motion traced through their relative proper motions indicate that they are excited by an energetic bipolar outflow. Our three-dimensional model fitting has shown that the maser kinematical structure in IRAS 06061+2151 is able to be explained by a biconical outflow with a large opening angle ($>$ 50\degr). The position angle of the flow major axis coincides very well with that of the large scale jet seen in 2.1$\:\mu\rmn{m}$ hydrogen emission. This maser geometry indicates the existence of dual structures composed of a collimated jet and a less collimated massive molecular flow. We have also detected a large velocity gradient in the southern maser group. This can be explained by a very small (on a scale of several tens of au) and clumpy (the density contrast by an order of magnitude or more) structure of the parental cloud. Such a structure may be formed by strong instability of shock front or splitting of high density core.

We discuss the star formation history of the Galaxy, based on the observations of extremely metal-poor stars (EMP) in the Galactic halo, to gain an insight into the evolution and structure formation in the early universe. The initialmass function (IMF) of EMP stars is derived from the observed fraction of carbon-enhanced EXP (CEMP) stars among the EMP survivors, which are thought to originate from the evolution in the close binary systems with mass transfer. Relying upon the theory of the evolution of EMP stars and of their binary evolution, we find that stars of metallicity [Fe/H]<-2.5 were formed at typical mass of ~10M_sun. The top heavy IMF thus obtained is applied to study the early chemical evolution of the Galaxy. We construct the merging history of our Galaxy semi-analytically and derive the metallicity distribution function (MDF) of low-mass EMP stars that survive to date with taking into account the contribution of binary systems. It is shown that the resultant MDF can well reproduce the observed distribution of EMP survivors, and, in particular, that they almost all stem from a less-mass companion in binary systems. We also investigate how first stars affect the MDF of EMP stars.

We investigate the evolution of dust that formed at Population III supernova (SN) explosions and its processing through the collisions with the reverse shocks resulting from the interaction of the SN ejecta with the ambient medium. In particular, we investigate the transport of the shocked dust within the SN remnant (SNR), and its effect on the chemical composition, the size distribution, and the total mass of dust surviving in SNR. We find that the evolution of the reverse shock, and hence its effect on the processing of the dust depends on the thickness of the envelope retained by the progenitor star. Furthermore, the transport and survival of the dust grains depend on their initial radius, a_{ini}, and composition: For Type II SNRs expanding into the interstellar medium (ISM) with a density of n_{H,0}=1 cm^{-3}, small grains with a_{ini} < ~ 0.05 micron are completely destroyed by sputtering in the postshock flow, while grains with a_{ini}= 0.05--0.2 micron are trapped into the dense shell behind the forward shock. Very large grains of a_{ini} > ~ 0.2 micron are ejected into the ISM without decreasing their sizes significantly. We find that the total mass fraction of dust that is destroyed by the reverse shock ranges from 0.2 to 1.0, depending on the energy of the explosion and the density of the ambient ISM. The results of our calculations have significant impact on the abundance pattern of subsequent generation of stars that form in the dense shell of primordial SNRs.

We performed phase-reference very long baseline interferometry (VLBI) observations on five radio-loud narrow-line Seyfert 1 galaxies (NLS1s) at 8.4 GHz with the Japanese VLBI Network (JVN). Each of the five targets (RXS J08066+7248, RXS J16290+4007, RXS J16333+4718, RXS J16446+2619, and B3 1702+457) in milli-Jansky levels were detected and unresolved in milli-arcsecond resolutions, i.e., with brightness temperatures higher than 10^7 K. The nonthermal processes of active galactic nuclei (AGN) activity, rather than starbursts, are predominantly responsible for the radio emissions from these NLS1s. Out of the nine known radio-loud NLS1s, including the ones chosen for this study, we found that the four most radio-loud objects exclusively have inverted spectra. This suggests a possibility that these NLS1s are radio-loud due to Doppler beaming, which can apparently enhance both the radio power and the spectral frequency.

We present the basic features and the activities of Japanese VLBI network (JVN), a newly-established VLBI network with baselines ranging from 50 to 2560 km spreading across the Japanese islands, and capable of observing at 6.7, 8.4, and 22 GHz. We show a number of results of JVN observations: 8.4-GHz continuum images of a Giga-hertz Peaked Spectrum (GPS) source and radio-loud Narrow-Line Seyfert 1 galaxies (NLS1s), the spatial and velocity structures of water masers in NML Cygni as well as methanol masers in Cep A, and demonstrative observations with the bigradient phase referencing.

We propose bigradient phase referencing (BPR), a new radio-observation technique, and report its performance using the Japanese very-long-baseline-interferometry network (JVN). In this method, a weak source is detected by phase-referencing using a primary calibrator, in order to play a role as a secondary calibrator for phase-referencing to a weak target. We will be given the opportunity to select a calibrator from lots of milli-Jansky sources, one of which may be located at the position closer to the target. With such a smaller separation, high-quality phase-referencing can be achieved. Furthermore, a subsequent more-sophisticated calibration can relocate array's focus to a hypothetical point much closer to the target; a higher quality of phase referencing is available. Our demonstrative observations with strong radio sources have proved the capabilities of BPR in terms of image dynamic ranges and astrometric reproducibility. The image dynamic range on a target has been improved with a factor of about six compared to that of normal phase-referencing; the resultant position difference of target's emission between two epochs was only 62+-50 micro-arcsecond, even with less than 2300-km baselines at 8.4 GHz and fast-switching of a target-calibrator pair of a 2.1-degree separation.

We compare the highly clustered populations of very high redshift galaxies with proto-clusters identified numerically in a standard $\Lambda$CDM universe ($\Omega_0=0.3, \lambda_0=0.7$) simulation. We evolve 256^3 dark matter particles in a comoving box of side 150h^{-1}Mpc. By the present day there are 63 cluster sized objects of mass in excess of 10^{14}h^{-1}Mo in this box. We trace these clusters back to higher redshift finding that their progenitors at z=4--5 are extended regions of typically 20--40 Mpc (comoving) in size, with dark halos of mass in excess of 10^{12}h^{-1}Mo and are overdense by typically 1.3--13 times the cosmological mean density. Comparison with the observation of Lyman alpha emitting (LAEs) galaxies at z=4.86 and at z=4.1 indicates that the observed excess clustering is consistent with that expected for a proto-cluster region if LAEs typically correspond to massive dark halos of more than 10^{12}h^{-1}Mo. We give a brief discussion on the relation between high redshift concentration of massive dark halos and present day rich clusters of galaxies.

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.

We investigate the destruction of dust grains by sputtering in the high-velocity interstellar shocks driven by supernovae (SNe) in the early universe to reveal the dependence of the time-scale of dust destruction on the gas density $n_{{\rm H}, 0}$ in the interstellar medium (ISM) as well as on the progenitor mass $M_{\rm pr}$ and explosion energy $E_{\rm 51}$ of SN. The sputtering yields for the combinations of dust and ion species of interest to us are evaluated by applying the so-called universal relation with a slight modification. The dynamics of dust grains and their destruction by sputtering in shock are calculated by taking into account the size distribution of each dust species, together with the time evolution of temperature and density of gas in spherically symmetric shocks. The results of calculations show that the efficiency of dust destruction depends not only on the sputtering yield but also on the initial size distribution of each grain species. The efficiency of dust destruction increases with increasing $E_{\rm 51}$ and/or increasing $n_{{\rm H}, 0}$, but is almost independent of $M_{\rm pr}$ as long as $E_{\rm 51}$ is the same. The mass of gas swept up by shock is the increasing function of $E_{\rm 51}$ and the decreasing function of $n_{{\rm H}, 0}$. Combining these results, we present the approximation formula for the time-scale of destruction for each grain species in the early universe as a function of $E_{\rm 51}$ and $n_{{\rm H}, 0}$. This formula is applicable for investigating the evolution of dust grains at the early epoch of the universe with the metallicity of $Z \la 10^{-3}$ $Z_\odot$. The effects of the cooling processes of gas on the destruction of dust are briefly discussed.

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.

We investigate how cluster morphology is affected by the cosmological constant in low-density universes. Using high-resolution cosmological N-body/SPH simulations of flat (\Omega_0 = 0.3, \lambda_0 = 0.7, \Lambda CDM) and open (\Omega_0 = 0.3, \lambda_0 = 0, OCDM) cold dark matter universes, we calculate statistical indicators to quantify the irregularity of the cluster morphologies. We study axial ratios, center shifts, cluster clumpiness, and multipole moment power ratios as indicators for the simulated clusters at z=0 and 0.5. Some of these indicators are calculated for both the X-ray surface brightness and projected mass distributions. In \Lambda CDM all these indicators tend to be larger than those in OCDM at z=0. This result is consistent with the analytical prediction of Richstone, Loeb, & Turner, that is, clusters in \Lambda CDM are formed later than in OCDM, and have more substructure at z=0. We make a Kolmogorov-Smirnov test on each indicator for these two models. We then find that the results for the multipole moment power ratios and the center shifts for the X-ray surface brightness are under the significance level (5%). We results also show that these two cosmological models can be distinguished more clearly at z=0 than z = 0.5 by these indicators.

The central kiloparsec of M81 has been observed in the CO(J=1-0) line and the 3 mm continuum at 100 pc resolution in an attempt to probe molecular gas, and to search for the nuclear inner Lindblad resonance (NILR), around the low-luminosity AGN M81*. We found the following. (1) Molecular gas in the central kpc is mainly on a ``pseudoring'' or a spiral arm at a radius of about 500 pc. (2) The region within ~300 pc from the nucleus is mostly devoid of molecular gas except for diffuse one; in particular, there is neither a giant molecular cloud that is now accreting on the nucleus nor a conspicuous gas feature that can be identified as an NILR. (3) The 3 mm continuum emission shows significant intraday variation, suggesting an emitting region of ~100 AU. (4) The 3-sigma upper limit for CO absorption toward the continuum source is tau_CO(0->1) dV < 0.1 for a linewidth of 10 km/s. The dearth of accreting molecular gas in the vicinity of the nucleus may explain the low luminosity of M81*.

We study the differential energy distribution of dark matter halos, carrying out cosmological N-body simulation.We give an analytical formula of the differential energy distribution of dark matter in the halos obtained by the numerical simulation.From the analytical formula we reconstruct the density profile described by the Navarro,Frenk, & White (NFW) profile.The NFW profile is consistent with the analytical formula of our fractional mass distribution. We find that a parameter in our analytical formula of differential energy distribution which is related with the slope of inner cusp of dark halo.We obtain the distribution function for the NFW profile which has sharp cut off at the high binding energy.We discuss physical reason of form of the analytical formula.

Using cosmological N-body simulations of critical (SCDM) and open ($\Omega = 0.3$, OCDM) cold dark matter models, we investigate evolution of cluster galaxies. From our numerical simulation, we construct merging history trees of the galaxies. By following their merging history, we show major merger fractions of the galaxies in cluster forming regions is roughly proportional to (1+z)^{4.5} at low redshifts (z < 2) and has a steep peak at $z \simeq 2.5$ and $z \simeq 3$ in SCDM and OCDM, respectively. We also show that the cluster galaxies are affected by tidal interaction after the clusters formed. Because the formation redshift of the cluster in SCDM, z_{form} = 0.15, is much more recent than that of the cluster in OCDM, z_{form} = 1.6, the cluster galaxies in SCDM show more rapid evolution by tidal interactions from z = 0.5 than those in OCDM.

We investigate effects of time evolution of a rich cluster of galaxies on its member galactic halos in the standard cold dark matter (SCDM) universe using high resolution N-body simulations. We identify several hundred galactic halos within virial radius of our simulated cluster. We also find that a large number of halos have been tidally disrupted at z = 0. Therefore we improve a method of deriving merging history trees of galaxies taking account of tidally stripped galaxies.

In our previous paper, we have shown that a gas disk in the nuclear region of a barred galaxy which contains a central supermassive black hole (SMBH) rapidly evolves into a nuclear gas ring by the effect of an additional inner Lindblad resonance caused by the SMBH. In this paper, we investigate the fate of the gas ring, involving self-gravity of gas, using two-dimensional hydrodynamical simulations. We find that the gas ring becomes gravitationally unstable for a gas surface density of gas above a critical value, and fragments into several gas clumps. Some denser clumps increase their mass via the accretion of the surrounding gas and collisions with other clumps, and finally a very massive gas clump (10^7 M_sun) is formed. Due to the torque from the most massive clump, a part of the gas in the ring loses its angular momentum and falls into the galactic center. As a result, a nuclear gas disk (50 pc) is formed around the SMBH. The accretion rate for $R<50$ pc attains about 1 M_sun/yr for 3.5*10^7 yr. At the final phase of the bar-driven fueling, self-gravity is crucial for the angular momentum transfer of the gas. This is a new mechanism for gas fueling to the vicinity of the SMBH.

We have studied the effects of Compton cooling on cooling flow by performing numerical hydrodynamic calculations of the time evolution of hot gas in clusters of galaxies with luminous quasars. We assumed various temperatures for the hot intracluster gas. We have shown that the Compton cooling due to very luminous quasar is effective in inducing cooling flow, before radiative cooling flow is realized. The mass flux due to the cooling flow increases with time, as the Compton cooled region expands. However, the mass flux of the Compton cooling flow is not large, less than 1 M_sun/yr in our model, since Compton cooled region is limited in an inner galactic region around a quasar. Even though the quasar active phase ceased, the cooling flow will continue for at least 10^9 yr. The accreted mass is enough to explain X-ray absorption lines in high red shift quasars, if the Compton cooled gas is compressed by high pressure intracluster gas.

When a supermassive black hole exists in the centre of a galaxy, an additional inner Lindblad resonance (ILR) exists inside the usual ILRs. We study gas dynamics in a weakly barred potential with a central supermassive black hole by using 2D numerical simulations, and we investigate the effect of the additional ILR on fuelling gas into nuclear starburst regions or AGNs. Our numerical results show that strong trailing spiral shocks are formed at the resonance region, and the gas in the shock region is rapidly fuelled into a central region and make a nuclear gas ring. As a result, a large amount of gas is concentrated in the nuclear region beyond the ILR in a dynamical time scale.

We have found an effective fueling process to a central region of galaxies with weak bar-like distortion by two dimensional hydrodynamical simulations. Gravitational instability of an elongated gas ring at the inner Lindblad resonance (ILR), which has been reported as an effective fueling mechanism, are not needed for this fueling process. A massive gaseous disk in a central region of galaxies sensitively responds to the weakly distorted potential, and a large amount of gas can be fed into within $1/20$ of a core radius of the potential in several $10^7$ yr. The ILRs, the dissipative nature of the gas, and self-gravity of the gas are essential for triggering this effective fueling. The accumulation process has not been ever known: the gas accumulates to form a dense `linear' structure inclined at about 45 degree with respect to the bar potential in a {\it leading sense}. We also found that a counter rotating gaseous core can be formed as a result of the fueling. The sense of the rotation of the core depends on a fraction of the gas mass to the background mass. Physical mechanism of the fueling process can be understood using a linear theory of gaseous orbits in a weak barred potential.

We have investigated the effects of primeval starbursts on the galactic structure, showing that the dynamical structure of young galaxies is changed by relaxation of a stellar system formed by a starburst-driven outflowing gas (superwind). When a superwind interacts with halo gas, an expanding dense supershell of shocked gas is formed. Using a similarity solution for the expanding shell and a condition for the gravitational instability of a gaseous shell, we show that an expanding gaseous shell with a mass of several $10^{10} M_\odot$ and a radius of several kpc becomes gravitationally unstable to form stars. A stellar shell is thus formed, relaxes and will evolve into a fat stellar system. \par In order to investigate the fate of the stellar shell and its dynamical influence on the host galaxy, we performed three dimensional $N$-body simulations of a stellar shell plus disk system. The evolution of the stellar shell was computed under two types of environment: an external flattened potential and a live stellar disk. We find that the relaxation process and resultant structure are significantly affected by the external disk potential, and also depend on the velocity dispersion of the initial shell stars and the mass of the shell. The final galactic structure ranges from a thick disk and a compact bulge with a high density core to a thin disk with a diffuse bulge. Isodensity maps of the computed disk plus bulge system resemble isophotes of S0 galaxies. We suggest that the effects of primeval starbursts would explain some observational properties of S0 galaxies.