We discovered a new growth mode of dust grains to km-sized bodies in protoplanetary disks that evolve by viscous accretion and magnetically driven disk winds (MDWs). We solved an approximate coagulation equation of dust grains with time-evolving disks that consist of both gas and solid components by a one dimensional model. With the collisional growth of dust grains, all solid particles initially drift inward toward the central star by the gas drag force. However, the radial profile of gas pressure, $P$, is modified by the MDW that disperses the gas in an inside-out manner. As a result, a local concentration of solid particles is created by the converging radial flux of drifting dust grains at the location with the convex upward profile of $P$. When the dimensionless stopping time, or the Stokes number, St, there exceeds unity, the solid particles spontaneously reach the growth dominated state because of the positive feedback between the suppressed radial drift and the enhanced accumulation of dust particles that drift from the outer part. These consecutive processes proceed even for spherical dust grains under the Epstein gas drag. Once the solid particles are in the drift limited state, the above mentioned condition of St $\gtrsim 1$ for the dust growth is equivalent with \begin{equation} \Sigma_{\rm d}/\Sigma_{\rm g} \gtrsim \eta, \nonumber \end{equation} where $\Sigma_{\rm d}/\Sigma_{\rm g}$ is the dust-to-gas surface-density ratio and $\eta$ is dimensionless radial pressure-gradient force. As a consequence of the successful growth of dust grains, a ring-like structure containing planetesimal-sized bodies is formed at the inner part of the protoplanetary disks. Such a ring-shaped concentration of planetesimals is expected to play a vital role in the subsequent planet formation.

We develop a framework for magnetohydrodynamical (MHD) simulations in a local cylindrical shearing box by extending the formulation of the Cartesian shearing box. We construct shearing-periodic conditions at the radial boundaries of a simulation box from the conservation relations of the basic MHD equations, taking into account the explicit radial dependence of physical quantities. We demonstrate quasi-steady mass accretion, which cannot be handled by the standard Cartesian shearing box model, with an ideal MHD simulation in a vertically unstratified cylindrical shearing box up to 200 rotations. In this demonstrative run we set up (i) net vertical magnetic flux, (ii) a locally isothermal equation of state, and (iii) a sub-Keplerian equilibrium rotation, whereas the sound velocity and the initial Alfven velocity have the same radial dependence as that of the Keplerian velocity. Inward mass accretion is induced to balance with the outward angular momentum flux of the MHD turbulence triggered by the magnetorotational instability in a self-consistent manner. We discuss detailed physical properties of the saturated magnetic field, in comparison to the results of a Cartesian shearing box simulation.

When a planet forms a deep gap in a protoplanetary disk, dust grains cannot pass through the gap. As a consequence, the density of the dust grains can increase up to the same level of the density of the gas at the outer edge. The feedback on the gas from the drifting dust grains is not negligible, in such a dusty region. We carried out two-dimensional two-fluid (gas and dust) hydrodynamic simulations. We found that when the radial flow of the dust grains across the gap is halted, a broad ring of the dust grains can be formed because of the dust feedback and the diffusion of the dust grains. The minimum mass of the planet to form the broad dust ring is consistent with the pebble-isolation mass, in the parameter range of our simulations. The broad ring of the dust grains is good environment for the formation of the protoplanetary solid core. If the ring is formed in the disk around the sun-like star at $\sim 2\ \mbox{AU}$, a massive solid core ($\sim 50M_{\oplus}$) can be formed within the ring, which may be connected to the formation of Hot Jupiters holding a massive solid core such as HD 149026b. In the disk of the dwarf star, a number of Earth-sized planets can be formed within the dust ring around $\sim 0.5\ \mbox{AU}$, which potentially explain the planet system made of multiple Earth-sized planets around the dwarf star such as TRAPPIST-1.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

We investigate the simultaneous evolution of dust and gas density profiles at a radial pressure bump located in a protoplanetary disk. If dust particles are treated as test particles, a radial pressure bump traps dust particles that drift radially inward. As the dust particles become more concentrated at the gas pressure bump, however, the drag force from dust to gas (back-reaction), which is ignored in a test-particle approach, deforms the pressure bump. We find that the pressure bump is completely deformed by the back-reaction when the dust-to-gas mass ratio reaches $\sim 1$ for a slower bump restoration. The direct gravitational instability of dust particles is inhibited by the bump destruction. In the dust-enriched region, the radial pressure support becomes $\sim 10-100$ times lower than the global value set initially. Although the pressure bump is a favorable place for streaming instability (SI), the flattened pressure gradient inhibits SI from forming large particle clumps corresponding to $100-1000$ km sized bodies, which has been previously proposed. If SI occurs there, the dust clumps formed would be $10-100$ times smaller, that is, of about $1 - 100$ km.