No Arabic abstract
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.
The mechanism of angular momentum transport in protoplanetary disks is fundamental to understand the distributions of gas and dust in the disks. The unprecedented, high spatial resolution ALMA observations taken toward HL Tau and subsequent radiative transfer modeling reveal that a high degree of dust settling is currently achieved at the outer part of the HL Tau disk. Previous observations however suggest a high disk accretion rate onto the central star. This configuration is not necessarily intuitive in the framework of the conventional viscous disk model, since efficient accretion generally requires a high level of turbulence, which can suppress dust settling considerably. We develop a simplified, semi-analytical disk model to examine under what condition these two properties can be realized in a single model. Recent, non-ideal MHD simulations are utilized to realistically model the angular momentum transport both radially via MHD turbulence and vertically via magnetically induced disk winds. We find that the HL Tau disk configuration can be reproduced well when disk winds are properly taken into account. While the resulting disk properties are likely consistent with other observational results, such an ideal situation can be established only if the plasma $beta$ at the disk midplane is $beta_0 simeq 2 times 10^4$ under the assumption of steady accretion. Equivalently, the vertical magnetic flux at 100 au is about 0.2 mG. More detailed modeling is needed to fully identify the origin of the disk accretion and quantitatively examine plausible mechanisms behind the observed gap structures in the HL Tau disk.
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 first long-baseline ALMA campaign resolved the disk around the young star HL Tau into a number of axisymmetric bright and dark rings. Despite the very young age of HL Tau these structures have been interpreted as signatures for the presence of (proto)planets. The ALMA images triggered numerous theoretical studies based on disk-planet interactions, magnetically driven disk structures, and grain evolution. Of special interest are the inner parts of disks, where terrestrial planets are expected to form. However, the emission from these regions in HL Tau turned out to be optically thick at all ALMA wavelengths, preventing the derivation of surface density profiles and grain size distributions. Here, we present the most sensitive images of HL Tau obtained to date with the Karl G. Jansky Very Large Array at 7.0 mm wavelength with a spatial resolution comparable to the ALMA images. At this long wavelength the dust emission from HL Tau is optically thin, allowing a comprehensive study of the inner disk. We obtain a total disk dust mass of 0.001 - 0.003 Msun, depending on the assumed opacity and disk temperature. Our optically thin data also indicate fast grain growth, fragmentation, and formation of dense clumps in the inner densest parts of the disk. Our results suggest that the HL Tau disk may be actually in a very early stage of planetary formation, with planets not already formed in the gaps but in the process of future formation in the bright rings.
Recent observations of HL Tau revealed remarkably detailed structure within the systems circumstellar disc. A range of hypotheses have been proposed to explain the morphology, including, e.g., planet-disc interactions, condensation fronts, and secular gravitational instabilities. While embedded planets seem to be able to explain some of the major structure in the disc through interactions with gas and dust, the substructure, such as low-contrast rings and bands, are not so easily reproduced. Here, we show that dynamical interactions between three planets (only two of which are modelled) and an initial population of large planetesimals can potentially explain both the major and minor banded features within the system. In this context, the small grains, which are coupled to the gas and reveal the disc morphology, are produced by the collisional evolution of the newly-formed planetesimals, which are ubiquitous in the system and are decoupled from the gas.
We conducted a detailed radiative transfer modeling of the dust emission from the circumstellar disk around HL Tau. The goal of our study is to derive the surface density profile of the inner disk and its structure. In addition to the Atacama Large Millimeter/submillimeter Array images at Band 3 (2.9mm), Band 6 (1.3mm), and Band 7 (0.87mm), the most recent Karl G. Jansky Very Large Array (VLA) observations at 7mm were included in the analysis. A simulated annealing algorithm was invoked to search for the optimum model. The radiative transfer analysis demonstrates that most radial components (i.e., >6AU) of the disk become optically thin at a wavelength of 7mm, which allows us to constrain, for the first time, the dust density distribution in the inner region of the disk. We found that a homogeneous grain size distribution is not sufficient to explain the observed images at different wavelengths simultaneously, while models with a shallower grain size distribution in the inner disk work well. We found clear evidence that larger grains are trapped in the first bright ring. Our results imply that dust evolution has already taken place in the disk at a relatively young (i.e., ~1Myr) age. We compared the midplane temperature distribution, optical depth, and properties of various dust rings with those reported previously. Using the Toomre parameter, we briefly discussed the gravitational instability as a potential mechanism for the origin of the dust clump detected in the first bright ring via the VLA observations.