A consensus view on the formation of planetesimals is now exposed to a threat, since recent numerical studies on the mechanical properties of dust aggregates tend to dispute the conceptual picture that submicrometer-sized grains conglomerate into planetesimals in protoplanetary disks. With the advent of precise laboratory experiments and extensive computer simulations on the interaction between elastic spheres comprising dust aggregates, we revisit a model for the tensile strength of dust aggregates consisting of small elastic grains. In the framework of contact mechanics and fracture mechanics, we examine outcomes of computer simulations and laboratory experiments on the tensile strength of dust aggregates. We provide a novel analytical formula that explicitly incorporates the volume effect on the tensile strength, namely, the dependence of tensile strength on the volume of dust aggregates. We find that our model for the tensile strength of dust aggregates well reproduces results of computer simulations and laboratory experiments, if appropriate values are adopted for the elastic parameters used in the model. Moreover, the model with dust aggregates of submicrometer-sized grains is in good harmony with the tensile strength of cometary dust and meteoroids derived from astronomical observations. Therefore, we reaffirm the commonly believed idea that the formation of planetesimals begins with conglomeration of submicrometer-sized grains condensed in protoplanetary disks.

We present the first astrometry catalog from the Japanese VLBI (very long baseline interferometer) project VERA (VLBI Exploration of Radio Astrometry). We have compiled all the astrometry results from VERA, providing accurate trigonometric annual parallax and proper motion measurements. In total, 99 maser sources are listed in the VERA catalog. Among them, 21 maser sources are newly reported while the rest of 78 sources are referred to previously published results or those in preparation for forthcoming papers. The accuracy in the VERA astrometry are revisited and compared with those from the other VLBI astrometry projects such as BeSSeL (The Bar and Spiral Structure Legacy) Survey and GOBELINS (the Gould's Belt Distances Survey) with the VLBA (Very Long Baseline Array). We have confirmed that most of the astrometry results are consistent with each other, and the largest error sources are due to source structure of the maser features and their rapid variation, along with the systematic calibration errors and different analysis methods. Combined with the BeSSeL results, we estimate the up-to-date fundamental Galactic parameter of $R_{0}=7.92\pm0.16_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~kpc and $\Omega_{\odot}=30.17\pm0.27_{\rm{stat.}}\pm0.3_{\rm{sys.}}$~km~s$^{-1}$~kpc$^{-1}$, where $R_{0}$ and $\Omega_{\odot}$ are the distance from the Sun to the Galactic center and the Sun's angular velocity of the Galactic circular rotation, respectively.

We studied the final phases of galactic mergers, focusing on interactions between supermassive black holes (SMBHs) and the interstellar medium in a central sub-kpc region, using an N-body/hydrodynamics code. This numerical experiment aims to understand the fate of the gas supplied by mergers of two or more galaxies with SMBHs, whose masses are $10^7 M_\odot$. We observed that the mass accretion rate to one SMBH exceeds the Eddington accretion rate when the distance between two black holes (BHs) rapidly decreases. However, this rapid accretion phase does not last for more than $10^7$ yrs, and it drops to $\sim$ 10% of the Eddington rate in the quasi-steady accretion phase. The rapid accretion is caused by the angular momentum transfer from the gas to the stellar component, and the moderate accretion in the quasi-steady phase is caused {by turbulent viscosity and gravitational torque in the disk. The second merger event enhances the mass accretion to the BHs; however, this phase takes place on a similar timescale to the first merger event. We also found that the AGN feedback and the mass accretion to BHs can coexist in the central region of merged galaxies, if the amount of feedback energy is given as $(2 \times 10^{-4} - 2 \times 10^{-3} )\dot{M} c^2$, where $\dot{M}$ is the accretion rate to $r= 1$ pc. The accretion rate is suppressed by $\sim$ 1/50 in the quasi-steady accretion phase for $0.02 \dot{M} c^2$. The fraction of the gas that finally falls to each BH is approximately 5-7% of the supplied total gas mass ($10^8 M_\odot$), and 15-20% of the gas forms a circumnuclear gas inside 100 pc. This remnant gas heavily obscures the luminous phase of the active galactic nuclei (AGN) during merger events, and the moderate AGN feedback does not alter this property.

We study the obscuring structure of circumnuclear disks (CNDs) by considering supernova (SN) feedbacks from nuclear starburst and the effect of anisotropic radiative pressure from AGNs. We suppose that the mass accretion onto a central supermassive black hole (SMBH) is triggered by SN-driven turbulence within CNDs, and we explore how the structures of CNDs depend on the BH mass ($M_{\rm BH}$) and AGN luminosity ($L_{\rm AGN}$). We find that the obscuring fraction ($f_{\rm obs}$) peaks at $\sim10\%$ of the Eddington luminosity ($L_{\rm Edd}$), and its maximal value is $f_\mathrm{obs} \sim 0.6$ for less massive SMBHs (e.g., $M_{\rm BH} < 10^{8}M_{\odot}$). This is because the scale height of CNDs is determined by the SN-driven accretion for a smaller $L_{\rm AGN}$, while the dusty molecular gas in CNDs is blown away by the radiation pressure from AGNs beyond the critical luminosity. On the other hand, for massive SMBHs (e.g., $M_{\rm BH} > 10^{8}M_{\odot}$), $f_{\rm obs}$ is always smaller than $0.2$, and it is almost independent of $L_{\rm AGN}$ because the scale height of CNDs is mainly controlled by the maximal star-formation efficiency ($C_{\rm *, max}$) in CNDs. By comparison with the obscuring fractions suggested from the mid-infrared observations of nearby AGNs, the SN plus radiative feedback model with $C_{\rm *, max} = 10^{-7}\, {\rm yr}^{-1}$ well reproduces the observations for $M_\mathrm{BH} = 10^8 M_\odot$. We also find that the intense starburst or the existence of dust-free absorbers inside CNDs are necessary, to explain X-ray observations.

As a natural consequence of the elementary processes of dust growth, we discovered that a new class of planets can be formed around supermassive black holes (SMBHs). We investigated a growth path from sub-micron sized icy dust monomers to Earth-sized bodies outside the "snow line", located several parsecs from SMBHs in low luminosity active galactic nuclei (AGNs). In contrast to protoplanetary disks, the "radial drift barrier" does not prevent the formation of planetesimals. In the early phase of the evolution, low collision velocity between dust particles promotes sticking; therefore, the internal density of the dust aggregates decreases with growth. When the porous aggregate's size reaches 0.1--1 cm, the collisional compression becomes effective, and the decrease in internal density stops. Once 10--100 m sized aggregates are formed, they are decoupled from gas turbulence, and the aggregate layer becomes gravitationally unstable, leading to the formation of planets by the fragmentation of the layer, with ten times the mass of the earth. The growth time scale depends on the turbulent strength of the circumnuclear disk and the black hole mass $M_{BH}$, and it is comparable to the AGN's lifetime ($\sim 10^8$ yr) for low mass ($M_{BH} \sim 10^6 M_\odot$) SMBHs.

Active galactic nuclei (AGNs) are believed to be obscured by an optical thick "torus" that covers a large fraction of solid angles for the nuclei. However, the physical origin of the tori and the differences in the tori among AGNs are not clear. In a previous paper based on three-dimensional radiation-hydorodynamic calculations, we proposed a physics-based mechanism for the obscuration, called "radiation-driven fountains," in which the circulation of the gas driven by central radiation naturally forms a thick disk that partially obscures the nuclear emission. Here, we expand this mechanism and conduct a series of simulations to explore how obscuration depends on the properties of AGNs. We found that the obscuring fraction f_obs for a given column density toward the AGNs changes depending on both the AGN luminosity and the black hole mass. In particular, f_obs for N_H \geq 10^22 cm^-2 increases from ~0.2 to ~0.6 as a function of the X-ray luminosity L_X in the 10^{42-44} ergs/s range, but f_obs becomes small (~0.4) above a luminosity (~10^{45} ergs/s). The behaviors of f_obs can be understood by a simple analytic model and provide insight into the redshift evolution of the obscuration. The simulations also show that for a given L_AGN, f_obs is always smaller (~0.2-0.3) for a larger column density (N_H \geq 10^23 cm^-2). We also found cases that more than 70% of the solid angles can be covered by the fountain flows.

It is widely believed that in active galactic nuclei (AGNs) a supermassive black hole with an accretion disk is surrounded by an optically and geometrically thick torus at sub-parsec scale. However, it is not clear how is the mass supply toward the central engine caused and how it is related with the internal structures of the tori. The magnetic field in the tori may contribute to the accretion process via the magneto-rotational instability (MRI). Using global three dimensional magnetohydrodynamic (MHD) simulations taking the effects of X-ray heating and radiative cooling into account studied the numerical resolution for azimuthal direction for MRI driving. We found that a strongly magnetized disk consisted of a cold ($< 10^3$ K) and warm ($10^4$ K) gas is developed in about 30 rotational periods. We also found in high resolution model that the mean azimuthal magnetic fields reverse their direction quasi-periodically. We confirmed that the typical wave length of the MRI should be resolved with a least 20 azimuthal grid cells.

We used the Atacama Large Millimeter/Submillimeter Array (ALMA) to map the CO(3-2) and [CI](1-0) lines, as well as their underlying continuum emission, from the central $\sim 200$ pc region of the Circinus galaxy that hosts the nearest type 2 Seyfert-class active galactic nucleus (AGN), with a spatial resolution of $\sim 6-15$ pc. The lines and continuum-emitting regions consist of a circumnuclear disk (CND; 74 pc $\times$ 34 pc) and spiral arms. The distribution of the continuum emission revealed a temperature-dependent dust geometry and possibly polar dust elongation in the torus region. The molecular mass of the CND is $M_{\rm H2} \sim 3 \times 10^6~M_\odot$ with a beam-averaged H$_2$ column density of $\sim 5 \times 10^{23}$ cm$^{-2}$ toward the AGN position, which contributes significantly to the nuclear obscuration. The [CI](1-0)/CO(3-2) ratio at the AGN position is unusually high, suggesting an X-ray dominated region-type chemistry. We decomposed the observed velocity fields into rotational and dispersion components, and revealed multi-phase dynamic nature in the $r \lesssim 10$ pc torus region, i.e., the diffuse atomic gas is more spatially extended along the vertical direction of the disk than the dense molecular gas. Through comparisons with our model predictions based on the radiation-driven fountain scheme, we indicate that atomic outflows are the driver of the geometrical thickness of the atomic disk. This supports the validity of the radiation-driven fountain scheme in the vicinity of this AGN, which would explain the long-lasting mystery, the physical origin of the AGN torus.

Inelastic collisions occur among regolith particles, such as those in the ejecta curtain from a crater, and may cause clustering or agglomeration of particles and thus produce discrete patterns of ejecta deposits around a crater. Previous studies have shown that clusters, and even agglomerates, are formed via mutual, inelastic collisions of spherical particles due to adhering forces between particles in granular streams. To investigate the condition of agglomerate formation in granular streams, we conducted laboratory experiments of granular streams using both spherical and irregular, non-spherical particles. Measurements of particle adhesion in this study were performed using a centrifugal separation method, in contrast to the previous study in which atomic force microscopy (AFM) was used. This enabled simultaneous measurements of multiple particles of various shapes for a statistical analysis of the results. With similar relative velocities and adhesion values, irregular particles were found to form agglomerates much more easily than spherical particles. The axial ratio of the agglomerates of spherical particles and irregular particles was similar and was in accordance with those observed in previous laboratory studies, whereas the size of the agglomerates of irregular particles was larger than the size of spherical particles. The degree of agglomeration and the size of agglomerates can be used as an indicator of the shape or adhesive force of the particles in granular stream. Our findings on agglomeration in granular streams could provide new insights into the origin of rays on airless bodies and grooves on Phobos.

We investigate the properties of the ionized gas irradiated by an active galactic nucleus (AGN) based on our "radiation-driven fountain" model for the nearest type-2 Seyfert galaxy, the Circinus galaxy (Wada et al. 2016). We conducted "quasi-three dimensional" spectral analysis using the spectral synthesis code CLOUDY and obtained the surface brightness distributions of lines, such as H$\alpha$, H$\beta$, [OIII], [NII], and [SII] for the central 16-parsec region. The ionized regions observed based on these lines show a conical morphology around the rotation axis, even if we do not phenomenologically postulate the presence of an optically thick "torus". This region also shows non-uniform internal structures, reflecting the inhomogeneous structure of fountain flows. Using ionization diagnostic diagrams, we investigated the spectral properties of the ionized gas. The diagrams based on the line ratios of [NII]/H$\alpha$ and [SII]/H$\alpha$ show that most regions of the cone have the same properties as those in the narrow line regions (NLRs) in AGNs, whereas using [OI]/H$\alpha$, the central 10-pc regions are rather LINER-like. The gas density, temperature, and ionizing parameter in regions identified as "NLR" are typically $n \sim 300-1500$ cm$^{-3}$, $T \sim 1-3 \times 10^4 $ K, and $ U \sim 0.01$, respectively. The morphology and [OIII] intensity are similar to the base of the observed [OIII] cone in the Circinus galaxy, implying some physical connections with the origin of the $\sim100$ parsec scale NLR.

Asteroid (162173) Ryugu is the target object of Hayabusa2, an asteroid exploration and sample return mission led by Japan Aerospace Exploration Agency (JAXA). Ground-based observations indicate that Ryugu is a C-type near-Earth asteroid with a diameter of less than 1 km, but the knowledge of its detailed properties is still very limited. This paper summarizes our best understanding of the physical and dynamical properties of Ryugu based on remote sensing and theoretical modeling. This information is used to construct a design reference model of the asteroid that is used for formulation of mission operations plans in advance of asteroid arrival. Particular attention is given to the surface properties of Ryugu that are relevant to sample acquisition. This reference model helps readers to appropriately interpret the data that will be directly obtained by Hayabusa2 and promotes scientific studies not only for Ryugu itself and other small bodies but also for the Solar System evolution that small bodies shed light on.

We present the results of our ALMA Cycle 4 high-spatial-resolution (0.04-0.07") observations, at HCN J=3-2 and HCO+ J=3-2 lines, of the nucleus of NGC 1068, the nearby prototypical type 2 active galactic nucleus (AGN). Our previous ALMA observations identified the compact emission of these lines at the putative location of the torus around a mass-accreting supermassive black hole. We now report that we have detected the rotation of this compact emission, with the eastern and western sides being redshifted and blueshifted, respectively. Unlike the previously reported CO J=6-5 emission, both the morphological and dynamical alignments of the HCN J=3-2 and HCO+ J=3-2 emission are roughly aligned along the east-west direction (i.e., the expected torus direction), suggesting that these molecular lines are better probes of a rotating dense molecular gas component in the torus. The western part of the torus exhibits larger velocity dispersion and stronger emission in the HCN J=3-2 and HCO+ J=3-2 lines than the eastern part, revealing a highly inhomogeneous molecular torus. The dense molecular gas in the torus and that of the host galaxy at 0.5-2.0" from the AGN along the torus direction are found to be counter-rotating, suggesting an external process happened in the past at the NGC 1068 nucleus.

In this study, we investigate the line emissions from cold molecular gas based on our previous "radiation-driven fountain model" (Wada 2016), which reliably explains the spectral energy distribution of the nearest type 2 Seyfert galaxy, the Circinus galaxy. Using a snapshot of the best-fit radiation-hydrodynamic model for the central r < 16 pc, in which non-equilibrium X-ray-dominated region chemistry is solved, we conduct post-processed, non-local thermodynamic equilibrium radiation transfer simulations for the CO lines. We obtain a spectral line energy distribution with a peak around J=6, and its distribution suggests that the lines are not thermalized. However, for a given line-of-sight, the optical depth distribution is highly non-uniform between $\tau \ll 1$ and $\tau \gg 1$. The CO-to-H2 conversion factor (X_CO), which can be directly obtained from the results, is not a constant and depends strongly on the integrated intensity, and it differs from the fiducial value for local objects. X_CO exhibits a large dispersion of more than one order of magnitude, reflecting the non-uniform internal structure of a "torus." We also found that the physical conditions differ between grid cells on a scale of a few parsecs along the observed lines of sight; therefore, a specific observed line ratio does not necessarily represent a single physical state of the ISM.

On both observational and theoretical grounds, the disk of our Galaxy should be accreting cool gas with temperature ~<10^5 K via the halo at a rate ~1 M_sun/yr. At least some of this accretion is mediated by high velocity clouds (HVCs), observed to be traveling in the halo with velocities of a few 100 km/s and occasionally impacting the disk at such velocities, especially in the outer regions of the Galaxy. We address the possibility of particle acceleration in shocks triggered by such HVC accretion events, and the detectability of consequent non-thermal emission in the radio to gamma-ray bands and high-energy neutrinos. For plausible shock velocities ~300 km/s and magnetic field strengths ~0.3-10 muG, electrons and protons may be accelerated up to ~1-10 TeV and ~30-10^3 TeV, respectively, in sufficiently strong adiabatic shocks during their lifetime of ~10^6 yr. The resultant pion decay and inverse Compton gamma-rays may be the origin of some unidentified Galactic GeV-TeV sources, particularly the "dark" source HESS J1503-582 that is spatially coincident with anomalous HI structure known as "forbidden-velocity wings". Correlation of their locations with star-forming regions may be weak, absent or even opposite. Non-thermal radio and X-ray emission from primary and/or secondary electrons may be detectable with deeper observations. The contribution of HVC accretion to Galactic cosmic rays is subdominant, but could be non-negligible in the outer Galaxy. As the thermal emission induced by HVC accretion is likely difficult to detect, observations of such phenomena may offer a unique perspective on probing gas accretion onto the Milky Way and other galaxies.

High-resolution observations from the sub-mm to the optical wavelength regime resolve the central few 100pc region of nearby galaxies in great detail. They reveal a large diversity of features: thick gas and stellar discs, nuclear starbursts, in- and outflows, central activity, jet interaction, etc. Concentrating on the role circumnuclear discs play in the life cycles of galactic nuclei, we employ 3D adaptive mesh refinement hydrodynamical simulations with the RAMSES code to self-consistently trace the evolution from a quasi-stable gas disc, undergoing gravitational (Toomre) instability, the formation of clumps and stars and the disc's subsequent, partial dispersal via stellar feedback. Our approach builds upon the observational finding that many nearby Seyfert galaxies have undergone intense nuclear starbursts in their recent past and in many nearby sources star formation is concentrated in a handful of clumps on a few 100pc distant from the galactic centre. We show that such observations can be understood as the result of gravitational instabilities in dense circumnuclear discs. By comparing these simulations to available integral field unit observations of a sample of nearby galactic nuclei, we find consistent gas and stellar masses, kinematics, star formation and outflow properties. Important ingredients in the simulations are the self-consistent treatment of star formation and the dynamical evolution of the stellar distribution as well as the modelling of a delay time distribution for the supernova feedback. The knowledge of the resulting simulated density structure and kinematics on pc scale is vital for understanding inflow and feedback processes towards galactic scales.

We used the Atacama Large Millimeter/submillimeter Array (ALMA) to map the CO(3-2) and the underlying continuum emissions around the type 1 low-luminosity active galactic nucleus (LLAGN; bolometric luminosity $\lesssim 10^{42}$ erg~s$^{-1}$) of NGC 1097 at $\sim 10$ pc resolution. These observations revealed a detailed cold gas distribution within a $\sim 100$ pc of this LLAGN. In contrast to the luminous Seyfert galaxy NGC 1068, where a $\sim 7$ pc cold molecular torus was recently revealed, a distinctively dense and compact torus is missing in our CO(3-2) integrated intensity map of NGC 1097. Based on the CO(3-2) flux, the gas mass of the torus of NGC 1097 would be a factor of $\gtrsim 2-3$ less than that found for NGC 1068 by using the same CO-to-H$_2$ conversion factor, which implies less active nuclear star formation and/or inflows in NGC 1097. Our dynamical modeling of the CO(3-2) velocity field implies that the cold molecular gas is concentrated in a thin layer as compared to the hot gas traced by the 2.12 $\mu$m H$_2$ emission in and around the torus. Furthermore, we suggest that NGC 1097 hosts a geometrically thinner torus than NGC 1068. Although the physical origin of the torus thickness remains unclear, our observations support a theoretical prediction that geometrically thick tori with high opacity will become deficient as AGNs evolve from luminous Seyferts to LLAGNs.

We have constructed the database of stars in the local group using the extended version of the SAGA (Stellar Abundances for Galactic Archaeology) database that contains stars in 24 dwarf spheroidal galaxies and ultra faint dwarfs. The new version of the database includes more than 4500 stars in the Milky Way, by removing the previous metallicity criterion of [Fe/H] <= -2.5, and more than 6000 stars in the local group galaxies. We examined a validity of using a combined data set for elemental abundances. We also checked a consistency between the derived distances to individual stars and those to galaxies in the literature values. Using the updated database, the characteristics of stars in dwarf galaxies are discussed. Our statistical analyses of alpha-element abundances show that the change of the slope of the [alpha/Fe] relative to [Fe/H] (so-called "knee") occurs at [Fe/H] = -1.0+-0.1 for the Milky Way. The knee positions for selected galaxies are derived by applying the same method. Star formation history of individual galaxies are explored using the slope of the cumulative metallicity distribution function. Radial gradients along the four directions are inspected in six galaxies where we find no direction dependence of metallicity gradients along the major and minor axes. The compilation of all the available data shows a lack of CEMP-s population in dwarf galaxies, while there may be some CEMP-no stars at [Fe/H] <~ -3 even in the very small sample. The inspection of the relationship between Eu and Ba abundances confirms an anomalously Ba-rich population in Fornax, which indicates a pre-enrichment of interstellar gas with r-process elements. We do not find any evidence of anti-correlations in O-Na and Mg-Al abundances, which characterises the abundance trends in the Galactic globular clusters.

High resolution submm observations are important in probing the morphology, column density and dynamics of obscured active galactic nuclei (AGNs). With high resolution (0.06 x 0.05) ALMA 690 GHz observations we have found bright (TB >80 K) and compact (FWHM 10x7 pc) CO 6-5 line emission in the nucleus of the extremely radio-quiet galaxy NGC1377. The integrated CO 6-5 intensity is aligned with the previously discovered jet/outflow of NGC1377 and is tracing the dense (n>1e4 cm-3), hot gas at the base of the outflow. The velocity structure is complex and shifts across the jet/outflow are discussed in terms of jet-rotation or separate, overlapping kinematical components. High velocity gas (deltaV +-145 km/s) is detected inside r<2-3 pc and we suggest that it is emerging from an inclined rotating disk or torus of position angle PA=140+-20 deg with a dynamical mass of approx 3e6 Msun. This mass is consistent with that of a supermassive black hole (SMBH), as inferred from the M-sigma relation. The gas mass of the proposed disk/torus constitutes <3% of the nuclear dynamical mass. In contrast to the intense CO 6-5 line emission, we do not detect dust continuum with an upper limit of S(690GHz)<2mJy. The corresponding, 5 pc, H2 column density is estimated to N(H2)<3e23 cm-2, which is inconsistent with a Compton Thick (CT) source. We discuss the possibility that CT obscuration may be occuring on small (subparsec) or larger scales. From SED fitting we suggest that half of the IR emission of NGC1377 is nuclear and the rest (mostly the far-infrared (FIR)) is more extended. The extreme radio quietness, and the lack of emission from other star formation tracers, raise questions on the origin of the FIR emission. We discuss the possibility that it is arising from the dissipation of shocks in the molecular jet/outflow or from irradiation by the nuclear source along the poles.

The structures and dynamics of molecular, atomic, and ionized gases are studied around a low-luminosity active galactic nucleus (AGN) with a small ($2\times 10^6 M_\odot$) black hole using 3D radiation hydrodynamic simulations. We studied, for the first time, the non-equilibrium chemistry for the X-ray dominated region in the "radiation-driven fountain" (Wada 2012) with supernova feedback. A double hollow cone structure is naturally formed without postulating a thick "torus" around a central source. The cone is occupied with an inhomogeneous, diffuse ionized gas and surrounded by a geometrically thick ($h/r \gtrsim 1$) atomic gas. Dense molecular gases are distributed near the equatorial plane, and energy feedback from supernovae enhances their scale height. Molecular hydrogen exists in a hot phase ( > 1000 K) as well as in a cold ( < 100 K), dense ( >$10^3$ cm$^{-3}$) phase. The velocity dispersion of H$_2$ in the vertical direction is comparable to the rotational velocity, which is consistent with near infrared observations of nearby Seyfert galaxies. Using 3D radiation transfer calculations for the dust emission, we find polar emission in the mid-infrared band (12$\mu m$), which is associated with bipolar outflows, as suggested in recent interferometric observations of nearby AGNs. If the viewing angle for the nucleus is larger than 75 deg, the spectral energy distribution (~ 2 -- 60 $\mu m$) of this model is consistent with that of the Circinus galaxy. The multi-phase interstellar medium observed in optical/infrared and X-ray observations is also discussed.

With high resolution (0."25 x 0."18) ALMA CO 3-2 observations of the nearby (D=21 Mpc, 1"=102 pc), extremely radio-quiet galaxy NGC1377, we have discovered a high-velocity, very collimated nuclear outflow which we interpret as a molecular jet with a projected length of +-150 pc. Along the jet axis we find strong velocity reversals where the projected velocity swings from -150 km/s to +150 km/s. A simple model of a molecular jet precessing around an axis close to the plane of the sky can reproduce the observations. The velocity of the outflowing gas is difficult to constrain due to the velocity reversals but we estimate it to be between 240 and 850 km/s and the jet to precess with a period P=0.3-1.1 Myr. The CO emission is clumpy along the jet and the total molecular mass in the high-velocity (+-(60 to 150 km/s)) gas lies between 2e6 Msun (light jet) and 2e7 Msun (massive jet). There is also CO emission extending along the minor axis of NGC1377. It holds >40% of the flux in NGC1377 and may be a slower, wide-angle molecular outflow which is partially entrained by the molecular jet. We discuss the driving mechanism of the molecular jet and suggest that it is either powered by a very faint radio jet or by an accretion disk-wind similar to those found towards protostars. The nucleus of NGC1377 harbours intense embedded activity and we detect emission from vibrationally excited HCN J=4-3 v_2=1f which is consistent with hot gas and dust. We find large columns of H2 in the centre of NGC1377 which may be a sign of a high rate of recent gas infall. The dynamical age of the molecular jet is short (<1 Myr), which could imply that it is young and consistent with the notion that NGC1377 is caught in a transient phase of its evolution. However, further studies are required to determine the age of the molecular jet, its mass and the role it is playing in the growth of the nucleus of NGC1377.

Adhesion forces between submicrometer-sized silicate grains play a crucial role in the formation of silicate dust agglomerates, rocky planetesimals, and terrestrial planets. The surface energy of silicate dust particles is the key to their adhesion and rolling forces in a theoretical model based on the contact mechanics. Here we revisit the cohesion of amorphous silica spheres by compiling available data on the surface energy for hydrophilic amorphous silica in various circumstances. It turned out that the surface energy for hydrophilic amorphous silica in a vacuum is a factor of 10 higher than previously assumed. Therefore, the previous theoretical models underestimated the critical velocity for the sticking of amorphous silica spheres, as well as the rolling friction forces between them. With the most plausible value of the surface energy for amorphous silica spheres, theoretical models based on the contact mechanics are in harmony with laboratory experiments. Consequently, we conclude that silicate grains with a radius of $0.1~\mu$m could grow to planetesimals via coagulation in a protoplanetary disk. We argue that the coagulation growth of silicate grains in a molecular cloud is advanced either by organic mantles rather than icy mantles or, if there are no mantles, by nanometer-sized grain radius.

Using two-dimensional radiation hydrodynamic simulations, we investigate origin of the ultra fast outflows (UFOs) that are often observed in luminous active galactic nuclei (AGNs). We found that the radiation force due to the spectral lines generates strong winds (line-driven disk winds) that are launched from the inner region of accretion disks (~30 Schwarzschild radii). A wide range of black hole masses ($M_{\rm BH}$) and Eddington ratios ($\varepsilon$) was investigated to study conditions for causing the line-driven winds. For $M_{\rm BH} = 10^6-10^9 M_\odot$ and $\varepsilon = 0.1-0.7$, funnel-shaped disk winds appear, in which dense matter is accelerated outward with an opening angle of 70-80 deg and with 10% of the light speed. If we observe the wind along its direction, the velocity, the column density, and the ionization state are consistent with those of the observed UFOs. As long as the obscuration by the torus does not affect the observations of X-ray bands, the UFOs could be statistically observed in about 13-28% of the luminous AGNs, which is not inconsistent with the observed ratio (~40%). We also found that the results are insensitive to the X-ray luminosity and the density of the disk surface. Thus, we can conclude that the UFOs could exist in any luminous AGNs, such as narrow-line Seyfert 1s (NLS1s) and quasars with $\varepsilon > 0.1$, in which fast line-driven winds are associated.

We reanalyze microlensing events in the published list of anomalous events that were observed from the OGLE lensing survey conducted during 2004-2008 period. In order to check the existence of possible degenerate solutions and extract extra information, we conduct analyses based on combined data from other survey and follow-up observation and consider higher-order effects. Among the analyzed events, we present analyses of 8 events for which either new solutions are identified or additional information is obtained. We find that the previous binary-source interpretations of 5 events are better interpreted by binary-lens models. These events include OGLE-2006-BLG-238, OGLE-2007-BLG-159, OGLE-2007-BLG-491, OGLE-2008-BLG-143, and OGLE-2008-BLG-210. With additional data covering caustic crossings, we detect finite-source effects for 6 events including OGLE-2006-BLG-215, OGLE-2006-BLG-238, OGLE-2006-BLG-450, OGLE-2008-BLG-143, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. Among them, we are able to measure the Einstein radii of 3 events for which multi-band data are available. These events are OGLE-2006-BLG-238, OGLE-2008-BLG-210, and OGLE-2008-BLG-513. For OGLE-2008-BLG-143, we detect higher-order effect induced by the changes of the observer's position caused by the orbital motion of the Earth around the Sun. In addition, we present degenerate solutions resulting from the known close/wide or ecliptic degeneracy. Finally, we note that the masses of the binary companions of the lenses of OGLE-2006-BLG-450 and OGLE-2008-BLG-210 are in the brown-dwarf regime.

We present the analysis of the gravitational microlensing event OGLE-2013-BLG-0102. The light curve of the event is characterized by a strong short-term anomaly superposed on a smoothly varying lensing curve with a moderate magnification $A_{\rm max}\sim 1.5$. It is found that the event was produced by a binary lens with a mass ratio between the components of $q = 0.13$ and the anomaly was caused by the passage of the source trajectory over a caustic located away from the barycenter of the binary. From the analysis of the effects on the light curve due to the finite size of the source and the parallactic motion of the Earth, the physical parameters of the lens system are determined. The measured masses of the lens components are $M_{1} = 0.096 \pm 0.013~M_{\odot}$ and $M_{2} = 0.012 \pm 0.002~M_{\odot}$, which correspond to near the hydrogen-burning and deuterium-burning mass limits, respectively. The distance to the lens is $3.04 \pm 0.31~{\rm kpc}$ and the projected separation between the lens components is $0.80 \pm 0.08~{\rm AU}$.

The central engines of Seyfert galaxies are thought to be enshrouded by geometrically thick gas and dust structures. In this article, we derive observable properties for a self-consistent model of such toroidal gas and dust distributions, where the geometrical thickness is achieved and maintained with the help of X-ray heating and radiation pressure due to the central engine. Spectral energy distributions (SEDs) and images are obtained with the help of dust continuum radiative transfer calculations with RADMC-3D. For the first time, we are able to present time-resolved SEDs and images for a physical model of the central obscurer. Temporal changes are mostly visible at shorter wavelengths, close to the combined peak of the dust opacity as well as the central source spectrum and are caused by variations in the column densities of the generated outflow. Due to the three-component morphology of the hydrodynamical models -- a thin disc with high density filaments, a surrounding fluffy component (the obscurer) and a low density outflow along the rotation axis -- we find dramatic differences depending on wavelength: whereas the mid-infrared images are dominated by the elongated appearance of the outflow cone, the long wavelength emission is mainly given by the cold and dense disc component. Overall, we find good agreement with observed characteristics, especially for those models, which show clear outflow cones in combination with a geometrically thick distribution of gas and dust, as well as a geometrically thin, but high column density disc in the equatorial plane.

We investigated the magnitude-phase relation of (162173) 1999 JU3, a target asteroid for the JAXA Hayabusa 2 sample return mission. We initially employed the international Astronomical Union's H-G formalism but found that it fits less well using a single set of parameters. To improve the inadequate fit, we employed two photometric functions, the Shevchenko and Hapke functions. With the Shevchenko function, we found that the magnitude-phase relation exhibits linear behavior in a wide phase angle range (alpha = 5-75 deg) and shows weak nonlinear opposition brightening at alpha< 5 deg, providing a more reliable absolute magnitude of Hv = 19.25 +- 0.03. The phase slope (0.039 +- 0.001 mag/deg) and opposition effect amplitude (parameterized by the ratio of intensity at alpha=0.3 deg to that at alpha=5 deg, I(0.3)/I(5)=1.31+-0.05) are consistent with those of typical C-type asteroids. We also attempted to determine the parameters for the Hapke model, which are applicable for constructing the surface reflectance map with the Hayabusa 2 onboard cameras. Although we could not constrain the full set of Hapke parameters, we obtained possible values, w=0.041, g=-0.38, B0=1.43, and h=0.050, assuming a surface roughness parameter theta=20 deg. By combining our photometric study with a thermal model of the asteroid (Mueller et al. in preparation), we obtained a geometric albedo of pv = 0.047 +- 0.003, phase integral q = 0.32 +- 0.03, and Bond albedo AB = 0.014 +- 0.002, which are commensurate with the values for common C-type asteroids.

We detect a cold, terrestrial planet in a binary-star system using gravitational microlensing. The planet has low mass (2 Earth masses) and lies projected at $a_{\perp,ph}$ ~ 0.8 astronomical units (AU) from its host star, similar to the Earth-Sun distance. However, the planet temperature is much lower, T<60 Kelvin, because the host star is only 0.10--0.15 solar masses and therefore more than 400 times less luminous than the Sun. The host is itself orbiting a slightly more massive companion with projected separation $a_{\perp,ch}=$10--15 AU. Straightforward modification of current microlensing search strategies could increase their sensitivity to planets in binary systems. With more detections, such binary-star/planetary systems could place constraints on models of planet formation and evolution. This detection is consistent with such systems being very common.

We report the discovery of MOA-2013-BLG-220Lb, which has a super-Jupiter mass ratio $q=3.01\pm 0.02\times 10^{-3}$ relative to its host. The proper motion, $\mu=12.5\pm 1\, {\rm mas}\,{\rm yr}^{-1}$, is one of the highest for microlensing planets yet discovered, implying that it will be possible to separately resolve the host within $\sim 7$ years. Two separate lines of evidence imply that the planet and host are in the Galactic disk. The planet could have been detected and characterized purely with follow-up data, which has important implications for microlensing surveys, both current and into the LSST era.

Characterizing a microlensing planet is done from modeling an observed lensing light curve. In this process, it is often confronted that solutions of different lensing parameters result in similar light curves, causing difficulties in uniquely interpreting the lens system, and thus understanding the causes of different types of degeneracy is important. In this work, we show that incomplete coverage of a planetary perturbation can result in degenerate solutions even for events where the planetary signal is detected with a high level of statistical significance. We demonstrate the degeneracy for an actually observed event OGLE-2012-BLG-0455/MOA-2012-BLG-206. The peak of this high-magnification event $(A_{\rm max}\sim400)$ exhibits very strong deviation from a point-lens model with $\Delta\chi^{2}\gtrsim4000$ for data sets with a total number of measurement 6963. From detailed modeling of the light curve, we find that the deviation can be explained by four distinct solutions, i.e., two very different sets of solutions, each with a two-fold degeneracy. While the two-fold (so-called "close/wide") degeneracy is well-understood, the degeneracy between the radically different solutions is not previously known. The model light curves of this degeneracy differ substantially in the parts that were not covered by observation, indicating that the degeneracy is caused by the incomplete coverage of the perturbation. It is expected that the frequency of the degeneracy introduced in this work will be greatly reduced with the improvement of the current lensing survey and follow-up experiments and the advent of new surveys.

We report the discovery of a massive planet OGLE-2008-BLG-355Lb. The light curve analysis indicates a planet:host mass ratio of q = 0.0118 +/- 0.0006 at a separation of 0.877 +/- 0.010 Einstein radii. We do not measure a significant microlensing parallax signal and do not have high angular resolution images that could detect the planetary host star. Therefore, we do not have a direct measurement of the host star mass. A Bayesian analysis, assuming that all host stars have equal probability to host a planet with the measured mass ratio implies a host star mass of M_h = 0.37_{-0.17}^{+0.30} M_Sun and a companion of mass M_P = 4.6^{+3.7}_{-2.2} M_Jup, at a projected separation of r_proj = 1.70^{+0.29}_{-0.30} AU. The implied distance to the planetary system is D_L = 6.8 +/- 1.1 kpc. A planetary system with the properties preferred by the Bayesian analysis would be a challenge to the core-accretion model of planet formation, as the core-accretion model predicts that massive planets are far more likely to form around more massive host stars. This core accretion model prediction is not consistent with our Bayesian prior of an equal probability of host stars of all masses to host a planet with the measured mass ratio. So, if the core accretion model prediction is right, we should expect that follow-up high angular resolution observations will detect a host star with a mass in the upper part of the range allowed by the Bayesian analysis. That is, the host would probably be a K or G dwarf.

The mass of the lenses giving rise to Galactic microlensing events can be constrained by measuring the relative lens-source proper motion and lens flux. The flux of the lens can be separated from that of the source, companions to the source, and unrelated nearby stars with high-resolution images taken when the lens and source are spatially resolved. For typical ground-based adaptive optics (AO) or space-based observations, this requires either inordinately long time baselines or high relative proper motions. We provide a list of microlensing events toward the Galactic Bulge with high relative lens-source proper motion that are therefore good candidates for constraining the lens mass with future high-resolution imaging. We investigate all events from 2004 -- 2013 that display detectable finite-source effects, a feature that allows us to measure the proper motion. In total, we present 20 events with mu >~ 8 mas/yr. Of these, 14 were culled from previous analyses while 6 are new, including OGLE-2004-BLG-368, MOA-2005-BLG-36, OGLE-2012-BLG-0211, OGLE-2012-BLG-0456, MOA-2012-BLG-532, and MOA-2013-BLG-029. In <~12 years the lens and source of each event will be sufficiently separated for ground-based telescopes with AO systems or space telescopes to resolve each component and further characterize the lens system. Furthermore, for the most recent events, comparison of the lens flux estimates from images taken immediately to those estimated from images taken when the lens and source are resolved can be used to empirically check the robustness of the single-epoch method currently being used to estimate lens masses for many events.

Global "second-generation" microlensing surveys aim to discover and characterize extrasolar planets and their frequency, by means of round-the-clock high-cadence monitoring of a large area of the Galactic bulge, in a controlled experiment. We report the discovery of a giant planet in microlensing event MOA-2011-BLG-322. This moderate-magnification event, which displays a clear anomaly induced by a second lensing mass, was inside the footprint of our second-generation microlensing survey, involving MOA, OGLE and the Wise Observatory. The event was observed by the survey groups, without prompting alerts that could have led to dedicated follow-up observations. Fitting a microlensing model to the data, we find that the timescale of the event was t_E=23.2 +/-0.8 days, and the mass ratio between the lens star and its companion is q=0.028 +/-0.001. Finite-source effects are marginally detected, and upper limits on them help break some of the degeneracy in the system parameters. Using a Bayesian analysis that incorporates a Galactic structure model, we estimate the mass of the lens at 0.39 +0.45/-0.19 M_sun, at a distance of 7.56 +/-0.91 kpc. Thus, the companion is likely a planet of mass 11.6 +13.4/-5.6 M_J, at a projected separation of 4.3 +1.5/-1.2 AU, rather far beyond the snow line. This is the first pure-survey planet reported from a second-generation microlensing survey, and shows that survey data alone can be sufficient to characterize a planetary model. With the detection of additional survey-only planets, we will be able to constrain the frequency of extrasolar planets near their systems' snow lines.

We report analysis of high microlensing event MOA-2008-BLG-379, which has a strong microlensing anomaly at its peak, due to a massive planet with a mass ratio of q = 6.9 x 10^{-3}. Because the faint source star crosses the large resonant caustic, the planetary signal dominates the light curve. This is unusual for planetary microlensing events, and as a result, the planetary nature of this light curve was not immediately noticed. The planetary nature of the event was found when the MOA Collaboration conducted a systematic study of binary microlensing events previously identified by the MOA alert system. We have conducted a Bayesian analysis based on a standard Galactic model to estimate the physical parameters of the lens system. This yields a host star mass of M_L = 0.66_{-0.33}^{+0.29} M_Sun orbited by a planet of mass m_P = 4.8_{-2.4}^{+2.1} M_Jup at an orbital separation of a = 4.1_{-1.5}^{+1.9} AU at a distance of D_L = 3.6 +/- 1.3 kpc. The faint source magnitude of I_S = 21.30 and relatively high lens-source relative proper motion of mu_rel = 7.6 +/- 1.6 mas/yr implies that high angular resolution adaptive optics or Hubble Space Telescope observations are likely to be able to detect the source star, which would determine the masses and distance of the planet and its host star.

We present the first microlensing candidate for a free-floating exoplanet-exomoon system, MOA-2011-BLG-262, with a primary lens mass of M_host ~ 4 Jupiter masses hosting a sub-Earth mass moon. The data are well fit by this exomoon model, but an alternate star+planet model fits the data almost as well. Nevertheless, these results indicate the potential of microlensing to detect exomoons, albeit ones that are different from the giant planet moons in our solar system. The argument for an exomoon hinges on the system being relatively close to the Sun. The data constrain the product M pi_rel, where M is the lens system mass and pi_rel is the lens-source relative parallax. If the lens system is nearby (large pi_rel), then M is small (a few Jupiter masses) and the companion is a sub-Earth-mass exomoon. The best-fit solution has a large lens-source relative proper motion, mu_rel = 19.6 +- 1.6 mas/yr, which would rule out a distant lens system unless the source star has an unusually high proper motion. However, data from the OGLE collaboration nearly rule out a high source proper motion, so the exoplanet+exomoon model is the favored interpretation for the best fit model. However, the alternate solution has a lower proper motion, which is compatible with a distant (so stellar) host. A Bayesian analysis does not favor the exoplanet+exomoon interpretation, so Occam's razor favors a lens system in the bulge with host and companion masses of M_host = 0.12 (+0.19 -0.06) M_solar and m_comp = 18 (+28 -100 M_earth, at a projected separation of a_perp ~ 0.84 AU. The existence of this degeneracy is an unlucky accident, so current microlensing experiments are in principle sensitive to exomoons. In some circumstances, it will be possible to definitively establish the low mass of such lens systems through the microlensing parallax effect. Future experiments will be sensitive to less extreme exomoons.

We present a detailed analysis of survey and follow-up observations of microlensing event OGLE-2012-BLG-0406 based on data obtained from 10 different observatories. Intensive coverage of the lightcurve, especially the perturbation part, allowed us to accurately measure the parallax effect and lens orbital motion. Combining our measurement of the lens parallax with the angular Einstein radius determined from finite-source effects, we estimate the physical parameters of the lens system. We find that the event was caused by a $2.73\pm 0.43\ M_{\rm J}$ planet orbiting a $0.44\pm 0.07\ M_{\odot}$ early M-type star. The distance to the lens is $4.97\pm 0.29$\ kpc and the projected separation between the host star and its planet at the time of the event is $3.45\pm 0.26$ AU. We find that the additional coverage provided by follow-up observations, especially during the planetary perturbation, leads to a more accurate determination of the physical parameters of the lens.

We analyze MOA-2010-BLG-311, a high magnification (A_max>600) microlensing event with complete data coverage over the peak, making it very sensitive to planetary signals. We fit this event with both a point lens and a 2-body lens model and find that the 2-body lens model is a better fit but with only Delta chi^2~80. The preferred mass ratio between the lens star and its companion is $q=10^(-3.7+/-0.1), placing the candidate companion in the planetary regime. Despite the formal significance of the planet, we show that because of systematics in the data the evidence for a planetary companion to the lens is too tenuous to claim a secure detection. When combined with analyses of other high-magnification events, this event helps empirically define the threshold for reliable planet detection in high-magnification events, which remains an open question.

We analyze the planetary microlensing event MOA-2010-BLG-328. The best fit yields host and planetary masses of Mh = 0.11+/-0.01 M_{sun} and Mp = 9.2+/-2.2M_Earth, corresponding to a very late M dwarf and sub-Neptune-mass planet, respectively. The system lies at DL = 0.81 +/- 0.10 kpc with projected separation r = 0.92 +/- 0.16 AU. Because of the host's a-priori-unlikely close distance, as well as the unusual nature of the system, we consider the possibility that the microlens parallax signal, which determines the host mass and distance, is actually due to xallarap (source orbital motion) that is being misinterpreted as parallax. We show a result that favors the parallax solution, even given its close host distance. We show that future high-resolution astrometric measurements could decisively resolve the remaining ambiguity of these solutions.

We present measurements of the microlensing optical depth and event rate toward the Galactic Bulge based on two years of the MOA-II survey. This sample contains ~1000 microlensing events, with an Einstein Radius crossing time of t_E < 200 days between -5 <l< 10 degree and -7 <b< -1 degree. Our event rate and optical depth analysis uses 474 events with well defined microlensing parameters. In the central fields with |l|< 5 degree, we find an event rates of \Gamma = [2.39+/-1.1]e^{[0.60\pm0.05](3-|b|)}x 10^{-5}/star/yr and an optical depth of \tau_{200} = [2.35+/-0.18]e^{[0.51+/-0.07](3-|b|)}x 10^{-6} for the 427 events using all sources brighter than I_s = 20 mag centered at (l,b)=(0.38, -3.72). We find that the event rate is maximized at low latitudes and a longitude of $l~1 degree. For the 111 events in 3.2 deg^2 of the central Galactic Bulge at |b| < 3.0 degree and 0.0 < l < 2.0, centered at (l,b)=(0.97, -2.26), we find \Gamma = 4.57_{-0.46}^{+0.51} x 10^{-5}/star/yr and \tau_{200} = 3.64_{-0.45}^{+ 0.51} x 10^{-6}. We also consider a Red Clump Giant (RCG) star sample with I_s<17.5 mag. Our results are consistent with previous optical depth measurements. We find that the previously observed difference in optical depth measurements between all-source and RCG samples may be largely due to statistical fluctuations. These event rate measurements towards the central galactic bulge are necessary to predict the microlensing event rate and to optimize the survey fields in the future space mission such as WFIRST.

Several lines of argument support the existence of a link between activity at the nuclei of galaxies, in the form of an accreting supermassive black hole, and star-formation activity in these galaxies. The exact nature of this link is still under debate. Radio jets have long been argued to be an ideal mechanism that allows AGN to interact with their host galaxy and regulate star-formation. In this context, we are using a sample of radio sources in the North Ecliptic Pole (NEP) field to study the nature of the putative link between AGN activity and star-formation. This is done by means of spectral energy distribution (SED) fitting. We use the excellent spectral coverage of the AKARI infrared space telescope together with the rich ancillary data available in the NEP to build SEDs extending from UV to far-IR wavelengths. Through SED fitting we constrain both the AGN and host galaxy components. We find a significant AGN component in our sample of relatively faint radio-sources ($<$mJy), that increases in power with increasing radio-luminosity. At the highest radio-luminosities, the presence of powerful jets dominates the radio emission of these sources. A positive correlation is found between the luminosity of the AGN component and that of star-formation in the host galaxy, independent of the radio luminosity. By contrast, for a given redshift and AGN luminosity, we find that increasing radio-luminosity leads to a decrease in the specific star-formation rate. The most radio-loud AGN are found to lie on the main sequence of star-formation for their respective redshifts. For the first time, such a two-sided feedback process is seen in the same sample. We conclude that radio jets do suppress star-formation in their host galaxies but appear not to totally quench it. Our results therefore support the maintenance nature of "radio-mode" feedback from radio-AGN jets.

A planetary microlensing signal is generally characterized by a short-term perturbation to the standard single lensing light curve. A subset of binary-source events can produce perturbations that mimic planetary signals, thereby introducing an ambiguity between the planetary and binary-source interpretations. In this paper, we present analysis of the microlensing event MOA-2012-BLG-486, for which the light curve exhibits a short-lived perturbation. Routine modeling not considering data taken in different passbands yields a best-fit planetary model that is slightly preferred over the best-fit binary-source model. However, when allowed for a change in the color during the perturbation, we find that the binary-source model yields a significantly better fit and thus the degeneracy is clearly resolved. This event not only signifies the importance of considering various interpretations of short-term anomalies, but also demonstrates the importance of multi-band data for checking the possibility of false-positive planetary signals.

Context: In planetesimal formation theory, several barriers have been proposed, which are bouncing, fragmentation, and radial drift problems. To understand the structure evolution of dust aggregates is a key in the planetesimal formation. Dust grains become fluffy by coagulation in protoplanetary disks. However, once they become fluffy, they are not sufficiently compressed by collisional compression to form compact planetesimals. Aims: We aim to reveal the pathway of the dust structure evolution from dust grains to compact planetesimals. Methods: Using the compressive strength formula, we analytically investigate how fluffy dust aggregates are compressed by static compression due to ram pressure of the disk gas and self gravity of the aggregates in protoplanetary disks. Results: We reveal the pathway of the porosity evolution from dust grains via fluffy aggregates to form planetesimals, circumventing the barriers in planetesimal formation. The aggregates are compressed by the disk gas to the density of 10^{-3} g/cm^3 in coagulation, which is more compact than the case with collisional compression. Then, they are compressed more by self gravity to 10^{-1} g/cm^3 when the radius is 10 km. Although the gas compression decelerate the growth, they grow enough rapidly to avoid the radial drift barrier when the orbital radius is < 6 AU in a typical disk. Conclusions: We propose fluffy dust growth scenario from grains to planetesimals. It enables the icy planetesimal formation in a wide range beyond the snowline in protoplanetary disks. This result proposes a concrete initial condition of planetesimals for the later stages of the planet formation.

Observations of accretion disks around young brown dwarfs have led to the speculation that they may form planetary systems similar to normal stars. While there have been several detections of planetary-mass objects around brown dwarfs (2MASS 1207-3932 and 2MASS 0441-2301), these companions have relatively large mass ratios and projected separations, suggesting that they formed in a manner analogous to stellar binaries. We present the discovery of a planetary-mass object orbiting a field brown dwarf via gravitational microlensing, OGLE-2012-BLG-0358Lb. The system is a low secondary/primary mass ratio (0.080 +- 0.001), relatively tightly-separated (~0.87 AU) binary composed of a planetary-mass object with 1.9 +- 0.2 Jupiter masses orbiting a brown dwarf with a mass 0.022 M_Sun. The relatively small mass ratio and separation suggest that the companion may have formed in a protoplanetary disk around the brown dwarf host, in a manner analogous to planets.

Although many models have been proposed, the physical mechanisms responsible for the formation of low-mass brown dwarfs are poorly understood. The multiplicity properties and minimum mass of the brown-dwarf mass function provide critical empirical diagnostics of these mechanisms. We present the discovery via gravitational microlensing of two very low-mass, very tight binary systems. These binaries have directly and precisely measured total system masses of 0.025 Msun and 0.034 Msun, and projected separations of 0.31 AU and 0.19 AU, making them the lowest-mass and tightest field brown-dwarf binaries known. The discovery of a population of such binaries indicates that brown dwarf binaries can robustly form at least down to masses of ~0.02 Msun. Future microlensing surveys will measure a mass-selected sample of brown-dwarf binary systems, which can then be directly compared to similar samples of stellar binaries.

Context: In protoplanetary disks, dust grains coagulate with each other and grow to form aggregates. As these aggregates grow by coagulation, their filling factor \phi decreases down to \phi << 1. However, comets, the remnants of these early planetesimals, have \phi ~ 0.1. Thus, static compression of porous dust aggregates is important in planetesimal formation. However, the static compression strength has been investigated only for relatively high density aggregates (\phi > 0.1). Aims: We investigate and find the compression strength of highly porous aggregates (\phi << 1). Methods: We perform three dimensional N-body simulations of aggregate compression with a particle-particle interaction model. We introduce a new method of static compression: the periodic boundary condition is adopted and the boundaries move with low speed to get closer. The dust aggregate is compressed uniformly and isotropically by themselves over the periodic boundaries. Results: We empirically derive a formula of the compression strength of highly porous aggregates (\phi << 1). We check the validity of the compression strength formula for wide ranges of numerical parameters, such as the size of initial aggregates, the boundary speed, the normal damping force, and material. We also compare our results to the previous studies of static compression in the relatively high density region (\phi > 0.1) and confirm that our results consistently connect to those in the high density region. The compression strength formula is also derived analytically.

We present the analysis of the gravitational microlensing event OGLE-2011-BLG-0251. This anomalous event was observed by several survey and follow-up collaborations conducting microlensing observations towards the Galactic Bulge. Based on detailed modelling of the observed light curve, we find that the lens is composed of two masses with a mass ratio q=1.9 x 10^-3. Thanks to our detection of higher-order effects on the light curve due to the Earth's orbital motion and the finite size of source, we are able to measure the mass and distance to the lens unambiguously. We find that the lens is made up of a planet of mass 0.53 +- 0.21,M_Jup orbiting an M dwarf host star with a mass of 0.26 +- 0.11 M_Sun. The planetary system is located at a distance of 2.57 +- 0.61 kpc towards the Galactic Centre. The projected separation of the planet from its host star is d=1.408 +- 0.019, in units of the Einstein radius, which corresponds to 2.72 +- 0.75 AU in physical units. We also identified a competitive model with similar planet and host star masses, but with a smaller orbital radius of 1.50 +- 0.50 AU. The planet is therefore located beyond the snow line of its host star, which we estimate to be around 1-1.5 AU.

Microlensing can provide an important tool to study binaries, especially those composed of faint or dark objects. However, accurate analysis of binary-lens light curves is often hampered by the well-known degeneracy between close (s<1) and wide (s>1) binaries, which can be very severe due to an intrinsic symmetry in the lens equation. Here s is the normalized projected binary separation. In this paper, we propose a method that can resolve the close/wide degeneracy using the effect of a lens orbital motion on lensing light curves. The method is based on the fact that the orbital effect tends to be important for close binaries while it is negligible for wide binaries. We demonstrate the usefulness of the method by applying it to an actually observed binary-lens event MOA-2011-BLG-040/OGLE-2011-BLG-0001, which suffers from severe close/wide degeneracy. From this, we are able to uniquely specify that the lens is composed of K and M-type dwarfs located at ~3.5 kpc from the Earth.

The disk wind, which is powered by the radiation force due to spectral lines (line force), is studied for broad absorption line (BAL) quasars. We investigate the structure of the disk wind based on the non-hydrodynamic method and compare with wind properties inferred from X-ray observations of BAL quasars. In this paper, we apply the stellar wind theory to the initial condition (the mass outflow rate at the base of the wind). We found the funnel-shaped winds with a half opening angle of 50^{circ} for the case of epsilon=0.3-0.9 and M_{BH}=10^{7-8.5}M_odot, where epsilon is the Eddington ratio and M_{BH} is the black hole mass. Thus, the absorption features are observed for an observer of which a viewing angle is around 50^{circ}. A probability of BAL quasars is 7-11%, which is roughly consistent the abundance ratio of BAL quasars, 10-15%. Here, the probability is estimated by the solid angle, that the absorbing features would be detected, divided by 4pi. In contrast, if the Eddington ratio is smaller than 0.01 or if the black hole is very massive, M_{BH} < 10^9M_{odot}, the disk wind is not launched due to the less effective line force. Then, the quasars are identified as non-BAL quasars independently of the observer's viewing angle.

We present an analysis of the anomalous microlensing event, MOA-2010-BLG-073, announced by the Microlensing Observations in Astrophysics survey on 2010-03-18. This event was remarkable because the source was previously known to be photometrically variable. Analyzing the pre-event source lightcurve, we demonstrate that it is an irregular variable over time scales >200d. Its dereddened color, $(V-I)_{S,0}$, is 1.221$\pm$0.051mag and from our lens model we derive a source radius of 14.7$\pm$1.3 $R_{\odot}$, suggesting that it is a red giant star. We initially explored a number of purely microlensing models for the event but found a residual gradient in the data taken prior to and after the event. This is likely to be due to the variability of the source rather than part of the lensing event, so we incorporated a slope parameter in our model in order to derive the true parameters of the lensing system. We find that the lensing system has a mass ratio of q=0.0654$\pm$0.0006. The Einstein crossing time of the event, $T_{\rm{E}}=44.3$\pm$0.1d, was sufficiently long that the lightcurve exhibited parallax effects. In addition, the source trajectory relative to the large caustic structure allowed the orbital motion of the lens system to be detected. Combining the parallax with the Einstein radius, we were able to derive the distance to the lens, $D_L$=2.8$\pm$0.4kpc, and the masses of the lensing objects. The primary of the lens is an M-dwarf with $M_{L,p}$=0.16$\pm0.03M_{\odot}$ while the companion has $M_{L,s}$=11.0$\pm2.0M_{\rm{J}}$ putting it in the boundary zone between planets and brown dwarfs.

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.

The Galactic bulge source MOA-2010-BLG-523S exhibited short-term deviations from a standard microlensing lightcurve near the peak of an Amax ~ 265 high-magnification microlensing event. The deviations originally seemed consistent with expectations for a planetary companion to the principal lens. We combine long-term photometric monitoring with a previously published high-resolution spectrum taken near peak to demonstrate that this is an RS CVn variable, so that planetary microlensing is not required to explain the lightcurve deviations. This is the first spectroscopically confirmed RS CVn star discovered in the Galactic bulge.

HI 1225+01 is an intergalactic gas cloud located on the outskirts of Virgo cluster. Its main components are two large clumps of comparable HI masses (M_HI ~ 10^9 Msun) separated by about 100 kpc. One of the clumps hosts a blue low-surface-brightness galaxy J1227+0136, while the other has no identified stellar emission and is sometimes referred to as a promising candidate of a "dark galaxy", an optically invisible massive intergalactic system. We present a deep optical image covering the whole HI 1225+01 structure for the first time, as well as a collection of archival data from ultraviolet to far-infrared (IR) spectral region of the brightest knot "R1" in J1227+0136. We find that R1 has a young stellar population of age 10-100 Myr and mass ~ 10^6 Msun, near-IR excess brightness which may point to the presence of hot dust with color temperature ~ 600 K, and relatively faint mid- to far-IR fluxes corresponding to the dust mass of up to ~ 100 Msun. Overall, it seems to share the general properties with low-metallicity blue compact dwarf galaxies. On the other hand, no optical counterpart to the other clump is found in our deepest-ever image. Now the limiting surface brightness reaches down to R_AB > 28 mag/arcsec2 for any emission extended over 10" (comparable to R1), which is more than one hundred times fainter than the brightest part of the companion galaxy J1227+0136.

Brown dwarfs are important objects because they may provide a missing link between stars and planets, two populations that have dramatically different formation history. In this paper, we present the candidate binaries with brown dwarf companions that are found by analyzing binary microlensing events discovered during 2004 - 2011 observation seasons. Based on the low mass ratio criterion of q < 0.2, we found 7 candidate events, including OGLE-2004-BLG-035, OGLE-2004-BLG-039, OGLE-2007-BLG-006, OGLE-2007-BLG-399/MOA-2007-BLG-334, MOA-2011-BLG-104/OGLE-2011-BLG-0172, MOA-2011-BLG-149, and MOA-201-BLG-278/OGLE-2011-BLG-012N. Among them, we are able to confirm that the companions of the lenses of MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149 are brown dwarfs by determining the mass of the lens based on the simultaneous measurement of the Einstein radius and the lens parallax. The measured mass of the brown dwarf companions are (0.02 +/- 0.01) M_Sun and (0.019 +/- 0.002) M_Sun for MOA-2011-BLG-104/OGLE-2011-BLG-0172 and MOA-2011-BLG-149, respectively, and both companions are orbiting low mass M dwarf host stars. More microlensing brown dwarfs are expected to be detected as the number of lensing events with well covered light curves increases with new generation searches.

We present the analysis of four candidate short duration binary microlensing events from the 2006-2007 MOA Project short event analysis. These events were discovered as a byproduct of an analysis designed to find short timescale single lens events that may be due to free-floating planets. Three of these events are determined to be microlensing events, while the fourth is most likely caused by stellar variability. For each of the three microlensing events, the signal is almost entirely due to a brief caustic feature with little or no lensing attributable mainly to the lens primary. One of these events, MOA-bin-1, is due to a planet, and it is the first example of a planetary event in which stellar host is only detected through binary microlensing effects. The mass ratio and separation are q = 4.9 +- 1.4 x 10^{-3} and s = 2.10 +- 0.05, respectively. A Bayesian analysis based on a standard Galactic model indicates that the planet, MOA-bin-1Lb, has a mass of m_p = 3.7 +- 2.1 M_{Jup}, and orbits a star of M_* = 0.75{+0.33 -0.41} M_solar at a semi-major axis of a = 8.3 {+4.5 -2.7} AU. This is one of the most massive and widest separation planets found by microlensing. The scarcity of such wide separation planets also has implications for interpretation of the isolated planetary mass objects found by this analysis. If we assume that we have been able to detect wide separation planets with a efficiency at least as high as that for isolated planets, then we can set limits on the distribution on planets in wide orbits. In particular, if the entire isolated planet sample found by Sumi et al. (2011) consists of planets bound in wide orbits around stars, we find that it is likely that the median orbital semi-major axis is > 30 AU.

We propose a plausible mechanism to explain the formation of the so-called "obscuring tori" around active galactic nuclei (AGNs) based on three-dimensional hydrodynamic simulations including radiative feedback from the central source. The X-ray heating and radiation pressure on the gas are explicitly calculated using a ray-tracing method. This radiation feedback drives a "fountain", that is, a vertical circulation of gas in the central a few to tens parsecs. Interaction between the non-steady outflows and inflows causes the formation of a geometrically thick torus with internal turbulent motion. As a result, the AGN is obscured for a wide range of solid angles. In a quasi-steady state, the opening angles for the column density toward a black hole < 10^23 cm^-2 are approximately +-30 deg and +-50 deg for AGNs with 10% and 1% Eddington luminosity, respectively. Mass inflows through the torus coexist with the outflow and internal turbulent motion, and the average mass accretion rate to the central parsec region is 2x10^-4 ~ 10^-3, M_sun/yr this is about ten times smaller than accretion rate required to maintain the AGN luminosity. This implies that relatively luminous AGN activity is intrinsically intermittent or that there are other mechanisms, such as stellar energy feedback, that enhance the mass accretion to the center.

Despite astrophysical importance of binary star systems, detections are limited to those located in small ranges of separations, distances, and masses and thus it is necessary to use a variety of observational techniques for a complete view of stellar multiplicity across a broad range of physical parameters. In this paper, we report the detections and measurements of 2 binaries discovered from observations of microlensing events MOA-2011-BLG-090 and OGLE-2011-BLG-0417. Determinations of the binary masses are possible by simultaneously measuring the Einstein radius and the lens parallax. The measured masses of the binary components are 0.43 $M_{\odot}$ and 0.39 $M_{\odot}$ for MOA-2011-BLG-090 and 0.57 $M_{\odot}$ and 0.17 $M_{\odot}$ for OGLE-2011-BLG-0417 and thus both lens components of MOA-2011-BLG-090 and one component of OGLE-2011-BLG-0417 are M dwarfs, demonstrating the usefulness of microlensing in detecting binaries composed of low-mass components. From modeling of the light curves considering full Keplerian motion of the lens, we also measure the orbital parameters of the binaries. The blended light of OGLE-2011-BLG-0417 comes very likely from the lens itself, making it possible to check the microlensing orbital solution by follow-up radial-velocity observation. For both events, the caustic-crossing parts of the light curves, which are critical for determining the physical lens parameters, were resolved by high-cadence survey observations and thus it is expected that the number of microlensing binaries with measured physical parameters will increase in the future.

Microlensing detections of cool planets are important for the construction of an unbiased sample to estimate the frequency of planets beyond the snow line, which is where giant planets are thought to form according to the core accretion theory of planet formation. In this paper, we report the discovery of a giant planet detected from the analysis of the light curve of a high-magnification microlensing event MOA-2010-BLG-477. The measured planet-star mass ratio is $q=(2.181\pm0.004)\times 10^{-3}$ and the projected separation is $s=1.1228\pm0.0006$ in units of the Einstein radius. The angular Einstein radius is unusually large $\theta_{\rm E}=1.38\pm 0.11$ mas. Combining this measurement with constraints on the "microlens parallax" and the lens flux, we can only limit the host mass to the range $0.13<M/M_\odot<1.0$. In this particular case, the strong degeneracy between microlensing parallax and planet orbital motion prevents us from measuring more accurate host and planet masses. However, we find that adding Bayesian priors from two effects (Galactic model and Keplerian orbit) each independently favors the upper end of this mass range, yielding star and planet masses of $M_*=0.67^{+0.33}_{-0.13}\ M_\odot$ and $m_p=1.5^{+0.8}_{-0.3}\ M_{\rm JUP}$ at a distance of $D=2.3\pm0.6$ kpc, and with a semi-major axis of $a=2^{+3}_{-1}$ AU. Finally, we show that the lens mass can be determined from future high-resolution near-IR adaptive optics observations independently from two effects, photometric and astrometric.

Rapid orbital drift of macroscopic dust particles is one of the major obstacles against planetesimal formation in protoplanetary disks. We reexamine this problem by considering porosity evolution of dust aggregates. We apply a porosity model based on recent N-body simulations of aggregate collisions, which allows us to study the porosity change upon collision for a wide range of impact energies. As a first step, we neglect collisional fragmentation and instead focus on dust evolution outside the snow line, where the fragmentation has been suggested to be less significant than inside the snow line because of a high sticking efficiency of icy particles. We show that dust particles can evolve into highly porous aggregates (with internal densities of much less than 0.1 g/cm^3) even if collisional compression is taken into account. We also show that the high porosity triggers significant acceleration in collisional growth. This acceleration is a natural consequence of particles' aerodynamical property at low Knudsen numbers, i.e., at particle radii larger than the mean free path of the gas molecules. Thanks to this rapid growth, the highly porous aggregates are found to overcome the radial drift barrier at orbital radii less than 10 AU (assuming the minimum-mass solar nebula model). This suggests that, if collisional fragmentation is truly insignificant, formation of icy planetesimals is possible via direct collisional growth of submicron-sized icy particles.

Geometrical cross sections of dust aggregates determine their coupling with disk gas, which governs their motions in protoplanetary disks. Collisional outcomes also depend on geometrical cross sections of initial aggregates. In the previous paper, we performed three-dimensional N-body simulations of sequential collisions of aggregates composed of a number of sub-micron-sized icy particles and examined radii of gyration (and bulk densities) of the obtained aggregates. We showed that collisional compression of aggregates is not efficient and that aggregates remain fluffy. In the present study, we examine geometrical cross sections of the aggregates. Their cross sections decreases due to the compression as well as their gyration radii. It is found that a relation between the cross section and the gyration radius proposed by Okuzumi et al. is valid for the compressed aggregates. We also refine the compression model proposed in our previous paper. The refined model enables us to calculate the evolution of both gyration radii and cross sections of growing aggregates and reproduces well our numerical results of sequential aggregate collisions. The refined model can describe non-equal-mass collisions as well as equal-mass case. Although we do not take into account oblique collisions in the present study, oblique collisions would further hinder compression of aggregates.

High-magnification microlensing events provide an important channel to detect planets. Perturbations near the peak of a high-magnification event can be produced either by a planet or a binary companion. It is known that central perturbations induced by both types of companions can be generally distinguished due to the basically different magnification pattern around caustics. In this paper, we present a case of central perturbations for which it is difficult to distinguish the planetary and binary interpretations. The peak of a lensing light curve affected by this perturbation appears to be blunt and flat. For a planetary case, this perturbation occurs when the source trajectory passes the negative perturbation region behind the back end of an arrowhead-shaped central caustic. For a binary case, a similar perturbation occurs for a source trajectory passing through the negative perturbation region between two cusps of an astroid-shaped caustic. We demonstrate the degeneracy for 2 high-magnification events of OGLE-2011-BLG-0526 and OGLE-2011-BLG-0950/MOA-2011-BLG-336. For OGLE-2011-BLG-0526, the $\chi^2$ difference between the planetary and binary model is $\sim$ 3, implying that the degeneracy is very severe. For OGLE-2011-BLG-0950/MOA-2011-BLG-336, the stellar binary model is formally excluded with $\Delta \chi^2 \sim$ 105 and the planetary model is preferred. However, it is difficult to claim a planet discovery because systematic residuals of data from the planetary model are larger than the difference between the planetary and binary models. Considering that 2 events observed during a single season suffer from such a degeneracy, it is expected that central perturbations experiencing this type of degeneracy is common.

We present the analysis of the light curves of 9 high-magnification single-lens gravitational microlensing events with lenses passing over source stars, including OGLE-2004-BLG-254, MOA-2007-BLG-176, MOA-2007-BLG-233/OGLE-2007-BLG-302, MOA-2009-BLG-174, MOA-2010-BLG-436, MOA-2011-BLG-093, MOA-2011-BLG-274, OGLE-2011-BLG-0990/MOA-2011-BLG-300, and OGLE-2011-BLG-1101/MOA-2011-BLG-325. For all events, we measure the linear limb-darkening coefficients of the surface brightness profile of source stars by measuring the deviation of the light curves near the peak affected by the finite-source effect. For 7 events, we measure the Einstein radii and the lens-source relative proper motions. Among them, 5 events are found to have Einstein radii less than 0.2 mas, making the lenses candidates of very low-mass stars or brown dwarfs. For MOA-2011-BLG-274, especially, the small Einstein radius of $\theta_{\rm E}\sim 0.08$ mas combined with the short time scale of $t_{\rm E}\sim 2.7$ days suggests the possibility that the lens is a free-floating planet. For MOA-2009-BLG-174, we measure the lens parallax and thus uniquely determine the physical parameters of the lens. We also find that the measured lens mass of $\sim 0.84\ M_\odot$ is consistent with that of a star blended with the source, suggesting that the blend is likely to be the lens. Although we find planetary signals for none of events, we provide exclusion diagrams showing the confidence levels excluding the existence of a planet as a function of the separation and mass ratio.

The microlensing event OGLE-2008-BLG-510 is characterised by an evident asymmetric shape of the peak, promptly detected by the ARTEMiS system in real time. The skewness of the light curve appears to be compatible both with binary-lens and binary-source models, including the possibility that the lens system consists of an M dwarf orbited by a brown dwarf. The detection of this microlensing anomaly and our analysis demonstrates that: 1) automated real-time detection of weak microlensing anomalies with immediate feedback is feasible, efficient, and sensitive, 2) rather common weak features intrinsically come with ambiguities that are not easily resolved from photometric light curves, 3) a modelling approach that finds all features of parameter space rather than just the `favourite model' is required, and 4) the data quality is most crucial, where systematics can be confused with real features, in particular small higher-order effects such as orbital motion signatures. It moreover becomes apparent that events with weak signatures are a silver mine for statistical studies, although not easy to exploit. Clues about the apparent paucity of both brown-dwarf companions and binary-source microlensing events might hide here.

Astrometry using H2O maser sources in star forming regions is expected to be a powerful tool to study the structures and dynamics of our Galaxy. Honma et al. (2007) (hereafter H2007) claimed that the annual parallax of S269 is determined within an error of 0.008 milliarcsec (mas), concluding that S269 is located at 5.3 kpc +- 0.2 kpc from the sun, and R= 13.1 kpc. They claimed that the rotational velocity of S269 is equal to that of the sun within a 3% error. This small error, however, is hardly understood when taking into account the results of other observations and theoretical studies of galactic dynamics. We here reanalyzed the VERA archival data using the self-calibration method (hybrid mapping), and found that clusters of maser features of S269 are distributed in much wider area than that investigated in H2007. We confirmed that, if we make a narrow region image without considering the presence of multiple maser spots, and only the phase calibration is applied, we can reproduce the same maser structures in H2007. The distribution extent of maser spots in the feature differs 0.2 mas from east to west between our results and H2007. Moreover, we found that change of relative positions of maser spots in the cluster reaches 0.1 mas or larger between observational epochs. This suggests that if one simply assumes the time-dependent, widely distributed maser sources as a stable single point source, it could cause errors of up to 0.1 mas in the annual parallax of S269. Taking into account the internal motions of maser spot clusters, the proper motion of S269 cannot be determined precisely. We estimated that the peculiar motion of S269 with respect to a Galactic circular rotation is ~20 km/s. These results imply that the observed kinematics of maser emissions in S269 cannot give a strong constraint on dynamics of the outer part of the Galaxy, in contrast to the claim by H2007.

Microlensing can provide a useful tool to probe binary distributions down to low-mass limits of binary companions. In this paper, we analyze the light curves of 8 binary lensing events detected through the channel of high-magnification events during the seasons from 2007 to 2010. The perturbations, which are confined near the peak of the light curves, can be easily distinguished from the central perturbations caused by planets. However, the degeneracy between close and wide binary solutions cannot be resolved with a $3\sigma$ confidence level for 3 events, implying that the degeneracy would be an important obstacle in studying binary distributions. The dependence of the degeneracy on the lensing parameters is consistent with a theoretic prediction that the degeneracy becomes severe as the binary separation and the mass ratio deviate from the values of resonant caustics. The measured mass ratio of the event OGLE-2008-BLG-510/MOA-2008-BLG-369 is $q\sim 0.1$, making the companion of the lens a strong brown-dwarf candidate.

Ultraluminous infrared galaxies (ULIRGs) with multiple ($\ge 3$) nuclei are frequently observed. It has been suggested that these nuclei are produced by multiple major mergers of galaxies. The expected rate of such mergers is, however, too low to reproduce the observed number of ULIRGs with multiple nuclei. We have performed high-resolution simulations of the merging of two gas-rich disk galaxies. We found that extremely massive and compact star clusters form from the strongly disturbed gas disks after the first or second encounter between the galaxies. The mass of such clusters reaches $\sim 10^8 M_{\odot}$, and their half-mass radii are $20-30 \rm{pc}$. Since these clusters consist of young stars, they appear to be several bright cores in the galactic central region ($\sim \rm{kpc}$). The peak luminosity of these clusters reaches $\sim 10%$ of the total luminosity of the merging galaxy. These massive and compact clusters are consistent with the characteristics of the observed multiple nuclei in ULIRGs. Multiple mergers are not necessary to explain multiple nuclei in ULIRGs.

[ABRIDGED] Based on high-resolution (R~42000 to 48000) and high signal-to-noise (S/N~50 to 150) spectra obtained with UVES/VLT, we present detailed elemental abundances (O, Na, Mg, Al, Si, Ca, Ti, Cr, Fe, Ni, Zn, Y, and Ba) and stellar ages for 26 microlensed dwarf and subgiant stars in the Galactic bulge. The analysis is based on equivalent width measurements and standard 1-D LTE MARCS model stellar atmospheres. We also present NLTE Li abundances based on line synthesis of the 7Li line at 670.8 nm. We show that the bulge metallicity distribution (MDF) is double-peaked; one peak at [Fe/H]= -0.6 and one at [Fe/H]=+0.3, and with a dearth of stars around solar metallicity. This is in contrast to the MDF derived from red giants in Baade's window, which peaks at this exact value. A simple significance test shows that it is extremely unlikely to have such a gap in the microlensed dwarf star MDF if the dwarf stars are drawn from the giant star MDF. To resolve this issue we discuss several possibilities, but we can not settle on a conclusive solution for the observed differences. We further find that the metal-poor bulge dwarf stars are predominantly old with ages greater than 10\,Gyr, while the metal-rich bulge dwarf stars show a wide range of ages. The metal-poor bulge sample is very similar to the Galactic thick disk in terms of average metallicity, elemental abundance trends, and stellar ages. Speculatively, the metal-rich bulge population might be the manifestation of the inner thin disk. If so, the two bulge populations could support the recent findings, based on kinematics, that there are no signatures of a classical bulge and that the Milky Way is a pure-disk galaxy. Also, recent claims of a flat IMF in the bulge based on the MDF of giant stars may have to be revised based on the MDF and abundance trends probed by our microlensed dwarf stars.

We report the discovery of two field brown dwarfs, ULAS J0128-0041 and ULAS J0321+0051, and the rediscovery of ULAS J0226+0051 (IfA 0230-Z1), in the Sloan Digital Sky Survey (SDSS) southern equatorial stripe. They are found in the course of our follow-up observation program of 1 um excess sources in the United Kingdom Infrared Telescope Infrared Deep Sky Survey. The Gemini Multi-Object Spectrographs spectra at red optical wavelengths (6500-10500 A) are presented, which reveal that they are early-T dwarfs. The classification is also supported by their optical to near-infrared colors. It is noted that ULAS J0321+0051 is one of the faintest currently known T dwarfs. The estimated distances to the three objects are 50-110 pc, thus they are among the most distant field T dwarfs known. Dense temporal coverage of the target fields achieved by the SDSS-II Supernova Survey allows us to perform a simple time-series analysis, which leads to the finding of significant proper motions of 150-290 mas/yr or the transverse velocities of 40-100 km/s for ULAS J0128-0041 and ULAS J0226+0051. We also find that there are no detectable, long-term (a-few-year) brightness variations above a few times 0.1 mag for the two brown dwarfs.

We present the discovery and mass measurement of the cold, low-mass planet MOA-2009-BLG-266Lb, made with the gravitational microlensing method. This planet has a mass of m_p = 10.4 +- 1.7 Earth masses and orbits a star of mass M_* = 0.56 +- 0.09 Solar masses at a semi-major axis of a = 3.2 (+1.9 -0.5) AU and an orbital period of P = 7.6 (+7.7 -1.5} yrs. The planet and host star mass measurements are enabled by the measurement of the microlensing parallax effect, which is seen primarily in the light curve distortion due to the orbital motion of the Earth. But, the analysis also demonstrates the capability to measure microlensing parallax with the Deep Impact (or EPOXI) spacecraft in a Heliocentric orbit. The planet mass and orbital distance are similar to predictions for the critical core mass needed to accrete a substantial gaseous envelope, and thus may indicate that this planet is a "failed" gas giant. This and future microlensing detections will test planet formation theory predictions regarding the prevalence and masses of such planets.

We present hydrodynamic simulations of gas clouds in the central kpc region of the Milky Way that is modeled with a three-dimensional bar potential. Our simulations consider realistic gas cooling and heating, star formation, and supernova feedback. A ring of dense gas clouds forms as a result of X1-X2 orbit transfer, and our potential model results in a ring radius of ~200 pc, which coincides with the extraordinary reservoir of dense molecular clouds in the inner bulge, the Central Molecular Zone (CMZ). The gas clouds accumulated in the CMZ can reach high enough densities to form stars, and with an appropriate choice of simulation parameters, we successfully reproduce the observed gas mass and the star formation rate (SFR) in the CMZ, ~2x10^7 Msun and ~0.1 Msun/yr. Star formation in our simulations takes place mostly in the outermost X2 orbits, and the SFR per unit surface area outside the CMZ is much lower. These facts suggest that the inner Galactic bulge may harbor a mild version of the nuclear star-forming rings seen in some external disk galaxies. Furthermore, from the relatively small size of the Milky Way's nuclear bulge, which is thought to be a result of sustained star formation in the CMZ, we infer that the Galactic inner bulge probably had a shallower density profile or stronger bar elongation in the past.

We propose a new dynamical picture of galactic stellar and gas spirals, based on hydrodynamic simulations in a `live' stellar disk. We focus especially on spiral structures excited in a isolated galactic disk without a stellar bar. Using high-resolution, 3-dimensional N-body/SPH simulations, we found that the spiral features of the gas in galactic disks are formed by essentially different mechanisms from the galactic shock in stellar density waves. The stellar spiral arms and the interstellar matter on average corotate in a galactic potential at any radii. Unlike the stream motions in the galactic shock, the interstellar matter flows into the local potential minima with irregular motions. The flows converge to form dense gas clouds/filaments near the bottom of the stellar spirals, whose global structures resemble dust-lanes seen in late-type spiral galaxies. The stellar arms are non-steady; they are wound and stretched by the galactic shear, and thus local densities of the arm change on a time scale of ~ 100 Myrs, due to bifurcating or merging with other arms. This makes the gas spirals associated with the stellar arms non-steady. The association of dense gas clouds are eventually dissolved into inter-arm regions with non-cirucular motions. Star clusters are formed from the cold, dense gases, whose ages are less than ~30 Myrs, and they are roughly associated with the background stellar arms without a clear spatial offset between gas spiral arms and distribution of young stars.

It has been believed that spirals in pure stellar disks, especially the ones spontaneously formed, decay in several galactic rotations due to the increase of stellar velocity dispersions. Therefore, some cooling mechanism, for example dissipational effects of the interstellar medium, was assumed to be necessary to keep the spiral arms. Here we show that stellar disks can maintain spiral features for several tens of rotations without the help of cooling, using a series of high-resolution three-dimensional $N$-body simulations of pure stellar disks. We found that if the number of particles is sufficiently large, e.g., $3\times 10^6$, multi-arm spirals developed in an isolated disk can survive for more than 10 Gyrs. We confirmed that there is a self-regulating mechanism that maintains the amplitude of the spiral arms. Spiral arms increase Toomre's $Q$ of the disk, and the heating rate correlates with the squared amplitude of the spirals. Since the amplitude itself is limited by the value of $Q$, this makes the dynamical heating less effective in the later phase of evolution. A simple analytical argument suggests that the heating is caused by gravitational scattering of stars by spiral arms, and that the self-regulating mechanism in pure-stellar disks can effectively maintain spiral arms on a cosmological timescale. In the case of a smaller number of particles, e.g., $3\times 10^5$, spiral arms grow faster in the beginning of the simulation (while $Q$ is small) and they cause a rapid increase of $Q$. As a result, the spiral arms become faint in several Gyrs.

Many efforts have been made to model the mass distribution and dynamical evolution of the circumnuclear gas in active galactic nuclei (AGNs). However, chemical evolution is not included in detail in three-dimensional (3-D) hydrodynamic simulations. The X-ray radiation from the AGN can drive the gas chemistry and affect the thermodynamics, as well as the excitation of the interstellar medium (ISM). Therefore, we estimate the effects (on chemical abundances and excitation) of X-ray irradiation by the AGN, for atomic and molecular gas in a 3-D hydrodynamic model of an AGN torus. We obtain the abundances of various species from an X-ray chemical model. A 3-D radiative transfer code estimates the level populations, which result in line intensity maps. Predictions for the CO J=1-0 to J=9-8 lines indicate that mid-J CO lines are excellent probes of density and dynamics in the central (<60 pc) region of the AGN, in contrast to the low-J CO lines. Analysis of the X_CO/\alpha conversion factors shows that only the higher-J CO lines can be used for gas mass determination in AGN tori. The [C II] 158 um emission traces mostly the hot (T_k>1000m K) central (<60 pc) region of the AGN torus. The [C II] 158 um line will be useful for ALMA observations of high redshift (z>1) AGNs. The spatial scales (>0.25 pc) probed with our simulations match the size of the structures that ALMA will resolve in nearby (<45 Mpc at 0.01") galaxies.

We studied the formation process of star clusters using high-resolution N-body/smoothed particle hydrodynamcs simulations of colliding galaxies. The total number of particles is 1.2x10^8 for our high resolution run. The gravitational softening is 5 pc and we allow gas to cool down to \sim 10 K. During the first encounter of the collision, a giant filament consists of cold and dense gas found between the progenitors by shock compression. A vigorous starburst took place in the filament, resulting in the formation of star clusters. The mass of these star clusters ranges from 10^{5-8} Msun. These star clusters formed hierarchically: at first small star clusters formed, and then they merged via gravity, resulting in larger star clusters.

We present the first example of binary microlensing for which the parameter measurements can be verified (or contradicted) by future Doppler observations. This test is made possible by a confluence of two relatively unusual circumstances. First, the binary lens is bright enough (I=15.6) to permit Doppler measurements. Second, we measure not only the usual 7 binary-lens parameters, but also the 'microlens parallax' (which yields the binary mass) and two components of the instantaneous orbital velocity. Thus we measure, effectively, 6 'Kepler+1' parameters (two instantaneous positions, two instantaneous velocities, the binary total mass, and the mass ratio). Since Doppler observations of the brighter binary component determine 5 Kepler parameters (period, velocity amplitude, eccentricity, phase, and position of periapsis), while the same spectroscopy yields the mass of the primary, the combined Doppler + microlensing observations would be overconstrained by 6 + (5 + 1) - (7 + 1) = 4 degrees of freedom. This makes possible an extremely strong test of the microlensing solution. We also introduce a uniform microlensing notation for single and binary lenses, we define conventions, summarize all known microlensing degeneracies and extend a set of parameters to describe full Keplerian motion of the binary lenses.

We report the gravitational microlensing discovery of a sub-Saturn mass planet, MOA-2009-BLG-319Lb, orbiting a K or M-dwarf star in the inner Galactic disk or Galactic bulge. The high cadence observations of the MOA-II survey discovered this microlensing event and enabled its identification as a high magnification event approximately 24 hours prior to peak magnification. As a result, the planetary signal at the peak of this light curve was observed by 20 different telescopes, which is the largest number of telescopes to contribute to a planetary discovery to date. The microlensing model for this event indicates a planet-star mass ratio of q = (3.95 +/- 0.02) x 10^{-4} and a separation of d = 0.97537 +/- 0.00007 in units of the Einstein radius. A Bayesian analysis based on the measured Einstein radius crossing time, t_E, and angular Einstein radius, \theta_E, along with a standard Galactic model indicates a host star mass of M_L = 0.38^{+0.34}_{-0.18} M_{Sun} and a planet mass of M_p = 50^{+44}_{-24} M_{Earth}, which is half the mass of Saturn. This analysis also yields a planet-star three-dimensional separation of a = 2.4^{+1.2}_{-0.6} AU and a distance to the planetary system of D_L = 6.1^{+1.1}_{-1.2} kpc. This separation is ~ 2 times the distance of the snow line, a separation similar to most of the other planets discovered by microlensing.

We have observed 7 new transits of the `hot Jupiter' WASP-5b using a 61 cm telescope located in New Zealand, in order to search for transit timing variations (TTVs) which can be induced by additional bodies existing in the system. When combined with other available photometric and radial velocity (RV) data, we find that its transit timings do not match a linear ephemeris; the best fit \chi^2 values is 32.2 with 9 degrees of freedom which corresponds to a confidence level of 99.982 % or 3.7 \sigma. This result indicates that excess variations of transit timings has been observed, due either to unknown systematic effects or possibly to real TTVs. The TTV amplitude is as large as 50 s, and if this is real, it cannot be explained by other effects than that due to an additional body or bodies. From the RV data, we put an upper limit on the RV amplitude caused by the possible secondary body (planet) as 21 m s^{-1}, which corresponds to its mass of 22-70 M_{Earth} over the orbital period ratio of the two planets from 0.2 to 5.0. From the TTVs data, using the numerical simulations, we place more stringent limits down to 2 M_{Earth} near 1:2 and 2:1 mean motion resonances (MMRs) with WASP-5b at the 3 \sigma level, assuming that the two planets are co-planer. We also put an upper limit on excess of Trojan mass as 43 M_{Earth} (3 \sigma) using both RV and photometric data. We also find that if the possible secondary planet has non- or a small eccentricity, its orbit would likely be near low-order MMRs. Further follow-up photometric and spectroscopic observations will be required to confirm the reality of the TTV signal, and results such as these will provide important information for the migration mechanisms of planetary systems.

We model the gas dynamics of barred galaxies using a three-dimensional, high-resolution, $N$-body+hydrodynamical simulation and apply it to the Milky Way in an attempt to reproduce both the large-scale structure and the clumpy morphology observed in Galactic H\emissiontype{I} and CO $l-v$ diagrams. Owing to including the multi-phase interstellar medium, self-gravity, star-formation and supernovae feedback, the clumpy morphology, as well as the large-scale features, in observed $l-v$ diagrams are naturally reproduced. We identify in our $l-v$ diagrams with a number of not only large-scale peculiar features such as the '3-kpc arm', '135-km s$^{-1}$ arm' and 'Connecting arm' but also clumpy features such as `Bania clumps', and then link these features in a face-on view of our model. We give suggestions on the real structure of the Milky Way and on the fate of gas clumps in the central region.

We report the result of the analysis of a dramatic repeating gravitational microlensing event OGLE-2009-BLG-092/MOA-2009-BLG-137, for which the light curve is characterized by two distinct peaks with perturbations near both peaks. We find that the event is produced by the passage of the source trajectory over the central perturbation regions associated with the individual components of a wide-separation binary. The event is special in the sense that the second perturbation, occurring $\sim 100$ days after the first, was predicted by the real-time analysis conducted after the first peak, demonstrating that real-time modeling can be routinely done for binary and planetary events. With the data obtained from follow-up observations covering the second peak, we are able to uniquely determine the physical parameters of the lens system. We find that the event occurred on a bulge clump giant and it was produced by a binary lens composed of a K and M-type main-sequence stars. The estimated masses of the binary components are $M_1=0.69 \pm 0.11\ M_\odot$ and $M_2=0.36\pm 0.06\ M_\odot$, respectively, and they are separated in projection by $r_\perp=10.9\pm 1.3\ {\rm AU}$. The measured distance to the lens is $D_{\rm L}=5.6 \pm 0.7\ {\rm kpc}$. We also detect the orbital motion of the lens system.

We report the result of the analysis of the light curve of the microlensing event MOA-2009-BLG-016. The light curve is characterized by a short-duration anomaly near the peak and an overall asymmetry. We find that the peak anomaly is due to a binary companion to the primary lens and the asymmetry of the light curve is explained by the parallax effect caused by the acceleration of the observer over the course of the event due to the orbital motion of the Earth around the Sun. In addition, we detect evidence for the effect of the finite size of the source near the peak of the event, which allows us to measure the angular Einstein radius of the lens system. The Einstein radius combined with the microlens parallax allows us to determine the total mass of the lens and the distance to the lens. We identify three distinct classes of degenerate solutions for the binary lens parameters, where two are manifestations of the previously identified degeneracies of close/wide binaries and positive/negative impact parameters, while the third class is caused by the symmetric cycloid shape of the caustic. We find that, for the best-fit solution, the estimated mass of the lower-mass component of the binary is (0.04 +- 0.01) M_sun, implying a brown-dwarf companion. However, there exists a solution that is worse only by \Delta\chi^2 ~ 3 for which the mass of the secondary is above the hydrogen-burning limit. Unfortunately, resolving these two degenerate solutions will be difficult as the relative lens-source proper motions for both are similar and small (~ 1 mas/yr) and thus the lens will remain blended with the source for the next several decades.

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.

Based on the physical model of a supermassive black hole (SMBH) growth via gas accretion in a circumnuclear disk (CND) proposed by Kawakatu & Wada (2008), we describe the formation of high-$z$ ($z > 6$) quasars (QSOs) whose BH masses are M_{BH}> 10^{9} M_{\odot}. We derive the necessary conditions to form QSOs at z > 6 by only gas accretion: (i) A large mass supply with M_{sup} > 10^{10}M_{\odot} from host galaxies to CNDs, because the final BH mass is only 1-10% of the total supplied mass from QSO hosts. (ii) High star formation efficiency for a rapid BH growth. We also find that if the BH growth is limited by the Eddington accretion, the final BH mass is greatly suppressed. Thus, the super-Eddington growth is required for the QSO formation. The evolution of the QSO luminosity depends on the redshift z_{i} at which accretion onto a seed BH is initiated. In other words, the brighter QSOs at z >6 favor the late growth of SMBHs (i.e., z_{i}=10) rather than early growth (i.e., z_{i}=30). Moreover, we predict the observable properties and the evolution of QSOs at z >6. In a QSO phase, there should exist a stellar rich massive CND, whose gas mass is about 10% of the dynamical mass inside 0.1-1 kpc}. On the other hand, in a phase where the BH grows (i.e., a proto-QSO phase), the proto-QSO has a gas rich massive CNDs whose gas mass is comparable to the dynamical mass (abridged).

We present new high resolution numerical simulations of the ISM in a central R ~32 parsecs region around a supermassive black hole at a galactic center. Three-dimensional hydrodynamic modeling of the ISM (Wada & Norman 2002) with the nuclear starburst now includes tracking of the formation of molecular hydrogen out of the neutral hydrogen phase. In a quasi equilibrium state, mass fraction of H2 is about 0.4 (total H2 mass is ~1.5 10^6 Msun) of the total gas mass for the uniform far UV (FUV) with G_0 = 10. The gas forms an inhomogeneous disk, whose scale-height becomes larger in the outer region. H2 forms a thin nuclear disk in the inner ~ 5 pc, which is surrounded by molecular clouds swelled up toward h < 10 pc. The velocity field of the disk is highly turbulent in the torus region, whose velocity dispersion is ~ 20 km/s on average. Average supernova rate (SNR) of ~ 5 10^-5/yr is large enough to energize these structures. Gas column densities toward the nucleus larger than 10^22 cm^-2 are observed if the viewing angle is smaller than \theta_v ~ 50 deg from the edge-on. However, the column densities are distributed over almost two orders of magnitude around the average for any given viewing angle due to the clumpy nature of the torus. For a stronger FUV (G_0 =100), the total H2 mass in an equillibrium is only slightly smaller (~ 0.35). Finally the morphology and kinematics of the circumnuclear molecular gas disks emerging from our models is similar to that revealed by recent near infrared observations using VLTI/Keck.

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

According to our common understandings, the original surface of a short-period comet nucleus has been lost by sublimation processes during its close approaches to the Sun. Sublimation results in the formation of a dust mantle on the retreated surface and in chemical differentiation of ices over tens or hundreds of meters below the mantle. In the course of NASA's Deep Impact mission, optical and infrared imaging observations of the ejecta plume were conducted by several researchers, but their interpretations of the data came as a big surprise: (1) The nucleus of comet 9P/Tempel 1 is free of a dust mantle, but maintains its pristine crust of submicron-sized carbonaceous grains; (2) Primordial materials are accessible already at a depth of several tens of cm with abundant silicate grains of submicrometer sizes. In this study, we demonstrate that a standard model of cometary nuclei explains well available observational data: (1) A dust mantle with a thickness of ~1-2 m builds up on the surface, where compact aggregates larger than tens of micrometers dominate; (2) Large fluffy aggregates are embedded in chemically differentiated layers as well as in the deepest part of the nucleus with primordial materials. We conclude that the Deep Impact results do not need any peculiar view of a comet nucleus.

We propose a new evolutionary model of a supermassive black hole (SMBH) and a circumnuclear disk (CND), taking into account the mass-supply from a host galaxy and the physical states of CND. In the model, two distinct accretion modes depending on gravitational stability of the CND play a key role on accreting gas to a SMBH. (i) If the CMD is gravitationally unstable, energy feedback from supernovae (SNe) supports a geometrically thick, turbulent gas disk. The accretion in this mode is dominated by turbulent viscosity, and it is significantly larger than that in the mode (ii), i.e., the CMD is supported by gas pressure. Once the gas supply from the host is stopped, the high accretion phase ($\sim 0.01- 0.1 M_{\odot} {\rm yr}^{-1}$) changes to the low one (mode (ii), $\sim 10^{-4} M_{\odot} {\rm yr}^{-1}$), but there is a delay with $\sim 10^{8}$ yr. Through this evolution, the gas-rich CND turns into the gas poor stellar disk. We found that not all the gas supplied from the host galaxy accrete onto the SMBH even in the high accretion phase (mode (i)), because the part of gas is used to form stars. As a result, the final SMBH mass ($M_{\rm BH,final}$) is not proportional to the total gas mass supplied from the host galaxy ($M_{\rm sup}$); $M_{\rm BH,final}/M_{\rm sup}$ decreases with $M_{\rm sup}$.This would indicate that it is difficult to form a SMBH with $\sim 10^{9} M_{\odot}$ observed at high-$z$ QSOs. The evolution of the SMBH and CND would be related to the evolutionary tracks of different type of AGNs.

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

The ``galactic shocks'' \citep{fujimoto68,roberts69} is investigated using a full three-dimensional hydrodynamic simulations, taking into account self-gravity of the ISM, radiative cooling, and star formation followed by energy feedback from supernovae. This is an essential progress from the previous numerical models, in which 2-D isothermal, non-self-gravitating gas is assumed. We find that the classic galactic shocks appears is unstable and transient, and it shifts to a globally quasi-steady, inhomogeneous pattern due to non-linear development of instabilities in the disk. The spiral patterns consists of many GMC-like dense condensations, but those local structures are not steady, and they evolves into irregular spurs in the inter-arm regions. Energy feedback from supernovae do not destroy the quasi-steady spiral arms, but it mainly contributes to vertical motion and structures of the ISM.

We explore millimeter line diagnostics of an obscuring molecular torus modeled by a hydrodynamic simulation with three-dimensional nonLTE radiative transfer calculations. Based on the results of high-resolution hydrodynamic simulation of the molecular torus around an AGN, we calculate intensities of HCN and HCO^{+} rotational lines as two representative high density tracers. The three-dimensional radiative transfer calculations shed light on a complicated excitation state in the inhomogeneous torus, even though a spatially uniform chemical structure is assumed. Our results suggest that HCN must be much more abundant than HCO^{+} in order to obtain a high ratio ($R_{HCN/HCO+}\sim 2$) observed in some of the nearby galaxies. There is a remarkable dispersion in the relation between integrated intensity and column density, indicative of possible shortcomings of HCN(1-0) and HCO^{+}(1-0) lines as high density tracers. The internal structures of the inhomogeneous molecular torus down to subparsec scale in external galaxies will be revealed by the forthcoming Atacama Large Millimeter/submillimeter Array (ALMA). The three-dimensional radiative transfer calculations of molecular lines with high-resolution hydrodynamic simulation prove to be a powerful tool to provide a physical basis for molecular line diagnostics of the central regions of external galaxies.

Possible origins of the molecular absorption discovered in some ULIRGs are investigated, based on a 3-D hydrodynamic model of star-forming interstellar gas in a galactic central region. The blue-shifted, warm ($\sim 200-300$ K), dense ($>10^6$ cm$^{-3}$) molecular gas suggested by CO absorption in IRAS 08752+3915 could be caused by the innermost region of the inhomogeneous inter-stellar medium (ISM) around a supermassive black hole. The infrequent observations of the dense gas with absorption in ULIRGs and Seyfert 2 galaxies could simply suggest that the high-density regions occupy only a very small volume fraction of the obscuring material. This is naturally expected if the inhomogeneous structure of the ISM is caused by non-linear development of instabilities. The model predicts a turbulent velocity field in the obscuring material, therefore blue- and red-shifted gases should be observable with nearly the same probability for the large enough statistical samples.

The probability distribution functions (PDF) of density of the ISM in galactic disks and global star formation rate are discussed. Three-dimensional hydrodynamic simulations show that the PDFs in globally stable, inhomogeneous ISM in galactic disks are well fitted by a single log-normal function over a wide density range. The dispersion of the log-normal PDF (LN-PDF) is larger for more gas-rich systems, and whereas the characteristic density of LN-PDF, for which the volume fraction becomes the maximum, does not significantly depend on the initial conditions. %At the characteristic density, thermal pressure and kinetic %pressure due to turbulent motion are comparable, and the material is %statistically stagnated in a turbulent flow. Supposing the galactic ISM is characterized by the LN-PDF, we give a global star formation rate (SFR) as a function of average gas density, a critical local density for star formation, and star formation efficiency. We find that the observed SFR is well-fitted by the theoretical SFR in a wide range of the global gas density ($10 - 10^4 M_\odot$ pc$^{-2}$). Star formation efficiency (SFE) for high density gas ($n > 10^3$ cm$^{-3}$) is SFE $= 0.001 - 0.01$ for normal spiral galaxies, and SFE $= 0.01 - 0.1$ for starburst galaxies. The LN-PDF and SFR proposed here could be applicable for modeling star formation on a kpc-scale in galaxies or numerical simulations of galaxy formation, in which the numerical resolution is not fine enough to describe the local star formation.

In order to examine the ``giant impact hypothesis'' for the Moon formation, we run the first grid-based, high-resolution hydrodynamic simulations for an impact between proto-Earth and a proto-planet. The spatial resolution for the impact-generated disk is greatly improved from previous particle-based simulations. This allows us to explore fine structures of a circumterrestrial debris disk and its long-term evolution. We find that in order to form a debris disk from which a lunar-sized satellite can be accumulated, the impact must result in a disk of mostly liquid or solid debris, where pressure is not effective, well before the accumulation process starts. If the debris is dominated by vapor gas, strong spiral shocks are generated, and therefore the circumterrestrial disk cannot survive more than several days. This suggests that there could be an appropriate mass range for terrestrial planets to harbor a large moon as a result of giant impacts, since vaporization during an impact depends on the impact energy.

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

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

This is the first attempt to construct detailed X-ray spectra of clusters of galaxies from the results of high-resolution hydrodynamic simulations and simulate X-ray observations in order to study velocity fields of the intracluster medium (ICM). The hydrodynamic simulations are based on the recently proposed tsunami model, in which cluster cores are affected by bulk motions of the ICM and turbulence is produced. We note that most other solutions of the cooling flow problem also involve the generation of turbulence in cluster cores. From the mock X-ray observations with Astro-E2 XRS, we find that turbulent motion of the ICM in cluster cores could be detected with the satellite. The Doppler shifts of the metal lines could be used to discriminate among turbulence models. The gas velocities measured through the mock observations are consistent with the line-emission weighted values inferred directly from hydrodynamic simulations.

We have performed three-dimensional, non-LTE radiative transfer calculations for 12CO and 13CO lines, applying them to our high-resolution hydrodynamic models of the `torus' around a SMBH in an active galactic nucleus. The hydrodynamic simulations reveal inhomogeneous and turbulent gas structure on a sub-pc scale in a circum-nuclear starburst region. Thick disks interlaced with filaments, clumps and holes are naturally formed due to the interplay among energy feedback from supernovae, self-gravity of the gas, galactic rotation, and radiative cooling. The intensity maps of the molecular lines for the circum-nuclear disks show a clumpy structure reflecting the intrinsic inhomogeneity and turbulent motion of the gas disk. The fine structure of the `torus' could be resolved in the nearby active galaxies using ALMA. We also found that the X-factor is not uniformly distributed in the central 100 pc region. The X-factor derived for 12CO(J=1-0) intensity depends strongly on the intensity, whereas the X for 12CO(3-2) is nearly constant over two orders of magnitude of the intensity. The suggested conversion factor for the molecular gas mass is X_CO (3-2) ~ 0.27 x 10^{20} cm^{-2} (K km s^{-1})^{-1}.

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

We report on infrared K- (2-2.5 micron) and L-band (2.8-4.1 micron) slit spectroscopy of 23 Seyfert 1 galaxies in the CfA and 12 micron samples. A polycyclic aromatic hydrocarbon (PAH) emission feature at 3.3 micron in the L band is primarily used to investigate nuclear star-forming activity in these galaxies. The 3.3 micron PAH emission is detected in 10 sources (=43%), demonstrating that detection of nuclear star-formation in a significant fraction of Seyfert 1 galaxies is now feasible. For the PAH-detected nuclei, the surface brightness values of the PAH emission are as high as those of typical starbursts, suggesting that the PAH emission probes the putative nuclear starbursts in the dusty tori around the central active galactic nuclei (AGNs). The magnitudes of the nuclear starbursts are quantitatively estimated from the observed 3.3 micron PAH emission luminosities. The estimated starburst luminosities relative to some indicators of AGN powers in these Seyfert 1s are compared with 32 Seyfert 2s in the same samples that we have previously observed. We find that there is no significant difference in nuclear starburst to AGN luminosity ratios between Seyfert 1 and 2 galaxies, and that nuclear starburst luminosity positively correlates with AGN power in both types of Seyferts. Our results favor a slightly-modified AGN unification model, which predicts that nuclear starbursts occurring in the dusty tori of Seyferts are physically connected to the central AGNs, rather than the classical unification paradigm, in which the dusty tori simply hide the central AGNs of Seyfert 2s and reprocess AGN radiation as infrared dust emission in Seyferts. No significant differences in nuclear star formation properties are recognizable between Seyfert 1s in the CfA and 12 micron samples.

Based on high-resolution two-dimensional hydrodynamic simulations, we show that the bulk gas motions in a cluster of galaxies, which are naturally expected during the process of hierarchical structure formation of the universe, have a serous impact on the core. We found that the bulk gas motions represented by acoustic-gravity waves create local but strong turbulence, which reproduces the complicated X-ray structures recently observed in cluster cores. Moreover, if the wave amplitude is large enough, they can suppress the radiative cooling of the cores. Contrary to the previous studies, the heating is operated by the turbulence, not weak shocks. The turbulence could be detected in near-future space X-ray missions such as ASTRO-E2.

We found that loosely wound spiral shocks in an isothermal gas disk caused by a non-axisymmetric potential are hydrodynamically unstable, if the shocks are strong enough. High resolution, global hydrodynamical simulations using three different numerical schemes, i.e. AUSM, CIP, and SPH, show similarly that trailing spiral shocks with the pitch angle of larger than ~10 deg wiggle, and clumps are developed in the shock-compressed layer. The numerical simulations also show clear wave crests that are associated with ripples of the spiral shocks. The spiral shocks tend to be more unstable in a rigidly rotating disk than in a flat rotation. This instability could be an origin of the secondary structures of spiral arms, i.e. the spurs/fins, observed in spiral galaxies. In spite of this local instability, the global spiral morphology of the gas is maintained over many rotational periods. The Kelvin-Helmholtz (K-H) instability in a shear layer behind the shock is a possible mechanism for the wiggle instability. The Richardson criterion for the K-H stability is expressed as a function of the Mach number, the pitch angle, and strength of the background spiral potential. The criterion suggests that spiral shocks with smaller pitch angles and smaller Mach numbers would be more stable, and this is consistent with the numerical results.

Using an analytical model and numerical simulations, we show that acoustic waves generated by turbulent motion in intracluster medium effectively heat the central region of a so-called ``cooling flow'' cluster. We assume that the turbulence is generated by substructure motion in a cluster or cluster mergers. Our analytical model can reproduce observed density and temperature profiles of a few clusters. We also show that waves can transfer more energy from the outer region of a cluster than thermal conduction alone. Numerical simulations generally support the results of the analytical study.

Supplying gas to the galactic central regions is one of key ingredients for AGN activity. I will review various fueling mechanisms for a R ~ 1 kpc region, determined mainly by numerical simulations over the last decade. I will also comment on the bars-within-bars mechanism. Observations suggest that the stellar bar is not a sufficient condition for gas fueling. Moreover, considering the various factors for the onset of gas accretion, stellar bars would not even be a necessary condition. I introduce recent progress obtained through our two- and three-dimensional, high resolution hydrodynamical simulations of the ISM in the central 0.1--1 kpc region of galaxies. Possible structure of the obscuring molecular tori around AGNs is also shown. The nuclear starburst is an important factor in determining the structure of the molecular tori and the mass accretion rate to the nucleus. It is natural that the ISM in the central 100 pc region is a highly inhomogeneous and turbulent structure. As a result, gas accretion to the central parsec region should be time dependent and stochastic. The conventional picture of gas fueling and the AGN unified model may be modified in many respects.

We investigate the chemo-dynamical effects of multiple supernova explosions in the central region of primordial galaxies using three-dimensional hydrodynamical simulations of the inhomogenous interstellar medium down to parsec-scales. We find that the final protogalactic structure and metal distribution depend strongly on the number of SNe. Specifically, 1) 1000 SNe after an instantaneous burst of star formation are sufficient to almost completely blow away the gas in these systems, whereas 2) 100 SN explosions trigger the collapse of the protogalactic cloud, leading to the formation of a cold, dense clumpy disk (n > 300 cm^-3) with metallicity, Z = 4 10^-4 Z_sun. These results imply that the metallicity of the ``second generation'' of stars could be Z ~ 10^-4 Z_sun, and that the environment to form metal-free stars in protogalaxies may be lost relatively quickly (< 10^7 yr) after the first burst of Z=0 star formation. The recently discovered ultra-metal-poor star (Christlieb et al. 2002) may represent surviving members of such second-generation star formation.

To interpret H_2 quasar absorption line observations in Damped Ly-alpha clouds (DLAs), we model the H_2 spatial distribution within a DLA. Based on numerical simulations of disk structures with parameters similar to those derived for such absorbers, we calculate the H_2 distribution as a function of ultraviolet background (UVB) intensity and dust-to-gas ratio. For typical values of these two quantities we find that the area in which the H_2 fraction exceeds 10^{-6} (typical observational detection limit) only covers $\la 10$% of the disk surface, i.e. H_2 has a very inhomogeneous, clumpy distribution even at these low abundance levels. This explains the relative paucity of H_2 detections in DLAs. We also show the dependence of the covering fraction of H_2 on dust-to-gas ratio and UVB intensity and we comment on the physics governing the H_2 chemical network at high redshift.

Rotation curves are often used to estimate the mass distribution of spiral galaxies, assuming circular rotation of disks. However, non-circular motions caused by a non-axisymmetric gravitational potential, such as a stellar bar, may disturb the velocity field, resulting in errors in mass estimation, especially in the central regions of galaxies, because the line-of-sight velocity depends on the viewing angles in a non-axisymmetric flow. Observing rotation curves of edge-on galaxies in time-dependent numerical simulations from different viewing angles, we obtain errors in the estimation of galactic mass from the rotation curves. In the most extreme case, the ellipticity of gas orbits is as high as 0.8 in the central regions, even if the bar potential is weak. When rotation curves are defined as the highest velocity envelope of position-velocity diagrams, the mass estimated from the rotation curves is larger than the true mass by a factor of five for 15% of the viewing angles, and the ratio between the apparent mass and true mass is less than six for any viewing angle. The overestimation in mass occurs more frequently than the underestimation.

High-resolution, 2-D hydrodynamical simulations with a large dynamic range are performed to study the turbulent nature of the interstellar medium (ISM) in galactic disks. The simulations are global, where the self-gravity of the ISM, realistic radiative cooling, and galactic rotation are taken into account. In the analysis undertaken here, feedback processes from stellar energy source are omitted. We find that the velocity field of the disk in a non-linear phase shows a steady power-law energy spectrum over three-orders of magnitude in wave number. This implies that the random velocity field can be modeled as fully-developed, stationary turbulence. Gravitational and thermal instabilities under the influence of galactic rotation contribute to form the turbulent velocity field. The Toomre effective Q value, in the non-linear phase, ranges over a wide range, and gravitationally stable and unstable regions are distributed patchily in the disk. These results suggest that large-scale galactic rotation coupled with the self-gravity of the gas can be the ultimate energy sources that maintain the turbulence in the local ISM. We find that our models of turbulent rotating disks are consistent with the velocity dispersion of an extended HI disk in the dwarf galaxy, NGC 2915, where there is no prominent active star formation. Numerical simulations show that the stellar bar in NGC 2915 enhances the velocity dispersion, and it also drives spiral arms as observed in the HI disk.

The structure of obscuring matter in the environment of active galactic nuclei with associated nuclear starbursts is investigated using 3-D hydrodynamical simulations. Simple analytical estimates suggest that the obscuring matter with energy feedback from supernovae has a torus-like structure with a radius of several tens of parsecs and a scale height of about 10 pc. These estimates are confirmed by the fully non-linear numerical simulations, in which the multi-phase inhomogeneous interstellar matter and its interaction with the supernovae are consistently followed. The globally stable, torus-like structure is highly inhomogeneous and turbulent. To achieve the high column densities (> 10^{24} cm^{-2}) as suggested by observations of some Seyfert 2 galaxies with nuclear starbursts, the viewing angle should be larger than about 70 degree from the pole-on for a 10^8 solar mass massive black hole. Due to the inhomogeneous internal structure of the torus, the observed column density is sensitive to the line-of-sight, and it fluctuates by a factor of order 100. The covering fraction for N > 10^{23} cm^{-2} is about 0.4. The average accretion rate toward R < 1 pc is 0.4 solar mass/yr, which is boosted to twice that in the model without the energy feedback.

The structure of the interstellar medium in the central kpc region of a galaxy with a weak bar-like potential is investigated taking into account realistic cooling and heating processes and the self-gravity of the gas. Using high resolution hydrodynamical simulations, it is revealed that the resonant structures (e.g. smooth spiral shocks and a nuclear ring) are very different from those seen in past numerical models where simple models of the ISM, i.e. non-self-gravitating, isothermal gas were assumed. We find that the pc-scale filaments and clumps form large scale spirals, which resemble those seen in real galaxies. The fine structures are different between the arms and in the nuclear region. The next generation millimeter interferometer (ALMA) may reveal the fine structures of the cold gas in nearby galaxies. We also find a large scale anisotropy in the gas temperature, which is caused due to non-circular velocity field of the gas.The damped orbit model based on the epicyclic approximation explains the distribution of the hot (> 10^4 K) and cold (< 100 K) gases appearing alternately around the galactic center. Because of the temperature anisotropy, cold gases observed by molecular lines do not necessarily represent the real gas distribution in galaxies. Position-Velocity diagrams depend strongly on the viewing angles. As a result, the rotational velocity inferred from the PV maps could be two times larger or smaller than the true circular velocity.

Using high-resolution, three-dimensional hydrodynamical simulations, we investigate the structure of the interstellar medium in the central hundred pc region in galaxies, taking into account self-gravity of the gas, radiative cooling from 10 K to $10^8$ K, and energy feedback from supernovae. Similar to the previous two-dimensional results produced by Wada and Norman, we find that a gravitationally and thermally unstable ISM evolves, in a self-stabilizing manner, into a quasi-stable thin disk, which is characterized by a network of cold ($T < 100$ K), dense clumps and filaments, and hot ($T > 10^6$ K), diffuse medium. Supernova explosions blow the diffuse gases from the disk, and as a result, a quasi-steady diffuse halo, which is not uniform but has a plume-like structure, is formed. The density probability distribution function (PDF) in a quasi-steady state is well fitted by a Log-Normal function over about seven orders of magnitude.

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

Based on our high resolution, two-dimensional hydrodynamical simulations, we propose that large cavities may be formed by the nonlinear development of the combined thermal and gravitational instabilities, without need for stellar energy injection in a galaxy modeling the Large Magellanic Clouds (LMC). Our numerical model of the star formation allows us to follow the evolution of the blastwaves due to supernovae in the inhomogeous, multi-phase, and turbulent-like media self-consistently. Formation of kpc-scale inhomogeneity, such as cavities, observed HI map of the LMC, is suppressed by frequent supernovae (average supernova rate for the whole disk is ~0.001/yr). However the supernova explosions are necessary for the hot component (T_g > 10^{6-7} K). Position-velocity maps show that kpc-scale shells/arcs formed through the nonlinear evolution in a model without stellar energy feedback has similar kinematics to explosional phenomena, such as supernovae. We also find that dense clumps and filamentary structure are formed due to a natural consequence of the non-linear evolution of the multi-phage ISM. Although the ISM in a small scale looks turbulent-like and transient, the global structure of the ISM is quasi-stable. We compare the observations of HI and molecular gas of the LMC with the numerically obtained HI and CO map. The morphology and statistical properties of the numerical HI and CO maps are discussed.

We investigate the origin of a unified scaling relation in spiral galaxies. Observed spiral galaxies are spread on a plane in the three-dimensionallogarithmic space of luminosity L, radius R and rotation velocity V. The plane is expressed as $L \propto (V R)^{\alpha}$ in I-passband, where $\alpha$ is a constant. On the plane, observed galaxies are distributed in an elongated region which looks like the shape of a surfboard. The well-known scaling relations, L-V (Tully-Fisher relation), V-R (also the Tully-Fisher relation) and R-L (Freeman's law), can be understood as oblique projections of the surfboard-like plane into 2-D spaces. This unified interpretation of the known scaling relations should be a clue to understand the physical origin of all the relations consistently. Furthermore, this interpretation can also explain why previous studies could not find any correlation between TF residuals and radius. In order to clarify the origin of this plane, we simulate formation and evolution of spiral galaxies with the N-body/SPH method, including cooling, star formation and stellar feedback. Initial conditions are set to isolated 14 spheres with two free parameters, such as mass and angular momentum. The CDM (h=0.5, $\Omega_0=1$) cosmology is considered as a test case. The simulations provide the following two conclusions: (a) The slope of the plane is well reproduced but the zero-point is not. This zero-point discrepancy could be solved in a low density ($\Omega_0<1$) and high expansion (h>0.5) cosmology. (b) The surfboard-shaped plane can be explained by the control of galactic mass and angular momentum.

A combination of Hubble Space Telescope (HST) WFPC2 and NICMOS images are used to investigate the gas/dust and stellar structure inside the central 300 pc of the nearby active galaxy NGC 4303. The NICMOS H-band (F160W) image reveals a bright core and a nuclear elongated bar-like structure of 250 pc in diameter. The bar is centered on the bright core, and its major axis is oriented in proyection along the spin axis of the nuclear gaseous rotating disk recently detected (Colina & Arribas 1999). The V-H (F606W - F160W) image reveals a complex gas/dust distribution with a two-arm spiral structure of about 225 pc in radius. The southwestern arm is traced by young star-forming knots while the northeastern arm is detected by the presence of dust lanes. These spirals do not have a smooth structure but rather they are made of smaller flocculent spirals or filament-like structures. The magnitudes and colors of the star-forming knots are typical of clusters of young stars with masses of 0.5 to 1 x $10^5 M_{solar}, and ages of 5 to 25 million years. The overall structure of the nuclear spirals as well as the size, number and masses of the star-forming knots are explained in the context of a massive gaseous nuclear disk subject to self-gravitational instabilities and to the gravitational field created by the nuclear bar. According to the model, the gaseous disk has a mass of about 5 x 10^7 M_{solar} inside a radius of 400 pc, the bar has a radius of 150 pc and a pattern speed of about 0.5 Myr^{-1}, and the average mass accretion rate into the core (R < 8 pc) is about 0.01 M_{solar}$ yr^{-1} for about 80 Myr.

We discuss the origin of the Tully-Fisher (TF) relation using the $N$-body/SPH method, which includes cooling, star formation and stellar feedback of energy, mass and metals. We consider initially rotating overdense spheres, and trace formation processes of disk galaxies from $z=25$ to $z=0$ in the Cold Dark Matter (CDM) cosmology. To clarify the origin of the TF relation, we simulate formation of 14 galaxies with different masses and spin parameters, and compute observable values, such as the total magnitude and the line-width. We find that the simulated galaxies reproduce the slope and scatter of the TF relation: the slope is originated in the difference of total galactic masses, and the scatter is produced by the difference of initial spin parameters. As well as the TF relation, observed features of spiral galaxies, such as the exponential light-profile and the flat rotation curve, are reproduced in our simulations, which were assumed {\it a priori} in past semi-analytical approaches.

Using high resolution, two-dimensional hydrodynamical simulations, we investigate the evolution of a self-gravitating multi-phase interstellar medium in the central kiloparsec region of a galactic disk. We find that a gravitationally and thermally unstable disk evolves, in a self-stabilizing manner, into a globally quasi-stable disk that consists of cold (T < 100 K), dense clumps and filaments surrounded by hot (T > 10^4 K), diffuse medium. The quasi-stationary, filamentary structure of the cold gas is remarkable. The hot gas, characterized by low-density holes and voids, is produced by shock heating. The shocks derive their energy from differential rotation and gravitational perturbations due to the formation of cold dense clumps. In the quasi-stable phase where cold and dense clouds are formed, the effective stability parameter, Q, has a value in the range 2-5. The dynamic range of our multi-phase calculations is 10^6 - 10^7 in both density and temperature. Phase diagrams for this turbulent medium are analyzed and discussed.

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.

Using N-body and SPH modeling we perform 3D numerical simulations of head-on collisions between gas rich disk galaxies, including collisions between counter-rotating disks and off-center collisions. Pure stellar intruders do not produce gaseous plumes similar to those seen in the Cartwheel and VII Zw466 complexes of interacting galaxies; the presence of gas in an intruder galaxy and radiative cooling are important for the formation of a gaseous plume extending from the disk of a target galaxy. A noticeable plume structure can be formed if the mass of an intruder is a few percent of the mass of the primary. The halo of the intruder is stripped in the collision, and dispersed particles form a broad stellar bridge connecting the two galaxies. The fraction of the intruder's halo dispersed in the collision depends on the total mass of the intruder, and low-mass intruders lose most of their mass.

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.