Source author record

Ruobing Dong

Ruobing Dong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

55works
6topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

55 published item(s)

preprint2024arXiv

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

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.

preprint2024arXiv

Shadowing in the protoplanetary disk of ZZ Tau IRS with HST

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2022arXiv

Investigating Protoplanetary Disk Cooling through Kinematics: Analytical GI Wiggle

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.

preprint2022arXiv

Observing planet-driven dust spirals with ALMA

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.

preprint2022arXiv

Substructures in protoplanetary disks imprinted by compact planetary systems

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.

preprint2022arXiv

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

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.

preprint2022arXiv

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

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.

preprint2021arXiv

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

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.

preprint2020arXiv

ALMA Observations of the Inner Cavity in the Protoplanetary Disk around Sz 84

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of a protoplanetary disk around the T Tauri star Sz~84 and analyses of the structures of the inner cavity in the central region of the dust disk. Sz~84's spectral energy distribution (SED) has been known to exhibit negligible infrared excess at $λ\lesssim$10~$μ$m due to the disk's cavity structure. Analyses of the observed visibilities of dust continuum at 1.3~mm and the SED indicate that the size of the cavity in the disk of large (millimeter size) dust grains is 8~au in radius and that in the disk of small (sub-micron size) dust grains is 60~au in radius. Furthermore, from the SED analyses, we estimate that the upper limit mass of small dust grains at $r<$60~au is less than $\sim$10$^{-3}$~$M_{\rm earth}$, which is $\lesssim$0.01~\% of the total (small~$+$~large) dust mass at $r<$60~au. These results suggest that large dust grains are dominant at $r<$60~au, implying that dust grains efficiently grow with less efficient fragmentation in this region, potentially due to weak turbulence and/or stickier dust grains. The balance of grain growth and dust fragmentation is an important factor for determining the size of large dust grains in protoplanetary disks, and thus Sz~84 could serve as a good testbed for investigations of grain growth in such disks.

preprint2020arXiv

Debris Disk Results from the Gemini Planet Imager Exoplanet Survey's Polarimetric Imaging Campaign

We report the results of a ${\sim}4$-year direct imaging survey of 104 stars to resolve and characterize circumstellar debris disks in scattered light as part of the Gemini Planet Imager Exoplanet Survey. We targeted nearby (${\lesssim}150$ pc), young (${\lesssim}500$ Myr) stars with high infrared excesses ($L_{\mathrm{IR}} / L_\star > 10^{-5}$), including 38 with previously resolved disks. Observations were made using the Gemini Planet Imager high-contrast integral field spectrograph in $H$-band (1.6 $μ$m) coronagraphic polarimetry mode to measure both polarized and total intensities. We resolved 26 debris disks and three protoplanetary/transitional disks. Seven debris disks were resolved in scattered light for the first time, including newly presented HD 117214 and HD 156623, and we quantified basic morphologies of five of them using radiative transfer models. All of our detected debris disks but HD 156623 have dust-poor inner holes, and their scattered-light radii are generally larger than corresponding radii measured from resolved thermal emission and those inferred from spectral energy distributions. To assess sensitivity, we report contrasts and consider causes of non-detections. Detections were strongly correlated with high IR excess and high inclination, although polarimetry outperformed total intensity angular differential imaging for detecting low inclination disks (${\lesssim} 70 °$). Based on post-survey statistics, we improved upon our pre-survey target prioritization metric predicting polarimetric disk detectability. We also examined scattered-light disks in the contexts of gas, far-IR, and millimeter detections. Comparing $H$-band and ALMA fluxes for two disks revealed tentative evidence for differing grain properties. Finally, we found no preference for debris disks to be detected in scattered light if wide-separation substellar companions were present.

preprint2020arXiv

Dynamical Evidence of a Spiral Arm--Driving Planet in the MWC 758 Protoplanetary Disk

More than a dozen young stars host spiral arms in their surrounding protoplanetary disks. The excitation mechanisms of such arms are under debate. The two leading hypotheses -- companion-disk interaction and gravitational instability (GI) -- predict distinct motion for spirals. By imaging the MWC 758 spiral arm system at two epochs spanning ${\sim}5$ yr using the SPHERE instrument on the Very Large Telescope (VLT), we test the two hypotheses for the first time. We find that the pattern speeds of the spirals are not consistent with the GI origin. Our measurements further evince the existence of a faint "missing planet" driving the disk arms. The average spiral pattern speed is $0.\!^\circ22\pm0.\!^\circ03$ yr$^{-1}$, pointing to a driver at $172_{-14}^{+18}$ au around a $1.9$ $M_\odot$ central star if it is on a circular orbit. In addition, we witness time varying shadowing effects on a global scale that are likely originated from an inner disk.

preprint2020arXiv

First Images of the Protoplanetary Disk Around PDS 201

Scattered light imaging has revealed nearly a dozen circumstellar disks around young Herbig Ae/Be stars$-$enabling studies of structures in the upper disk layers as potential signs of on-going planet formation. We present the first images of the disk around the variable Herbig Ae star PDS 201 (V* V351 Ori), and an analysis of the images and spectral energy distribution through 3D Monte-Carlo radiative transfer simulations and forward modelling. The disk is detected in three datasets with LBTI/LMIRCam at the LBT, including direct observations in the $Ks$ and $L'$ filters, and an $L'$ observation with the 360$^\circ$ vector apodizing phase plate coronagraph. The scattered light disk extends to a very large radius of $\sim$250 au, which places it among the largest of such disks. Exterior to the disk, we establish detection limits on substellar companions down to $\sim$5 M$_{Jup}$ at $\gtrsim$1.5" ($\gtrsim$500 au), assuming the Baraffe et al. (2015) models. The images show a radial gap extending to $\sim$0.4" ($\sim$140 au at a distance of 340 pc) that is also evident in the spectral energy distribution. The large gap is a possible signpost of multiple high-mass giant planets at orbital distances ($\sim$60-100 au) that are unusually massive and widely-separated compared to those of planet populations previously inferred from protoplanetary disk substructures.

preprint2020arXiv

GW Ori: Interactions Between a Triple-star System and its Circumtriple Disk in Action

GW Ori is a hierarchical triple system which has a rare circumtriple disk. We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of 1.3 mm dust continuum and 12CO J=2-1 molecular gas emission of the disk. For the first time, we identify three dust rings in the disk at ~46, 188, and 338 AU, with estimated dust mass of ~70-250 Earth masses, respectively. To our knowledge, the outer ring in GW Ori is the largest dust ring ever found in protoplanetary disks. We use visibility modelling of dust continuum to show that the disk has misaligned parts and the innermost dust ring is eccentric. The disk misalignment is also suggested by the CO kinematics modelling. We interpret these substructures as evidence of ongoing dynamical interactions between the triple stars and the circumtriple disk.

preprint2020arXiv

Long Baseline Observations of the HD100546 Protoplanetary Disk with ALMA

Using the Atacama Large Millimeter/submillimeter Array (ALMA), we observed the young Herbig star HD 100546, host to a prominent disk with a deep, wide gap in the dust. The high-resolution 1.3 mm continuum observation reveals fine radial and azimuthal substructures in the form of a complex maze of ridges and trenches sculpting a dust ring. The $^{12}$CO(2-1) channel maps are modulated by wiggles or kinks that deviate from Keplerian kinematics particularly over the continuum ring, where deviations span 90$^\circ$ in azimuth, covering 5 km s$^{-1}$. The most pronounced wiggle resembles the imprint of an embedded massive planet of at least 5 M$_{\rm Jup}$ predicted from previous hydrodynamical simulations (Perez, Casassus, & Benitez-Llambay 2018). Such planet is expected to open a deep gap in both gas and dust density fields within a few orbital timescales, yet the kinematic wiggles lie near ridges in the continuum. The lesser strength of the wiggles in the $^{13}$CO and C$^{18}$O isotopologues show that the kinematic signature weakens at lower disk heights, and suggests qualitatively that it is due to vertical flows in the disk surface. Within the gap, the velocity field transitions from Keplerian to strongly non-Keplerian via a twist in position angle, suggesting the presence of another perturber and/or an inner warp. We also present VLT/SPHERE sparse aperture masking data which recovers scattered light emission from the gap's edges but shows no evidence for signal within the gap, discarding a stellar binary origin for its opening.

preprint2020arXiv

Near-Infrared Imaging of a Spiral in the CQ Tau Disk

We present $L^\prime$-band Keck/NIRC2 imaging and $H$-band Subaru/AO188+HiCIAO polarimetric observations of CQ Tau disk with a new spiral arm. Apart from the spiral feature our observations could not detect any companion candidates. We traced the spiral feature from the $r^2$-scaled HiCIAO polarimetric intensity image and the fitted result is used for forward modeling to reproduce the ADI-reduced NIRC2 image. We estimated the original surface brightness after throughput correction in $L^\prime$-band to be $\sim126$ mJy/arcsec$^2$ at most. We suggest that the grain temperature of the spiral may be heated up to $\sim$200 K in order to explain both of the $H$- and $L^{\prime}$-bands results. The $H$-band emission at the location of the spiral originates from the scattering from the disk surface while both scattering and thermal emission may contribute to the $L^{\prime}$-band emission. If the central star is only the light source of scattered light, the spiral emission at $L^\prime$-band should be thermal emission. If an inner disk also acts as the light source, the scattered light and the thermal emission may equally contribute to the $L^\prime$-band spiral structure.

preprint2020arXiv

The Gemini Planet Imager view of the HD 32297 debris disk

We present new $H$-band scattered light images of the HD 32297 edge-on debris disk obtained with the Gemini Planet Imager (GPI). The disk is detected in total and polarized intensity down to a projected angular separation of 0.15", or 20au. On the other hand, the large scale swept-back halo remains undetected, likely a consequence of its markedly blue color relative to the parent body belt. We analyze the curvature of the disk spine and estimate a radius of $\approx$100au for the parent body belt, smaller than past scattered light studies but consistent with thermal emission maps of the system. We employ three different flux-preserving post-processing methods to suppress the residual starlight and evaluate the surface brightness and polarization profile along the disk spine. Unlike past studies of the system, our high fidelity images reveal the disk to be highly symmetric and devoid of morphological and surface brightness perturbations. We find the dust scattering properties of the system to be consistent with those observed in other debris disks, with the exception of HR 4796. Finally, we find no direct evidence for the presence of a planetary-mass object in the system.

preprint2020arXiv

The Planetary Luminosity Problem: "Missing Planets" and the Observational Consequences of Episodic Accretion

The high occurrence rates of spiral arms and large central clearings in protoplanetary disks, if interpreted as signposts of giant planets, indicate that gas giants form commonly as companions to young stars ($<$ few Myr) at orbital separations of 10--300\,au. However, attempts to directly image this giant planet population as companions to more mature stars ($> 10$\, Myr) have yielded few successes. This discrepancy could be explained if most giant planets form "cold start," i.e., by radiating away much of their formation energy as they assemble their mass, rendering them faint enough to elude detection at later times. In that case, giant planets should be bright at early times, during their accretion phase, and yet forming planets are detected only rarely through direct imaging techniques. Here we explore the possibility that the low detection rate of accreting planets is the result of episodic accretion through a circumplanetary disk. We also explore the possibility that the companion orbiting the Herbig Ae star HD~142527 may be a giant planet undergoing such an accretion outburst.

preprint2019arXiv

Development Plans for the Atacama Large Millimeter/submillimeter Array (ALMA)

(abridged) The Atacama Large Millimeter/submillimeter Array (ALMA) was the top-ranked priority for a new ground-based facility in the 2000 Canadian Long Range Plan. Ten years later, at the time of LRP2010, ALMA construction was well underway, with first science observations anticipated for 2011. In the past 8 years, ALMA has proved itself to be a high-impact, high-demand observatory, with record numbers of proposals submitted to the annual calls and large numbers of highly cited scientific papers across fields from protoplanetary disks to high-redshift galaxies and quasars. The LRP2010 ALMA white paper laid out 8 specific metrics that could be used to judge the success of Canada's participation in ALMA. Among these metrics were publications (number; impact), collaborations (international; multi-wavelength), and student training. To call out one particular metric, Canadians are making excellent use of ALMA in training graduate students and postdocs: as of June 2018, 12 of 23 Canadian first-author papers were led by a graduate student, and a further 4 papers were led by postdocs. All 8 metrics argue for Canada's involvement in ALMA over the past decade to be judged a success. The successful achievement of these wide-ranging goals argues strongly for Canada's continuing participation in ALMA over the next decade and beyond. Looking forward, our community needs to: (1) maintain Canadian access to ALMA and our competitiveness in using ALMA; (2) preserve full Canadian funding for our share of ALMA operations; (3) identify components of ALMA development in which Canada can play a significant role, including stimulating expertise in submillimetre instrumentation to capitalize on future opportunities; and (4) keep Canadians fully trained and engaged in ALMA, as new capabilities become available, reaching the widest possible community of potential users.

preprint2019arXiv

High-Resolution Near-Infrared Polarimetry and Sub-Millimeter Imaging of FS Tau A: Possible Streamers in Misaligned Circumbinary Disk System

We analyzed the young (2.8-Myr-old) binary system FS Tau A using near-infrared (H-band) high-contrast polarimetry data from Subaru/HiCIAO and sub-millimeter CO (J=2-1) line emission data from ALMA. Both the near-infrared and sub-millimeter observations reveal several clear structures extending to $\sim$240 AU from the stars. Based on these observations at different wavelengths, we report the following discoveries. One arm-like structure detected in the near-infrared band initially extends from the south of the binary with a subsequent turn to the northeast, corresponding to two bar-like structures detected in ALMA observations with an LSRK velocity of 1.19-5.64 km/s. Another feature detected in the near-infrared band extends initially from the north of the binary, relating to an arm-like structure detected in ALMA observations with an LSRK velocity of 8.17-16.43 km/s. From their shapes and velocities, we suggest that these structures can mostly be explained by two streamers that connect the outer circumbinary disk and the central binary components. These discoveries will be helpful for understanding the evolution of streamers and circumstellar disks in young binary systems.

preprint2019arXiv

LRP2020: Signposts of planet formation in protoplanetary disks

Successful exoplanet surveys in the last decade have revealed that planets are ubiquitous throughout the Milky Way, and show a large diversity in mass, location and composition. At the same time, new facilities such as the Atacama Large Millimeter/submillimeter Array (ALMA) and optical/infrared facilities including Gemini/GPI have provided us with sharper images than ever before of protoplanetary disks around young stars, the birth cradles of planets. The high spatial resolution has revealed astonishing structures in disks, such as rings, gaps, asymmetries and spiral arms, and the enormous jump in sensitivity has provided the tools for both large, statistically relevant surveys and deep, sensitive molecular line studies. These observations have revolutionized our view of planet formation, disk formation and disk evolution, bringing model simulations and observations closer to the same level of detail, with many contributions from Canadian researchers on theoretical, observational and technological sides. The new results have inevitably led to a range of new questions, which require next generation instruments such as the Next Generation Very Large Array (ngVLA) and large scale optical infrared facilities. In this white paper we will discuss the current transformation in our understanding of planet formation and the next steps and challenges in connecting theory with exoplanet demographics and protoplanetary disk observations for Canadian research.

preprint2019arXiv

The detection of dust gap-ring structure in the outer region of the CR Cha protoplanetary disk

We observe the dust continuum at 225 GHz and CO isotopologue (12CO, 13CO, and C18O) J=2-1 emission lines toward the CR Cha protoplanetary disk using the Atacama Large Millimeter/Submillimeter Array (ALMA). The dust continuum image shows a dust gap-ring structure in the outer region of the dust disk. A faint dust ring is also detected around 120 au beyond the dust gap. The CO isotopologue lines indicate that the gas disk is more extended than the dust disk. The peak brightness temperature of the 13CO line shows a small bump around 130 au while 12CO and C18O lines do not show. We investigate two possible mechanisms for reproducing the observed dust gap-ring structure and a gas temperature bump. First, the observed gap structure can be opened by a Jupiter mass planet using the relation between the planet mass and the gap depth and width. Meanwhile, the radiative transfer calculations based on the observed dust surface density profile show that the observed dust ring could be formed by dust accumulation at the gas temperature bump, that is, the gas pressure bump produced beyond the outer edge of the dust disk.

preprint2019arXiv

The Formation of Stars -- From Filaments to Cores to Protostars and Protoplanetry Disks

Star formation involves the flow of gas and dust within molecular clouds into protostars and young stellar objects (YSOs) due to gravity. Along the way, these flows are shaped significantly by many other mechanisms, including pressure, turbulent motions, magnetic fields, stellar feedback, jets, and angular momentum. How all these mechanisms interact nonlinearly with each other on various length scales leads to the formation and evolution of substructures within clouds, including filaments, clumps, cores, disks, outflows, the protostars/YSOs themselves, and planets. In this white paper, prepared for the 2020 Long Range Plan panel which will recommend Canada's future directions for astronomy, we describe the observational and theoretical leadership in the star formation field that Canada's vibrant community has demonstrated over the past decade. Drawing from this extensive background, we identify five key questions that must be addressed for further progress to be made in understanding star formation in the next decade. Addressing these questions will improve our understanding of the dynamics of the dense gas and the role of the magnetic field in star formation, the optical properties of the dust used to trace mass and magnetic fields, the sources of variability in star-forming objects on short timescales, and the physical processes that specifically promote the clustering of stars. We further highlight key facilities in which Canada should become involved to continue making progress in this field. Single-dish facilities we recommend include LSST, trans-atmospheric far-infrared telescopes like BLAST-TNG and SPICA, and ground-based telescopes like JCMT, GBT, and CCAT-p. Interferometric facilities we recommend include ALMA, ngVLA, and SKA1.

preprint2016arXiv

Architecture design study and technology roadmap for the Planet Formation Imager (PFI)

The Planet Formation Imager (PFI) Project has formed a Technical Working Group (TWG) to explore possible facility architectures to meet the primary PFI science goal of imaging planet formation in situ in nearby star- forming regions. The goals of being sensitive to dust emission on solar system scales and resolving the Hill-sphere around forming giant planets can best be accomplished through sub-milliarcsecond imaging in the thermal infrared. Exploiting the 8-13 micron atmospheric window, a ground-based long-baseline interferometer with approximately 20 apertures including 10km baselines will have the necessary resolution to image structure down 0.1 milliarcseconds (0.014 AU) for T Tauri disks in Taurus. Even with large telescopes, this array will not have the sensitivity to directly track fringes in the mid-infrared for our prime targets and a fringe tracking system will be necessary in the near-infrared. While a heterodyne architecture using modern mid-IR laser comb technology remains a competitive option (especially for the intriguing 24 and 40μm atmospheric windows), the prioritization of 3-5μm observations of CO/H2O vibrotational levels by the PFI-Science Working Group (SWG) pushes the TWG to require vacuum pipe beam transport with potentially cooled optics. We present here a preliminary study of simulated L- and N-band PFI observations of a realistic 4-planet disk simulation, finding 21x2.5m PFI can easily detect the accreting protoplanets in both L and N-band but can see non-accreting planets only in L band. (abridged -- see PDF for full abstract)

preprint2016arXiv

Circumstellar Disks of the Most Vigorously Accreting Young Stars

Young stellar objects (YSOs) may not accumulate their mass steadily, as was previously thought, but in a series of violent events manifesting themselves as sharp stellar brightening. These events can be caused by fragmentation due to gravitational instabilities in massive gaseous disks surrounding young stars, followed by migration of dense gaseous clumps onto the star. We report our high angular resolution, coronagraphic near-infrared polarization imaging observations using the High Contrast Instrument for the Subaru Next Generation Adaptive Optics (HiCIAO) of the Subaru 8.2 m Telescope, towards four YSOs which are undergoing luminous accretion outbursts. The obtained infrared images have verified the presence of several hundred AUs scale arms and arcs surrounding these YSOs. In addition, our hydrodynamics simulations and radiative transfer models further demonstrate that these observed structures can indeed be explained by strong gravitational instabilities occurring at the beginning of the disk formation phase. The effect of those tempestuous episodes of disk evolution on star and planet formation remains to be understood.

preprint2016arXiv

How Spirals and Gaps Driven by Companions in Protoplanetary Disks Appear in Scattered Light at Arbitrary Viewing Angles

Direct imaging observations of protoplanetary disks at near-infrared (NIR) wavelengths have revealed structures of potentially planetary origin. Investigations of observational signatures from planet-induced features have so far focused on disks viewed face-on. Combining 3D hydrodynamics and radiative transfer simulations, we study how the appearance of the spiral arms and the gap produced in a disk by a companion varies with inclination and position angle in NIR scattered light. We compare the cases of a $3M_{\rm J}$ and a $0.1M_{\odot}$ companion, and make predictions suitable for testing with Gemini/GPI, VLT/NACO/SPHERE, and Subaru/HiCIAO/SCExAO. We find that the two trailing arms produced by an external perturber can have a variety of morphologies in inclined systems - they may appear as one trailing arm; two trailing arms on the same side of the disk; or two arms winding in opposite directions. The disk ring outside a planetary gap may also mimic spiral arms when viewed at high inclinations. We suggest potential explanations for the features observed in HH 30, HD 141569 A, AK Sco, HD 100546, and AB Aur. We emphasize that inclined views of companion-induced features cannot be converted into face-on views using simple and commonly practiced image deprojections.

preprint2016arXiv

Signatures of Gravitational Instability in Resolved Images of Protostellar Disks

Protostellar (class 0/I) disks, having masses comparable to those of their nascent host stars, and fed continuously from their natal infalling envelopes, are prone to gravitational instability (GI). Motivated by advances in near-infrared (NIR) adaptive optics imaging and mm-wave interferometry, we explore the observational signatures of GI in disks, using hydrodynamical and Monte Carlo radiative transfer simulations to synthesize NIR scattered light images and mm dust continuum maps. Spiral arms induced by GI, located at disk radii of hundreds of AUs, are local overdensities and have their photospheres displaced to higher altitudes above the disk midplane, arms therefore scatter more NIR light from their central stars than inter-arm regions, and are detectable at distances up to 1 kpc by Gemini/GPI, VLT/SPHERE, and Subaru/HiCIAO/SCExAO. By contrast, collapsed clumps formed by disk fragmentation have such strong local gravitational fields that their scattering photospheres are at lower altitudes, such fragments appear fainter than their surroundings in NIR scattered light. Spiral arms and streamers recently imaged in four FU Ori systems at NIR wavelengths resemble GI-induced structures and support the interpretation that FUors are gravitationally unstable protostellar disks. At mm wavelengths, both spirals and clumps appear brighter in thermal emission than the ambient disk and can be detected by ALMA at distances up to 0.4 kpc with one-hour integration times at ~0.1 arcsec resolution. Collapsed fragments having masses >=1 MJ can be detected by ALMA within ~10 minutes.

preprint2016arXiv

Stability and Occurrence Rate Constraints on the Planetary Sculpting Hypothesis for "Transitional" Disks

Transitional disks, protoplanetary disks with deep and wide central gaps, may be the result of planetary sculpting. By comparing numerical planet-opening-gap models with observed gaps, we find systems of 3-6 giant planets are needed in order to open gaps with the observed depths and widths. We explore the dynamical stability of such multi-planet systems using N-body simulations that incorporate prescriptions for gas effects. We find they can be stable over a typical disk lifetime, with the help of eccentricity damping from the residual gap gas that facilitates planets locking into mean motion resonances. However, in order to account for the occurrence rate of transitional disks, the planet sculpting scenario demands gap-opening-friendly disk conditions, in particular, a disk viscosity $α\lesssim0.001$. In addition, the demography of giant planets at $\sim 3-30$ AU separations, poorly constrained by current data, has to largely follow occurrence rates extrapolated outward from radial velocity surveys, not the lower occurrence rates extrapolated inward from direct imaging surveys. Even with the most optimistic occurrence rates, transitional disks cannot be a common phase that most gas disks experience at the end of their life, as popularly assumed, simply because there are not enough planets to open these gaps. Finally, as consequences of demanding almost all giant planets at large separations participate in transitional disk sculpting, the majority of such planets must form early and end up in a chain of mean motion resonances at the end of disk lifetime.

preprint2016arXiv

The PDS 66 Circumstellar Disk as seen in Polarized Light with the Gemini Planet Imager

We present H and K band imaging polarimetry for the PDS 66 circumstellar disk obtained during the commissioning of the Gemini Planet Imager (GPI). Polarization images reveal a clear detection of the disk in to the 0.12'' inner working angle (IWA) in H band, almost 3 times as close to the star as the previous HST observations with NICMOS and STIS (0.35'' effective IWA). The centro-symmetric polarization vectors confirm that the bright inner disk detection is due to circumstellar scattered light. A more diffuse disk extends to a bright outer ring centered at 80 AU. We discuss several physical mechanisms capable of producing the observed ring + gap structure. GPI data confirm enhanced scattering on the East side of the disk which is inferred to be nearer to us. We also detect a lateral asymmetry in the South possibly due to shadowing from material within the inner working angle. This likely corresponds to a temporally variable azimuthal asymmetry observed in HST/STIS coronagraphic imaging.

preprint2015arXiv

Absence of Significant Cool Disks in Young Stellar Objects Exhibiting Repetitive Optical Outbursts

We report Submillimeter Array (SMA) 1.3 mm high angular resolution observations towards the four EXor type outbursting young stellar objects (YSOs) VY Tau, V1118 Ori, V1143 Ori, and NY Ori. The data mostly show low dust masses $M_{dust}$ in the associated circumstellar disks. Among the sources, NY Ori possesses a relatively massive disk with $M_{dust} \sim 9 \times 10^{-4}$ $M_{\odot}$. V1118 Ori has a marginal detection equivalent to $M_{dust} \sim 6 \times 10^{-5}$ $M_{\odot}$. V1143 Ori has a non-detection also equivalent to $M_{dust} < 6 \times 10^{-5}$ $M_{\odot}$. For the nearest source VY Tau, we get a surprising non-detection which provides a stringent upper limit $M_{dust} < 6 \times 10^{-6}$ $M_{\odot}$. We interpret our findings as suggesting that the gas and dust reservoirs that feed the short duration, repetitive optical outbursts seen in some EXors may be limited to the small scale, innermost region of their circumstellar disks. This hot dust may have escaped our detection limits. Follow-up, more sensitive millimeter observations are needed to improve our understanding of the triggering mechanisms of EXor type outbursts.

preprint2015arXiv

An M Dwarf Companion and Its Induced Spiral Arms in the HD 100453 Protoplanetary Disk

Recent VLT/SPHERE near-infrared imaging observations revealed two spiral arms with a near m=2 rotational symmetry in the protoplanetary disk around the ~1.7 solar mass Herbig star HD 100453. A ~0.3 solar mass M dwarf companion, HD 100453 B, was also identified at a projected separation of 120 AU from the primary. In this Letter, we carry out hydrodynamic and radiative transfer simulations to examine the scattered light morphology of the HD 100453 disk as perturbed by the companion on a circular and coplanar orbit. We find that the companion truncates the disk at ~45 AU in scattered light images, and excites two spiral arms in the remaining (circumprimary) disk with a near m=2 rotational symmetry. Both the truncated disk size and the morphology of the spirals are in excellent agreement with the SPHERE observations at Y , J, H, and K1-bands, suggesting that the M dwarf companion is indeed responsible for the observed double-spiral-arm pattern. Our model suggests that the disk is close to face on (inclination angle ~5 degree), and that the entire disk-companion system rotates counterclockwise on the sky. The HD 100453 observations, along with our modeling work, demonstrate that double spiral arm patterns in near-infrared scattered light images can be generically produced by companions, and support future observations to identify the companions responsible for the arms observed in the MWC 758 and SAO 206462 systems.

preprint2015arXiv

Direct imaging of an asymmetric debris disk in the HD 106906 planetary system

We present the first scattered light detections of the HD 106906 debris disk using Gemini/GPI in the infrared and HST/ACS in the optical. HD 106906 is a 13 Myr old F5V star in the Sco-Cen association, with a previously detected planet-mass candidate HD 106906b projected 650 AU from the host star. Our observations reveal a near edge-on debris disk that has a central cleared region with radius $\sim$50 AU, and an outer extent $>$500 AU. The HST data show the outer regions are highly asymmetric, resembling the ''needle'' morphology seen for the HD 15115 debris disk. The planet candidate is oriented $\sim$21$°$ away from the position angle of the primary's debris disk, strongly suggesting non-coplanarity with the system. We hypothesize that HD 106906b could be dynamically involved in the perturbation of the primary's disk, and investigate whether or not there is evidence for a circumplanetary dust disk or cloud that is either primordial or captured from the primary. We show that both the existing optical properties and near-infrared colors of HD 106906b are weakly consistent with this possibility, motivating future work to test for the observational signatures of dust surrounding the planet.

preprint2015arXiv

Inferring Planet Mass from Spiral Structures in Protoplanetary Disks

Recent observations of protoplanetary disk have reported spiral structures that are potential signatures of embedded planets, and modeling efforts have shown that a single planet can excite multiple spiral arms, in contrast to conventional disk-planet interaction theory. Using two and three-dimensional hydrodynamics simulations to perform a systematic parameter survey, we confirm the existence of multiple spiral arms in disks with a single planet, and discover a scaling relation between the azimuthal separation of the primary and secondary arm, $ϕ_{\rm sep}$, and the planet-to-star mass ratio $q$: $ϕ_{\rm sep} = 102^{\circ} (q/0.001)^{0.2}$ for companions between Neptune mass and 16 Jupiter masses around a 1 solar mass star, and $ϕ_{\rm sep} = 180^{\circ}$ for brown dwarf mass companions. This relation is independent of the disk's temperature, and can be used to infer a planet's mass to within an accuracy of about 30% given only the morphology of a face-on disk. Combining hydrodynamics and Monte-Carlo radiative transfer calculations, we verify that our numerical measurements of $ϕ_{\rm sep}$ are accurate representations of what would be measured in near-infrared scattered light images, such as those expected to be taken by Gemini/GPI, VLT/SPHERE, or Subaru/SCExAO in the future. Finally, we are able to infer, using our scaling relation, that the planet responsible for the spiral structure in SAO 206462 has a mass of about 6 Jupiter masses.

preprint2015arXiv

Observational Signatures of Planets in Protoplanetary Disks I: Gaps Opened by Single and Multiple Young Planets in Disks

It has been suggested that the gaps and cavities recently discovered in transitional disks are opened by planets. To explore this scenario, we combine two-dimensional two fluid (gas + particle) hydrodynamical calculations with three-dimensional Monte Carlo Radiative Transfer simulations, and study the observational signatures of gaps opened by one or several planets, making qualitative comparisons with observations. We find that a single planet as small as 0.2 MJ can produce a deep gap at millimeter (mm) wavelengths and almost no features at near-infrared (NIR) wavelengths, while multiple planets can open up a few *10 AU wide common gap at both wavelengths. Both the contrast ratio of the gaps and the wavelength dependence of the gap sizes are broadly consistent with data. We also confirm previous results that NIR gap sizes may be smaller than mm gap sizes due to dust-gas coupling and radiative transfer effects. When viewed at a moderate inclination angle, a physically circular on-centered gap could appear to be off-centered from the star due to shadowing. Planet-induced spiral arms are more apparent at NIR than at mm wavelengths. Overall, our results suggest that the planet-opening-gap scenario is a promising way to explain the origin of the transitional disks. Finally, inspired by the recent ALMA release of the image of the HL Tau disk, we show that multiple narrow gaps, well separated by bright rings, can be opened by 0.2 MJ planets soon after their formation in a relatively massive disk.

preprint2015arXiv

Observational Signatures of Planets in Protoplanetary Disks II: Spiral Arms Observed in Scattered Light Imaging Can be Induced by Planets

Using 3D global hydro simulations coupled with radiative transfer calculations, we study the appearance of density waves induced by giant planets in direct imaging observations at near infrared wavelengths. We find that a 6 MJ planet in a typical disk around a 1 M_sun star can produce prominent and detectable spiral arms both interior and exterior to its orbit. The inner arms have (1) two well separated arms in roughly m=2 symmetry, (2) exhibit ~10-15 degrees pitch angles, (3) ~180-270 degrees extension in the azimuthal direction, and (4) ~150% surface brightness enhancement, all broadly consistent with observed spiral arms in the SAO 206462 and MWC 758 systems. The outer arms cannot explain observations as they are too tightly wound given typical disk scale height. We confirm previous results that the outer density waves excited by a 1 MJ planet exhibit low contrast in the IR and are practically not detectable. We also find that 3D effects of the waves are important. Compared to isothermal models, density waves in adiabatic disks exhibit weaker contrast in surface density but stronger contrast in scattered light images, due to a more pronounced vertical structure in the former caused by shock heating. To drive observed pairs of arms with an external companion on a circular orbit, a massive planet, possibly a brown dwarf, is needed at around [r~0.7", PA~10 degrees] (position angle PA from north to east) in SAO 206462 and [r~0.6", PA~10 degrees] in MWC 758. Their existence may be confirmed by direct imaging planet searches.

preprint2015arXiv

Peering into the Giant Planet Forming Region of the TW Hydrae Disk with the Gemini Planet Imager

We present Gemini Planet Imager (GPI) adaptive optics near-infrared images of the giant planet-forming regions of the protoplanetary disk orbiting the nearby (D = 54 pc), pre-main sequence (classical T Tauri) star TW Hydrae. The GPI images, which were obtained in coronagraphic/polarimetric mode, exploit starlight scattered off small dust grains to elucidate the surface density structure of the TW Hya disk from 80 AU to within 10 AU of the star at 1.5 AU resolution. The GPI polarized intensity images unambiguously confirm the presence of a gap in the radial surface brightness distribution of the inner disk. The gap is centered near 23 AU, with a width of 5 AU and a depth of 50%. In the context of recent simulations of giant planet formation in gaseous, dusty disks orbiting pre-main sequence stars, these results indicate that at least one young planet with a mass 0.2 M_J could be present in the TW Hya disk at an orbital semi-major axis similar to that of Uranus. If this (proto)planet is actively accreting gas from the disk, it may be readily detectable by GPI or a similarly sensitive, high-resolution infrared imaging system.

preprint2015arXiv

Shallow Cavities in Multiple-Planet Systems

Large cavities are often observed in protoplanetary disks, which might suggest the presence of planets opening gaps in the disk. Multiple planets are necessary to produce a wide cavity in the gas. However, multiple planets may also be a burden to the carving out of very deep gaps. When additional planets are added to the system, the time-dependent perturbations from these additional satellites can stir up gas in the gap, suppressing cavity opening. In this study, we perform two-dimensional numerical hydro calculations of gap opening for single and multiple planets, showing the effect that additional planets have on the gap depths. We show that multiple planets produce much shallower cavities than single planets, so that more massive planets are needed in the multiple-planet case to produce an equivalent gap depth as in the single-planet case. To deplete a gap by a factor of 100 for the parameters chosen in this study, one only requires $M_p \approx 3.5M_J$ in the single-planet case, but much more massive planets, $M_p \approx 7M_J$ are required in the multiple-planet case. This requirement of high-mass planets implies that such planets may be detectable in the next generation of direct imaging projects, in gaps whose depths are constrained to be sufficiently deep by ALMA.

preprint2015arXiv

Spiral Arms in Gravitationally Unstable Protoplanetary Disks as Imaged in Scattered Light

Combining 3D smoothed-particle hydrodynamics and Monte Carlo radiative transfer calculations, we examine the morphology of spiral density waves induced by gravitational instability (GI) in protoplanetary disks, as they would appear in direct images at near-infrared (NIR) wavelengths. We find that systems with disk-to-star-mass ratios q=M_disk/M_star that are ~0.25 or more may produce prominent spiral arms in NIR imaging, remarkably resembling features observed in the MWC 758 and SAO 206462 systems. The contrast of GI-induced arms at NIR wavelengths can reach a factor of ~3, and their pitch angles are about 10-15 degree. The dominant azimuthal wavenumber of GI-induced spiral arms roughly obeys m~1/q in the range 2<~1/q<~8. In particular, a massive disk with q~0.5 can exhibit grand-design m=2 spirals. GI-induced arms are in approximate corotation with the local disk, and may therefore trap dust particles by pressure drag. Although GI can produce NIR spiral arms with morphologies, contrasts, and pitch angles similar to those reported in recent observations, it also makes other demands that may or may not be satisfied in any given system. A GI origin requires that the spirals be relatively compact, on scales <~100 AU; that the disk be massive, q>~0.25; and that the accretion rate Mdot be high, on the order of 1e-6 solar mass per year.

preprint2015arXiv

The Effects of Self-Shadowing by a Puffed up Inner Rim in Scattered Light Images of Protoplanetary Disks

We explore whether protoplanetary disks with self-shadowing from puffed up inner rims exhibit observable features in scattered light images. We use both self-consistent hydrostatic equilibrium calculations and parameterized models to produce the vertically puffed up inner rims. We find that, in general, the transition between the shadowed and flared regions occurs in a smooth manner over a broad radius range, and no sudden jump exists at the outer edge of the shadow in either the disk temperature or density structures. As a result, a puffed up rim cannot create sharp ring/arc/spiral-arm-like features in the outer disk as have been detected in recent direct NIR imaging of disks. On the other hand, if the puffed up rim has a sharp edge in the vertical direction, the shadowing effect can produce a distinct 3-stage broken power law in the radial intensity profile of the scattered light, with 2 steep surface brightness radial profiles in the inner and outer disk joined by a shallow transition region around the shadow edge. These types of scattered light profiles may have already been observed, such as in the recent Subaru direct imaging of the TW Hydrae system.

preprint2015arXiv

The Structure of Spiral Shocks Excited by Planetary-mass Companions

Direct imaging observations have revealed spiral structures in protoplanetary disks. Previous studies have suggested that planet-induced spiral arms cannot explain some of these spiral patterns, due to the large pitch angle and high contrast of the spiral arms in observations. We have carried out three dimensional (3-D) hydrodynamical simulations to study spiral wakes/shocks excited by young planets. We find that, in contrast with linear theory, the pitch angle of spiral arms does depend on the planet mass, which can be explained by the non-linear density wave theory. A secondary (or even a tertiary) spiral arm, especially for inner arms, is also excited by a massive planet. With a more massive planet in the disk, the excited spiral arms have larger pitch angle and the separation between the primary and secondary arms in the azimuthal direction is also larger. We also find that although the arms in the outer disk do not exhibit much vertical motion, the inner arms have significant vertical motion, which boosts the density perturbation at the disk atmosphere. Combining hydrodynamical models with Monte-Carlo radiative transfer calculations, we find that the inner spiral arms are considerably more prominent in synthetic near-IR images using full 3-D hydrodynamical models than images based on 2-D models assuming vertical hydrostatic equilibrium, indicating the need to model observations with full 3-D hydrodynamics. Overall, companion-induced spiral arms not only pinpoint the companion's position but also provide three independent ways (pitch angle, separation between two arms, and contrast of arms) to constrain the companion's mass.

preprint2014arXiv

Surface Geometry of Protoplanetary Disks Inferred From Near-Infrared Imaging Polarimetry

We present a new method of analysis for determining the surface geometry of five protoplanetary disks observed with near-infrared imaging polarimetry using Subaru-HiCIAO. Using as inputs the observed distribution of polarized intensity (PI), disk inclination, assumed properties for dust scattering, and other reasonable approximations, we calculate a differential equation to derive the surface geometry. This equation is numerically integrated along the distance from the star at a given position angle. We show that, using these approximations, the local maxima in the PI distribution of spiral arms (SAO 206462, MWC 758) and rings (2MASS J16042165-2130284, PDS 70) are associated with local concave-up structures on the disk surface. We also show that the observed presence of an inner gap in scattered light still allows the possibility of a disk surface that is parallel to the light path from the star, or a disk that is shadowed by structures in the inner radii. Our analysis for rings does not show the presence of a vertical inner wall as often assumed in studies of disks with an inner gap. Finally, we summarize the implications of spiral and ring structures as potential signatures of ongoing planet formation.

preprint2014arXiv

The Science Case for the Planet Formation Imager (PFI)

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

preprint2013arXiv

Mapping H-band Scattered Light Emission in the Mysterious SR21 Transitional Disk

We present the first near infrared (NIR) spatially resolved images of the circumstellar transitional disk around SR21. These images were obtained with the Subaru HiCIAO camera, adaptive optics and the polarized differential imaging (PDI) technique. We resolve the disk in scattered light at H-band for stellocentric 0.1"<r<0.6" (12<r<75AU). We compare our results with previously published spatially-resolved 880 micron continuum Submillimeter Array (SMA) images that show an inner r<36AU cavity in SR21. Radiative transfer models reveal that the large disk depletion factor invoked to explain SR21's sub-mm cavity cannot be "universal" for all grain sizes. Even significantly more moderate depletions (delta=0.1, 0.01 relative to an undepleted disk) than those that reproduce the sub-mm cavity (delta~10^-6) are inconsistent with our H-band images when they are assumed to carry over to small grains, suggesting that surface grains scattering in the NIR either survive or are generated by whatever mechanism is clearing the disk midplane. In fact, the radial polarized intensity profile of our H-band observations is smooth and steeply inwardly-increasing (r^-3), with no evidence of a break at the 36AU sub-mm cavity wall. We hypothesize that this profile is dominated by an optically thin disk envelope or atmosphere component. We also discuss the compatibility of our data with the previously postulated existence of a sub-stellar companion to SR21 at r~10-20AU, and find that we can neither exclude nor verify this scenario. This study demonstrates the power of multiwavelength imaging of transitional disks to inform modeling efforts, including the debate over precisely what physical mechanism is responsible for clearing these disks of their large midplane grains.

preprint2012arXiv

Dust Filtration by Planet-Induced Gap Edges: Implications for Transitional Disks

By carrying out two-dimensional two-fluid global simulations, we have studied the response of dust to gap formation by a single planet in the gaseous component of a protoplanetary disk - the so-called "dust filtration" mechanism. We have found that a gap opened by a giant planet at 20 AU in a α=0.01, \dot{M}=10^{-8} Msun/yr disk can effectively stop dust particles larger than 0.1 mm drifting inwards, leaving a sub-millimeter dust cavity/hole. However, smaller particles are difficult to filter by a planet-induced gap due to 1) dust diffusion, and 2) a high gas accretion velocity at the gap edge. An analytic model is also derived to understand what size particles can be filtered by the gap edge. Finally, with our updated understanding of dust filtration, we have computed Monte-Carlo radiative transfer models with variable dust size distributions to generate the spectral energy distributions (SEDs) of disks with gaps. By comparing with transitional disk observations (e.g. GM Aur), we have found that dust filtration alone has difficulties to deplete small particles sufficiently to explain the near-IR deficit of transitional disks, except under some extreme circumstances. The scenario of gap opening by multiple planets studied previously suffers the same difficulty. One possible solution is by invoking both dust filtration and dust growth in the inner disk. In this scenario, a planet induced gap filters large dust particles in the disk, and the remaining small dust particles passing to the inner disk can grow efficiently without replenishment from fragmentation of large grains. Predictions for ALMA have also been made based on all these scenarios. We conclude that dust filtration with planet(s) in the disk is a promising mechanism to explain submm observations of transitional disks but it may need to be combined with other processes (e.g. dust growth) to explain the near-IR deficit.

preprint2012arXiv

Polarimetric Imaging of Large Cavity Structures in the Pre-transitional Protoplanetary Disk around PDS 70: Observations of the disk

We present high resolution H-band polarized intensity (PI; FWHM = 0."1: 14 AU) and L'-band imaging data (FWHM = 0."11: 15 AU) of the circumstellar disk around the weak-lined T Tauri star PDS 70 in Centaurus at a radial distance of 28 AU (0."2) up to 210 AU (1."5). In both images, a giant inner gap is clearly resolved for the first time, and the radius of the gap is ~70 AU. Our data show that the geometric center of the disk shifts by ~6 AU toward the minor axis. We confirm that the brown dwarf companion candidate to the north of PDS 70 is a background star based on its proper motion. As a result of SED fitting by Monte Carlo radiative transfer modeling, we infer the existence of an optically thick inner disk at a few AU. Combining our observations and modeling, we classify the disk of PDS 70 as a pre-transitional disk. Furthermore, based on the analysis of L'-band imaging data, we put an upper limit mass of companions at ~30 to ~50MJ within the gap. Taking account of the presence of the large and sharp gap, we suggest that the gap could be formed by dynamical interactions of sub-stellar companions or multiple unseen giant planets in the gap.

preprint2012arXiv

The Structure of Pre-transitional Protoplanetary Disks I: Radiative Transfer Modeling of the Disk+Cavity in the PDS 70 system

Through detailed radiative transfer modeling, we present a disk+cavity model to simultaneously explain both the SED and Subaru H-band polarized light imaging for the pre-transitional protoplanetary disk PDS 70. Particularly, we are able to match not only the radial dependence, but also the absolute scale, of the surface brightness of the scattered light. Our disk model has a cavity 65 AU in radius, which is heavily depleted of sub-micron-sized dust grains, and a small residual inner disk which produces a weak but still optically thick NIR excess in the SED. To explain the contrast of the cavity edge in the Subaru image, a factor of ~1000 depletion for the sub-micron-sized dust inside the cavity is required. The total dust mass of the disk may be on the order of 1e-4 M_sun, only weakly constrained due to the lack of long wavelength observations and the uncertainties in the dust model. The scale height of the sub-micron-sized dust is ~6 AU at the cavity edge, and the cavity wall is optically thick in the vertical direction at H-band. PDS 70 is not a member of the class of (pre-)transitional disks identified by Dong et al. (2012), whose members only show evidence of the cavity in the millimeter-sized dust but not the sub-micron-sized dust in resolved images. The two classes of (pre-)transitional disks may form through different mechanisms, or they may just be at different evolution stages in the disk clearing process.

preprint2012arXiv

X-ray Properties of Intermediate-mass Black Holes in Active Galaxies. III. Spectral Energy Distribution and Possible Evidence for Intrinsically X-ray-weak AGNs

We present a systematic X-ray study, the third in a series, of 49 active galactic nuclei with intermediate-mass black holes (IMBH; ~10^5-10^6 M_sun) using Chandra observations. We detect 42 out of 49 targets with a 0.5-2 keV X-ray luminosity 10^41-10^43 erg/s. We perform spectral fitting for the 10 objects with enough counts (>200), and they are all well fit by a simple power-law model modified by Galactic absorption, with no sign of significant intrinsic absorption. While we cannot fit the X-ray spectral slope directly for the rest of the sample, we estimate it from the hardness ratio and find a range of photon indices consistent with those seen in more luminous and massive objects. The X-ray-to-optical spectral slope (alphaox) of our IMBH sample is systematically flatter than in active galaxies with more massive black holes, consistent with the well-known correlation between alphaox and UV luminosity. Thanks to the wide dynamic range of our sample, we find evidence that alphaox increases with decreasing M_BH as expected from accretion disk models, where the UV emission systematically decreases as M_BH decreases and the disk temperature increases. We also find a long tail toward low alphaox values. While some of these sources may be obscured, given the high L_bol/L_Eddington values in the sample, we argue that some may be intrinsically X-ray-weak, perhaps owing to a rare state that radiates very little coronal emission.

preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks I: Linear Regime

Density waves excited by planets embedded in protoplanetary disks play a central role in planetary migration and gap opening processes. We carry out 2D shearing sheet simulations to study the linear regime of wave evolution with the grid-based code Athena, and provide detailed comparisons with the theoretical predictions. Low mass planets (down to ~0.03 Earth mass at 1 AU) and high spatial resolution (256 grid points per scale height) are chosen to mitigate the effects of wave nonlinearity. To complement the existing numerical studies, we focus on the primary physical variables such as the spatial profile of the wave, torque density, and the angular momentum flux carried by the wave, instead of secondary quantities such as the planetary migration rate. Our results show percent level agreement with theory in both physical and Fourier space. New phenomena such as the change of the toque density sign far from the planet are discovered and discussed. Also, we explore the effect of the numerical algorithms, and find that a high order of accuracy, high resolution, and an accurate planetary potential are crucial to achieve good agreement with the theory. We find that the use of a too large time-step without properly resolving the dynamical time scale around the planet produces incorrect results, and may lead to spurious gap opening. Global simulations of planet migration and gap opening violating this requirement may be affected by spurious effects resulting in e.g. the incorrect planetary migration rate and gap opening mass.

preprint2011arXiv

Density Waves Excited by Low-Mass Planets in Protoplanetary Disks II: High-Resolution Simulations of the Nonlinear Regime

We investigate numerically the propagation of density waves excited by a low-mass planet in a protoplanetary disk in the nonlinear regime, using 2D local shearing box simulations with the grid-based code Athena at high spatial resolution (256 grid points per scale height h). The nonlinear evolution results in the wave steepening into a shock, causing damping and angular momentum transfer to the disk. On long timescales this leads to spatial redistribution of the disk density, causing migration feedback and potentially resulting in gap opening. Previous numerical studies concentrated on exploring these secondary phenomena as probes of the nonlinear wave evolution. Here we focus on exploring the evolution of the basic wave properties, such as its density profile evolution, shock formation, post-shock wave behavior, and provide comparison with analytical theory. The generation of potential vorticity at the shock is computed analytically and is subsequently verified by simulations and used to pinpoint the shock location. We confirm the theoretical relation between the shocking length and the planet mass (including the effect of the equation of state), and the post-shock decay of the angular momentum flux carried by the wave. The post-shock evolution of the wave profile is explored, and we quantitatively confirm its convergence to the theoretically expected N-wave shape. The accuracy of various numerical algorithms used to compute the nonlinear wave evolution is also investigated: we find that higher order spatial reconstruction and high resolution are crucial for capturing the shock formation correctly.

preprint2011arXiv

Investigation of the errors in SDSS proper-motion measurements using samples of quasars

We investigate in detail the probability distribution function (pdf) of the proper-motion measurement errors in the SDSS+USNO-B proper-motion catalog of \citet{mun04} using clean quasar samples. The pdf of the errors is well-represented by a Gaussian core with extended wings, plus a very small fraction ($<0.1%$) of "outliers". We find while formally the pdf could be well-fit by a five-parameter fitting function, for many purposes it is also adequately to represent the pdf with a one-parameter approximation to this function. We apply this pdf to the calculation of the confidence intervals on the true proper motion for a SDSS+USNO-B proper motion measurement, and discuss several scientific applications of the SDSS proper motion catalogue. Our results have various applications in studies of the galactic structure and stellar kinematics. Specifically, they are crucial for searching hyper-velocity stars in the Galaxy.

preprint2010arXiv

A Systematic Search for X-ray Cavities in the Hot Gas of Galaxy Groups

We have performed a systematic search for X-ray cavities in the hot gas of 51 galaxy groups with Chandra archival data. The cavities are identified based on two methods: subtracting an elliptical beta model fitted to the X-ray surface brightness, and performing unsharp masking. 13 groups in the sample 25% are identified as clearly containing cavities, with another 13 systems showing tentative evidence for such structures. We find tight correlations between the radial and tangential radii of the cavities, and between their size and projected distance from the group center, in quantitative agreement with the case for more massive clusters. This suggests that similar physical processes are responsible for cavity evolution and disruption in systems covering a large range in total mass. We see no clear association between the detection of cavities and the current 1.4 GHz radio luminosity of the central brightest group galaxy, but there is a clear tendency for systems with a cool core to be more likely to harbor detectable cavities. To test the efficiency of the adopted cavity detection procedures, we employ a set of mock images designed to mimic typical Chandra data of our sample, and find that the model-fitting approach is generally more reliable than unsharp masking for recovering cavity properties. Finally, we find that the detectability of cavities is strongly influenced by a few factors, particularly the signal-to-noise ratio of the data, and that the real fraction of X-ray groups with prominent cavities could be substantially larger than the 25--50% suggested by our analysis.

preprint2010arXiv

Dusty Disks around White Dwarfs I: Origin of Debris Disks

A significant fraction of the mature FGK stars have cool dusty disks at least an orders of magnitudes brighter than the solar system's outer zodiacal light. Since such dusts must be continually replenished, they are generally assumed to be the collisional fragments of residual planetesimals analogous to the Kuiper Belt objects. At least 10% of solar type stars also bear gas giant planets. The fraction of stars with known gas giants or detectable debris disks (or both) appears to increase with the stellar mass. Here, we examine the dynamical evolution of systems of long-period gas giant planets and residual planetesimals as their host stars evolve off the main sequence, lose mass, and form planetary nebula around remnant white dwarf cores. The orbits of distant gas giant planets and super-km-size planetesimals expand adiabatically. During the most intense AGB mass loss phase, sub-meter-size particles migrate toward their host stars due to the strong hydrodynamical drag by the intense stellar wind. Along their migration paths, gas giant planets capture and sweep up sub-km-size planetesimals onto their mean-motion resonances. These planetesimals also acquire modest eccentricities which are determined by the mass of the perturbing planets, the rate and speed of stellar mass loss. The swept-up planetesimals undergo disruptive collisions which lead to the production of grains with an extended size range. The radiation drag on these particles is ineffective against the planets' resonant barrier and they form 30-to-150-AU-sizes rings which can effective reprocess the stellar irradiation in the form of FIR continuum. We identify the recently discovered dust ring around the white dwarf WD 2226-210 at the center of the Helix nebula as a prototype of such disks and suggest such rings may be common.

preprint2010arXiv

H-alpha and Free-Free Emission from the WIM

Recent observations have found the ratio of H-alpha to free-free radio continuum to be surprisingly high in the diffuse ionized ISM (the so-called WIM), corresponding to an electron temperature of only ~3000K. Such low temperatures were unexpected in gas that was presumed to be photoionized. We consider a 3-component model for the observed diffuse emission, consisting of a mix of (1) photoionized gas, (2) gas that is recombining and cooling, and (3) cool H I gas. This model can successfully reproduce the observed intensities of free-free continuum, H-alpha, and collisionally-excited lines such as NII 6583. To reproduce the low observed value of free-free to H-alpha, the PAH abundance in the photoionized regions must be lowered by a factor ~3, and ~20% of the diffuse H-alpha must be reflected from dust grains, as suggested by Wood & Reynolds (1999).

preprint2009arXiv

Buoyant Bubbles in Intracluster Gas: Effects of Magnetic Fields and Anisotropic Viscosity

Recent observations by Chandra and XMM-Newton indicate there are complex structures at the cores of galaxy clusters, such as cavities and filaments. One plausible model for the formation of such structures is the interaction of radio jets with the intracluster medium (ICM). To investigate this idea, we use three-dimensional magnetohydrodynamic simulations including anisotropic (Braginskii) viscosity to study the effect of magnetic fields on the evolution and morphology of buoyant bubbles in the ICM. We investigate a range of different initial magnetic field geometries and strengths, and study the resulting x-ray surface brightness distribution for comparison to observed clusters. Magnetic tension forces and viscous transport along field lines tend to suppress instabilities parallel, but not perpendicular, to field lines. Thus, the evolution of the bubble depends strongly on the initial field geometry. We find toroidal field loops initially confined to the interior of the bubble are best able reproduce the observed cavity structures.