[en] The disk surrounding SAO 206462, an 8 Myr old Herbig Ae star, has recently been reported to exhibit spiral arms, an asymmetric dust continuum, and a dust-depleted inner cavity. By carrying out two-dimensional, two-fluid hydrodynamic calculations, we find that a planetary-mass companion located at the outer disk could be responsible for these observed structures. In this model, the planet excites primary and secondary arms interior to its orbit. It also carves a gap and generates a local pressure bump at the inner gap edge where a vortex forms through Rossby wave instability. The vortex traps radially drifting dust particles, forming a dust-depleted cavity in the inner disk. We propose that the vortex is responsible for the brightest southwestern peak seen in infrared scattered light and sub-millimeter dust continuum emission. In particular, it is possible that the scattered light is boosted as one of the spiral arms passes through the high density vortex region, although the vortex alone may be able to explain the peak. We suggest that a planetary companion with a mass of 10–15 $M_J$ is orbiting SAO 206462 at 100–120 au. Monitoring of the brightest peak over the next few years will help reveal its origin because the spiral arms and vortex will show distinguishable displacement.

[en] We carry out two-fluid, two-dimensional global hydrodynamic simulations to test whether protostellar infall can trigger the Rossby wave instability (RWI) in protoplanetry disks. Our results show that infall can trigger the RWI and generate vortices near the outer edge of the mass landing on the disk (i.e., centrifugal radius). We find that the RWI is triggered under a variety of conditions, although the details depend on the disk parameters and the infall pattern. The common key feature of triggering the RWI is the steep radial gradient of the azimuthal velocity induced by the local increase in density at the outer edge of the infall region. Vortices form when the instability enters the nonlinear regime. In our standard model where self-gravity is neglected, vortices merge together to a single vortex within ∼20 local orbital times, and the merged vortex survives for the remaining duration of the calculation (>170 local orbital times). The vortex takes part in outward angular momentum transport, with a Reynolds stress of ≲10"−"2. Our two-fluid calculations show that vortices efficiently trap dust particles with stopping times of the order of the orbital time, locally enhancing the dust to gas ratio for particles of the appropriate size by a factor of ∼40 in our standard model. When self-gravity is considered, however, vortices tend to be impeded from merging and may eventually dissipate. We conclude it may well be that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation

[en] We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3R_{\mathrm{p}}⩽<!--\; ???\; -->R⩽<!--\; ???\; -->0.7R_{\mathrm{p}}$, where R _p denotes the semimajor axis of the planetary orbit (assumed to be 5 au), the estimated vertical diffusion rate associated with the turbulence is characterized by ${\mathrm{\alpha <!--\; ??\; -->}}_{\mathrm{d}\mathrm{i}\mathrm{f}\mathrm{f}}\sim <!--\; ???\; -->(0.2\text{--}1.2)\times <!--\; ??\; -->{10}^{-<!--\; ???\; -->2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr timescales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

[en] We present results from a suite of three-dimensional global hydrodynamic simulations that shows that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport and vertical mixing at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter α ∼ 5 × 10⁻⁴ in our reference model). The instability is found to operate in a wide range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity ν ≤ 10⁻⁵. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own solar system) and FU Orionis outburst models. Further work is required to determine the nature of the instability and to evaluate its observational consequences in physically more complete disk models than we have considered in this paper.

[en] Besides the spirals induced by the Lindblad resonances, planets can generate a family of tightly wound spirals through buoyancy resonances. The excitation of buoyancy resonances depends on the thermal relaxation timescale of the gas. By computing timescales of various processes associated with thermal relaxation, namely, radiation, diffusion, and gas–dust collision, we show that the thermal relaxation in protoplanetary disks’ surface layers (Z/R ≳ 0.1) and outer disks (R ≳ 100 au) is limited by infrequent gas–dust collisions. The use of the isothermal equation of state or rapid cooling, common in protoplanetary disk simulations, is therefore not justified. Using three-dimensional hydrodynamic simulations, we show that the collision-limited slow thermal relaxation provides favorable conditions for buoyancy resonances to develop. Buoyancy resonances produce predominantly vertical motions, whose magnitude at the ¹²CO emission surface is of the order of 100 m s⁻¹ for Jovian-mass planets, sufficiently large to detect using molecular line observations with ALMA. We generate synthetic observations and describe characteristic features of buoyancy resonances in Keplerian-subtracted moment maps and velocity channel maps. Based on the morphology and magnitude of the perturbation, we propose that the tightly wound spirals observed in TW Hya could be driven by a (sub-)Jovian-mass planet at 90 au. We discuss how non-Keplerian motions driven by buoyancy resonances can be distinguished from those driven by other origins. We argue that observations of multiple lines tracing different heights, with sufficiently high spatial/spectral resolution and sensitivity to separate the emission arising from the near and far sides of the disk, will help constrain the origin of non-Keplerian motions.

[en] We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting for ∼< 20% of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that the initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the M-dot-R_wall plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of the order of 10 AU

[en] We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the ^dead zone⁾. We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models, which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, R ∼< 0.5 AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

[en] We improve on our previous treatments of the long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e., PdV work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead zone. The GI-induced spiral waves propagate well inside of the gravitationally unstable region before they trigger outbursts at R ≲ 1 AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local α treatments for GI. In our standard model where zero dead-zone residual viscosity (α_rd) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as 0.01 in terms of the generic α parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite α_rd triggers thermally driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when α_rd = 0. The inclusion of non-zero residual viscosity in the dead zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a 'quiescent' accretion phase with M-dot _acc∼10⁻⁸--10⁻⁷ M_⊙ yr⁻¹ over 60% of the time and spend less than 15% of the infall phase in accretion outbursts. While our models indicate that episodic mass accretion during protostellar evolution can qualitatively help explain the low accretion luminosities seen in most low-mass protostars, detailed tests of the mechanism will require model calculations for a range of protostellar masses with some constraint on the initial core angular momentum, which affects the length of time spent in a quasi-steady disk accretion phase.

[en] As host to two accreting planets, PDS 70 provides a unique opportunity to probe the chemical complexity of atmosphere-forming material. We present ALMA Band 6 observations of the PDS 70 disk and report the first chemical inventory of the system. With a spatial resolution of 0.″4–0.″5 (∼50 au), 12 species are detected, including CO isotopologs and formaldehyde, small hydrocarbons, HCN and HCO⁺ isotopologs, and S-bearing molecules. SO and CH₃OH are not detected. All lines show a large cavity at the center of the disk, indicative of the deep gap carved by the massive planets. The radial profiles of the line emission are compared to the (sub)millimeter continuum and infrared scattered light intensity profiles. Different molecular transitions peak at different radii, revealing the complex interplay between density, temperature, and chemistry in setting molecular abundances. Column densities and optical depth profiles are derived for all detected molecules, and upper limits obtained for the nondetections. Excitation temperature is obtained for H₂CO. Deuteration and nitrogen fractionation profiles from the hydrocyanide lines show radially increasing fractionation levels. Comparison of the disk chemical inventory to grids of chemical models from the literature strongly suggests a disk molecular layer hosting a carbon-to-oxygen ratio C/O > 1, thus providing for the first time compelling evidence of planets actively accreting high C/O ratio gas at present time.

[en] We have carried out the first very long baseline interferometry (VLBI) imaging of a 44 GHz class I methanol maser (7₀-6₁ A ⁺) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151–1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds × 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than ≈ 650 km corresponding to 100 Mλ in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ∼5 mas × 2 mas, which corresponds to the linear size of ∼15 AU × 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ∼3.5 × 10⁸ to 1.0 × 10¹⁰ K, which are higher than the estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ∼50 mas. The 44 GHz class I methanol maser in IRAS 18151–1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser