Search

The traditional paradigm of viscosity-dominated evolution of protoplanetary discs has been recently challenged by magnetized disc winds. However, distinguishing wind-driven and turbulence-driven accretion through observations has been difficult. In this study, we present a novel approach to identifying their separate contribution to angular momentum transport by studying the gap and ring morphology of planet-forming discs in the ALMA continuum. We model the gap-opening process of planets in discs with both viscous evolution and wind-driven accretion by 2D multi-fluid hydrodynamical simulations. Our results show that gap-opening planets in wind-driven accreting discs generate characteristic substructures that differ from those in purely viscous discs. Specifically, we demonstrate that discs, where wind-driven accretion dominates the production of substructures, exhibit significant asymmetries. Based on the diverse outputs of mock images in the ALMA continuum, we roughly divide the planet-induced features into four regimes (moderate-viscosity dominated, moderate-wind dominated, strong-wind dominated, inviscid). The classification of these regimes sets up a potential method to constrain the strength of magnetized disc wind and viscosity based on the observed gap and ring morphology. We discuss the asymmetry feature in our mock images and its potential manifestation in ALMA observations., Comment: 11 pages, 6 figures, accepted by MNRAS. Welcome any comments and suggestions!

The azimuthal asymmetries of dust rings in protoplanetary disks such as a crescent around young stars are often interpreted as dust traps, and thus as ideal locations for planetesimal and planet formations. Whether such dust traps effectively promote planetesimal formation in disks around very-low-mass stars (VLM; a mass of $\lesssim$0.2~$M_\odot$) is debatable, as the dynamical and grain growth timescales in such systems are long. To investigate grain growth in such systems, we studied the dust ring with crescent around the VLM star ZZ~Tau~IRS using the Karl G. Jansky Very Large Array (JVLA) at centimeter wavelengths. Significant signals were detected around ZZ~Tau~IRS. To estimate the maximum grain size ($a_{\rm max}$) in the crescent, we compared the observed spectral energy distribution (SED) with SEDs for various $a_{\rm max}$ values predicted by radiative transfer calculations. We found $a_{\rm max} \gtrsim$~1~mm and $\lesssim$~60~$\mu$m in the crescent and ring, respectively, though our modeling efforts rely on uncertain dust properties. Our results suggest that grain growth occurred in the ZZ~Tau~IRS disk, relative to sub-micron-sized interstellar medium. Planet formation in crescent with mm-sized pebbles might proceed more efficiently than in other regions with sub-millimeter-sized pebbles via pebble accretion scenarios., Comment: 11 pages, 3 figures, accepted in ApJ

Planet-disc interactions build up local pressure maxima that may halt the radial drift of protoplanetary dust, and pile it up in rings and crescents. ALMA observations of the HD135344B disc revealed two rings in the thermal continuum stemming from ~mm-sized dust. At higher frequencies the inner ring is brighter relative to the outer ring, which is also shaped as a crescent rather than a full ring. In near-IR scattered light images, the disc is modulated by a 2-armed grand-design spiral originating inside the ALMA inner ring. Such structures may be induced by a massive companion evacuating the central cavity, and by a giant planet in the gap separating both rings, that channels the accretion of small dust and gas through its filamentary wakes while stopping the larger dust from crossing the gap. Here we present ALMA observations in the J=(2-1)CO isotopologue lines and in the adjacent continuum, with up to 12km baselines. Angular resolutions of 0.03" reveal the tentative detection of a filament connecting both rings, and which coincides with a local discontinuity in the pitch angle of the IR spiral, proposed previously as the location of the protoplanet driving this spiral. Line diagnostics suggest that turbulence, or superposed velocity components, is particularly strong in the spirals. The 12CO(2-1) 3-D rotation curve points at stellocentric accretion at radii within the inner dust ring, with a radial velocity of up to ~6%+-0.5% Keplerian, which corresponds to an excessively large accretion rate of ~2E-6M_sun/yr if all of the CO layer follows the 12CO(2-1) kinematics. This suggests that only the surface layers of the disc are undergoing accretion, and that the line broadening is due to superposed laminar flows., accepted to MNRAS

Direct images of protoplanets embedded in disks around infant stars provide the key to understanding the formation of gas giant planets like Jupiter. Using the Subaru Telescope and Hubble Space Telescope, we find evidence for a jovian protoplanet around AB Aurigae orbiting at a wide projected separation (93 au), likely responsible for multiple planet-induced features in the disk. Its emission is reproducible as reprocessed radiation from an embedded protoplanet. We also identify two structures located at 430-580 au that are candidate sites of planet formation. These data reveal planet formation in the embedded phase and a protoplanet discovery at wide, > 50 au separations characteristic of most imaged exoplanets. With at least one clump-like protoplanet and multiple spiral arms, the AB Aur system may also provide the evidence for a long-considered alternative to the canonical model for Jupiter's formation: disk (gravitational) instability., Author's personal version: 19 pages, 5 Figures, 1 Table; 32 Supplementary pages, 18 Supplementary Figures, 1 Supplementary Table; Accepted for Publication in Nature Astronomy. Published version: https://www.nature.com/articles/s41550-022-01634-x

Azimuthally asymmetric structures have been discovered in millimeter continuum emission from many protoplanetary disks. One hypothesis is that they are vortices produced by the Rossby wave instability, for example at edges of planet-opened gaps or deadzones. Confirming the vortex nature of these structures will have profound implications to planet formation. One way to test the hypothesis is to compare the observed morphology of vortex candidates in near-infrared scattered light with theoretical expectations. To this end, we synthesize the appearance of vortices in H-band polarized light by combining hydrodynamic and radiative transfer simulations of the Rossby wave instability at a deadzone edge. In a disk at 140 pc, at the peak in its evolution a vortex at 65 au may appear as a radially narrow arc 50% - 70% brighter compared with an axisymmetric disk model. The contrast depends on the inclination of the disk and the position angle of the vortex only weakly. Such contrast levels are well detectable in imaging observations of bright disks using instruments such as VLT/SPHERE, Subaru/SCExAO, and Gemini/GPI. A vortex also casts a shadow in the outer disk, which may aid its identification. Finally, at modest to high inclinations (e.g., 60 degrees) a vortex may mimic a one-armed spiral. In the HD 34282 disk, such a one-armed spiral with a shadowed region on the outside has been found in scattered light. This feature roughly coincides with an azimuthal asymmetry in mm continuum emission, signifying the presence of a vortex., 10 pages, 7 figures, ApJ accepted

The substructures observed in protoplanetary disks may be the signposts of embedded planets carving gaps or creating vortices. The inferred masses of these planets often fall in the Jovian regime despite their low abundance compared to lower-mass planets, partly because previous works often assume that a single substructure (a gap or vortex) is caused by a single planet. In this work, we study the possible imprints of compact systems composed of Neptune-like planets ($\sim10-30\;M_\oplus$) and show that long-standing vortices are a prevalent outcome when their inter-planetary separation ($\Delta a$) falls below $\sim8$ times $H_{{\rm p}}$ -- the average disk's scale height at the planets locations. In simulations where a single planet is unable to produce long-lived vortices, two-planet systems can preserve them for at least $5,000$ orbits in two regimes: i) fully-shared density gaps with elongated vortices around the stable Lagrange points $L_4$ and $L_5$ for the most compact planet pairs ($\Delta a \lesssim 4.6\; H_{{\rm p}}$); ii) partially-shared gaps for more widely spaced planets ($\Delta a \sim 4.6 - 8\;H_{{\rm p}}$) forming vortices in a density ring between the planets through the Rossby wave instability. The latter case can produce vortices with a wide range of aspect ratios down to $\sim3$ and can occur for planets captured into the 3:2 (2:1) mean-motion resonances for disk's aspects ratios of $h\gtrsim 0.033$ ($h\gtrsim 0.057$). We suggest that their long lifetimes are sustained by the interaction of spiral density waves launched by the neighboring planets. Overall, our results show that distinguishing imprint of compact systems with Neptune-mass planets are long-lived vortices inside the density gaps, which in turn are shallower than single-planet gaps for a fixed gap width., Comment: 17 pages, 7 figures, 1 table, submitted to The Astrophysical Journal

Sparse aperture masking interferometry (SAM) is a high resolution observing technique that allows for imaging at and beyond a telescope's diffraction limit. The technique is ideal for searching for stellar companions at small separations from their host star; however, previous analysis of SAM observations of young stars surrounded by dusty disks have had difficulties disentangling planet and extended disk emission. We analyse VLT/SPHERE-IRDIS SAM observations of the transition disk LkCa\,15, model the extended disk emission, probe for planets at small separations, and improve contrast limits for planets. We fit geometrical models directly to the interferometric observables and recover previously observed extended disk emission. We use dynamic nested sampling to estimate uncertainties on our model parameters and to calculate evidences to perform model comparison. We compare our extended disk emission models against point source models to robustly conclude that the system is dominated by extended emission within 50 au. We report detections of two previously observed asymmetric rings at $\sim$17 au and $\sim$45 au. The peak brightness location of each model ring is consistent with the previous observations. We also, for the first time with imaging, robustly recover an elliptical Gaussian inner disk, previously inferred via SED fitting. This inner disk has a FWHM of ~5 au and a similar inclination and orientation as the outer rings. Finally, we recover no clear evidence for candidate planets. By modelling the extended disk emission, we are able to place a lower limit on the near infrared companion contrast of at least 1000., Comment: Accepted for publication in ApJ

Spiral density waves generated by an embedded planet are understood to cause ``kinks'' in observed velocity channel maps of CO surface emission, by perturbing the gas motion within the spiral arms. If velocity kinks are a reliable probe of embedded planets, we should expect to see the planet-driven spiral arms in other observational tracers. We test this prediction by searching the dust continuum for the midplane counterparts of the spirals responsible for all of the velocity kink planet candidates reported to date, whose orbits lie inside the dust continuum disk. We find no clear detection of any spiral structure in current continuum observations for 6 of the 10 velocity kink planet candidates in our sample (DoAr 25, GW Lup, Sz 129, HD 163296 #2, P94, and HD 143006), despite the high planet masses inferred from the kink amplitude. The remaining 4 cases include 3 clear detections of two-armed dust spirals (Elias 27, IM Lup and WaOph 6) wherein neither spiral arm aligns with a wake originating from reported planet location, suggesting that under the planetary-origin hypothesis, an accurate method for inferring the location of the planet in the midplane may need to encompass vertical effects. The 10th case, HD 97048, is inconclusive with current knowledge of the disk geometry., Comment: 5 figures, 13 pages. Accepted for publication in ApJ Letters. Supplementary material available at https://doi.org/10.6084/m9.figshare.21330426

As one of the most commonly observed disk substructures, dust rings from high-resolution disk surveys appear to have different radial widths. Recent observations on PDS 70 and AB Aur reveal not only planets in the disk, but also the accompanying wide dust rings. We use three-dimensional dust-and-gas disk simulations to study whether gap-opening planets are responsible for the large ring width in disk observations. We find that gap-opening planets can widen rings of dust trapped at the pressure bump via planetary perturbations, even with the mid-plane dust-to-gas ratio approaching order unity and with the dust back-reaction accounted for. We show that the planet-related widening effect of dust rings can be quantified using diffusion-advection theory, and provide a generalized criterion for an equilibrated dust ring width in three-dimensional disk models. We also suggest that the ring width can be estimated using the gas turbulent viscosity $\alpha_{\rm turb}$, but with cautions about the Schmidt number greater than order unity., Comment: 17 pages, 6+3 figures, 2 tables, accepted in ApJ

Large inner dust gaps in transition disks are frequently posited as evidence of giant planets sculpting gas and dust in the disk, or the opening of a gap by photoevaporative winds. Although the former hypothesis is strongly supported by the observations of planets and deep depletions in gas within the gap some disks, many T Tauri stars hosting transition disks accrete at rates typical for an undepleted disk, raising the question of how gap opening occurs in these objects. We thus present an analysis of the structure of the transition disk around the T Tauri star DM Tau, which is strongly accreting ($\sim 10^{-8.3}~\mathrm{M}_\odot~ \mathrm{yr}^{-1}$) and turbulent ($\alpha=0.078 \pm 0.02$). Using the DALI thermochemical code, we fit disk models to simultaneously reproduce the accretion rate, high level of turbulence, the gas traced by ALMA band 6 observations of $^{12}$CO, $^{13}$CO, and C$^{18}$O J=2--1 lines, and the observed dust emission from the mm continuum and spectral energy distribution. We find a shallow depletion in gas surface density of $\sim 10$ relative to the outer disk and a gas rich inner disk is consistent with the observations. The planet mass of $, Comment: Accepted for publication in the Astronomical Journal. 36 pages, 28 figures

Close encounters between young stellar objects in star forming clusters are expected to dramatically perturb circumstellar disks. Such events are witnessed in numerical simulations of star formation, but few direct observations of ongoing encounters have been made. Here we report sub-0".1 resolution Atacama Large Millimeter Array (ALMA) and Jansky Very Large Array (JVLA) observations towards the million year old binary protostar Z CMa in dust continuum and molecular line emission. A point source ~4700 au from the binary has been discovered at both millimeter and centimeter wavelengths. It is located along the extension of a ~2000 au streamer structure previously found in scattered light imaging, whose counterpart in dust and gas emission is also newly identified. Comparison with simulations shows signposts of a rare flyby event in action. Z CMa is a "double burster", as both binary components undergo accretion outbursts, which may be facilitated by perturbations to the host disk by flybys., Comment: Published in Nature Astronomy. Here is the authors' version with the Supplementary Information integrated into the Methods section

Scattered light imaging studies have detected nearly two dozen spiral arm systems in circumstellar disks, yet the formation mechanisms for most of them are still under debate. Although existing studies can use motion measurements to distinguish leading mechanisms such as planet-disk interaction and disk self-gravity, close-in stellar flybys can induce short-lived spirals and even excite arm-driving planets into highly eccentric orbits. With unprecedented stellar location and proper motion measurements from Gaia DR3, here we study for known spiral arm systems their flyby history with their stellar neighbours by formulating an analytical on-sky flyby framework. For stellar neighbors currently located within 10 pc from the spiral hosts, we restrict the flyby time to be within the past $10^4$ yr and the flyby distance to be within $10$ times the disk extent in scattered light. Among a total of $12570$ neighbors that are identified in Gaia DR3 for $20$ spiral systems, we do not identify credible flyby candidates for isolated systems. Our analysis suggests that close-in recent flyby is not the dominant formation mechanism for isolated spiral systems in scattered light., 14 pages, 4 figures, 2 tables, ApJS accepted. Data files for Table 2 in the ancillary folder. Code framework available at https://github.com/slinling/afm-spirals

Formed in protoplanetary disks around young stars, giant planets can leave observational features such as spirals and gaps in their natal disks through planet-disk interactions. Although such features can indicate the existence of giant planets, protoplanetary disk signals can overwhelm the innate luminosity of planets. Therefore, in order to image planets that are embedded in disks, it is necessary to remove the contamination from the disks to reveal the planets possibly hiding within their natal environments. We observe and directly model the detected disk in the Keck/NIRC2 vortex coronagraph $L'$-band observations of the single-armed protoplanetary disk around HD 34282. Despite a non-detection of companions for HD 34282, this direct disk modeling improves planet detection sensitivity by up to a factor of 2 in flux ratio and ${\sim}10 M_{\rm Jupiter}$ in mass. This suggests that performing disk modeling can improve directly imaged planet detection limits in systems with visible scattered light disks, and can help to better constrain the occurrence rates of self-luminous planets in these systems., 8 pages, 3 figures, 1 table, ApJ Letters accepted. Data files available in the ancillary folder

We present a high-contrast imaging search for Pa$\beta$ line emission from protoplanets in the PDS~70 system with Keck/OSIRIS integral field spectroscopy. We applied the high-resolution spectral differential imaging technique to the OSIRIS $J$-band data but did not detect the Pa$\beta$ line at the level predicted using the parameters of \cite{Hashimoto2020}. This lack of Pa$\beta$ emission suggests the MUSE-based study may have overestimated the line width of H$\alpha$. We compared our Pa$\beta$ detection limits with the previous H$\alpha$ flux and H$\beta$ limits and estimated $A_{\rm V}$ to be $\sim0.9$ and 2.0 for PDS~70~b and c respectively. In particular, PDS~70~b's $A_{\rm V}$ is much smaller than implied by high-contrast near-infrared studies, which suggests the infrared-continuum photosphere and the hydrogen-emitting regions exist at different heights above the forming planet., Comment: 8pages, 5 figs, accepted for publication in AJ

FU Orionis-type objects (FUors) are low-mass pre-main sequence stars undergoing a temporary, but significant increase of mass accretion rate from the circumstellar disk onto the protostar. It is not yet clear what triggers the accretion bursts and whether the disks of FUors are in any way different from disks of non-bursting young stellar objects. Motivated by this, we conducted a 1.3 mm continuum survey of ten FUors and FUor-like objects with ALMA, using both the 7 m array and the 12 m array in two different configurations to recover emission at the widest possible range of spatial scales. We detected all targeted sources and several nearby objects as well. To constrain the disk structure, we fit the data with models of increasing complexity from 2D Gaussian to radiative transfer, enabling comparison with other samples modeled in a similar way. The radiative transfer modeling gives disk masses that are significantly larger than what is obtained from the measured millimeter fluxes assuming optically thin emission, suggesting that the FUor disks are optically thick at this wavelength. In comparison with samples of regular Class II and Class I objects, the disks of FUors are typically a factor of 2.9-4.4 more massive and a factor of 1.5-4.7 smaller in size. A significant fraction of them (65-70%) may be gravitationally unstable., Comment: 34 pages, 8 tables, 29 figures, accepted for publication in the Astrophysical Journal Supplement Series

Tidal encounters in star clusters perturb discs around young protostars. In Cuello et al. (2019a, Paper I) we detailed the dynamical signatures of a stellar flyby in both gas and dust. Flybys produce warped discs, spirals with evolving pitch angles, increasing accretion rates, and disc truncation. Here we present the corresponding observational signatures of these features in optical/near-infrared scattered light and (sub-) millimeter continuum and CO line emission. Using representative prograde and retrograde encounters for direct comparison, we post-process hydrodynamical simulations with radiative transfer methods to generate a catalogue of multi-wavelength observations. This provides a reference to identify flybys in recent near-infrared and sub-millimetre observations (e.g., RW Aur, AS 205, HV Tau & DO Tau, FU Ori, V2775 Ori, and Z CMa)., 11 pages, 5 figures, accepted for publication in MNRAS

Dust gaps and rings appear ubiquitous in bright protoplanetary disks. Disk-planet interaction with dust-trapping at the edges of planet-induced gaps is one plausible explanation. However, the sharpness of some observed dust rings indicate that sub-mm-sized dust grains have settled to a thin layer in some systems. We test whether or not such dust around gas gaps opened by planets can remain settled by performing three-dimensional, dust-plus-gas simulations of protoplanetary disks with an embedded planet. We find planets massive enough to open gas gaps stir small, sub-mm-sized dust grains to high disk elevations at the gap edges, where the dust scale-height can reach ~70% of the gas scale-height. We attribute this dust 'puff-up' to the planet-induced meridional gas flows previously identified by Fung & Chiang and others. We thus emphasize the importance of explicit 3D simulations to obtain the vertical distribution of sub-mm-sized grains around gas gaps opened by massive planets. We caution that the gas-gap-opening planet interpretation of well-defined dust rings is only self-consistent with large grains exceeding mm in size., 15 pages, 7+2 figures, accepted in ApJ

The detection of a wide range of substructures such as rings, cavities and spirals has become a common outcome of high spatial resolution imaging of protoplanetary disks, both in the near-infrared scattered light and in the thermal millimetre continuum emission. The most frequent interpretation of their origin is the presence of planetary-mass companions perturbing the gas and dust distribution in the disk (perturbers), but so far the only bona-fide detection has been the two giant planets around PDS 70. Here, we collect a sample of 15 protoplanetary disks showing substructures in SPHERE scattered light images and present a homogeneous derivation of planet detection limits in these systems. We also estimate the mass of these perturbers through a Hill radius prescription and a comparison to ALMA data. Assuming that one single planet carves each substructure in scattered light, we find that more massive perturbers are needed to create gaps within cavities than rings, and that we might be close to a detection in the cavities of RX J1604, RX J1615, Sz Cha, HD 135344B and HD 34282. We reach typical mass limits in these cavities of 3-10 Mjup. For planets in the gaps between rings, we find that the detection limits of SPHERE are about an order of magnitude away in mass, and that the gaps of PDS 66 and HD 97048 seem to be the most promising structures for planet searches. The proposed presence of massive planets causing spiral features in HD 135344B and HD 36112 are also within SPHERE's reach assuming hot-start models.These results suggest that current detection limits are able to detect hot-start planets in cavities, under the assumption that they are formed by a single perturber located at the centre of the cavity. More realistic planet mass constraints would help to clarify whether this is actually the case, which might point to perturbers not being the only way of creating substructures., 26 pages, 33 figures, accepted for publication in A&A

We report an analysis of the dust disk around DM~Tau, newly observed with the Atacama Large Millimeter/submillimeter Array (ALMA) at 1.3 mm. The ALMA observations with high sensitivity (8.4~$\mu$Jy/beam) and high angular resolution (35~mas, 5.1~au) detect two asymmetries on the ring at $r\sim$20~au. They could be two vortices in early evolution, the destruction of a large scale vortex, or double continuum emission peaks with different dust sizes. We also found millimeter emissions with $\sim$50~$\mu$Jy (a lower limit dust mass of 0.3~$M_{\rm Moon}$) inside the 3-au ring. To characterize these emissions, we modeled the spectral energy distribution (SED) of DM~Tau using a Monte Carlo radiative transfer code. We found that an additional ring at $r=$ 1~au could explain both the DM~Tau SED and the central point source. The disk midplane temperature at the 1-au ring calculated in our modeling is less than the typical water sublimation temperature of 150~K, prompting the possibility of forming small icy planets there., Comment: 26 pages, 17 figures

We present Atacama Large Millimeter/submillimeter Array (ALMA) Band 6 observations of dust continuum emission of the disk around WW Cha. The dust continuum image shows a smooth disk structure with a faint (low-contrast) dust ring, extending from $\sim 40$ au to $\sim 70$ au, not accompanied by any gap. We constructed the simple model to fit the visibility of the observed data by using MCMC method and found that the bump (we call the ring without the gap the bump) has two peaks at $40$ au and $70$ au. The residual map between the model and observation indicates asymmetric structures at the center and the outer region of the disk. These asymmetric structures are also confirmed by model-independent analysis of the imaginary part of the visibility. The asymmetric structure at the outer region is consistent with a spiral observed by SPHERE. To constrain physical quantities of the disk (dust density and temperature), we carried out radiative transfer simulations. We found that the midplane temperature around the outer peak is close to the freezeout temperature of CO on water ice ($\sim 30$ K). The temperature around the inner peak is about $50$ K, which is close to the freezeout temperature of H$_2$S and also close to the sintering temperature of several species. We also discuss the size distribution of the dust grains using the spectral index map obtained within the Band 6 data., Comment: 24 pages, 25 figures, accepted for publication in The Astrophysical Journal