When two planets are born in a protoplanetary disk, they may enter into a mean-motion resonance as a consequence of the convergent planetary migration. The formation of mean-motion resonances is important for understanding how the planetary systems are shaped in the disk environments. Motivated by recent progress in the comprehension of the migration of partial gap-opening planets, we have investigated the orbital evolution of the planet pairs in a wide range of masses and disk properties with the aim to find out when the resonance capture is likely to happen. Using the formula for the migration timescale of a gap-opening planet developed in our previous work, we have derived a simple criterion that allows us to predict when the migration will be convergent (divergent). Further, we have verified the criterion using two-dimensional hydrodynamic simulations. We have found that the resonant pair of planets formed at the early phase of evolution can depart from the resonance at later times because the migration speed of the outer planet slows down due to the gap formation. Moreover, adopting our formula of the migration timescale, we have also carried out three-body simulations, which confirm the results of hydrodynamic simulations. Finally, we have compared our predictions with the observations, selecting a sample of known two-planet systems.

In the conventional view of type II migration, a giant planet migrates inward in the viscous velocity of the accretion disc in the so-call disc-dominate case. Recent hydrodynamic simulations, however, showed that planets migrate with velocities much faster than the viscous one in massive discs. Such fast migration cannot be explained by the conventional picture. Scardoni et al. (2020) has recently argued this new picture. By carrying out similar hydrodynamic simulations, they found that the migration velocity slows down with time and eventually reaches the prediction by the conventional theory. They interpreted the fast migration as an initial transient one and concluded that the conventional type II migration is realised after the transient phase. We show that the migration velocities obtained by Scardoni et al. (2020) are consistent with the previous simulations even in the transient phase that they proposed. We also find that the transient fast migration proposed by Scardoni et al. (2020) is well described by a new model of Kanagawa et al. (2018). The new model can appropriately describe significant inward migration during the initial transient phase that Scardoni et al. (2020) termed. Hence, we conclude that the time-variation of the transient migration velocity is due to the changes of the orbital radius of the planet and its background surface density during the migration.

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

A planet is formed within a protoplanetary disk. Recent observations have revealed substructures such as gaps and rings, which may indicate forming planets within the disk. Due to disk--planet interaction, the planet migrates within the disk, which can affect a shape of the planet-induced gap. In this paper, we investigate effects of fast inward migration of the planet on the gap shape, by carrying out hydrodynamic simulations. We found that when the migration timescale is shorter than the timescale of the gap-opening, the orbital radius is shifted inward as compared to the radial location of the gap. We also found a scaling relation between the radial shift of the locations of the planet and the gap as a function of the ratio of the timescale of the migration and gap-opening. Our scaling relation also enables us to constrain the gas surface density and the viscosity when the gap and the planet are observed. Moreover, we also found the scaling relation between the location of the secondary gap and the aspect ratio. By combining the radial shift and the secondary gap, we may constrain the physical condition of the planet formation and how the planet evolves in the protoplanetary disk, from the observational morphology.

A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolution of planets within the protoplanetary disk and after the disk dispersal. Before the disk dispersal, our N-body simulations include both migration and mass-growth of the planet coupled with evolution of the disk. By varying the disk parameters, we produce a variety of widely-separated planetary systems consisting of three super-Jupiters at the end of disk dispersal. We found the outer planet is more massive than the inner one, and the migration of the innermost planet is inefficient due to the accretion of outer planet(s). We also showed how the final configuration and the final planetary mass depend on disk parameters. The migration is found to be convergent and no planet-pair has a period ratio less than 2. After the disk dispersal, we switch to pure gravitational N-body simulations and integrate the orbits up to 10 Gyr. Most simulated systems remain stable for at least 10 Gyr. We discuss implications of our result in terms of the observed widely-separated planetary systems HR 8799 and PDS 70.

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

We report the detection of an excess in dust continuum emission at 233~GHz (1.3~mm in wavelength) in the protoplanetary disk around TW~Hya revealed through high-sensitivity observations at $\sim$3~au resolution with the Atacama Large Millimeter/submillimeter Array (ALMA). The sensitivity of the 233~GHz image has been improved by a factor of 3 with regard to that of our previous cycle 3 observations. The overall structure is mostly axisymmetric, and there are apparent gaps at 25 and 41 au as previously reported. The most remarkable new finding is a few au-scale excess emission in the south-west part of the protoplanetary disk. The excess emission is located at 52 au from the disk center and is 1.5 times brighter than the surrounding protoplanetary disk at a significance of 12$\sigma$. We performed a visibility fitting to the extracted emission after subtracting the axisymmetric protoplanetary disk emission and found that the inferred size and the total flux density of the excess emission are 4.4$\times$1.0~au and 250~$\mu$Jy, respectively. The dust mass of the excess emission corresponds to 0.03~$M_\oplus$ if a dust temperature of 18~K is assumed. Since the excess emission can also be marginally identified in the Band 7 image at almost the same position, the feature is unlikely to be a background source. The excess emission can be explained by a dust clump accumulated in a small elongated vortex or a massive circumplanetary disk around a Neptune mass forming-planet.

The planet migration due to the disk--planet interaction is one of the most important processes to determine the architecture of planetary systems. A sufficiently massive planet forms a density gap and migrates together with the gap. By carrying out two-dimensional and two-fluid (gas and dust grains) hydrodynamic simulations, we investigated the effects of the dust feedback on the migration of the gap-opening planet, which was not considered in previous studies. We found that the gas surface density at the outer edge of the gap becomes smaller due to the dust feedback, and thus the torque exerted from the outer disk decreases. This mechanism becomes effective as the gap becomes wider and deeper. In particular, when the mass of the planet is Jupiter-size and turbulent viscosity is $\alpha = 3\times 10^{-4}$, the planet can migrate outward due to the reduction of the torque exerted from the outer disk. Even for a smaller planet, the migration becomes significantly slow down. This termination of the inward migration triggered by the dust feedback may explain why ring and gap structures can be frequently observed within the protoplanetary disks.

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.

The mass and semimajor axis distribution of gas giants in exoplanetary systems obtained by radial velocity surveys shows that super-jupiter-mass planets are piled up at > 1 au, while jupiter/sub-jupiter-mass planets are broadly distributed from ~0.03 au to beyond 1 au. This feature has not been explained by theoretical predictions. In order to reconcile this inconsistency, we investigate evolution of gas giants with a new type II migration formula by Kanagawa et al. (2018), by comparing the migration, growth timescales of gas giants, and disk lifetime and by population synthesis simulation. While the classical migration model assumes that a gas giant opens up a clear gap in the protoplanetary disk and the planet migration tied to the disk gas accretion, recent high-resolution simulations show that the migration of gap-opening planets is decoupled from the disk gas accretion and Kanagawa et al. (2018) proposed that type II migration speed is no other than type I migration speed with the reduced disk gas surface density in the gap. We show that with this new formula, type II migration is significantly reduced for super-jupiter-mass planets, if the disk accretion is driven by the disk wind as suggested by recent MHD simulations. Population synthesis simulations show that super-jupiter-mass planets remain at > 1 au without any additional ingredient such as disk photoevaporation. Therefore, the mystery of the pile-up of gas giants at > 1 au would be theoretically solved, if the new formula is confirmed and wind-driven disk accretion dominates.

A large planet orbiting a star in a protoplanetary disk opens a density gap along its orbit due to the strong disk-planet interaction and migrates with the gap in the disk. It is expected that in the ideal case, a gap-opening planet migrates at the viscous drift speed, which is referred to as type II migration. However, recent hydrodynamic simulations have shown that in general, the gap-opening planet is not locked to the viscous disk evolution. A new physical model is required to explain the migration speed of such a planet. For this reason, we re-examined the migration of a planet in the disk, by carrying out the two-dimensional hydrodynamic simulations in a wide parameter range. We have found that the torque exerted on the gap-opening planet depends on the surface density at the bottom of the gap. The planet migration slows down as the surface density of the bottom of the gap decreases. Using the gap model developed in our previous studies, we have constructed an empirical formula of the migration speed of the gap-opening planets, which is consistent with the results given by the hydrodynamic simulations performed by us and other researchers. Our model easily explains why the migration speed of the gap-opening planets can be faster than the viscous gas drift speed. It can also predict the planet mass at which the type I migration is no longer adequate due to the gap development in the disk, providing a gap formation criterion based on planetary migration.

A giant planet embedded in a protoplanetary disk creates a gap. This process is important for both theory and observations. Using results of a survey for a wide parameter range with two-dimensional hydrodynamic simulations, we constructed an empirical formula for the gap structure (i.e., the radial surface density distribution), which can reproduce the gap width and depth obtained by two-dimensional simulations. This formula enables us to judge whether an observed gap is likely to be caused by an embedded planet or not. The propagation of waves launched by the planet is closely connected to the gap structure. It makes the gap wider and shallower as compared with the case where an instantaneous wave damping is assumed. The hydrodynamic simulations shows that the waves do not decay immediately at the launching point of waves, even when the planet is as massive as Jupiter. Based on the results of hydrodynamic simulations, we also obtained an empirical model of wave propagation and damping for the cases of deep gaps. The one-dimensional gap model with our wave propagation model is able to well reproduce the gap structures in hydrodynamic simulations. In the case of a Jupiter-mass planet, we also found that the waves with smaller wavenumber (e.g., $m=2$) are excited and transport the angular momentum to the location far away from the planet. The wave with $m=2$ is closely related with a secondary wave launched by the site opposite from the planet.

The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8" (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2" (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_\odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_\odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.

The total amount of dust (or "metallicity") and the dust distribution in protoplanetary disks are crucial for planet formation. Dust grains radially drift due to gas--dust friction, and the gas is affected by the feedback from dust grains. We investigate the effects of the feedback from dust grains on the viscous evolution of the gas, taking into account the vertical dust settling. The feedback from the grains pushes the gas outward. When the grains are small and the dust-to-gas mass ratio is much smaller than unity, the radial drift velocity is reduced by the feedback effect but the gas still drifts inward. When the grains are sufficiently large or piled-up, the feedback is so effective that forces the gas flows outward. Although the dust feedback is affected by dust settling, we found that the 2D approximation reasonably reproduces the vertical averaged flux of gas and dust. We also performed the 2D two-fluid hydrodynamic simulations to examine the effect of the feedback from the grains on the evolution of the gas disk. We show that when the feedback is effective, the gas flows outward and the gas density at the region within $\sim 10\ \mbox{AU}$ is significantly depleted. As a result, the dust-to-gas mass ratio at the inner radii may significantly excess unity, providing the environment where planetesimals are easily formed via, e.g., streaming instability. We also show that a simplified 1D model well reproduces the results of the 2D two-fluid simulations, which would be useful for future studies.

We report $\sim$3 au resolution imaging observations of the protoplanetary disk around TW Hya at 145 and 233 GHz with the Atacama Large Millimeter/Submillimeter Array. Our observations revealed two deep gaps ($\sim$25--50 \%) at 22 and 37~au and shallower gaps (a few \%) at 6, 28, and 44~au, as recently reported by Andrews et al. (2016). The central hole with a radius of $\sim3$~au was also marginally resolved. The most remarkable finding is that the spectral index $\alpha (R)$ between bands 4 and 6 peaks at the 22~au gap. The derived power-law index of the dust opacity $\beta (R)$ is $\sim1.7$ at the 22~au gap and decreases toward the disk center to $\sim0$. The most prominent gap at 22~au could be caused by the gravitational interaction between the disk and an unseen planet with a mass of $\lesssim$1.5 $M_\mathrm{Neptune}$ although other origins may be possible. The planet-induced gap is supported by the fact that $\beta (R)$ is enhanced at the 22~au gap, indicating a deficit of $\sim$mm-sized grains within the gap due to dust filtration by a planet.

A giant planet creates a gap in a protoplanetary disk, which might explain the observed gaps in protoplanetary disks. The width and depth of the gaps depend on the planet mass and disk properties. We have performed two--dimensional hydrodynamic simulations for various planet masses, disk aspect ratios and viscosities, to obtain an empirical formula for the gap width. The gap width is proportional to the square root of the planet mass, -3/4 power of the disk aspect ratio and -1/4 power of the viscosity. This empirical formula enables us to estimate the mass of a planet embedded in the disk from the width of an observed gap. We have applied the empirical formula for the gap width to the disk around HL~Tau, assuming that each gap observed by ALMA observations is produced by planets, and discussed the planet masses within the gaps. The estimate of planet masses from the gap widths is less affected by the observational resolution and dust filtration than that from the gap depth.

We report the first detection of a gap and a ring in 336 GHz dust continuum emission from the protoplanetary disk around TW Hya, using the Atacama Large Millimeter/Submillimeter Array (ALMA). The gap and ring are located at around 25 and 41 au from the central star, respectively, and are associated with the CO snow line at ~30 au. The gap has a radial width of less than 15 au and a mass deficit of more than 23%, taking into account that the observations are limited to an angular resolution of ~15 au. In addition, the 13CO and C18O J = 3 - 2 lines show a decrement in CO line emission throughout the disk, down to ~10 au, indicating a freeze-out of gas-phase CO onto grain surfaces and possible subsequent surface reactions to form larger molecules. The observed gap could be caused by gravitational interaction between the disk gas and a planet with a mass less than super-Neptune (2M_{Neptune}), or could be the result of the destruction of large dust aggregates due to the sintering of CO ice.

Coronagraphic imagery of the circumstellar disk around HD 169142 in H-band polarized intensity (PI) with Subaru/HiCIAO is presented. The emission scattered by dust particles at the disk surface in 0.2" <= r <= 1.2", or 29 <= r <= 174 AU, is successfully detected. The azimuthally-averaged radial profile of the PI shows a double power-law distribution, in which the PIs in r=29-52 AU and r=81.2-145 AU respectively show r^{-3}-dependence. These two power-law regions are connected smoothly with a transition zone (TZ), exhibiting an apparent gap in r=40-70 AU. The PI in the inner power-law region shows a deep minimum whose location seems to coincide with the point source at \lambda = 7 mm. This can be regarded as another sign of a protoplanet in TZ. The observed radial profile of the PI is reproduced by a minimally flaring disk with an irregular surface density distribution or with an irregular temperature distribution or with the combination of both. The depletion factor of surface density in the inner power-law region (r< 50 AU) is derived to be <= 0.16 from a simple model calculation. The obtained PI image also shows small scale asymmetries in the outer power-law region. Possible origins for these asymmetries include corrugation of the scattering surface in the outer region, and shadowing effect by a puffed up structure in the inner power-law region.

A giant planet embedded in a protoplanetary disk forms a gap. An analytic relationship among the gap depth, planet mass $M_{p}$, disk aspect ratio $h_p$, and viscosity $\alpha$ has been found recently, and the gap depth can be written in terms of a single parameter $K= (M_{p}/M_{\ast})^2 h_p^{-5} \alpha^{-1}$. We discuss how observed gap features can be used to constrain the disk and/or planet parameters based on the analytic formula for the gap depth. The constraint on the disk aspect ratio is critical in determining the planet mass so the combination of the observations of the temperature and the image can provide a constraint on the planet mass. We apply the formula for the gap depth to observations of HL~Tau and HD~169142. In the case of HL~Tau, we propose that a planet with $\gtrsim 0.3$ is responsible for the observed gap at $30$~AU from the central star based on the estimate that the gap depth is $\lesssim 1/3$. In the case of HD~169142, the planet mass that causes the gap structure recently found by VLA is $\gtrsim 0.4 M_J$. We also argue that the spiral structure, if observed, can be used to estimate the lower limit of the disk aspect ratio and the planet mass.

The gap formation induced by a giant planet is important in the evolution of the planet and the protoplanetary disc. We examine the gap formation by a planet with a new formulation of one-dimensional viscous discs which takes into account the deviation from Keplerian disc rotation due to the steep gradient of the surface density. This formulation enables us to naturally include the Rayleigh stable condition for the disc rotation. It is found that the derivation from Keplerian disc rotation promotes the radial angular momentum transfer and makes the gap shallower than in the Keplerian case. For deep gaps, this shallowing effect becomes significant due to the Rayleigh condition. In our model, we also take into account the propagation of the density waves excited by the planet, which widens the range of the angular momentum deposition to the disc. The effect of the wave propagation makes the gap wider and shallower than the case with instantaneous wave damping. With these shallowing effects, our one-dimensional gap model is consistent with the recent hydrodynamic simulations.

We examine the gas circulation near a gap opened by a giant planet in a protoplanetary disk. We show with high resolution 3D simulations that the gas flows into the gap at high altitude over the mid-plane, at a rate dependent on viscosity. We explain this observation with a simple conceptual model. From this model we derive an estimate of the amount of gas flowing into a gap opened by a planet with Hill radius comparable to the scale-height of a layered disk (i. e. a disk with viscous upper layer and inviscid midplane). Our estimate agrees with modern MRI simulations(Gressel et al., 2013). We conclude that gap opening in a layered disk can not slow down significantly the runaway gas accretion of Saturn to Jupiter-mass planets.

In the core accretion model for the formation of gas giant planets, runaway gas accretion onto a core is the primary requisite, triggered when the core mass reaches a critical value. The recently revealed wide diversity of the extrasolar giant planets suggests the necessity to further the understanding of the conditions resulting in the critical core mass that initiates runaway accretion. We study the internal structure of protoplanets under hydrostatic and thermal equilibria represented in terms of a polytropic equation of state to investigate what factors determine and affect the critical core mass. We find that the protoplanets, embedded in protoplanetary disks, have the same configuration as red giants, characterized by the envelope of the centrally-condensed type solution. Applying the theory of stellar structure with homology invariants, we demonstrate that there are three types of criteria for the critical core mass depending on the stiffness of polytrope and the nature of outer boundary condition. For the stiff polytropes of index $N \le 3$ with the Bondi radius as the outer boundary, the criterion governing the critical core mass occurs at the surface. For stiff polytropes with the Hill outer boundary and for soft polytropes of $N>3$, this criterion acts at the bottom of gaseous envelope. Further, we elucidate the roles and effects of coexistent radiative and convective zones in the envelope of critical core mass. Based on the results, we discuss the relevance of Bondi and Hill surface conditions and explore the parameter dependences of critical core mass.