Showing 1–50 of 151 results for author: Dong, R

First JVLA Radio Observation on PDS70

Authors: Hauyu Baobab Liu, Simon Casassus, Ruobing Dong, Kiyoaki Doi, Jun Hashimoto, Takayuki Muto

Abstract: PDS~70 is a protoplanetary system that hosts two actively accreting gas giants, namely PDS~70b and PDS~70c. The system has a $\sim$60--100 au dusty ring that has been resolved by ALMA, along with circumplanetary disks around the two gas giants. Here we report the first JVLA Q (40--48 GHz), Ka (29--37 GHz), K (18--26 GHz), and X (8--12 GHz) bands continuum observations, and the complementary ALMA B… ▽ More PDS~70 is a protoplanetary system that hosts two actively accreting gas giants, namely PDS~70b and PDS~70c. The system has a $\sim$60--100 au dusty ring that has been resolved by ALMA, along with circumplanetary disks around the two gas giants. Here we report the first JVLA Q (40--48 GHz), Ka (29--37 GHz), K (18--26 GHz), and X (8--12 GHz) bands continuum observations, and the complementary ALMA Bands 3 ($\sim$98 GHz) and 4 ($\sim$145 GHz) observations towards PDS~70. The dusty ring appears azimuthally asymmetric in our ALMA images. We obtained firm detections at Ka and K bands without spatially resolving the source; we obtained a marginal detection at Q band, and no detection at X band. The spectral indices ($α$) are 5$\pm$1 at 33--44 GHz and 0.6$\pm$0.2 at 22--33 GHz. At 10--22 GHz, the conservative lower limit of $α$ is 1.7. The 33--44 GHz flux density is likely dominated by the optically thin thermal emission of grown dust with $\gtrsim$1 mm maximum grain sizes, which may be associated with the azimuthally asymmetric substructure induced by planet-disk interaction. Since PDS~70 was not detected at X band, we found it hard to explain the low spectral index at 22--33 GHz only with free-free emission. Hence, we attribute the dominant emission at 22--33 GHz to the emission of spinning nanometer-sized dust particles, while free-free emission may partly contribute to emission at this frequency range. In some protoplanetary disks, the emission of spinning nanometer-sized dust particles may resemble the 20--50 GHz excess in the spectra of millimeter-sized dust. The finding of strong continuum emission of spinning nanometer-sized particles can complicate the procedure of constraining the properties of grown dust. Future high-resolution, multi-frequency JVLA/ngVLA and SKA observations may shed light on this issue. △ Less

Submitted 28 June, 2024; originally announced June 2024.

Observational characteristics of circum-planetary-mass-object disks in the era of James Webb Space Telescope

Authors: Xilei Sun, Pinghui Huang, Ruobing Dong, Shang-Fei Liu

Abstract: Recent observations have confirmed circumplanetary disks (CPDs) embedded in parental protoplanetary disks (PPDs). On the other hand, planetary-mass companions (PMCs) and planetary-mass objects (PMOs) are likely to harbor their own accretion disks. Unlike PPDs, CPDs and other disks around planet analogues are generally too compact to be spatially resolved by current instrumentation. In this study,… ▽ More Recent observations have confirmed circumplanetary disks (CPDs) embedded in parental protoplanetary disks (PPDs). On the other hand, planetary-mass companions (PMCs) and planetary-mass objects (PMOs) are likely to harbor their own accretion disks. Unlike PPDs, CPDs and other disks around planet analogues are generally too compact to be spatially resolved by current instrumentation. In this study, we generate over 4,000 spectral energy distributions (SEDs) of circum-PMO-disks (CPMODs) with various host temperature and disk properties, which can be categorized into four prototypes, i.e., full, pre-transitional, transitional and evolved CPMODs. We propose a classification scheme based on their near-to-mid-infrared colors. Using those CPMOD models, we synthesize JWST (NIRCam and MIRI) photometry for F444W, F1000W and F2550W wide filters. We show F444W - F1000W and F444 - F2550W colors can be applied to distinguish different types of CPMODs, especially for those around hot hosts. Our results indicate that the ongoing and future JWST observations are promising to unveil structures and properties of CPMODs. △ Less

Submitted 13 June, 2024; originally announced June 2024.

Comments: 18 pages, 7 figures, accepted for publication in ApJ

The James Webb Interferometer: Space-based interferometric detections of PDS 70 b and c at 4.8 $μ$m

Authors: Dori Blakely, Doug Johnstone, Gabriele Cugno, Anand Sivaramakrishnan, Peter Tuthill, Ruobing Dong, Benjamin J. S. Pope, Loïc Albert, Max Charles, Rachel A. Cooper, Matthew De Furio, Louis Desdoigts, René Doyon, Logan Francis, Alexandra Z. Greenbaum, David Lafrenière, James P. Lloyd, Michael R. Meyer, Laurent Pueyo, Shrishmoy Ray, Joel Sánchez-Bermúdez, Anthony Soulain, Deepashri Thatte, Thomas Vandal

Abstract: We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of… ▽ More We observed the planet-hosting system PDS 70 with the James Webb Interferometer, JWST's Aperture Masking Interferometric (AMI) mode within NIRISS. Observing with the F480M filter centered at 4.8 $μ$m, we simultaneously fit a geometric model to the outer disk and the two known planetary companions. We re-detect the protoplanets PDS 70 b and c at an SNR of 21 and 11, respectively. Our photometry of both PDS 70 b and c provide evidence for circumplanetary disk emission through fitting SED models to these new measurements and those found in the literature. We also newly detect emission within the disk gap at an SNR of $\sim$4, at a position angle of $207^{+11}_{-10}$ degrees, and an unconstrained separation within $\sim$200 mas. Follow-up observations will be needed to determine the nature of this emission. We place a 5$σ$ upper limit of $Δ$mag = 7.56 on the contrast of the candidate PDS 70 d at 4.8 $μ$m, which indicates that if the previously observed emission at shorter wavelengths is due to a planet, this putative planet has a different atmospheric composition than PDS 70 b or c. Finally, we place upper limits on emission from any additional planets in the disk gap. We find an azimuthally averaged 5$σ$ upper limit of $Δ$mag $\approx$ 7.5 at separations greater than 125 mas. These are the deepest limits to date within $\sim$250 mas at 4.8 $μ$m and the first space-based interferometric observations of this system. △ Less

Submitted 19 April, 2024; originally announced April 2024.

Comments: Submitted to ApJ

JWST/NIRCam Imaging of Young Stellar Objects III: Detailed Imaging of the Nebular Environment Around the HL Tau Disk

Authors: Camryn Mullin, Ruobing Dong, Jarron Leisenring, Gabriele Cugno, Thomas Greene, Doug Johnstone, Michael R. Meyer, Kevin R. Wagner, Schuyler G. Wolff, Martha Boyer, Scott Horner, Klaus Hodapp, Don McCarthy, George Rieke, Marcia Rieke, Erick Young

Abstract: As part of the James Webb Space Telescope (JWST) Guaranteed Time Observation (GTO) program "Direct Imaging of YSOs" (program ID 1179), we use JWST NIRCam's direct imaging mode in F187N, F200W, F405N, and F410M to perform high contrast observations of the circumstellar structures surrounding the protostar HL Tau. The data reveal the known stellar envelope, outflow cavity, and streamers, but do not… ▽ More As part of the James Webb Space Telescope (JWST) Guaranteed Time Observation (GTO) program "Direct Imaging of YSOs" (program ID 1179), we use JWST NIRCam's direct imaging mode in F187N, F200W, F405N, and F410M to perform high contrast observations of the circumstellar structures surrounding the protostar HL Tau. The data reveal the known stellar envelope, outflow cavity, and streamers, but do not detect any companion candidates. We detect scattered light from an in-flowing spiral streamer previously detected in $\textrm{HCO}^+$ by ALMA, and part of the structure connected to the c-shaped outflow cavity. For detection limits in planet mass we use BEX evolutionary tracks when $M_\textrm{p}<2M_\textrm{J}$ and AMES-COND evolutionary tracks otherwise, assuming a planet age of 1 Myr (youngest available age). Inside the disk region, due to extended envelope emission, our point-source sensitivities are $\sim5$ mJy ($37~M_{\rm J}$) at 40 AU in F187N, and $\sim0.37$ mJy ($5.2~M_{\rm J}$) at 140 AU in F405N. Outside the disk region, the deepest limits we can reach are $\sim0.01$ mJy ($0.75~M_{\rm J}$) at a projected separation of $\sim525$ AU. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: 13 pages, 6 figures, 2 tables, accepted to AAS Astronomical Journal

Forming localized dust concentrations in a dust ring: DM Tau case study

Authors: Hauyu Baobab Liu, Takayuki Muto, Mihoko Konishi, Chia-Ying Chung, Jun Hashimoto, Kiyoaki Doi, Ruobing Dong, Tomoyuki Kudo, Yasuhiro Hasegawa, Yuka Terada, Akimasa Kataoka

Abstract: The previous, high angular resolution 225 GHz ($\sim$1.3 mm) continuum observations on the transitional disk DM Tau have resolved an outer ring at 20-120 au radii that is weakly azimuthally asymmetric. We aimed to examine dust growth and filtration in the outer ring. We performed the $\sim$0$''$.06 ($\sim$8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40-48 GHz ($\sim$7 mm; Q band) cont… ▽ More The previous, high angular resolution 225 GHz ($\sim$1.3 mm) continuum observations on the transitional disk DM Tau have resolved an outer ring at 20-120 au radii that is weakly azimuthally asymmetric. We aimed to examine dust growth and filtration in the outer ring. We performed the $\sim$0$''$.06 ($\sim$8.7 au) resolution Karl G. Jansky Very Large Array (JVLA) 40-48 GHz ($\sim$7 mm; Q band) continuum observations and the complementary observations at lower frequencies. In addition, we analyzed the archival JVLA observations that were taken since 2010. Intriguingly, the Q band image resolved the azimuthally highly asymmetric, knotty dust emission sources close to the inner edge of the outer ring. Fitting the 8-700 GHz spectral energy distribution (SED) with two dust components indicates that the maximum grain size in these knotty dust emission sources is likely $\gtrsim$300 $μ$m while it is $\lesssim$50 $μ$m in the rest of the ring. These results may be explained by trapping of inward migrating grown dust close to the ring inner edge. The exact mechanism for developing the azimuthal asymmetry has not yet been identified, which may be due to planet-disk interaction that might also be responsible for the creation of the dust cavity and pressure bump, or the fluid instabilities and vortex formation due to shear motions. Finally, we remark that the asymmetries in DM Tau are hard to diagnose from the $\gtrsim$225 GHz observations owing to a high optical depth at the ring. In other words, the apparent symmetric or asymmetric morphology of the transitional disks may be related to the optical depths of those disks at the observing frequency. △ Less

Submitted 13 February, 2024; v1 submitted 5 February, 2024; originally announced February 2024.

Comments: 16 pages, 7 figures; accepted to A&A

JWST/NIRCam Imaging of Young Stellar Objects. II. Deep Constraints on Giant Planets and a Planet Candidate Outside of the Spiral Disk Around SAO 206462

Authors: Gabriele Cugno, Jarron Leisenring, Kevin R. Wagner, Camryn Mullin, Roubing Dong, Thomas Greene, Doug Johnstone, Michael R. Meyer, Schuyler G. Wolff, Charles Beichman, Martha Boyer, Scott Horner, Klaus Hodapp, Doug Kelly, Don McCarthy, Thomas Roellig, George Rieke, Marcia Rieke, John Stansberry, Erick Young

Abstract: We present JWST/NIRCam F187N, F200W, F405N and F410M direct imaging data of the disk surrounding SAO 206462. Previous images show a very structured disk, with a pair of spiral arms thought to be launched by one or more external perturbers. The spiral features are visible in three of the four filters, with the non-detection in F410M due to the large detector saturation radius. We detect with a sign… ▽ More We present JWST/NIRCam F187N, F200W, F405N and F410M direct imaging data of the disk surrounding SAO 206462. Previous images show a very structured disk, with a pair of spiral arms thought to be launched by one or more external perturbers. The spiral features are visible in three of the four filters, with the non-detection in F410M due to the large detector saturation radius. We detect with a signal-to-noise ratio of 4.4 a companion candidate (CC1) that, if on a coplanar circular orbit, would orbit SAO 206462 at a separation of $\sim300$ au, $2.25σ$ away from the predicted separation for the driver of the eastern spiral. According to the BEX models, CC1 has a mass of $M_\mathrm{CC1}=0.8\pm0.3~M_\mathrm{J}$. No other companion candidates were detected. At the location predicted by simulations of both spirals generated by a single massive companion, the NIRCam data exclude objects more massive than $\sim2.2~M_\mathrm{J}$ assuming the BEX evolutionary models. In terms of temperatures, the data are sensitive to objects with $T_{\text{eff}}\sim650-850$ K, when assuming planets emit like blackbodies ($R_\mathrm{p}$ between 1 and $3 R_\mathrm{J}$). From these results, we conclude that if the spirals are driven by gas giants, these must be either cold or embedded in circumplanetary material. In addition, the NIRCam data provide tight constraints on ongoing accretion processes. In the low extinction scenario we are sensitive to mass accretion rates of the order $\dot{M}\sim10^{-9} M_\mathrm{J}$ yr$^{-1}$. Thanks to the longer wavelengths used to search for emission lines, we reach unprecedented sensitivities to processes with $\dot{M}\sim10^{-7} M_\mathrm{J}$ yr$^{-1}$ even towards highly extincted environments ($A_\mathrm{V}\approx50$~mag). △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: 18 pages, 8 figures, 3 tables

JWST/NIRCam Imaging of Young Stellar Objects. I. Constraints on Planets Exterior to The Spiral Disk Around MWC 758

Authors: Kevin Wagner, Jarron Leisenring, Gabriele Cugno, Camryn Mullin, Ruobing Dong, Schuyler G. Wolff, Thomas Greene, Doug Johnstone, Michael R. Meyer, Charles Beichman, Martha Boyer, Scott Horner, Klaus Hodapp, Doug Kelly, Don McCarthy, Tom Roellig, George Rieke, Marcia Rieke, Michael Sitko, John Stansberry, Erick Young

Abstract: MWC 758 is a young star hosting a spiral protoplanetary disk. The spirals are likely companion-driven, and two previously-identified candidate companions have been identified -- one at the end the Southern spiral arm at ~0.6 arcsec, and one interior to the gap at ~0.1 arcsec. With JWST/NIRCam, we provide new images of the disk and constraints on planets exterior to ~1". We detect the two-armed spi… ▽ More MWC 758 is a young star hosting a spiral protoplanetary disk. The spirals are likely companion-driven, and two previously-identified candidate companions have been identified -- one at the end the Southern spiral arm at ~0.6 arcsec, and one interior to the gap at ~0.1 arcsec. With JWST/NIRCam, we provide new images of the disk and constraints on planets exterior to ~1". We detect the two-armed spiral disk, a known background star, and a spatially resolved background galaxy, but no clear companions. The candidates that have been reported are at separations that are not probed by our data with sensitivity sufficient to detect them -- nevertheless, these observations place new limits on companions down to ~2 Jupiter-masses at ~150 au and ~0.5 Jupiter masses at ~600 au. Owing to the unprecedented sensitivity of JWST and youth of the target, these are among the deepest mass-detection limits yet obtained through direct imaging observations, and provide new insights into the system's dynamical nature. △ Less

Submitted 5 January, 2024; originally announced January 2024.

Comments: Accepted for publication in AJ

Shadowing in the protoplanetary disk of ZZ Tau IRS with HST

Authors: Jun Hashimoto, Ruobing Dong, Takayuki Muto, Hauyu Baobab Liu, Yuka Terada

Abstract: An inner component misaligned from an outer component in a protoplanetary disk can result in the former casting shadows on the latter. We present a new instance of shadowing on the outer disk around a very low mass star, ZZ~Tau~IRS. Through the analysis of near-infrared (NIR) archival data at $λ=1.6$~$μ$m acquired with the Wide Field Camera 3 on the Hubble Space Telescope, we identified brightness… ▽ More An inner component misaligned from an outer component in a protoplanetary disk can result in the former casting shadows on the latter. We present a new instance of shadowing on the outer disk around a very low mass star, ZZ~Tau~IRS. Through the analysis of near-infrared (NIR) archival data at $λ=1.6$~$μ$m acquired with the Wide Field Camera 3 on the Hubble Space Telescope, we identified brightness asymmetries in the top and bottom halves of the highly inclined outer disk, separated by a dark lane. The brighter sides in the top and bottom halves are on the opposite sides, which we attributed to shadows cast by a misaligned inner disk. Radiative transfer modeling of the system with a misaligned angle of 15~deg between the inner and outer disks well reproduced the observations. Additionally, we found an elevated brightness temperature of $^{12}$CO~(3-2) at $r\sim30$~au on the brighter side in NIR wavelengths in the top half disk, which can be explained by the shadowing effect too. While the origin of the misaligned inner disk remains unclear, future monitoring observations to search for temporal variations in brightness asymmetries will likely provide useful clues. △ Less

Submitted 3 January, 2024; originally announced January 2024.

Comments: 12 pages, 8 figures, accepted in AJ

A Uniform Analysis of Debris Disks with the Gemini Planet Imager I: An Empirical Search for Perturbations from Planetary Companions in Polarized Light Images

Authors: Katie A. Crotts, Brenda C. Matthews, Gaspard Duchêne, Thomas M. Esposito, Ruobing Dong, Justin Hom, Rebecca Oppenheimer, Malena Rice, Schuyler G. Wolff, Christine H. Chen, Clarissa R. Do Ó, Paul Kalas, Briley L. Lewis, Alycia J. Weinberger, David J. Wilner, Mark Ammons, Pauline Arriaga, Robert J. De Rosa, John H. Debes, Michael P. Fitzgerald, Eileen C. Gonzales, Dean C. Hines, Sasha Hinkley, A. Meredith Hughes, Ludmilla Kolokolova , et al. (15 additional authors not shown)

Abstract: The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, p… ▽ More The Gemini Planet Imager (GPI) has excelled in imaging debris disks in the near-infrared. The GPI Exoplanet Survey (GPIES) imaged twenty-four debris disks in polarized $H$-band light, while other programs observed half of these disks in polarized $J$- and/or $K1$-bands. Using these data, we present a uniform analysis of the morphology of each disk to find asymmetries suggestive of perturbations, particularly those due to planet-disk interactions. The multi-wavelength surface brightness, the disk color and geometry permit identification of any asymmetries such as warps or disk offsets from the central star. We find that nineteen of the disks in this sample exhibit asymmetries in surface brightness, disk color, disk geometry, or a combination of the three, suggesting that for this sample, perturbations, as seen in scattered light, are common. The relationship between these perturbations and potential planets in the system are discussed. We also explore correlations among stellar temperatures, ages, disk properties, and observed perturbations. We find significant trends between the vertical aspect ratio and the stellar temperature, disk radial extent, and the dust grain size distribution power-law, $q$. We also confirm a trend between the disk color and stellar effective temperature, where the disk becomes increasingly red/neutral with increasing temperature. Such results have important implications on the evolution of debris disk systems around stars of various spectral types. △ Less

Submitted 24 November, 2023; originally announced November 2023.

Comments: 46 pages, 20 figures, 6 tables, accepted for publication in ApJ

A Companion in V1247 Ori Supported by Spiral Arm Pattern Motion

Authors: Bin B. Ren, Chen Xie, Myriam Benisty, Ruobing Dong, Jaehan Bae, Tomas Stolker, Rob G. van Holstein, John H. Debes, Antonio Garufi, Christian Ginski, Stefan Kraus

Abstract: While there have been nearly two dozen of spiral arms detected from planet-forming disks in near-infrared scattered light, none of their substellar drivers have been confirmed. By observing spiral systems in at least two epochs spanning multiple years, and measuring the motion of the spirals, we can distinguish the cause of the spirals, and locate the orbits of the driving planets if they trigger… ▽ More While there have been nearly two dozen of spiral arms detected from planet-forming disks in near-infrared scattered light, none of their substellar drivers have been confirmed. By observing spiral systems in at least two epochs spanning multiple years, and measuring the motion of the spirals, we can distinguish the cause of the spirals, and locate the orbits of the driving planets if they trigger the spirals. Upon a recent validation of this approach using the co-motion between a stellar companion and a spiral, we obtained a second epoch observation for the spiral system in the disk of V1247 Ori in the $H$-band polarized scattered light using VLT/SPHERE/IRDIS. Combining our observations with archival IRDIS data, we established a $4.8$ yr timeline to constrain the V1247 Ori spiral motion. We obtained a pattern speed of $0.40^{\circ} \pm 0.09^{\circ}$ yr$^{-1}$ for the north-east spiral. This corresponds to an orbital period of $900\pm200$ yr, and thus the semi-major axis of the hidden planetary driver is $118\pm19$ au for a 2.0 $\pm$ 0.1 M$_\odot$ central star. The location agrees with the gap in ALMA dust continuum observations, providing joint support for the existence of a companion driving the scattered-light spirals while carving a millimeter gap. With an angular separation of 0.29" $\pm$ 0.05", this hidden companion is an ideal target for JWST imaging. △ Less

Submitted 7 December, 2023; v1 submitted 23 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Astronomy and Astrophysics; 6 pages, 5 figures

The dependence of the structure of planet-opened gaps in protoplanetary disks on radiative cooling

Authors: Minghao Zhang, Pinghui Huang, Ruobing Dong

Abstract: Planets can excite density waves and open annular gas gaps in protoplanetary disks. The depth of gaps is influenced by the evolving angular momentum carried by density waves. While the impact of radiative cooling on the evolution of density waves has been studied, a quantitative correlation to connect gap depth with the cooling timescale is lacking. To address this gap in knowledge, we employ the… ▽ More Planets can excite density waves and open annular gas gaps in protoplanetary disks. The depth of gaps is influenced by the evolving angular momentum carried by density waves. While the impact of radiative cooling on the evolution of density waves has been studied, a quantitative correlation to connect gap depth with the cooling timescale is lacking. To address this gap in knowledge, we employ the grid-based code Athena++ to simulate disk-planet interactions, treating cooling as a thermal relaxation process. We establish quantitative dependences of steady-state gap depth (Eq. 36) and width (Eq. 41) on planetary mass, Shakura-Sunyaev viscosity, disk scale height, and thermal relaxation timescale $(β)$. We confirm previous results that gap opening is the weakest when thermal relaxation timescale is comparable to local dynamical timescale. Significant variations in gap depth, up to an order of magnitude, are found with different $β$. In terms of width, a gap is at its narrowest around $β=1$, approximately $10\%$ to $20\%$ narrower compared to the isothermal case. When $β\sim100$, it can be $\sim20\%$ wider, and higher viscosity enhances this effect. We derive possible masses of the gas gap-opening planets in AS 209, HD 163296, MWC 480, and HL Tau, accounting for the uncertainties in local thermal relaxation timescale. △ Less

Submitted 18 October, 2023; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: 19 pages, 16 figures, 4 tables, accepted for publication in ApJ

Formation of misaligned second-generation discs through flyby encounters

Authors: Jeremy L. Smallwood, Rebecca Nealon, Nicolás Cuello, Ruobing Dong, Richard A. Booth

Abstract: Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary di… ▽ More Observations reveal protoplanetary discs being perturbed by flyby candidates. We simulate a scenario where an unbound perturber, i.e., a flyby, undergoes an inclined grazing encounter, capturing material and forming a second-generation protoplanetary disc. We run $N$--body and three-dimensional hydrodynamical simulations of a parabolic flyby grazing a particle disc and a gas-rich protoplanetary disc, respectively. In both our $N$--body and hydrodynamic simulations, we find that the captured, second-generation disc forms at a tilt twice the initial flyby tilt. This relationship is robust to variations in the flyby's tilt, position angle, periastron, and mass. We extend this concept by also simulating the case where the flyby has a disc of material prior to the encounter but we do not find the same trend. An inclined disc with respect to the primary disc around a misaligned flyby is tilted by a few degrees, remaining close to its initial disc tilt. Therefore, if a disc is present around the flyby before the encounter, the disc may not tilt up to twice the perturber tilt depending on the balance between the angular momentum of the circumsecondary disc and captured particles. In the case where the perturber has no initial disc, analyzing the orientation of these second-generation discs can give information about the orbital properties of the flyby encounter. △ Less

Submitted 5 October, 2023; originally announced October 2023.

Comments: 17 pages, 15 figures, accepted to MNRAS

Evolution of the Planetary Obliquity: The Eccentric Kozai-Lidov Mechanism Coupled with Tide

Authors: Xiumin Huang, Jianghui Ji, Shangfei Liu, Ruobing Dong, Su Wang

Abstract: The planetary obliquity plays a significant role in determining physical properties of planetary surfaces and climate. As direct detection is constrained due to the present observation accuracy, kinetic theories are helpful to predict the evolution of the planetary obliquity. Here the coupling effect between the eccentric Kozai-Lidov (EKL) effect and the equilibrium tide is extensively investigate… ▽ More The planetary obliquity plays a significant role in determining physical properties of planetary surfaces and climate. As direct detection is constrained due to the present observation accuracy, kinetic theories are helpful to predict the evolution of the planetary obliquity. Here the coupling effect between the eccentric Kozai-Lidov (EKL) effect and the equilibrium tide is extensively investigated, the planetary obliquity performs to follow two kinds of secular evolution paths, based on the conservation of total angular momentum. The equilibrium timescale of the planetary obliquity $t_{\mathrm{eq}}$ varies along with $r_{t}$, which is defined as the initial timescale ratio of the tidal dissipation and secular perturbation. We numerically derive the linear relationship between $t_{\mathrm{eq}}$ and $r_{t}$ with the maximum likelihood method. The spin-axis orientation of S-type terrestrials orbiting M-dwarfs reverses over $90^\circ$ when $r_{t} > 100$, then enter the quasi-equilibrium state between $40^\circ$ and $60^\circ$, while the maximum obliquity can reach $130^\circ$ when $r_{t} > 10^4 $. Numerical simulations show that the maximum obliquity increases with the semi-major axis ratio $a_1$/$a_2$, but is not so sensitive to the eccentricity $e_2$. The likelihood of obliquity flip for S-type terrestrials in general systems with $a_2 < 45$ AU is closely related to $m_1$. The observed potential oblique S-type planets HD 42936 b, GJ 86 Ab and $τ$ Boot Ab are explored to have a great possibility to be head-down over the secular evolution of spin. △ Less

Submitted 26 August, 2023; originally announced August 2023.

Comments: 18 pages, 12 figures, accepted for publication in ApJ

Journal ref: The Astrophysical Journal, Volume 956, Number 1, article id.45 (2023)

Centimeter-sized Grains in the Compact Dust Ring around Very Low Mass Star CIDA 1

Authors: Jun Hashimoto, Hauyu Baobab Liu, Ruobing Dong, Beibei Liu, Takayuki Muto, Yuka Terada

Abstract: We examined the grain size in the dust ring encircling the 0.19~$M_\sun$ T Tauri star CIDA 1 using the Karl G. Jansky Very Large Array (JVLA) at multiple centimeter wavelengths, with a spatial resolution of 0$\farcs$2--0$\farcs$9. We detected distinct partial-ring structures at these wavelengths around CIDA~1. Based on spatial distributions and spectral indexes, we determined that these centimeter… ▽ More We examined the grain size in the dust ring encircling the 0.19~$M_\sun$ T Tauri star CIDA 1 using the Karl G. Jansky Very Large Array (JVLA) at multiple centimeter wavelengths, with a spatial resolution of 0$\farcs$2--0$\farcs$9. We detected distinct partial-ring structures at these wavelengths around CIDA~1. Based on spatial distributions and spectral indexes, we determined that these centimeter emissions originated from dust, rather than free-free or synchrotron emissions. To estimate the maximum grain size ($a_{\rm max}$) within the ring, we compared the observed spectral energy distribution (SED) with SEDs calculated for different $a_{\rm max}$ values using radiative transfer calculations. Our findings indicate an $a_{\rm max}$ value of approximately 2.5~cm in the ring, assuming the dust opacity can be approximated by the DSHARP models. These results suggest that grain growth took place within the CIDA~1 ring, potentially facilitating more efficient planet formation through pebble accretion scenarios involving centimeter-sized pebbles. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 14 pages, 5 figures, accepted in AJ

Chemical footprints of giant planet formation. Role of planet accretion in shaping the C/O ratio of protoplanetary disks

Authors: Haochang Jiang, Yu Wang, Chris W. Ormel, Sebastiaan Krijt, Ruobing Dong

Abstract: Protoplanetary disks, the birthplaces of planets, commonly feature bright rings and dark gaps in both continuum and line emission maps. Accreting planets are interacting with the disk, not only through gravity, but also by changing the local irradiation and elemental abundances, which are essential ingredients for disk chemistry. We propose that giant planet accretion can leave chemical footprints… ▽ More Protoplanetary disks, the birthplaces of planets, commonly feature bright rings and dark gaps in both continuum and line emission maps. Accreting planets are interacting with the disk, not only through gravity, but also by changing the local irradiation and elemental abundances, which are essential ingredients for disk chemistry. We propose that giant planet accretion can leave chemical footprints in the gas local to the planet, which potentially leads to the spatial coincidence of molecular emissions with the planet in ALMA observation. Through 2D multi-fluid hydrodynamical simulations in Athena++ with built-in sublimation, we simulate the process of an accreting planet locally heating up its vicinity, opening a gas gap in the disk, and creating the conditions for C-photochemistry. An accreting planet located outside the methane snowline can render the surrounding gas hot enough to sublimate the C-rich organics off pebbles before they are accreted by the planet. This locally elevates the disk gas-phase C/O ratio, providing a potential explanation for the C$_2$H line-emission rings observed with ALMA. In particular, our findings provide an explanation for the MWC480 disk, where previous work has identified a statistically significant spatial coincidence of line-emission rings inside a continuum gap. Our findings present a novel view of linking the gas accretion of giant planets and their natal disks through the chemistry signals. This model demonstrates that giant planets can actively shape their forming chemical environment, moving beyond the traditional understanding of the direct mapping of primordial disk chemistry onto planets. △ Less

Submitted 17 July, 2023; originally announced July 2023.

Comments: 18 pages, 13 figures, Accepted for publication in A&A

Journal ref: A&A 678, A33 (2023)

Direct images and spectroscopy of a giant protoplanet driving spiral arms in MWC 758

Authors: Kevin Wagner, Jordan Stone, Andrew Skemer, Steve Ertel, Ruobing Dong, Dániel Apai, Eckhart Spalding, Jarron Leisenring, Michael Sitko, Kaitlin Kratter, Travis Barman, Mark Marley, Brittany Miles, Anthony Boccaletti, Korash Assani, Ammar Bayyari, Taichi Uyama, Charles E. Woodward, Phil Hinz, Zackery Briesemeister, Kellen Lawson, François Ménard, Eric Pantin, Ray W. Russell, Michael Skrutskie , et al. (1 additional authors not shown)

Abstract: Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets t… ▽ More Understanding the driving forces behind spiral arms in protoplanetary disks remains a challenge due to the faintness of young giant planets. MWC 758 hosts such a protoplanetary disk with a two-armed spiral pattern that is suggested to be driven by an external giant planet. We present new thermal infrared observations that are uniquely sensitive to redder (i.e., colder or more attenuated) planets than past observations at shorter wavelengths. We detect a giant protoplanet, MWC 758c, at a projected separation of ~100 au from the star. The spectrum of MWC 758c is distinct from the rest of the disk and consistent with emission from a planetary atmosphere with Teff = 500 +/- 100 K for a low level of extinction (AV<30), or a hotter object with a higher level of extinction. Both scenarios are commensurate with the predicted properties of the companion responsible for driving the spiral arms. MWC 758c provides evidence that spiral arms in protoplanetary disks can be caused by cold giant planets or by those whose optical emission is highly attenuated. MWC 758c stands out both as one of the youngest giant planets known, and also as one of the coldest and/or most attenuated. Furthermore, MWC 758c is among the first planets to be observed within a system hosting a protoplanetary disk. △ Less

Submitted 8 July, 2023; originally announced July 2023.

Comments: Published in Nature Astronomy

Dynamical detection of a companion driving a spiral arm in a protoplanetary disk

Authors: Chen Xie, Bin B. Ren, Ruobing Dong, Élodie Choquet, Arthur Vigan, Jean-François Gonzalez, Kevin Wagner, Taotao Fang, Maria Giulia Ubeira-Gabellini

Abstract: Radio and near-infrared observations have observed dozens of protoplanetary disks that host spiral arm features. Numerical simulations have shown that companions may excite spiral density waves in protoplanetary disks via companion-disk interaction. However, the lack of direct observational evidence for spiral-driving companions poses challenges to current theories of companion-disk interaction. H… ▽ More Radio and near-infrared observations have observed dozens of protoplanetary disks that host spiral arm features. Numerical simulations have shown that companions may excite spiral density waves in protoplanetary disks via companion-disk interaction. However, the lack of direct observational evidence for spiral-driving companions poses challenges to current theories of companion-disk interaction. Here we report multi-epoch observations of the binary system HD 100453 with the Spectro-Polarimetric High-contrast Exoplanet REsearch (SPHERE) facility at the Very Large Telescope. By recovering the spiral features via robustly removing starlight contamination, we measure spiral motion across 4 yr to perform dynamical motion analyses. The spiral pattern motion is consistent with the orbital motion of the eccentric companion. With this first observational evidence of a companion driving a spiral arm among protoplanetary disks, we directly and dynamically confirm the long-standing theory on the origin of spiral features in protoplanetary disks. With the pattern motion of companion-driven spirals being independent of companion mass, here we establish a feasible way of searching for hidden spiral-arm-driving planets that are beyond the detection of existing ground-based high-contrast imagers. △ Less

Submitted 17 June, 2023; v1 submitted 15 June, 2023; originally announced June 2023.

Comments: Accepted for publication in Astronomy and Astrophysics; 12 pages, 9 figures

Journal ref: A&A 675, L1 (2023)

PPDONet: Deep Operator Networks for Fast Prediction of Steady-State Solutions in Disk-Planet Systems

Authors: Shunyuan Mao, Ruobing Dong, Lu Lu, Kwang Moo Yi, Sifan Wang, Paris Perdikaris

Abstract: We develop a tool, which we name Protoplanetary Disk Operator Network (PPDONet), that can predict the solution of disk-planet interactions in protoplanetary disks in real-time. We base our tool on Deep Operator Networks (DeepONets), a class of neural networks capable of learning non-linear operators to represent deterministic and stochastic differential equations. With PPDONet we map three scalar… ▽ More We develop a tool, which we name Protoplanetary Disk Operator Network (PPDONet), that can predict the solution of disk-planet interactions in protoplanetary disks in real-time. We base our tool on Deep Operator Networks (DeepONets), a class of neural networks capable of learning non-linear operators to represent deterministic and stochastic differential equations. With PPDONet we map three scalar parameters in a disk-planet system -- the Shakura \& Sunyaev viscosity $α$, the disk aspect ratio $h_\mathrm{0}$, and the planet-star mass ratio $q$ -- to steady-state solutions of the disk surface density, radial velocity, and azimuthal velocity. We demonstrate the accuracy of the PPDONet solutions using a comprehensive set of tests. Our tool is able to predict the outcome of disk-planet interaction for one system in less than a second on a laptop. A public implementation of PPDONet is available at \url{https://github.com/smao-astro/PPDONet}. △ Less

Submitted 18 May, 2023; originally announced May 2023.

Comments: 10 pages, 6 figures, 2 tables; ApJL accepted

Distinguishing Magnetized Disc Winds from Turbulent Viscosity through Substructure Morphology in Planet-forming Discs

Authors: Yinhao Wu, Yi-Xian Chen, Haochang Jiang, Ruobing Dong, Enrique Macías, Min-Kai Lin, Giovanni P. Rosotti, Vardan Elbakyan

Abstract: 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 morpholog… ▽ More 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. △ Less

Submitted 27 May, 2023; v1 submitted 14 May, 2023; originally announced May 2023.

Comments: 11 pages, 6 figures, accepted by MNRAS. Welcome any comments and suggestions!

Exciting spiral arms in protoplanetary discs from flybys

Authors: Jeremy L. Smallwood, Chao-Chin Yang, Zhaohuan Zhu, Rebecca G. Martin, Ruobing Dong, Nicolás Cuello, Andrea Isella

Abstract: Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges f… ▽ More Spiral arms are observed in numerous protoplanetary discs. These spiral arms can be excited by companions, either on bound or unbound orbits. We simulate a scenario where an unbound perturber, i.e. a flyby, excites spiral arms during a periastron passage. We run three-dimensional hydrodynamical simulations of a parabolic flyby encountering a gaseous protoplanetary disc. The perturber mass ranges from $10\, \rm M_J$ to $1\, \rm M_{\odot}$. The perturber excites a two-armed spiral structure, with a more prominent spiral feature for higher mass perturbers. The two arms evolve over time, eventually winding up, consistent with previous works. We focus on analysing the pattern speed and pitch angle of these spirals during the whole process. The initial pattern speed of the two arms are close to the angular velocity of the perturber at periastron, and then it decreases over time. The pitch angle also decreases over time as the spiral winds up. The spirals disappear after several local orbital times. An inclined prograde orbit flyby induces similar disc substructures as a coplanar flyby. A solar-mass flyby event causes increased eccentricity growth in the protoplanetary disc, leading to an eccentric disc structure which dampens over time. The spirals' morphology and the disc eccentricity can be used to search for potential unbound stars or planets around discs where a flyby is suspected. Future disc observations at high resolution and dedicated surveys will help to constrain the frequency of such stellar encounters in nearby star-forming regions. △ Less

Submitted 10 March, 2023; originally announced March 2023.

Comments: 17 pages, 18 figures, accepted to MNRAS

Multiple Rings and Asymmetric Structures in the Disk of SR 21

Authors: Yi Yang, Hauyu Baobab Liu, Takayuki Muto, Jun Hashimoto, Ruobing Dong, Kazuhiro Kanagawa, Munetake Momose, Eiji Akiyama, Yasuhiro Hasegawa, Takashi Tsukagoshi, Mihoko Konishi, Motohide Tamura

Abstract: Crescent-like asymmetric dust structures discovered in protoplanetary disks indicate dust aggregations. Thus, the research on them helps us understand the planet formation process. Here we analyze the ALMA data of the protoplanetary disk around the T-Tauri star SR 21, which has asymmetric structures detected in previous sub-millimeter observations. Imaged at ALMA Band 6 (1.3 mm) with a spatial res… ▽ More Crescent-like asymmetric dust structures discovered in protoplanetary disks indicate dust aggregations. Thus, the research on them helps us understand the planet formation process. Here we analyze the ALMA data of the protoplanetary disk around the T-Tauri star SR 21, which has asymmetric structures detected in previous sub-millimeter observations. Imaged at ALMA Band 6 (1.3 mm) with a spatial resolution of about 0.$\arcsec$04, the disk is found to consist of two rings and three asymmetric structures, with two of the asymmetric structures being in the same ring. Compared to the Band 6 image, the Band 3 (2.7 mm) image also shows the three asymmetric structures but with some clumps. The elongated asymmetric structures in the outer ring could be due to the interactions of a growing planet. By fitting the Band 3 and Band 6 dust continuum data, two branches of solutions of maximum dust size in the disk are suggested: one is larger than 1 mm, and the other is smaller than 300 $μm$. High-resolution continuum observations at longer wavelengths as well as polarization observations can help break the degeneracy. We also suggest that the prominent spiral previously identified in VLT/SPHERE observations to the south of the star at 0.$\arcsec$25 may be the scattered light counterpart of the Inner Arc, and the structure is a dust-trapping vortex in nature. The discovered features in SR 21 make it a good target for studying the evolution of asymmetric structures and planet formation. △ Less

Submitted 9 March, 2023; originally announced March 2023.

Comments: Accepted by The Astrophysical Journal

The ALMA view of MP Mus (PDS 66): a protoplanetary disk with no visible gaps down to 4 au scales

Authors: Á. Ribas, E. Macías, P. Weber, S. Pérez, N. Cuello, R. Dong, A. Aguayo, C. Cáceres, J. Carpenter, W. R. F. Dent, I. de Gregorio-Monsalvo, G. Duchêne, C. C. Espaillat, P. Riviere-Marichalar, M. Villenave

Abstract: We present ALMA multiwavelength observations of the protoplanetary disk around the nearby (d$\sim$100 pc) young solar analog MP Mus (PDS 66). These observations at 0.89 mm, 1.3 mm, and 2.2 mm have angular resolutions of $\sim$ 1", 0.05", and 0.25", respectively, and probe the dust and gas in the system with unprecedented detail and sensitivity. The disk appears smooth down to the 4 au resolution o… ▽ More We present ALMA multiwavelength observations of the protoplanetary disk around the nearby (d$\sim$100 pc) young solar analog MP Mus (PDS 66). These observations at 0.89 mm, 1.3 mm, and 2.2 mm have angular resolutions of $\sim$ 1", 0.05", and 0.25", respectively, and probe the dust and gas in the system with unprecedented detail and sensitivity. The disk appears smooth down to the 4 au resolution of the 1.3 mm observations, in contrast with most disks observed at comparable spatial scales. The dust disk has a radius of 60$\pm$5 au, a dust mass of $0.14_{-0.06}^{+0.11} M_{\rm Jup}$, and a mm spectral index $<2$ in the inner 30 au, suggesting optically thick emission from grains with high albedo in this region. Several molecular gas lines are also detected extending up to 130$\pm$15 au, similar to small grains traced by scattered light observations. Comparing the fluxes of different CO isotopologues with previous models yields a gas mass of $0.1-1 M_{\rm Jup}$, implying a gas to dust ratio of 1-10. We also measure a dynamical stellar mass of $M_{\rm dyn}$=1.30$\pm$0.08 $M_\odot$ and derive an age of 7-10 Myr for the system. The survival of large grains in an evolved disk without gaps/rings is surprising, and it is possible that existing substructures remain undetected due to optically thick emission at 1.3 mm. Alternatively, small structures may still remain unresolved with the current observations. Based on simple scaling relations for gap-opening planets and gap widths, this lack of substructures places upper limits to the masses of planets in the disk as low as 2 $M_\oplus$-0.06 $M_{\rm Jup}$ at $r > 40$ au. The lack of mm emission at radii $r > 60$ au also suggests that the gap in scattered light between 30-80 au is likely not a gap in the disk density, but a shadow cast by a puffed-up inner disk. △ Less

Submitted 22 February, 2023; originally announced February 2023.

Comments: 17 pages + 3 appendices. Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 673, A77 (2023)

ALMA Band 6 high-resolution observations of the transitional disk around SY Cha

Authors: Ryuta Orihara, Munetake Momose, Takayuki Muto, Jun Hashimoto, Hauyu Baobab Liu, Takashi Tsukagoshi, Tomoyuki Kudo, Sanemichi Takahashi, Yi Yang, Yasuhiro Hasegawa, Ruobing Dong, Mihoko Konishi, Eiji Akiyama

Abstract: In this study, we reported the results of high-resolution (0.14 arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 225 GHz dust continuum and CO molecular emission lines from the transitional disk around SY Cha. Our high-resolution observations clearly revealed the inner cavity and the central point source for the first time. The radial profile of the ring can be appro… ▽ More In this study, we reported the results of high-resolution (0.14 arcsec) Atacama Large Millimeter/submillimeter Array (ALMA) observations of the 225 GHz dust continuum and CO molecular emission lines from the transitional disk around SY Cha. Our high-resolution observations clearly revealed the inner cavity and the central point source for the first time. The radial profile of the ring can be approximated by a bright narrow ring superimposed on a fainter wide ring. Furthermore, we found that there is a weak azimuthal asymmetry in dust continuum emission. For gas emissions, we detected $\rm{}^{12}CO$(2$-$1), $\rm{}^{13}CO$(2$-$1) and $\rm{}C^{18}O$(2$-$1), from which we estimated the total gas mass of the disk to be $2.2\times10^{-4}M_\odot$, assuming a CO/H$_2$ ratio of $10^{-4}$. The observations showed that the gas is present inside the dust cavity. The analysis of the velocity structure of the $\rm{}^{12}CO$(2$-$1) emission line revealed that the velocity is distorted at the location of the dust inner disk, which may be owing to warping of the disk or radial gas flow within the cavity of the dust disk. High-resolution observations of SY Cha showed that this system is composed of a ring and a distorted inner disk, which may be common, as indicated by the survey of transitional disk systems at a resolution of $\sim$0.1~arcsec. △ Less

Submitted 11 February, 2023; originally announced February 2023.

Comments: 33 pages, 22 figures

Testing velocity kinks as a planet-detection method: Do velocity kinks in surface gas emission trace planetary spiral wakes in the midplane continuum?

Authors: Jessica Speedie, Ruobing Dong

Abstract: 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… ▽ More 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. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 5 figures, 13 pages. Accepted for publication in ApJ Letters. Supplementary material available at https://doi.org/10.6084/m9.figshare.21330426

Grain Growth in the Dust Ring with Crescent around Very Low Mass Star ZZ Tau IRS with JVLA

Authors: Jun Hashimoto, Hauyu Baobab Liu, Ruobing Dong, Beibei Liu, Takayuki Muto

Abstract: 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 t… ▽ More 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~$μ$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. △ Less

Submitted 2 November, 2022; originally announced November 2022.

Comments: 11 pages, 3 figures, accepted in ApJ

Gap-opening Planets Make Dust Rings Wider

Authors: Jiaqing Bi, Min-Kai Lin, Ruobing Dong

Abstract: 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 d… ▽ More 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 $α_{\rm turb}$, but with cautions about the Schmidt number greater than order unity. △ Less

Submitted 20 October, 2022; originally announced October 2022.

Comments: 17 pages, 6+3 figures, 2 tables, accepted in ApJ

Stellar Flyby Analysis for Spiral Arm Hosts with Gaia DR3

Authors: Linling Shuai, Bin B. Ren, Ruobing Dong, Xingyu Zhou, Laurent Pueyo, Robert J. De Rosa, Taotao Fang, Dimitri Mawet

Abstract: 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-drivin… ▽ More 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. △ Less

Submitted 7 October, 2022; originally announced October 2022.

Comments: 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

Gap Opening and Inner Disk Structure in the Strongly Accreting Transition Disk of DM Tau

Authors: Logan Francis, Nienke van der Marel, Doug Johnstone, Eiji Akiyama, Simon Bruderer, Ruobing Dong, Jun Hashimoto, Hauyu Baobab Liu, Takayuki Muto, Yi Yang

Abstract: 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 undep… ▽ More 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 ($α=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 $<1$ M$_\mathrm{Jup}$ implied by the gap depth is in tension with predictions for dust trapping in a highly viscous disk, which requires a more massive planet of of $\sim10$M$_\mathrm{Jup}$. Photoevaporative models including a dead zone can qualitatively reproduce some features of the DM Tau disk, but still struggle to explain the high accretion rates and the observed mm continuum flux. △ Less

Submitted 2 August, 2022; originally announced August 2022.

Comments: Accepted for publication in the Astronomical Journal. 36 pages, 28 figures

Distributions of gas and small and large grains in the LkH$α\,330$ disk trace a young planetary system

Authors: P. Pinilla, M. Benisty, N. T. Kurtovic, J. Bae, R. Dong, Z. Zhu, S. Andrews, J. Carpenter, C. Ginski, J. Huang, A. Isella, L. Pérez, L. Ricci, G. Rosotti, M. Villenave, D. Wilner

Abstract: [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''… ▽ More [abridged] We present new scattered light and millimeter observations of the protoplanetary disk around LkH$α\,330$, using SPHERE/VLT and ALMA, respectively. The scattered-light SPHERE observations reveal an asymmetric ring at around 45au from the star in addition to two spiral arms with similar radial launching points at around 90au. The millimeter observations from ALMA (resolution of 0.06''$\times$0.04'') mainly show an asymmetric ring located at 110au from the star. In addition to this asymmetry, there are two faint symmetric rings at 60au and 200au. The $^{12}$CO, $^{13}$CO, and C$^{18}$O lines seem to be less abundant in the inner disk (these observations have a resolution of 0.16''$\times$0.11''). The $^{13}$CO peaks at a location similar to the inner ring observed with SPHERE, suggesting that this line is optically thick and traces variations of disk temperature instead of gas surface-density variations, while the C$^{18}$O peaks slightly further away at around 60au. We compare our observations with hydrodynamical simulations that include gas and dust evolution, and conclude that a 10$M_{\rm{Jup}}$ mass planet at 60au and in an eccentric orbit ($e=0.1$) can qualitatively explain most of the observed structures. A planet in a circular orbit leads to a much narrower concentration in the millimeter emission, while a planet in a more eccentric orbit leads to a very eccentric cavity as well. In addition, the outer spiral arm launched by the planet changes its pitch angle along the spiral due to the eccentricity and when it interacts with the vortex, potentially appearing in observations as two distinct spirals. Our observations and models show that LkH$α\,330$ is an interesting target to search for (eccentric-) planets while they are still embedded in their parental disk, making it an excellent candidate for studies on planet-disk interaction. △ Less

Submitted 25 July, 2022; v1 submitted 20 June, 2022; originally announced June 2022.

Comments: Accepted for publication in A&A, minor changes after language edition

Journal ref: A&A 665, A128 (2022)

Substructures in protoplanetary disks imprinted by compact planetary systems

Authors: Juan Garrido-Deutelmoser, Cristobal Petrovich, Leonardo Krapp, Kaitlin M. Kratter, Ruobing Dong

Abstract: 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… ▽ More 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 ($Δ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 ($Δa \lesssim 4.6\; H_{\rm p}$); ii) partially-shared gaps for more widely spaced planets ($Δ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. △ Less

Submitted 19 April, 2022; originally announced April 2022.

Comments: 17 pages, 7 figures, 1 table, submitted to The Astrophysical Journal

Two Rings and a Marginally Resolved, 5 AU, Disk Around LkCa 15 Identified Via Near Infrared Sparse Aperture Masking Interferometry

Authors: Dori Blakely, Logan Francis, Doug Johnstone, Anthony Soulain, Peter Tuthill, Anthony Cheetham, Joel Sanchez-Bermudez, Anand Sivaramakrishnan, Ruobing Dong, Nienke van der Marel, Rachel Cooper, Arthur Vigan, Faustine Cantalloube

Abstract: 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 exten… ▽ More 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. △ Less

Submitted 14 April, 2022; originally announced April 2022.

Comments: Accepted for publication in ApJ

Images of Embedded Jovian Planet Formation At A Wide Separation Around AB Aurigae

Authors: Thayne Currie, Kellen Lawson, Glenn Schneider, Wladimir Lyra, John Wisniewski, Carol Grady, Olivier Guyon, Motohide Tamura, Takayuki Kotani, Hajime Kawahara, Timothy Brandt, Taichi Uyama, Takayuki Muto, Ruobing Dong, Tomoyuki Kudo, Jun Hashimoto, Misato Fukagawa, Kevin Wagner, Julien Lozi, Jeffrey Chilcote, Taylor Tobin, Tyler Groff, Kimberly Ward-Duong, William Januszewski, Barnaby Norris , et al. (8 additional authors not shown)

Abstract: 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 r… ▽ More 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. △ Less

Submitted 1 April, 2022; originally announced April 2022.

Comments: 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

The Appearance of Vortices in Protoplanetary Disks in Near-Infrared Scattered Light

Authors: Metea Marr, Ruobing Dong

Abstract: 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 co… ▽ More 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. △ Less

Submitted 22 March, 2022; originally announced March 2022.

Comments: 10 pages, 7 figures, ApJ accepted

Planet-Disk Interactions

Authors: Sijme-Jan Paardekooper, Ruobing Dong, Paul Duffell, Jeffrey Fung, Frederic S. Masset, Gordon Ogilvie, Hidekazu Tanaka

Abstract: Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments.… ▽ More Planet-disk interactions, where an embedded massive body interacts gravitationally with the protoplanetary disk it was formed in, can play an important role in reshaping both the disk and the orbit of the planet. Spiral density waves are launched into the disk by the planet, which, if they are strong enough, can lead to the formation of a gap. Both effects are observable with current instruments. The back-reaction of perturbations induced in the disk, both wave-like and non-wavelike, is a change in orbital elements of the planet. The efficiency of orbital migration is a long-standing problem in planet formation theory. We discuss recent progress in planet-disk interactions for different planet masses and disk parameters, in particular the level of turbulence, and progress in modeling observational signatures of embedded planets. △ Less

Submitted 8 November, 2022; v1 submitted 17 March, 2022; originally announced March 2022.

Comments: 28 pages, 12 figures, chapter to appear in Protostars and Planets VII, eds. Shu-ichiro Inutsuka, Yuri Aikawa, Takayuki Muto, Kengo Tomida, and Motohide Tamura. Updated version including community input

Observing planet-driven dust spirals with ALMA

Authors: Jessica Speedie, Richard A. Booth, Ruobing Dong

Abstract: ALMA continuum observations of thermal emission from the dust component of protoplanetary disks have revealed an abundance of substructures that may be interpreted as evidence for embedded planets, but planet-driven spiral arms -- perhaps one of the most compelling lines of evidence -- have proven comparatively elusive. In this work, we test the capabilities of ALMA to detect the planet-driven spi… ▽ More ALMA continuum observations of thermal emission from the dust component of protoplanetary disks have revealed an abundance of substructures that may be interpreted as evidence for embedded planets, but planet-driven spiral arms -- perhaps one of the most compelling lines of evidence -- have proven comparatively elusive. In this work, we test the capabilities of ALMA to detect the planet-driven spiral signal in continuum emission. Carrying out hydrodynamic simulations and radiative transfer calculations, we present synthetic Band 7 continuum images for a wide range of disk and observing conditions. We show that thermal mass planets at tens of au typically drive spirals detectable within a few hours of integration time, and the detectable planet mass may be as low as $\sim$Neptune mass ($0.3 \, M_{\rm th}$). The grains probed by ALMA form spirals morphologically identical to the underlying gas spiral. The temperature of the dust spiral is crucial in determining its contrast, and spirals are easier to detect in disks with an adiabatic equation of state and longer cooling times. Resolving the spiral is not necessary for its detection; with the help of residual maps, the optimal beam size is a few times the spiral width at a constant noise level. Finally, we show how the presence of gaps and rings can impair our ability to recognize co-located spirals. Our work demonstrates the planet-finding potential of the current design specification of ALMA, and suggests that observing capability is not the bottleneck in searching for spirals induced by thermal mass planets. △ Less

Submitted 1 March, 2022; originally announced March 2022.

Comments: 11 figures, 20 pages. Accepted for publication in The Astrophysical Journal. Image repository available at https://doi.org/10.6084/m9.figshare.19148912

A likely flyby of binary protostar Z CMa caught in action

Authors: Ruobing Dong, Hauyu Baobab Liu, Nicolas Cuello, Christophe Pinte, Peter Abraham, Eduard Vorobyov, Jun Hashimoto, Agnes Kospal, Eugene Chiang, Michihiro Takami, Lei Chen, Michael Dunham, Misato Fukagawa, Joel Green, Yasuhiro Hasegawa, Thomas Henning, Yaroslav Pavlyuchenkov, Tae-Soo Pyo, Motohide Tamura

Abstract: 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 m… ▽ More 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. △ Less

Submitted 20 February, 2022; v1 submitted 14 January, 2022; originally announced January 2022.

Comments: Published in Nature Astronomy. Here is the authors' version with the Supplementary Information integrated into the Methods section

Improving Planet Detection with Disk Modeling: Keck/NIRC2 Imaging of the HD 34282 Single-armed Protoplanetary Disk

Authors: Juan Quiroz, Nicole L. Wallack, Bin Ren, Ruobing Dong, Jerry W. Xuan, Dimitri Mawet, Maxwell A. Millar-Blanchaer, Garreth Ruane

Abstract: 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 necessa… ▽ More 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. △ Less

Submitted 24 November, 2021; originally announced November 2021.

Comments: 8 pages, 3 figures, 1 table, ApJ Letters accepted. Data files available in the ancillary folder

Journal ref: ApJL 924 (2022) L4

UV Spectropolarimetry with Polstar: Protoplanetary Disks

Authors: John P. Wisniewski, Andrei V. Berdyugin, Svetlana V. Berdyugina, William C. Danchi, Ruobing Dong, Rene D. Oudmaijer, Vladimir S. Airapetian, Sean D. Brittain, Ken Gayley, Richard Ignace, Maud Langlois, Kellen D. Lawson, Jamie R. Lomax, Motohide Tamura, Jorick S. Vink, Paul A. Scowen

Abstract: Polstar is a proposed NASA MIDEX mission that would feature a high resolution UV spectropolarimeter capable of measure all four Stokes parameters onboard a 60cm telescope. The mission would pioneer the field of time-domain UV spectropolarimetry. Time domain UV spectropolarimetry offers the best resource to determine the geometry and physical conditions of protoplanetary disks from the stellar surf… ▽ More Polstar is a proposed NASA MIDEX mission that would feature a high resolution UV spectropolarimeter capable of measure all four Stokes parameters onboard a 60cm telescope. The mission would pioneer the field of time-domain UV spectropolarimetry. Time domain UV spectropolarimetry offers the best resource to determine the geometry and physical conditions of protoplanetary disks from the stellar surface to <5 AU. We detail two key objectives that a dedicated time domain UV spectropolarimetry survey, such as that enabled by Polstar, could achieve: 1) Test the hypothesis that magneto-accretion operating in young planet-forming disks around lower-mass stars transitions to boundary layer accretion in planet-forming disks around higher mass stars; and 2) Discriminate whether transient events in the innermost regions of planet-forming disks of intermediate mass stars are caused by inner disk mis-alignments or from stellar or disk emissions. △ Less

Submitted 9 December, 2021; v1 submitted 12 November, 2021; originally announced November 2021.

Comments: 16 pages, 7 figures

CI Tau: A controlled experiment in disk-planet interaction

Authors: Dhruv Muley, Ruobing Dong

Abstract: CI Tau is a young (~2 Myr) T Tauri system with a substantial near-infrared excess in its SED, indicating that the protoplanetary disk extends very close to its star. This is seemingly at odds with the radial-velocity discovery of CI Tau b, a ~12 MJ planet at ~0.1 au, which would be expected to carve a wide, deep cavity in the innermost disk. To investigate this apparent contradiction, we run 2D hy… ▽ More CI Tau is a young (~2 Myr) T Tauri system with a substantial near-infrared excess in its SED, indicating that the protoplanetary disk extends very close to its star. This is seemingly at odds with the radial-velocity discovery of CI Tau b, a ~12 MJ planet at ~0.1 au, which would be expected to carve a wide, deep cavity in the innermost disk. To investigate this apparent contradiction, we run 2D hydrodynamics simulations to study the effect of the planet on the disk, then post-process the results with radiative transfer to obtain an SED. We find that at ~0.1 au, even such a massive companion has little impact on the near-infrared excess, a result that holds regardless of planetary eccentricity and dust size distribution. Conversely, the observed full-disk signature in CI Tau's SED is consistent with the existence of the hot super-Jupiter CI Tau b. As our simulations uncover, clear transition-disk signatures in SEDs are more likely to be signposts of nascent "warm" Jupiters, located at around 1 AU in the future habitable zones of their host stars. △ Less

Submitted 27 October, 2021; v1 submitted 25 October, 2021; originally announced October 2021.

Comments: Accepted to ApJL. Comments and questions welcome. v2: fixed abstract to reflect accepted version

GW Ori: circumtriple rings and planets

Authors: Jeremy L. Smallwood, Rebecca Nealon, Cheng Chen, Rebecca G. Martin, Jiaqing Bi, Ruobing Dong, Christophe Pinte

Abstract: GW Ori is a hierarchical triple star system with a misaligned circumtriple protoplanetary disc. Recent ALMA observations have identified three dust rings with a prominent gap at $100\, \rm au$ and misalignments between each of the rings. A break in the gas disc may be driven either by the torque from the triple star system or a planet that is massive enough to carve a gap in the disc. Once the dis… ▽ More GW Ori is a hierarchical triple star system with a misaligned circumtriple protoplanetary disc. Recent ALMA observations have identified three dust rings with a prominent gap at $100\, \rm au$ and misalignments between each of the rings. A break in the gas disc may be driven either by the torque from the triple star system or a planet that is massive enough to carve a gap in the disc. Once the disc is broken, the rings nodally precess on different timescales and become misaligned. We investigate the origins of the dust rings by means of $N$-body integrations and 3-dimensional hydrodynamic simulations. We find that for observationally-motivated parameters of protoplanetary discs, the disc does not break due to the torque from the star system. We suggest that the presence of a massive planet (or planets) in the disc separates the inner and outer disc. We conclude that the disc breaking in GW Ori is likely caused by undetected planets -- the first planet(s) in a circumtriple orbit. △ Less

Submitted 20 September, 2021; originally announced September 2021.

Comments: 17 pages, 16 figures, accepted to MNRAS

Keck/OSIRIS Pa$β$ high-contrast imaging and updated constraints on PDS~70b

Authors: Taichi Uyama, Chen Xie, Yuhiko Aoyama, Charles A. Beichman, Jun Hashimoto, Ruobing Dong, Yasuhiro Hasegawa, Masahiro Ikoma, Dimitri Mawet, Michael W. McElwain, Jean-Baptiste Ruffio, Kevin R. Wagner, Jason J. Wang, Yifan Zhou

Abstract: We present a high-contrast imaging search for Pa$β$ 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$β$ line at the level predicted using the parameters of \cite{Hashimoto2020}. This lack of Pa$β$ emission suggests the MUSE-b… ▽ More We present a high-contrast imaging search for Pa$β$ 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$β$ line at the level predicted using the parameters of \cite{Hashimoto2020}. This lack of Pa$β$ emission suggests the MUSE-based study may have overestimated the line width of H$α$. We compared our Pa$β$ detection limits with the previous H$α$ flux and H$β$ 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. △ Less

Submitted 14 September, 2021; originally announced September 2021.

Comments: 8pages, 5 figs, accepted for publication in AJ

Investigating Protoplanetary Disk Cooling through Kinematics: Analytical GI Wiggle

Authors: Cristiano Longarini, Giuseppe Lodato, Claudia Toci, Benedetta Veronesi, Cassandra Hall, Ruobing Dong, Jason Patrick Terry

Abstract: It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic… ▽ More It is likely that young protostellar disks undergo a self-gravitating phase. Such systems are characterized by the presence of a spiral pattern that can be either in a quasi-steady state or in a nonlinear unstable condition. This spiral wave affects both the gas dynamics and kinematics, resulting in deviations from the Keplerian rotation. Recently, a lot of attention has been devoted to kinematic studies of planet-forming environments, and we are now able to measure even small perturbations of velocity field ($\lesssim 1 \%$ of the Keplerian speed) thanks to high spatial and spectral resolution observations of protostellar disks. In this work, we investigate the kinematic signatures of gravitational instability: we perform an analytical study of the linear response of a self-gravitating disk to a spiral-like perturbation, focusing our attention on the velocity field perturbations. We show that unstable disks have clear kinematic imprints into the gas component across the entire disk extent, due to the GI spiral wave perturbation, resulting in deviations from Keplerian rotation. The shape of these signatures depends on several parameters, but they are significantly affected by the cooling factor: by detecting these features, we can put constraints on protoplanetary disk cooling. △ Less

Submitted 11 January, 2022; v1 submitted 25 August, 2021; originally announced August 2021.

Comments: Accepted for publication in APJL, 13 pages, 4 figures

Observational Signatures of Planets in Protoplanetary Disks: Temperature structures in spiral arms

Authors: Dhruv Muley, Ruobing Dong, Jeffrey Fung

Abstract: High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with $β$-cooling to a vertically-stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a c… ▽ More High-resolution imaging of protoplanetary disks has unveiled a rich diversity of spiral structure, some of which may arise from disk-planet interaction. Using 3D hydrodynamics with $β$-cooling to a vertically-stratified background, as well as radiative-transfer modeling, we investigate the temperature rise in planet-driven spirals. In rapidly cooling disks, the temperature rise is dominated by a contribution from stellar irradiation, 0.3-3% inside the planet radius but always <0.5% outside. When cooling time equals or exceeds dynamical time, however, this is overwhelmed by hydrodynamic PdV work, which introduces a 10-20% perturbation within a factor of 2 from the planet's orbital radius. We devise an empirical fit of the spiral amplitude $Δ(T)$ to take into account both effects. Where cooling is slow, we find also that temperature perturbations from buoyancy spirals -- a strictly 3D, non-isothermal phenomenon -- become nearly as strong as those from Lindblad spirals, which are amenable to 2D and isothermal studies. Our findings may help explain observed thermal features in disks like TW Hydrae and CQ Tauri, and underscore that 3D effects have a qualitatively important effect on disk structure. △ Less

Submitted 13 July, 2021; originally announced July 2021.

Comments: Accepted to AJ. Comments and questions welcome

Massive compact disks around FU Orionis-type young eruptive stars revealed by ALMA

Authors: Á. Kóspál, F. Cruz-Sáenz de Miera, J. A. White, P. Ábrahám, L. Chen, T. Csengeri, R. Dong, M. M. Dunham, O. Fehér, J. D. Green, J. Hashimoto, Th. Henning, M. Hogerheijde, T. Kudo, H. B. Liu, M. Takami, E. I. Vorobyov

Abstract: 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… ▽ More 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. △ Less

Submitted 28 June, 2021; originally announced June 2021.

Comments: 34 pages, 8 tables, 29 figures, accepted for publication in the Astrophysical Journal Supplement Series

A dusty filament and turbulent CO spirals in HD135344B-SAO206462

Authors: Simon Casassus, Valentin Christiaens, Miguel Carcamo, Sebastian Perez, Philipp Weber, Barbara Ercolano, Nienke van der Marel, Christophe Pinte, Ruobing Dong, Clement Baruteau, Lucas Cieza, Ewine van Dishoeck, Andres Jordan, Daniel Price, Olivier Absil, Carla Arce-Tord, Virginie Faramaz, Christian Flores, Maddalena Reggiani

Abstract: 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.… ▽ More 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. △ Less

Submitted 12 August, 2021; v1 submitted 16 April, 2021; originally announced April 2021.

Comments: accepted to MNRAS

An asymmetric dust ring around a very low mass star ZZ Tau IRS

Authors: Jun Hashimoto, Ruobing Dong, Takayuki Muto

Abstract: We present Atacama Large Millimeter/submillimeter Array (ALMA) gas and dust observations at band 7 (339~GHz: 0.89~mm) of the protoplanetary disk around a very low mass star ZZ~Tau~IRS with a spatial resolution of 0\farcs25. The $^{12}$CO~$J=3\rightarrow2$ position--velocity diagram suggests a dynamical mass of ZZ~Tau~IRS of $\sim$0.1--0.3~$M_{\sun}$. The disk has a total flux density of 273.9 mJy,… ▽ More We present Atacama Large Millimeter/submillimeter Array (ALMA) gas and dust observations at band 7 (339~GHz: 0.89~mm) of the protoplanetary disk around a very low mass star ZZ~Tau~IRS with a spatial resolution of 0\farcs25. The $^{12}$CO~$J=3\rightarrow2$ position--velocity diagram suggests a dynamical mass of ZZ~Tau~IRS of $\sim$0.1--0.3~$M_{\sun}$. The disk has a total flux density of 273.9 mJy, corresponding to an estimated mass of 24--50~$M_\oplus$ in dust. The dust emission map shows a ring at $r=$ 58~au and an azimuthal asymmetry at $r=$ \jh{45}~au with a position angle of 135\degr. The properties of the asymmetry, including radial width, aspect ratio, contrast, and contribution to the total flux, were found to be similar to the asymmetries around intermediate mass stars ($\sim$2~$M_{\sun}$) such as MWC~758 and IRS~48. This implies that the asymmetry in the ZZ~Tau~IRS disk shares a similar origin with others, despite the star being $\sim$10 times less massive. Our observations also suggest that the inner and outer parts of the disk may be misaligned. Overall, the ZZ~Tau~IRS disk shows evidence of giant planet formation at $\sim$10 au scale at a few Myr. If confirmed, it will challenge existing core accretion models, in which such planets have been predicted to be extremely hard to form around very low mass stars. △ Less

Submitted 30 March, 2021; originally announced March 2021.

Comments: 19 pages, 7 figures, accepted in AJ

Puffed up Edges of Planet-opened Gaps in Protoplanetary Disks. I. hydrodynamic simulations

Authors: Jiaqing Bi, Min-Kai Lin, Ruobing Dong

Abstract: 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… ▽ More 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. △ Less

Submitted 16 March, 2021; originally announced March 2021.

Comments: 15 pages, 7+2 figures, accepted in ApJ

Perturbers: SPHERE detection limits to planetary-mass companions in protoplanetary disks

Authors: R. Asensio-Torres, Th. Henning, F. Cantalloube, P. Pinilla, D. Mesa, A. Garufi, S. Jorquera, R. Gratton, G. Chauvin, J. Szulagyi, R. van Boekel, R. Dong, G. -D. Marleau, M. Benisty, M. Villenave, C. Bergez-Casalou, C. Desgrange, M. Janson, M. Keppler, M. Langlois, F. Menard, E. Rickman, T. Stolker, M. Feldt, T. Fusco , et al. (3 additional authors not shown)

Abstract: 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… ▽ More 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. △ Less

Submitted 9 March, 2021; originally announced March 2021.

Comments: 26 pages, 33 figures, accepted for publication in A&A

Journal ref: A&A 652, A101 (2021)

ALMA Observations of the Asymmetric Dust Disk around DM Tau

Authors: Jun Hashimoto, Takayuki Muto, Ruobing Dong, Hauyu Baobab Liu, Nienke van der Marel, Logan Francis, Yasuhiro Hasegawa, Takashi Tsukagoshi

Abstract: 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~$μ$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 con… ▽ More 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~$μ$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~$μ$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. △ Less

Submitted 11 February, 2021; originally announced February 2021.

Comments: 26 pages, 17 figures

ALMA observation of the protoplanetary disk around WW Cha: faint double-peaked ring and asymmetric structure

Authors: Kazuhiro D. Kanagawa, Jun Hashimoto, Takayuki Muto, Takashi Tsukagoshi, Sanemichi Z. Takahashi, Yasuhiro Hasegawa, Mihoko Konishi, Hideko Nomura, Hauyu Baobab Liu, Ruobing Dong, Akimasa Kataoka, Munetake Momose, Tomohiro Ono, Michael Sitko, Michihiro Takami, Kengo Tomida

Abstract: 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… ▽ More 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. △ Less

Submitted 17 March, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: 24 pages, 25 figures, accepted for publication in The Astrophysical Journal