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.

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 $\mu$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 \deg$). 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.

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.

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.

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.

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.

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.

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.

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.

The SR24 multi-star system hosts both circumprimary and circumsecondary disks, which are strongly misaligned from each other. The circumsecondary disk is circumbinary in nature. Interestingly, both disks are interacting, and they possibly rotate in opposite directions. To investigate the nature of this unique twin disk system, we present 0.''1 resolution near-infrared polarized intensity images of the circumstellar structures around SR24, obtained with HiCIAO mounted on the Subaru 8.2 m telescope. Both the circumprimary disk and the circumsecondary disk are resolved and have elongated features. While the position angle of the major axis and radius of the NIR polarization disk around SR24S are 55$^{\circ}$ and 137 au, respectively, those around SR24N are 110$^{\circ}$ and 34 au, respectively. With regard to overall morphology, the circumprimary disk around SR24S shows strong asymmetry, whereas the circumsecondary disk around SR24N shows relatively strong symmetry. Our NIR observations confirm the previous claim that the circumprimary and circumsecondary disks are misaligned from each other. Both the circumprimary and circumsecondary disks show similar structures in $^{12}$CO observations in terms of its size and elongation direction. This consistency is because both NIR and $^{12}$CO are tracing surface layers of the flared disks. As the radius of the polarization disk around SR24N is roughly consistent with the size of the outer Roche lobe, it is natural to interpret the polarization disk around SR24N as a circumbinary disk surrounding the SR24Nb-Nc system.

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.

We present the first polarimetric detection of the inner disk component around the pre-main sequence B9.5 star HD 141569A. Gemini Planet Imager H-band (1.65 micron) polarimetric differential imaging reveals the highest signal-to-noise ratio detection of this ring yet attained and traces structure inwards to 0.25" (28 AU at a distance of 111 pc). The radial polarized intensity image shows the east side of the disk, peaking in intensity at 0.40" (44 AU) and extending out to 0.9" (100 AU). There is a spiral arm-like enhancement to the south, reminiscent of the known spiral structures on the outer rings of the disk. The location of the spiral arm is coincident with 12CO J=3-2 emission detected by ALMA, and hints at a dynamically active inner circumstellar region. Our observations also show a portion of the middle dusty ring at ~220 AU known from previous observations of this system. We fit the polarized H-band emission with a continuum radiative transfer Mie model. Our best-fit model favors an optically thin disk with a minimum dust grain size close to the blow-out size for this system: evidence of on-going dust production in the inner reaches of the disk. The thermal emission from this model accounts for virtually all of the far-infrared and millimeter flux from the entire HD 141569A disk, in agreement with the lack of ALMA continuum and CO emission beyond ~100 AU. A remaining 8-30 micron thermal excess a factor of ~2 above our model argues for a yet-unresolved warm innermost 5-15 AU component of the disk.

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.

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

(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.

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.

We present ALMA observations of 12CO, 13CO, and C18O J=2--1 lines and the 230 GHz continuum for the FU Ori-type object (FUor) V900 Mon (d~1.5 kpc), for which the accretion burst was triggered between 1953 and 2009. We identified CO emission associated with a molecular bipolar outflow extending up to a ~10^4 au scale and a rotating molecular envelope extending over >10^4 au. The interaction with the hot energetic FUor wind, which was observed using optical spectroscopy, appears limited to a region within ~400 au of the star. The envelope mass and the collimation of the extended CO outflow suggest that the progenitor of this FUor is a low-mass Class I young stellar object (YSO). These parameters for V900 Mon, another FUor, and a few FUor-like stars are consistent with the idea that FUor outbursts are associated with normal YSOs. The continuum emission is marginally resolved in our observations with a 0."2x0."15 (~300x225 au) beam, and a Gaussian model provides a deconvolved FWHM of ~90 au. The emission is presumably associated with a dusty circumstellar disk, plus a possible contribution from a wind or a wind cavity close to the star. The warm compact nature of the disk continuum emission could be explained with viscous heating of the disk, while gravitational fragmentation in the outer disk and/or a combination of grain growth and their inward drift may also contribute to its compact nature.

Some circumstellar disks are observed to show prominent spiral arms in infrared scattered light or (sub-)millimeter dust continuum. The spirals might be formed from self-gravity, shadows, or planet-disk interactions. Recently, it was hypothesized that massive vortices can drive spiral arms in protoplanetary disks in a way analogous to planets. In this paper, we study the basic properties of vortex-driven spirals by the Rossby Wave Instability in 2D hydrodynamics simulations. We study how the surface density contrast, the number, and the shape of vortex-driven spirals depend on the properties of the vortex. We also compare vortex-driven spirals with those induced by planets. The surface density contrast of vortex-driven spirals in our simulations are comparable to those driven by a sub-thermal mass planet, typically a few to a few tens of Earth masses. In addition, different from the latter, the former is not sensitive to the mass of the vortex. Vortex-driven spiral arms are not expected to be detectable in current scattered light observations, and the prominent spirals observed in scattered light in a few protoplanetary disks, such as SAO 206462 (HD 135344B), MWC 758, and LkHa 330, are unlikely to be induced by the candidate vortices in them.

Theoretical studies suggest that a giant planet around the young star MWC 758 could be responsible for driving the spiral features in its circumstellar disk. Here, we present a deep imaging campaign with the Large Binocular Telescope with the primary goal of imaging the predicted planet. We present images of the disk in two epochs in the $L^{\prime}$ filter (3.8 $\mu m$) and a third epoch in the $M^{\prime}$ filter (4.8 $\mu m$). The two prominent spiral arms are detected in each observation, which constitute the first images of the disk at $M^\prime$, and the deepest yet in $L^\prime$ ($\Delta L^\prime=$12.1 exterior to the disk at 5$\sigma$ significance). We report the detection of a S/N$\sim$3.9 source near the end of the Sourthern arm, and, from the source's detection at a consistent position and brightness during multiple epochs, we establish a $\sim$90% confidence-level that the source is of astrophysical origin. We discuss the possibilities that this feature may be a) an unresolved disk feature, and b) a giant planet responsible for the spiral arms, with several arguments pointing in favor of the latter scenario. We present additional detection limits on companions exterior to the spiral arms, which suggest that a $\lesssim$4 M$_{Jup}$ planet exterior to the spiral arms could have escaped detection. Finally, we do not detect the companion candidate interior to the spiral arms reported recently by Reggiani et al. (2018), although forward modelling suggests that such a source would have likely been detected.

We investigate how the detectability of signatures of self-gravity in a protoplanetary disc depends on its temporal evolution. We run a one-dimensional model for secular timescales to follow the disc mass as a function of time. We then combine this with three-dimensional global hydrodynamics simulations that employ a hybrid radiative transfer method to approximate realistic heating and cooling. We simulate ALMA continuum observations of these systems, and find that structures induced by the gravitational instability (GI) are readily detectable when $q=M_\mathrm{disc}/M_*\gtrsim 0.25$ and $R_\mathrm{outer}\lesssim 100$ au. The high accretion rate generated by gravito-turbulence in such a massive disc drains its mass to below the detection threshold in $\sim10^4$ years, or approximately 1 % of the typical disc lifetime. Therefore, discs with spiral arms detected in ALMA dust observations, if generated by self-gravity, must either be still receiving infall to maintain a high $q$ value, or have just emerged from their natal envelope. Detection of substructure in systems with lower $q$ is possible, but would require a specialist integration with the most extended configuration over several days. This disfavours the possibility of GI-caused spiral structure in systems with $q<0.25$ being detected in relatively short integration times, such as those found in the DSHARP ALMA survey (Andrews et al. 2018; Huang et al. 2018). We find no temporal dependence of detectability on dynamical timescales

Protoplanets can produce structures in protoplanetary disks via gravitational disk-planet interactions. Once detected, such structures serve as signposts of planet formation. Here we investigate the kinematic signatures in disks produced by multi-Jupiter mass ($M_{\rm J}$) planets using 3D hydrodynamics and radiative transfer simulations. Such a planet opens a deep gap, and drives transonic vertical motions inside. Such motions include both a bulk motion of the entire half-disk column, and turbulence on scales comparable to and smaller than the scale height. They significantly broaden molecular lines from the gap, producing double-peaked line profiles at certain locations, and a kinematic velocity dispersion comparable to thermal after azimuthal averaging. The same planet does not drive fast vertical motions outside the gap, except at the inner spiral arms and the disk surface. Searching for line broadening induced by multi-$M_{\rm J}$ planets inside gaps requires an angular resolution comparable to the gap width, an assessment of the gap gas temperature to within a factor of 2, and a high sensitivity needed to detect line emission from the gap.

Over the past three decades instruments on the ground and in space have discovered thousands of planets outside the solar system. These observations have given rise to an astonishingly detailed picture of the demographics of short-period planets, but are incomplete at longer periods where both the sensitivity of transit surveys and radial velocity signals plummet. Even more glaring is that the spectra of planets discovered with these indirect methods are either inaccessible (radial velocity detections) or only available for a small subclass of transiting planets with thick, clear atmospheres. Direct detection can be used to discover and characterize the atmospheres of planets at intermediate and wide separations, including non-transiting exoplanets. Today, a small number of exoplanets have been directly imaged, but they represent only a rare class of young, self-luminous super-Jovian-mass objects orbiting tens to hundreds of AU from their host stars. Atmospheric characterization of planets in the <5 AU regime, where radial velocity (RV) surveys have revealed an abundance of other worlds, is technically feasible with 30-m class apertures in combination with an advanced AO system, coronagraph, and suite of spectrometers and imagers. There is a vast range of unexplored science accessible through astrometry, photometry, and spectroscopy of rocky planets, ice giants, and gas giants. In this whitepaper we will focus on one of the most ambitious science goals --- detecting for the first time habitable-zone rocky (<1.6 R_Earth) exoplanets in reflected light around nearby M-dwarfs

Rings and radial gaps are ubiquitous in protoplanetary disks, yet their possible connection to planet formation is currently subject to intense debates. In principle, giant planet formation leads to wide gaps which separate the gas and dust mass reservoir in the outer disk, while lower mass planets lead to shallow gaps which are manifested mainly on the dust component. We used the Atacama Large Millimeter/submillimeter Array (ALMA) to observe the star HD169142, host to a prominent disk with deep wide gaps that sever the disk into inner and outer regions. The new ALMA high resolution images allow for the outer ring to be resolved as three narrow rings. The HD169142 disk thus hosts both the wide gaps trait of transition disks and a narrow ring system similar to those observed in the TW Hya and HL Tau systems. The mass reservoir beyond a deep gap can thus host ring systems. The observed rings are narrow in radial extent (width/radius of 1.5/57.3, 1.8/64.2 and 3.4/76.0, in au) and have asymmetric mutual separations: the first and middle ring are separated by 7 au while the middle and outermost ring are distanced by ~12 au. Using hydrodynamical modeling we found that a simple explanation, involving a single migrating low mass planet (10 M$_\oplus$), entirely accounts for such an apparently complex phenomenon. Inward migration of the planet naturally explains the ring's asymmetric mutual separation. The isolation of HD169142's outer rings thus allows a proof of concept to interpret the detailed architecture of the outer region of protoplanetary disks with low mass planet formation of mini-Neptune's size, i.e. as in the protosolar nebula.

Two studies utilizing sparse aperture masking (SAM) interferometry and $H_{\rm \alpha}$ differential imaging have reported multiple jovian companions around the young solar-mass star, LkCa 15 (LkCa 15 bcd): the first claimed direct detection of infant, newly-formed planets ("protoplanets"). We present new near-infrared direct imaging/spectroscopy from the SCExAO system coupled with the CHARIS integral field spectrograph and multi-epoch thermal infrared imaging from Keck/NIRC2 of LkCa 15 at high Strehl ratios. These data provide the first direct imaging look at the same wavelengths and in the same locations where previous studies identified the LkCa 15 protoplanets and thus offer the first decisive test of their existence. The data do not reveal these planets. Instead, we resolve extended emission tracing a dust disk with a brightness and location comparable to that claimed for LkCa 15 bcd. Forward-models attributing this signal to orbiting planets are inconsistent with the combined SCExAO/CHARIS and Keck/NIRC2 data. An inner disk provides a more compelling explanation for the SAM detections and perhaps also the claimed $H_{\alpha}$ detection of LkCa 15 b. We conclude that there is currently no clear, direct evidence for multiple protoplanets orbiting LkCa 15, although the system likely contains at least one unseen jovian companion. To identify jovian companions around LkCa 15 from future observations, the inner disk should be detected and its effect modeled, removed, and shown to be distinguishable from planets. Protoplanet candidates identified from similar systems should likewise be clearly distinguished from disk emission through modeling.

The most dramatic phases of terrestrial planet formation are thought to be oligarchic and chaotic growth, on timescales of up to 100-200 Myr, when violent impacts occur between large planetesimals of sizes up to proto-planets. Such events are marked by the production of large amounts of debris as has been observed in some exceptionally bright and young debris disks (termed extreme debris disks). Here we report five years of Spitzer measurements of such systems around two young solar-type stars: ID8 and P1121. The short-term (weekly to monthly) and long-term (yearly) disk variability is consistent with the aftermaths of large impacts involving large asteroid-size bodies. We demonstrate that an impact-produced clump of optically thick dust, under the influence of the dynamical and viewing geometry effects, can produce short-term modulation in the disk light curves. The long-term disk flux variation is related to the collisional evolution within the impact-produced fragments once released into a circumstellar orbit. The time-variable behavior observed in the P1121 system is consistent with a hypervelocity impact prior to 2012 that produced vapor condensates as the dominant impact product. Two distinct short-term modulations in the ID8 system argue for two violent impacts at different times and locations. Its long-term variation is consistent with the collisional evolution of two different populations of impact-produced debris dominated by either vapor condensates or escaping boulders. The bright, variable emission from the dust produced in large impacts from extreme debris disks provides a unique opportunity to study violent events during the era of terrestrial planet formation.

We analyze the highest-resolution millimeter continuum and near-infrared (NIR) scattered-light images presented to date of the circumbinary disk orbiting V4046 Sgr, a ~20 Myr old actively accreting, close binary T Tauri star system located a mere 72.4 pc from Earth. We observed the disk with the Atacama Large Millimeter/submillimeter Array (ALMA) at 870-micron during Cycle 4, and we analyze these data in conjunction with archival NIR (H band) polarimetric images obtained with SPHERE/IRDIS on the ESO Very Large Telescope. At 0.3'' (20 au) resolution, the 870-micron image reveals a marginally resolved ring that peaks at ~32 au and has an extension of ~ 90 au. We infer a lower limit on dust mass of ~ 60.0 M_earth within the 870-micron ring, and confirm that the ring is well aligned with the larger-scale gaseous disk. A second, inner dust ring is also tentatively detected in the ALMA observations; its position appears coincident with the inner (~14 au radius) ring detected in scattered light. Using synthetic 870 micron and H-band images obtained from disk-planet interaction simulations, we attempt to constrain the mass of the putative planet orbiting at 20 au. Our trials suggest that a circumbinary Jovian-mass planet may be responsible for generating the dust ring and gap structures detected within the disk. We discuss the longevity of the gas-rich disk orbiting V4046 Sgr in the context of the binary nature of the system.

In recent years, the disk populations in a number of young star-forming regions have been surveyed with ALMA. Understanding the disk properties and their correlation with those of the central star is critical to understand planet formation. In particular, a decrease of the average measured disk dust mass with the age of the region has been observed. We conducted high-sensitivity continuum ALMA observations of 43 Class II young stellar objects in CrA at 1.3 mm (230 GHz). The typical spatial resolution is 0.3". The continuum fluxes are used to estimate the dust masses of the disks, and a survival analysis is performed to estimate the average dust mass. We also obtained new VLT/X-Shooter spectra for 12 of the objects in our sample. 24 disks are detected, and stringent limits have been put on the average dust mass of the non-detections. Accounting for the upper limits, the average disk mass in CrA is $6\pm3\,\rm M_\oplus$, significantly lower than that of disks in other young (1-3 Myr) star forming regions (e.g. Lupus) and appears consistent with the 5-10 Myr old Upper Sco. The position of the stars in our sample on the HR diagram, however, seems to confirm that that CrA has age similar to Lupus. Neither external photoevaporation nor a lower than usual stellar mass distribution can explain the low disk masses. On the other hand, a low-mass disk population could be explained if the disks are small, which could happen if the parent cloud has a low temperature or intrinsic angular momentum, or if the the angular momentum of the cloud is removed by some physical mechanism such as magnetic braking. In order to fully explain and understand the dust mass distribution of protoplanetary disks and their evolution, it may also be necessary to take into consideration the initial conditions of star and disk formation process, which may vary from region to region, and affect planet formation.

We present H-band polarized scattered light imagery and JHK high-contrast spectroscopy of the protoplanetary disk around HD 163296 observed with the HiCIAO and SCExAO/CHARIS instruments at Subaru Observatory. The polarimetric imagery resolve a broken ring structure surrounding HD 163296 that peaks at a distance along the major axis of 0.65 (66 AU) and extends out to 0.98 (100 AU) along the major axis. Our 2011 H-band data exhibit clear axisymmetry, with the NW- and SE- side of the disk exhibiting similar intensities. Our data are clearly different than 2016 epoch H-band observations from VLT/SPHERE that found a strong 2.7x asymmetry between the NW- and SE-side of the disk. Collectively, these results indicate the presence of time variable, non-azimuthally symmetric illumination of the outer disk. Based on our 3D-MCRT modeling of contemporaneous IR spectroscopic and H-band polarized intensity imagery of the system, we suggest that while the system could plausibly host an inclined inner disk component, such a component is unlikely to be responsible for producing the observed time-dependent azimuthal variations in the outer scattered light disk of the system. While our SCExAO/CHARIS data are sensitive enough to recover the planet candidate identified from NIRC2 in the thermal IR, we fail to detect an object with a corresponding JHK brightness estimated from the atmospheric models of Baraffe et al. 2003. This suggests that the candidate is either fainter in JHK bands than model predictions, possibly due to extinction from the disk or atmospheric dust/clouds, or that it is an artifact of the dataset/data processing. Our SCExAO/CHARIS data lower the IR mass limits for planets inferred at larger stellocentric separations; however, these ALMA-predicted protoplanet candidates are currently still consistent with direct imaging constraints.

Resolved ALMA and VLA observations indicate the existence of two dust traps in the protoplanetary disc MWC 758. By means of 2D gas+dust hydrodynamical simulations post-processed with 3D dust radiative transfer calculations, we show that the spirals in scattered light, the eccentric, asymmetric ring and the crescent-shaped structure in the (sub)millimetre can all be caused by two giant planets: a 1.5-Jupiter mass planet at 35 au (inside the spirals) and a 5-Jupiter mass planet at 140 au (outside the spirals). The outer planet forms a dust-trapping vortex at the inner edge of its gap (at ~85 au), and the continuum emission of this dust trap reproduces the ALMA and VLA observations well. The outer planet triggers several spiral arms which are similar to those observed in polarised scattered light. The inner planet also forms a vortex at the outer edge of its gap (at ~50 au), but it decays faster than the vortex induced by the outer planet, as a result of the disc's turbulent viscosity. The vortex decay can explain the eccentric inner ring seen with ALMA as well as the low signal and larger azimuthal spread of this dust trap in VLA observations. Finding the thermal and kinematic signatures of both giant planets could verify the proposed scenario.

We have now accumulated a wealth of observations of the planet-formation environment and of mature planetary systems. These data allow us to test and refine theories of gas-giant planet formation by placing constraints on the conditions and timescale of this process. Yet a number of fundamental questions remain unanswered about how protoplanets accumulate material, their photospheric properties and compositions, and how they interact with protoplanetary disks. While we have begun to detect protoplanet candidates during the last several years, we are presently only sensitive to the widest separation, highest mass / accretion rate cases. Current observing facilities lack the angular resolution and inner working angle to probe the few-AU orbital separations where giant planet formation is thought to be most efficient. They also lack the contrast to detect accretion rates that would form lower mass gas giants and ice giants. Instruments and telescopes coming online over the next decade will provide high contrast in the inner giant-planet-forming regions around young stars, allowing us to build a protoplanet census and to characterize planet formation in detail for the first time.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the CO ($J$=2--1) line emission from the protoplanetary disk around T-Tauri star SU Aurigae (hereafter SU Aur). Previous observations in optical and near infrared wavelengths find a unique structure in SU Aur. One of the highlights of the observational results is that an extended tail-like structure is associated with the disk, indicating mass transfer from or into the disk. Here we report the discovery of the counterpart of the tail-like structure in CO gas extending more than 1000 au long. Based on geometric and kinematic perspectives, both of the disk and the tail-like structure components physically connect to each other. Several theoretical studies predicted the observed tail-like structure via the following possible scenarios, 1) a gaseous stream from the molecular cloud remnant, 2) collision with a (sub)stellar intruder or a gaseous blob from the ambient cloud, and 3) ejection of a planetary or brown dwarf mass object due to gravitational instability via multi-body gravitational interaction. Since the tail-like structures associated with the SU Aur disk is a new example following RW Aurigae, some disks may experience the internal or external interaction and drastically lose mass during disk evolution.

Since the discovery of the multi-ring structure of the HL Tau disk, ALMA data suggest that the dust continuum emission of many, if not all, protoplanetary disks consists of rings and gaps, no matter their spectral type or age. The origin of these gaps so far remains unclear. We present a sample study of 16 disks with multiple ring-like structures in the continuum, using published ALMA archival data, to compare their morphologies and gap locations in a systematic way. The 16 targets range from early to late type stars, from <0.5 Myr to >10 Myr, from ~0.2 to 40 L_Sun and include both full and transitional disks with cleared inner dust cavities. Stellar ages are revised using new Gaia distances. Gap locations are derived using a simple radial fit to the intensity profiles. Using a radiative transfer model, the temperature profiles are computed. The gap radii generally do not correspond to the orbital radii of snow lines of the most common molecules. A snow line model can likely be discarded as a common origin of multi-ring systems. In addition, there are no systematic trends in the gap locations that could be related to resonances of planets. Finally, the outer radius of the disks decreases for the oldest disks in the sample, indicating that if multi-ring disks evolve in a similar way, outer dust rings either dissipate with the gas or grow into planetesimal belts.

The large crescents imaged by ALMA in transition disks suggest that azimuthal dust trapping concentrates the larger grains, but centimetre-wavelengths continuum observations are required to map the distribution of the largest observable grains. A previous detection at ~1cm of an unresolved clump along the outer ring of MWC758 (Clump1), and buried inside more extended sub-mm continuum, motivates followup VLA observations. Deep multiconfiguration integrations reveal the morphology of Clump 1 and additional cm-wave components which we characterize via comparison with a deconvolution of recent 342GHz data (~1mm). Clump1, which concentrates ~1/3 of the whole disk flux density at ~1cm, is resolved as a narrow arc with a deprojected aspect ratio Chi>5.6, and with half the azimuthal width than at 342 GHz. The spectral trends in the morphology of Clump1 are quantitatively consistent with the Lyra-Lin prescriptions for dust trapping in an anticyclonic vortex, provided with porous grains (f~0.2+-0.2) in a very elongated (Chi~14+-3) and cold (T~23+-2K) vortex. The same prescriptions constrain the turbulence parameter alpha and the gas surface density Sigma_g through log10( alpha x Sigma_g /g/cm2)~-2.3+-0.4, thus requiring values for Sigma_g larger than a factor of a few compared to that reported in the literature from the CO isotopologues, if alpha <~ 1E-3. Such physical conditions imply an appreciably optically thick continuum even at cm-wavelengths (tau(33GHz)~0.2). A secondary and shallower peak at 342GHz is about twice fainter relative to Clump1 at 33GHz. Clump2 appears to be less efficient at trapping large grains.

We present Atacama Large Millimeter/submillimeter Array (ALMA) observations of the dust continuum emission at 1.3 mm and 12CO J=2-1 line emission of the transitional disk around DM Tau. DM Tau's disk is thought to possess a dust-free inner cavity inside a few au, from the absence of near-infrared excess on its spectral energy distribution (SED). Previous submillimeter observations were, however, unable to detect the cavity; instead, a dust ring ~20 au in radius was seen. The excellent angular resolution achieved in the new ALMA observations, 43 x 31 mas, allows discovery of a 4 au radius inner dust ring, confirming previous SED modeling results. This inner ring is symmetric in continuum emission, but asymmetric in 12CO emission. The known (outer) dust ring at ~20 au is recovered and shows azimuthal asymmetry with a strong-weak side contrast of ~1.3. The gap between these two rings is depleted by a factor of ~40 in dust emission relative to the outer ring. An extended outer dust disk is revealed, separated from the outer ring by another gap. The location of the inner ring is comparable to that of the main asteroid belt in the solar system. As a disk with a "proto-asteroid belt," the DM Tau system offers valuable clues to disk evolution and planet formation in the terrestrial planet forming region.

ALMA has revolutionized our view of protoplanetary disks, revealing structures such as gaps, rings and asymmetries that indicate dust trapping as an important mechanism in the planet formation process. However, the high resolution images have also shown that the optically thin assumption for millimeter continuum emission may not be valid and the low values of the spectral index may be related to optical depth rather than dust growth. Longer wavelength observations are essential to properly disentangle these effects. The high sensitivity and spatial resolution of the next-generation Very Large Array (ngVLA) will open up the possibilities to spatially resolve disk continuum emission at centimeter wavelengths and beyond, which allows the study of dust growth in disks in the optically thin regime and further constrain models of planet formation.

Pairs of azimuthal intensity decrements at near symmetric locations have been seen in a number of protoplanetary disks. They are most commonly interpreted as the two shadows cast by a highly misaligned inner disk. Direct evidence of such an inner disk, however, remain largely illusive, except in rare cases. In 2012, a pair of such shadows were discovered in scattered light observations of the near face-on disk around 2MASS J16042165-2130284, a transitional object with a cavity $\sim$60 AU in radius. The star itself is a `dipper', with quasi-periodic dimming events on its light curve, commonly hypothesized as caused by extinctions by transiting dusty structures in the inner disk. Here, we report the detection of a gas disk inside the cavity using ALMA observations with $\sim0$\farcs2 angular resolution. A twisted butterfly pattern is found in the moment 1 map of CO (3-2) emission line towards the center, which is the key signature of a high misalignment between the inner and outer disks. In addition, the counterparts of the shadows are seen in both dust continuum emission and gas emission maps, consistent with these regions being cooler than their surroundings. Our findings strongly support the hypothesized misaligned-inner-disk origin of the shadows in the J1604-2130 disk. Finally, the inclination of inner disk would be close to -45 $^{\circ}$ in contrast with 45 $^{\circ}$; it is possible that its internal asymmetric structures cause the variations on the light curve of the host star.

ALMA has found multiple dust gaps and rings in a number of protoplanetary disks in continuum emission at millimeter wavelengths. The origin of such structures is in debate. Recently, we documented how one super-Earth planet can open multiple (up to five) dust gaps in a disk with low viscosity ($\alpha\lesssim10^{-4}$). In this paper, we examine how the positions, depths, and total number of gaps opened by one planet depend on input parameters, and apply our results to real systems. Gap locations (equivalently, spacings) are the easiest metric to use when making comparisons between theory and observations, as positions can be robustly measured. We fit the locations of gaps empirically as functions of planet mass and disk aspect ratio. We find that the locations of the double gaps in HL Tau and TW Hya, and of all three gaps in HD 163296, are consistent with being opened by a sub-Saturn mass planet. This scenario predicts the locations of other gaps in HL Tau and TW Hya, some of which appear consistent with current observations. We also show how the Rossby wave instability may develop at the edges of several gaps and result in multiple dusty vortices, all caused by one planet. A planet as low in mass as Mars may produce multiple dust gaps in the terrestrial planet forming region.

We present near-IR imaging polarimetry of five classical FU Ori-type objects (FU Ori, V1057 Cyg, V1515 Cyg, V1735 Cyg, Z CMa) with a $\sim$0\farcs1 resolution observed using HiCIAO+AO188 at Subaru Telescope. We observed scattered light associated with circumstellar dust around four of them (i.e., all but V1515 Cyg). Their polarized intensity distribution shows a variety of morphologies with arms, tails or streams, spikes and fragmented distributions, many of which were reported in our previous paper. The morphologies of these reflection nebulae significantly differ from many other normal young stellar objects (Class I-II objects). These structures are attributed to gravitationally unstable disks, trails of clump ejections, dust blown by a wind or a jet, and a stellar companion. We can consistently explain our results with the scenario that their accretion outbursts (FUor outbursts) are triggered by gravitationally fragmenting disks, and with the hypothesis that many low-mass young stellar objects experience such outbursts.

The Planet Formation Imager (PFI, www.planetformationimager.org) is a next-generation infrared interferometer array with the primary goal of imaging the active phases of planet formation in nearby star forming regions. PFI will be sensitive to warm dust emission using mid-infrared capabilities made possible by precise fringe tracking in the near-infrared. An L/M band combiner will be especially sensitive to thermal emission from young exoplanets (and their disks) with a high spectral resolution mode to probe the kinematics of CO and H2O gas. In this paper, we give an overview of the main science goals of PFI, define a baseline PFI architecture that can achieve those goals, point at remaining technical challenges, and suggest activities today that will help make the Planet Formation Imager facility a reality.

Spiral arm structures seen in scattered light observations of protoplanetary disks can potentially serve as signposts of planetary companions. They can also lend unique insights into disk masses, which are critical in setting the mass budget for planet formation but are difficult to determine directly. A surprisingly high fraction of disks that have been well-studied in scattered light have spiral arms of some kind (8/29), as do a high fraction (6/11) of well-studied Herbig intermediate mass stars (i.e., Herbig stars $> 1.5M_\odot$). Here we explore the origin of spiral arms in Herbig systems by studying their occurrence rates, disk properties, and stellar accretion rates. We find that two-arm spirals are more common in disks surrounding Herbig intermediate mass stars than are directly imaged giant planet companions to mature A and B stars. If two-arm spirals are produced by such giant planets, this discrepancy suggests that giant planets are much fainter than predicted by hot start models. In addition, the high stellar accretion rates of Herbig stars, if sustained over a reasonable fraction of their lifetimes, suggest that disk masses are much larger than inferred from their submillimeter continuum emission.As a result, gravitational instability is a possible explanation for multi-arm spirals. Future observations can lend insights into the issues raised here.

Large-scale spiral arms have been revealed in scattered light images of a few protoplanetary disks. Theoretical models suggest that such arms may be driven by and co-rotate with giant planets, which has called for remarkable observational efforts to look for them. By examining the rotation of the spiral arms for the MWC 758 system over a 10-yr timescale, we are able to provide dynamical constraints on the locations of their perturbers. We present reprocessed Hubble Space Telescope (HST)/NICMOS F110W observations of the target in 2005, and the new Keck/NIRC2 $L'$-band observations in 2017. MWC~758's two well-known spiral arms are revealed in the NICMOS archive at the earliest observational epoch. With additional Very Large Telescope (VLT)/SPHERE data, our joint analysis leads to a pattern speed of ${0.6^\circ}^{+3.3^\circ}_{-0.6^\circ}\, \mathrm{yr}^{-1}$ at $3\sigma$ for the two major spiral arms. If the two arms are induced by a perturber on a near-circular orbit, its best fit orbit is at $89$ au ($0.59"$), with a $3\sigma$ lower limit of 30 au ($0.20"$). This finding is consistent with the simulation prediction of the location of an arm-driving planet for the two major arms in the system.

We present $H$- and $K_{\rm s}$-bands observations of the LkH$\alpha$ 330 disk with a multi-band detection of the large gap and spiral-like structures. The morphology of the outer disk ($r\sim$$0\farcs3$) at PA=0--45$^\circ$ and PA=180--290$^\circ$ are likely density wave-induced spirals and comparison between our observational results and simulations suggests a planet formation. We have also investigated the azimuthal profiles at the ring and the outer-disk regions as well as radial profiles in the directions of the spiral-like structures and semi-major axis. Azimuthal analysis shows a large variety in wavelength and implies that the disk has non-axisymmetric dust distributions. The radial profiles in the major-axis direction (PA=$271^\circ$) suggest that the outer region ($r\geq0\farcs25$) may be influenced by shadows of the inner region of the disk. The spiral-like directions (PA=10$^\circ$ and 230$^\circ$) show different radial profiles, which suggests that the surfaces of the spiral-like structures are highly flared and/or have different dust properties. Finally, a color-map of the disk shows a lack of an outer eastern region in the $H$-band disk, which may hint the presence of an inner object that casts a directional shadow onto the disk.

Spatially resolved structures in protoplanetary disks hint at unseen planets. Previous imaging observations of the transitional disk around MWC 758 revealed an inner cavity, a ring-like outer disk, emission clumps, and spiral arms, all possibly generated by companions. We present ALMA dust continuum observations of MWC 758 at 0.87 millimeter (mm) wavelength with 43$\times$39 mas angular resolution (6.9$\times$6.2 AU) and 20 $\mu$Jy beam$^{-1}$ rms. The central sub-mm emission cavity is revealed to be eccentric; once deprojected, its outer edge can be well-fitted by an ellipse with an eccentricity of 0.1 and one focus on the star. The broad ring-like outer disk is resolved into three narrow rings with two gaps in between. The outer two rings tentatively show the same eccentricity and orientation as the innermost ring bounding the inner cavity. The two previously known dust emission clumps are resolved in both the radial and azimuthal directions, with radial widths equal to $\sim$4$\times$ the local scale height. Only one of the two spiral arms previously imaged in near-infrared (NIR) scattered light is revealed in ALMA dust emission, at a slightly larger stellocentric distance owing to projection effects. We also submit evidence of disk truncation at $\sim$100 AU based on comparing NIR imaging observations with models. The spirals, the north clump, and the truncated disk edge are all broadly consistent with the presence of one companion exterior to the spirals at roughly 100 AU.

We present ALMA 0.87 mm continuum, HCO+ J=4--3 emission line, and CO J=3--2 emission line data of the disk of material around the young, Sun-like star PDS 70. These data reveal the existence of a possible two component transitional disk system with a radial dust gap of 0."2 +/- 0."05, an azimuthal gap in the HCO+ J=4--3 moment zero map, as well as two bridge-like features in the gas data. Interestingly these features in the gas disk have no analogue in the dust disk making them of particular interest. We modeled the dust disk using the Monte Carlo radiative transfer code HOCHUNK3D (Whitney et al. 2013) using a two disk components. We find that there is a radial gap that extends from 15-60 au in all grain sizes which differs from previous work.

The discovery of thousands of exoplanets over the last couple of decades has shown that the birth of planets is a very efficient process in nature. Theories invoke a multitude of mechanisms to describe the assembly of planets in the disks around pre-main-sequence stars, but observational constraints have been sparse on account of insufficient sensitivity and resolution. Understanding how planets form and interact with their parental disk is crucial also to illuminate the main characteristics of a large portion of the full population of planets that is inaccessible to current and near-future observations. This White Paper describes some of the main issues for our current understanding of the formation and evolution of planets, and the critical contribution expected in this field by the Next Generation Very Large Array.

HD 100453AB is a 10+/-2 Myr old binary whose protoplanetary disk was recently revealed to host a global two-armed spiral structure. Given the relatively small projected separation of the binary (1.05", or ~108 au), gravitational perturbations by the binary seemed to be a likely driving force behind the formation of the spiral arms. However, the orbit of these stars remained poorly understood, which prevented a proper treatment of the dynamical influence of the companion on the disk. We observed HD 100453AB between 2015-2017 utilizing extreme adaptive optics systems on the Very Large Telescope and Magellan Clay Telescope. We combined the astrometry from these observations with published data to constrain the parameters of the binary's orbit to a=1.06"+/-0.09", e=0.17+/-0.07, and i=32.5+/- 6.5 degrees. We utilized publicly available ALMA CO data to constrain the inclination of the disk to i~28 degrees, which is relatively co-planar with the orbit of the companion and consistent with previous estimates from scattered light images. Finally, we input these constraints into hydrodynamical and radiative transfer simulations to model the structural evolution of the disk. We find that the spiral structure and truncation of the circumprimary disk in HD 100453 are consistent with a companion-dirven origin. Furthermore, we find that the primary star's rotation, its outer disk, and the companion exhibit roughly the same direction of angular momentum, and thus the system likely formed from the same parent body of material.

We present high resolution ($0.11"\times0.06"$) 3mm ALMA observations of the highly inclined transition disk around the star T Cha. Our continuum image reveals multiple dust structures: an inner disk, a spatially resolved dust gap, and an outer ring. When fitting sky-brightness models to the real component of the 3mm visibilities, we infer that the inner emission is compact ($\le1$au in radius), the gap width is between 18-28 au, and the emission from the outer ring peaks at $\sim36$ au. We compare our ALMA image with previously published 1.6$\mu$m VLT/SPHERE imagery. This comparison reveals that the location of the outer ring is wavelength dependent. More specifically, the peak emission of the 3mm ring is at a larger radial distance than that of the 1.6$\mu$m ring, suggesting that millimeter-sized grains in the outer disk are located further away from the central star than micron-sized grains. We discuss different scenarios to explain our findings, including dead zones, star-driven photoevaporation, and planet-disk interactions. We find that the most likely origin of the dust gap is from an embedded planet, and estimate --- for a single planet scenario --- that T Cha's gap is carved by a $1.2M_{jup}$ planet.

We study the quasi-normal modes of asymptotically anti-de Sitter black holes in a class of shift-symmetric Horndeski theories where a gravitational scalar is derivatively coupled to the Einstein tensor. The space-time differs from exact Schwarzschild-anti-de Sitter, resulting in a different effective potential for the quasi-normal modes and a different spectrum. We numerically compute this spectrum for a massless test scalar coupled both minimally to the metric, and non-minimally to the gravitational scalar. We find interesting differences from the Schwarzschild-anti-de Sitter black hole found in general relativity.

High contrast imaging instruments such as GPI and SPHERE are discovering gap structures in protoplanetary disks at an ever faster pace. Some of these gaps may be opened by planets forming in the disks. In order to constrain planet formation models using disk observations, it is crucial to find a robust way to quantitatively back out the properties of the gap-opening planets, in particular their masses, from the observed gap properties, such as their depths and widths. Combing 2D and 3D hydrodynamics simulations with 3D radiative transfer simulations, we investigate the morphology of planet-opened gaps in near-infrared scattered light images. Quantitatively, we obtain correlations that directly link intrinsic gap depths and widths in the gas surface density to observed depths and widths in images of disks at modest inclinations under finite angular resolution. Subsequently, the properties of the surface density gaps enable us to derive the disk scale height at the location of the gap $h$, and to constrain the quantity $M_{\rm p}^2/\alpha$, where $M_{\rm p}$ is the mass of the gap-opening planet and $\alpha$ characterizes the viscosity in the gap. As examples, we examine the gaps recently imaged by VLT/SPHERE, Gemini/GPI, and Subaru/HiCIAO in HD 97048, TW Hya, HD 169142, LkCa 15, and RX J1615.3-3255. Scale heights of the disks and possible masses of the gap-opening planets are derived assuming each gap is opened by a single planet. Assuming $\alpha=10^{-3}$, the derived planet mass in all cases are roughly between 0.1-1 $M_{\rm J}$.

Recently, high angular resolution imaging instruments such as SPHERE and GPI have discovered many spiral-arm-like features in near-infrared scattered light images of protoplanetary disks. Theory and simulations have suggested that these arms are most likely excited by planets forming in the disks; however, a quantitative relation between the arm-to-disk brightness contrast and planet mass is still missing. Using 3D hydrodynamics and radiative transfer simulations, we examine the morphology and contrast of planet-induced arms in disks. We find a power-law relation for the face-on arm contrast ($\delta_{\rm max}$) as a function of planet mass ($M_{\rm p}$) and disk aspect ratio ($h/r$): $\delta_{\rm max}\approx((M_{\rm p}/M_{\rm J})/(h/r)^{1.38})^{0.22}$. With current observational capability, at a 30 AU separation, the minimum planet mass for driving detectable arms in a disk around a 1 Myr 1$M_\odot$ star at 140 pc at low inclinations is around Saturn mass. For planets more massive than Neptune masses, they typically drive multiple arms. Therefore in observed disks with spirals, it is unlikely that each spiral arm originates from a different planet. We also find only massive perturbers with at least multi-Jupiter masses are capable of driving bright arms with $\delta_{\rm max}\gtrsim2$ as found in SAO 206462, MWC 758, and LkH$\alpha$~330, and these arms do not follow linear wave propagation theory. Additionally, we find the morphology and contrast of the primary and secondary arms are largely unaffected by a modest level of viscosity with $\alpha\lesssim0.01$. Finally, the contrast of the arms in the SAO 206462 disk suggests that the perturber SAO 206462 b at $\sim100$ AU is about $5-10M_{\rm J}$ in mass.

We investigate the observational signatures of super-Earths (i.e., Earth-to-Neptune mass planets) in their natal disks of gas and dust. Combining two-fluid global hydrodynamics simulations with a radiative transfer code, we calculate the distributions of gas and of sub-mm-sized dust in a disk perturbed by a super-Earth, synthesizing images in near-infrared scattered light and the mm-wave thermal continuum for direct comparison with observations. In low viscosity gas ($\alpha\lesssim10^{-4}$), a super-Earth opens two annular gaps to either side of its orbit by the action of Lindblad torques. This double gap and its associated gas pressure gradients cause dust particles to be dragged by gas into three rings: one ring sandwiched between the two gaps, and two rings located at the gap edges farthest from the planet. Depending on system parameters, additional rings may manifest for a single planet. A double gap located at tens of AUs from a host star in Taurus can be detected in the dust continuum by the Atacama Large Millimeter Array (ALMA) at an angular resolution of ~0".03 after two hours of integration. Ring and gap features persist in a variety of background disk profiles, last for thousands of orbits, and change their relative positions and dimensions depending on the speed and direction of planet migration. Candidate double gaps have been observed by ALMA in systems like HL Tau (D5 and D6) and TW Hya (at 37 and 43 AU); we submit that each double gap is carved by one super-Earth in nearly inviscid gas.

We present optical and near-infrared high contrast images of the transitional disk HD 100546 taken with the Magellan Adaptive Optics system (MagAO) and the Gemini Planet Imager (GPI). GPI data include both polarized intensity and total intensity imagery, and MagAO data are taken in Simultaneous Differential Imaging mode at H{\alpha}. The new GPI H -band total intensity data represent a significant enhancement in sensitivity and field rotation compared to previous data sets and enable a detailed exploration of substructure in the disk. The data are processed with a variety of differential imaging techniques (polarized, angular, reference, and simultaneous differential imaging) in an attempt to identify the disk structures that are most consistent across wavelengths, processing techniques, and algorithmic parameters. The inner disk cavity at 15 au is clearly resolved in multiple datasets, as are a variety of spiral features. While the cavity and spiral structures are identified at levels significantly distinct from the neighboring regions of the disk under several algorithms and with a range of algorithmic parameters, emission at the location of HD 100546 c varies from point-like under aggressive algorithmic parameters to a smooth continuous structure with conservative parameters, and is consistent with disk emission. Features identified in the HD100546 disk bear qualitative similarity to computational models of a moderately inclined two-armed spiral disk, where projection effects and wrapping of the spiral arms around the star result in a number of truncated spiral features in forward-modeled images.

We report ALMA Cycle 2 observations of 230 GHz (1.3 mm) dust continuum emission, and $^{12}$CO, $^{13}$CO, and C$^{18}$O J = 2-1 line emission, from the Upper Scorpius transitional disk [PZ99] J160421.7-213028, with an angular resolution of ~0".25 (35 AU). Armed with these data and existing H-band scattered light observations, we measure the size and depth of the disk's central cavity, and the sharpness of its outer edge, in three components: sub-$\mu$m-sized "small" dust traced by scattered light, millimeter-sized "big" dust traced by the millimeter continuum, and gas traced by line emission. Both dust populations feature a cavity of radius $\sim$70 AU that is depleted by factors of at least 1000 relative to the dust density just outside. The millimeter continuum data are well explained by a cavity with a sharp edge. Scattered light observations can be fitted with a cavity in small dust that has either a sharp edge at 60 AU, or an edge that transitions smoothly over an annular width of 10 AU near 60 AU. In gas, the data are consistent with a cavity that is smaller, about 15 AU in radius, and whose surface density at 15 AU is $10^{3\pm1}$ times smaller than the surface density at 70 AU; the gas density grades smoothly between these two radii. The CO isotopologue observations rule out a sharp drop in gas surface density at 30 AU or a double-drop model as found by previous modeling. Future observations are needed to assess the nature of these gas and dust cavities, e.g., whether they are opened by multiple as-yet-unseen planets or photoevaporation.

We resolved FU Ori at 29-37 GHz using the JVLA with $\sim$0$''$.07 resolution, and performed the complementary JVLA 8-10 GHz observations, the SMA 224 GHz and 272 GHz observations, and compared with archival ALMA 346 GHz observations to obtain the SEDs. Our 8-10 GHz observations do not find evidence for the presence of thermal radio jets, and constrain the radio jet/wind flux to at least 90 times lower than the expected value from the previously reported bolometric luminosity-radio luminosity correlation. The emission at $>$29 GHz may be dominated by the two spatially unresolved sources, which are located immediately around FU Ori and its companion FU Ori S, respectively. Their deconvolved radii at 33 GHz are only a few au. The 8-346 GHz SEDs of FU Ori and FU Ori S cannot be fit by constant spectral indices (over frequency). The more sophisticated models for SEDs suggest that the $>$29 GHz emission is contributed by a combination of free-free emission from ionized gas, and thermal emission from optically thick and optically thin dust components. We hypothesize that dust in the innermost parts of the disks ($\lesssim$0.1 au) has been sublimated, and thus the disks are no more well shielded against the ionizing photons. The estimated overall gas and dust mass based on SED modeling, can be as high as a fraction of a solar mass, which is adequate for developing disk gravitational instability. Our present explanation for the observational data is that the massive inflow of gas and dust due to disk gravitational instability or interaction with a companion/intruder, was piled up at the few au scale due to the development of a deadzone with negligible ionization. The piled up material subsequently triggered the thermal and the MRI instabilities when the ionization fraction in the inner sub-au scale region exceeded a threshold value, leading to the high protostellar accretion rate.

In this paper we analyze in detail a rarely discussed question of gravity wave production from evaporating primordial black holes. These black holes emit gravitons which are, at classical level, registered as gravity waves. We use the latest constraints on their abundance, and calculate the power emitted in gravitons at the time of their evaporation. We then solve the coupled system of equations that gives us the evolution of the frequency and amplitude of gravity waves during the expansion of the universe. The spectrum of gravitational waves that can be detected today depends on multiple factors: fraction of the total energy density which was occupied by primordial black holes, the epoch in which they were formed, and quantities like their mass and angular momentum. We conclude that very small primordial black holes which evaporate before the big-bang nucleosynthesis emit gravitons whose spectral energy fraction today can be as large as $10^{-7.5}$. On the other hand, those which are massive enough so that they still exist now can yield a signal as high as $10^{-6.5}$. However, typical frequencies of the gravity waves from primordial black holes are still too high to be observed with the current and near future gravity wave observations.

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{\mu}m atmospheric windows), the prioritization of 3-5{\mu}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)

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 $\alpha\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.

Dust trapping accelerates the coagulation of dust particles, and thus it represents an initial step toward the formation of planetesimals. We report $H$-band (1.6 um) linear polarimetric observations and 0.87 mm interferometric continuum observations toward a transitional disk around LkHa 330. As results, a pair of spiral arms were detected in the $H$-band emission and an asymmetric (potentially arm-like) structure was detected in the 0.87 mm continuum emission. We discuss the origin of the spiral arm and the asymmetric structure, and suggest that a massive unseen planet is the most plausible explanation. The possibility of dust trapping and grain growth causing the asymmetric structure was also investigated through the opacity index (beta) by plotting the observed SED slope between 0.87 mm from our SMA observation and 1.3 mm from literature. The results imply that grains are indistinguishable from ISM-like dust in the east side ($beta = 2.0 pm 0.5$), but much smaller in the west side $beta = 0.7^{+0.5}_{-0.4}$, indicating differential dust size distribution between the two sides of the disk. Combining the results of near-infrared and submillimeter observations, we conjecture that the spiral arms exist at the upper surface and an asymmetric structure resides in the disk interior. Future observations at centimeter wavelengths and differential polarization imaging in other bands (Y to K) with extreme AO imagers are required to understand how large dust grains form and to further explore the dust distribution in the disk.

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.

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.

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.

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.

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.

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.

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, $\phi_{\rm sep}$, and the planet-to-star mass ratio $q$: $\phi_{\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 $\phi_{\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 $\phi_{\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.

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.

Directly detecting thermal emission from young extrasolar planets allows measurement of their atmospheric composition and luminosity, which is influenced by their formation mechanism. Using the Gemini Planet Imager, we discovered a planet orbiting the \$sim$20 Myr-old star 51 Eridani at a projected separation of 13 astronomical units. Near-infrared observations show a spectrum with strong methane and water vapor absorption. Modeling of the spectra and photometry yields a luminosity of L/LS=1.6-4.0 x 10-6 and an effective temperature of 600-750 K. For this age and luminosity, "hot-start" formation models indicate a mass twice that of Jupiter. This planet also has a sufficiently low luminosity to be consistent with the "cold- start" core accretion process that may have formed Jupiter.

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$\deg$ 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.

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.

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.

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.

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.

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.

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.

The formation scenario of a gapped disk, i.e., transitional disk, and its asymmetry is still under debate. Proposed scenarios such as disk-planet interaction, photoevaporation, grain growth, anticyclonic vortex, eccentricity, and their combinations would result in different radial distributions of the gas and the small (sub-$\mu$m size) and large (millimeter size) dust grains as well as asymmetric structures in a disk. Optical/near-infrared (NIR) imaging observations and (sub-)millimeter interferometry can trace small and large dust grains, respectively; therefore multi-wavelength observations could help elucidate the origin of complicated structures of a disk. Here we report SMA observations of the dust continuum at 1.3~mm and $^{12}$CO~$J=2\rightarrow1$ line emission of the pre-transitional protoplanetary disk around the solar-mass star PDS~70. PDS~70, a weak-lined T Tauri star, exhibits a gap in the scattered light from its disk with a radius of $\sim$65~AU at NIR wavelengths. However, we found a larger gap in the disk with a radius of $\sim$80~AU at 1.3~mm. Emission from all three disk components (the gas and the small and large dust grains) in images exhibits a deficit in brightness in the central region of the disk, in particular, the dust-disk in small and large dust grains has asymmetric brightness. The contrast ratio of the flux density in the dust continuum between the peak position to the opposite side of the disk reaches 1.4. We suggest the asymmetries and different gap-radii of the disk around PDS~70 are potentially formed by several (unseen) accreting planets inducing dust filtration.

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.

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.

We define a new inflationary scenario in which inflation starts naturally after the Big Bang when the energy density drops below some critical value. As a model, we use recently proposed symmetron field whose effective potential depends on the energy density of the environment. At high densities, right after the Big Bang, the potential for the symmetron is trivial, and the field sits in equilibrium at the bottom of the potential. When the density drops below some critical value, the potential changes its shape into a symmetry breaking potential, and the field starts rolling down. This scenario does not require any special initial conditions for inflation to start. We also construct a concrete model with two fields, i.e. with symmetron as an inflaton and an additional scalar field which describes the matter content in the early universe. For the simplest coupling, the amplitude and shape of the power spectrum are the same as in the single field slow-roll inflation.

We have updated our publicly available dust radiative transfer code (HOCHUNK3D) to include new emission processes and various 3-D geometries appropriate for forming stars. The 3-D geometries include warps and spirals in disks, accretion hotspots on the central star, fractal clumping density enhancements, and misaligned inner disks. Additional axisymmetric (2-D) features include gaps in disks and envelopes, "puffed-up inner rims" in disks, multiple bipolar cavity walls, and iteration of disk vertical structure assuming hydrostatic equilibrium. We include the option for simple power-law envelope geometry, which combined with fractal clumping, and bipolar cavities, can be used to model evolved stars as well as protostars. We include non-thermal emission from PAHs and very small grains, and external illumination from the interstellar radiation field. The grid structure was modified to allow multiple dust species in each cell; based on this, a simple prescription is implemented to model dust stratification. We describe these features in detail, and show example calculations of each. Some of the more interesting results include the following: 1) Outflow cavities may be more clumpy than infalling envelopes. 2) PAH emission in high-mass stars may be a better indicator of evolutionary stage than the broadband SED slope; and related to this, 3) externally illuminated clumps and high-mass stars in optically thin clouds can masquerade as YSOs. 4) Our hydrostatic equilibrium models suggest that dust settling is likely ubiquitous in T Tauri disks, in agreement with previous observations.

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.

We present the first near-IR scattered light detection of the transitional disk associated with the Herbig Ae star MWC 758 using data obtained as part of the Strategic Exploration of Exoplanets and Disks with Subaru, and 1.1 micron HST/NICMOS data. While sub-millimeter studies suggested there is a dust-depleted cavity with r=0.35, we find scattered light as close as 0.1 (20-28 AU) from the star, with no visible cavity at H, K', or Ks. We find two small-scaled spiral structures which asymmetrically shadow the outer disk. We model one of the spirals using spiral density wave theory, and derive a disk aspect ratio of h ~ 0.18, indicating a dynamically warm disk. If the spiral pattern is excited by a perturber, we estimate its mass to be 5+3,-4 Mj, in the range where planet filtration models predict accretion continuing onto the star. Using a combination of non-redundant aperture masking data at L' and angular differential imaging with Locally Optimized Combination of Images at K' and Ks, we exclude stellar or massive brown dwarf companions within 300 mas of the Herbig Ae star, and all but planetary mass companions exterior to 0.5. We reach 5-sigma contrasts limiting companions to planetary masses, 3-4 MJ at 1.0 and 2 MJ at 1.55 using the COND models. Collectively, these data strengthen the case for MWC 758 already being a young planetary system.

We report high-resolution (0.07 arcsec) near-infrared polarized intensity images of the circumstellar disk around the star 2MASS J16042165-2130284 obtained with HiCIAO mounted on the Subaru 8.2 m telescope. We present our $H$-band data, which clearly exhibits a resolved, face-on disk with a large inner hole for the first time at infrared wavelengths. We detect the centrosymmetric polarization pattern in the circumstellar material as has been observed in other disks. Elliptical fitting gives the semimajor axis, semiminor axis, and position angle (P.A.) of the disk as 63 AU, 62 AU, and -14 $^{\circ}$, respectively. The disk is asymmetric, with one dip located at P.A.s of $\sim85^{\circ}$. Our observed disk size agrees well with a previous study of dust and CO emission at submillimeter wavelength with Submillimeter Array. Hence, the near-infrared light is interpreted as scattered light reflected from the inner edge of the disk. Our observations also detect an elongated arc (50 AU) extending over the disk inner hole. It emanates at the inner edge of the western side of the disk, extending inward first, then curving to the northeast. We discuss the possibility that the inner hole, the dip, and the arc that we have observed may be related to the existence of unseen bodies within the disk.

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.

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.

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.

We present H-band polarimetric imagery of UX Tau A taken with HiCIAO/AO188 on the Subaru Telescope. UX Tau A has been classified as a pre-transitional disk object, with a gap structure separating its inner and outer disks. Our imagery taken with the 0.15 (21 AU) radius coronagraphic mask has revealed a strongly polarized circumstellar disk surrounding UX Tau A which extends to 120 AU, at a spatial resolution of 0.1 (14 AU). It is inclined by 46 \pm 2 degree as the west side is nearest. Although SED modeling and sub-millimeter imagery suggested the presence of a gap in the disk, with the inner edge of the outer disk estimated to be located at 25 - 30 AU, we detect no evidence of a gap at the limit of our inner working angle (23 AU) at the near-infrared wavelength. We attribute the observed strong polarization (up to 66 %) to light scattering by dust grains in the disk. However, neither polarization models of the circumstellar disk based on Rayleigh scattering nor Mie scattering approximations were consistent with the observed azimuthal profile of the polarization degrees of the disk. Instead, a geometric optics model of the disk with nonspherical grains with the radii of 30 micron meter is consistent with the observed profile. We suggest that the dust grains have experienced frequent collisional coagulations and have grown in the circumstellar disk of UX Tau A.

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 \alpha=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.

We present high-resolution, H-band, imaging observations, collected with Subaru/HiCIAO, of the scattered light from the transitional disk around SAO 206462 (HD 135344B). Although previous sub-mm imagery suggested the existence of the dust-depleted cavity at r~46AU, our observations reveal the presence of scattered light components as close as 0.2" (~28AU) from the star. Moreover, we have discovered two small-scale spiral structures lying within 0.5" (~70AU). We present models for the spiral structures using the spiral density wave theory, and derive a disk aspect ratio of h~0.1, which is consistent with previous sub-mm observations. This model can potentially give estimates of the temperature and rotation profiles of the disk based on dynamical processes, independently from sub-mm observations. It also predicts the evolution of the spiral structures, which can be observable on timescales of 10-20 years, providing conclusive tests of the model. While we cannot uniquely identify the origin of these spirals, planets embedded in the disk may be capable of exciting the observed morphology. Assuming that this is the case, we can make predictions on the locations and, possibly, the masses of the unseen planets. Such planets may be detected by future multi-wavelengths observations.

Transitional circumstellar disks around young stellar objects have a distinctive infrared deficit around 10 microns in their Spectral Energy Distributions (SED), recently measured by the Spitzer Infrared Spectrograph (IRS), suggesting dust depletion in the inner regions. These disks have been confirmed to have giant central cavities by imaging of the submillimeter (sub-mm) continuum emission using the Submillimeter Array (SMA). However, the polarized near-infrared scattered light images for most objects in a systematic IRS/SMA cross sample, obtained by HiCIAO on the Subaru telescope, show no evidence for the cavity, in clear contrast with SMA and Spitzer observations. Radiative transfer modeling indicates that many of these scattered light images are consistent with a smooth spatial distribution for micron-sized grains, with little discontinuity in the surface density of the micron-sized grains at the cavity edge. Here we present a generic disk model that can simultaneously account for the general features in IRS, SMA, and Subaru observations. Particularly, the scattered light images for this model are computed, which agree with the general trend seen in Subaru data. Decoupling between the spatial distributions of the micron-sized dust and mm-sized dust inside the cavity is suggested by the model, which, if confirmed, necessitates a mechanism, such as dust filtration, for differentiating the small and big dust in the cavity clearing process. Our model also suggests an inwardly increasing gas-to-dust-ratio in the inner disk, and different spatial distributions for the small dust inside and outside the cavity, echoing the predictions in grain coagulation and growth models.

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.

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.

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.

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).

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.

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.

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.