Lopsided dust rings in transition disks

Context. Particle trapping in local or global pressure maxima in protoplanetary disks is one of the new paradigms in the theory of the first stages of planet formation. However, finding observational evidence for this effect is not easy. Recent work suggests that the large ring-shaped outer disks observed in transition disk sources may in fact be lopsided and constitute large banana-shaped vortices. Aims. We wish to investigate how effective dust can accumulate along the azimuthal direction. We also want to find out if the size- sorting resulting from this can produce a detectable signatures at millimeter wavelengths. Methods. To keep the numerical cost under control we develop a 1+1D method in which the azimuthal variations are treated sepa- rately from the radial ones. The azimuthal structure is calculated analytically for a steady-state between mixing and azimuthal drift. We derive equilibration time scales and compare the analytical solutions to time-dependent numerical simulations. Results. We find that weak, but long-lived azimuthal density gradients in the gas can induce very strong azimuthal accumulations of dust. The strength of the accumulations depends on the P\'eclet number, which is the relative importance of advection and diffusion. We apply our model to transition disks and our simulated observations show that this effect would be easily observable with ALMA and in principle allows to put constraints on the strength of turbulence and the local gas density.Comment: 4 pages, 4 figures, accepted for publication in A&A Letter

Trapping dust particles in the outer regions of protoplanetary disks

Aims. We attempt to explain grain growth to mm sized particles and their retention in the outer regions of protoplanetary disks, as observed at sub-mm and mm wavelengths, by investigating whether strong inhomogeneities in the gas density profiles can decelerate excessive radial drift and help the dust particles to grow. Methods. We use coagulation/fragmentation and disk-structure models, to simulate the evolution of dust in a bumpy surface density profile, which we mimic with a sinusoidal disturbance. For different values of the amplitude and length scale of the bumps, we investigate the ability of this model to produce and retain large particles on million-year timescales. In addition, we compare the pressure inhomogeneities considered in this work with the pressure profiles that come from magnetorotational instability. Using the Common Astronomy Software Applications ALMA simulator, we study whether there are observational signatures of these pressure inhomogeneities that can be seen with ALMA. Results. We present the conditions required to trap dust particles and the corresponding calculations predicting the spectral slope in the mm-wavelength range, to compare with current observations. Finally, we present simulated images using different antenna configurations of ALMA at different frequencies, to show that the ring structures will be detectable at the distances of either the Taurus Auriga or Ophiucus star-forming regions

Scattered light images of spiral arms in marginally gravitationally unstable discs with an embedded planet

Scattered light images of transition discs in the near-infrared often show non-axisymmetric structures in the form of wide-open spiral arms in addition to their characteristic low-opacity inner gap region. We study self-gravitating discs and investigate the influence of gravitational instability on the shape and contrast of spiral arms induced by planet-disc interactions. Two-dimensional non-isothermal hydrodynamical simulations including viscous heating and a cooling prescription are combined with three-dimensional dust continuum radiative transfer models for direct comparison to observations. We find that the resulting contrast between the spirals and the surrounding disc in scattered light is by far higher for pressure scale height variations, i.e. thermal perturbations, than for pure surface density variations. Self-gravity effects suppress any vortex modes and tend to reduce the opening angle of planet-induced spirals, making them more tightly wound. If the disc is only marginally gravitationally stable with a Toomre parameter around unity, an embedded massive planet (planet-to-star mass ratio of $10^{-2}$) can trigger gravitational instability in the outer disc. The spirals created by this instability and the density waves launched by the planet can overlap resulting in large-scale, more open spiral arms in the outer disc. The contrast of these spirals is well above the detection limit of current telescopes.Comment: Accepted for publication in MNRAS; 13 pages, 8 figure

Millimetre spectral indices of transition disks and their relation to the cavity radius

Transition disks are protoplanetary disks with inner depleted dust cavities and excellent candidates to investigate the dust evolution under the existence of a pressure bump. A pressure bump at the outer edge of the cavity allows dust grains from the outer regions to stop their rapid inward migration towards the star and efficiently grow to millimetre sizes. Dynamical interactions with planet(s) have been one of the most exciting theories to explain the clearing of the inner disk. We look for evidence of the presence of millimetre dust particles in transition disks by measuring their spectral index with new and available photometric data. We investigate the influence of the size of the dust depleted cavity on the disk integrated millimetre spectral index. We present the 3mm photometric observations carried out with PdBI of four transition disks: LkHa330, UXTauA, LRLL31, and LRLL67. We use available values of their fluxes at 345GHz to calculate their spectral index, as well as the spectral index for a sample of twenty transition disks. We compare the observations with two kind of models. In the first set of models, we consider coagulation and fragmentation of dust in a disk in which a cavity is formed by a massive planet located at different positions. The second set of models assumes disks with truncated inner parts at different radius and with power-law dust size distributions, where the maximum size of grains is calculated considering turbulence as the source of destructive collisions. We show that the integrated spectral index is higher for transition disks than for regular protoplanetary disks. For transition disks, the probability that the measured spectral index is positively correlated with the cavity radius is 95%. High angular resolution imaging of transition disks is needed to distinguish between the dust trapping scenario and the truncated disk case.Comment: Accepted for publication in A&A, including language editio

Testing particle trapping in transition disks with ALMA

We present new Atacama Large Millimeter/submillimeter Array (ALMA) continuum observations at 336GHz of two transition disks, SR21 and HD135344B. In combination with previous ALMA observations from Cycle 0 at 689GHz, we compare the visibility profiles at the two frequencies and calculate the spectral index ($\alpha_{\rm{mm}}$). The observations of SR21 show a clear shift in the visibility nulls, indicating radial variations of the inner edge of the cavity at the two wavelengths. Notable radial variations of the spectral index are also detected for SR21 with values of $\alpha_{\rm{mm}}{\sim}3.8-4.2$ in the inner region ($r<35$ AU) and $\alpha_{\rm{mm}}{\sim}2.6-3.0$ outside. An axisymmetric ring (which we call the ring model) or a ring with the addition of an azimuthal Gaussian profile, for mimicking a vortex structure (which we call the vortex model), is assumed for fitting the disk morphology. For SR21, the ring model better fits the emission at 336GHz, conversely the vortex model better fits the 689GHz emission. For HD135344B, neither a significant shift in the null of the visibilities nor radial variations of $\alpha_{\rm{mm}}$ are detected. Furthermore, for HD135344B, the vortex model fits both frequencies better than the ring model. However, the azimuthal extent of the vortex increases with wavelength, contrary to model predictions for particle trapping by anticyclonic vortices. For both disks, the azimuthal variations of $\alpha_{\rm{mm}}$ remain uncertain to confirm azimuthal trapping. The comparison of the current data with a generic model of dust evolution that includes planet-disk interaction suggests that particles in the outer disk of SR21 have grown to millimetre sizes and have accumulated in a radial pressure bump, whereas with the current resolution there is not clear evidence of radial trapping in HD135344B, although it cannot be excluded either.Comment: Minor changes after language edition. Accepted for publication in A&A (abstract slightly shortened for arXiv

A Multi-Wavelength Analysis of Dust and Gas in the SR 24S Transition Disk

We present new Atacama Large Millimeter/sub-millimeter Array (ALMA) 1.3 mm continuum observations of the SR 24S transition disk with an angular resolution $\lesssim0.18"$ (12 au radius). We perform a multi-wavelength investigation by combining new data with previous ALMA data at 0.45 mm. The visibilities and images of the continuum emission at the two wavelengths are well characterized by a ring-like emission. Visibility modeling finds that the ring-like emission is narrower at longer wavelengths, in good agreement with models of dust trapping in pressure bumps, although there are complex residuals that suggest potentially asymmetric structures. The 0.45 mm emission has a shallower profile inside the central cavity than the 1.3 mm emission. In addition, we find that the $^{13}$CO and C$^{18}$O (J=2-1) emission peaks at the center of the continuum cavity. We do not detect either continuum or gas emission from the northern companion to this system (SR 24N), which is itself a binary system. The upper limit for the dust disk mass of SR 24N is $\lesssim 0.12\,M_{\bigoplus}$, which gives a disk mass ratio in dust between the two components of $M_{\mathrm{dust, SR\,24S}}/M_{\mathrm{dust, SR\,24N}}\gtrsim840$. The current ALMA observations may imply that either planets have already formed in the SR 24N disk or that dust growth to mm-sizes is inhibited there and that only warm gas, as seen by ro-vibrational CO emission inside the truncation radii of the binary, is present.Comment: Accepted for publication in Ap

A giant planet shaping the disk around the very low-mass star CIDA 1

Context. Exoplanetary research has provided us with exciting discoveries of planets around very low-mass (VLM) stars (0.08 M⊙ ≲ M* ≲ 0.3 M⊙; e.g., TRAPPIST-1 and Proxima Centauri). However, current theoretical models still strive to explain planet formation in these conditions and do not predict the development of giant planets. Recent high-resolution observations from the Atacama Large Millimeter/submillimeter Array (ALMA) of the disk around CIDA 1, a VLM star in Taurus, show substructures that hint at the presence of a massive planet. Aims. We aim to reproduce the dust ring of CIDA 1, observed in the dust continuum emission in ALMA Band 7 (0.9 mm) and Band 4 (2.1 mm), along with its 12CO (J = 3−2) and 13CO (J = 3−2) channel maps, assuming the structures are shaped by the interaction of the disk with a massive planet. We seek to retrieve the mass and position of the putative planet, through a global simulation that assesses planet-disk interactions to quantitatively reproduce protoplanetary disk observations of both dust and gas emission in a self-consistent way. Methods. Using a set of hydrodynamical simulations, we model a protoplanetary disk that hosts an embedded planet with a starting mass of between 0.1 and 4.0 MJup and initially located at a distance of between 9 and 11 au from the central star. We compute the dust and gas emission using radiative transfer simulations, and, finally, we obtain the synthetic observations, treating the images as the actual ALMA observations. Results. Our models indicate that a planet with a minimum mass of ~1.4 MJup orbiting at a distance of ~9−10 au can explain the morphology and location of the observed dust ring in Band 7 and Band 4. We match the flux of the dust emission observation with a dust-to-gas mass ratio in the disk of ~10−2. We are able to reproduce the low spectral index (~2) observed where the dust ring is detected, with a ~40−50% fraction of optically thick emission. Assuming a 12CO abundance of 5 × 10−5 and a 13CO abundance 70 times lower, our synthetic images reproduce the morphology of the 12CO (J = 3−2) and 13CO (J = 3−2) observed channel maps where the cloud absorption allowed a detection. From our simulations, we estimate that a stellar mass M* = 0.2 M⊙ and a systemic velocity vsys = 6.25 km s−1 are needed to reproduce the gas rotation as retrieved from molecular line observations. Applying an empirical relation between planet mass and gap width in the dust, we predict a maximum planet mass of ~4−8 MJup. Conclusions. Our results suggest the presence of a massive planet orbiting CIDA 1, thus challenging our understanding of planet formation around VLM stars

Ring shaped dust accumulation in transition disks

Context.Transition disks are believed to be the final stages of protoplanetary disks, during which a forming planetary system or photoevaporation processes open a gap in the inner disk, drastically changing the disk structure. From theoretical arguments it is expected that dust growth, fragmentation and radial drift are strongly influenced by gas disk structure, and pressure bumps in disks have been suggested as key features that may allow grains to converge and grow efficiently. Aims. We want to study how the presence of a large planet in a disk influences the growth and radial distribution of dust grains, and how observable properties are linked to the mass of the planet. Methods. We combine two-dimensional hydrodynamical disk simulations of disk-planet interactions with state-of-the-art coagulation/fragmentation models to simulate the evolution of dust in a disk which has a gap created by a massive planet. We compute images at different wavelengths and illustrate our results using the example of the transition disk LkCa15. Results. The gap opened by a planet and the long-range interaction between the planet and the outer disk create a single large pressure bump outside the planetary orbit. Millimeter-sized particles form and accumulate at the pressure maximum and naturally produce ring-shaped sub-millimeter emission that is long-lived because radial drift no longer depletes the large grain population of the disk. For large planet masses around 9 $M_{\mathrm{Jup}}$, the pressure maximum and, therefore, the ring of millimeter particles is located at distances that can be more than twice the star-planet separation, creating a large spatial separation between the gas inner edge of the outer disk and the peak millimeter emission. Smaller grains do get closer to the gap and we predict how the surface brightness varies at different wavelengths.Comment: Accepted for publication in Astronomy and Astrophysic