We present high resolution millimeter continuum ALMA observations of the disks around the T Tauri stars LkCa 15 and J1610. These disks host dust-depleted inner regions, possibly carved by massive planets, and are of prime interest to study the imprints of planet-disk interactions. While at moderate angular resolution they appear as a broad ring surrounding a cavity, the continuum emission resolves into multiple rings at a resolution of ~60$\times$40 mas (~7.5 au for LkCa 15, ~6 au for J1610) and ~$7\,\mu$Jy beam$^{-1}$ rms at 1.3 mm. In addition to a broad extended component, LkCa 15 and J1610 host 3 and 2 narrow rings, respectively, with two bright rings in LkCa 15 being radially resolved. The rings look marginally optically thick, with peak optical depths of ~0.5 (neglecting scattering), in agreement with high angular resolution observations of full disks. We perform hydrodynamical simulations with an embedded, sub-Jovian-mass planet and show that the observed multi-ringed substructure can be qualitatively explained as the outcome of the planet-disk interaction. We note however that the choice of the disk cooling timescale alone can significantly impact the resulting gas and dust distributions around the planet, leading to different numbers of rings and gaps and different spacings between them. We propose that the massive outer disk regions of transition disks are favorable places for planetesimals and possibly second generation planet formation of objects with a lower mass than the planets carving the inner cavity (typically few $M_{\rm Jup}$), and that the annular substructures observed in LkCa 15 and J1610 may be indicative of planetary core formation within dust-rich pressure traps. Current observations are compatible with other mechanisms being at the origin of the observed substructures, in particular with narrow rings generated at the edge of the CO and N$_2$ snowlines.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

We present recent results obtained with the VLT/MUSE Integral Field Spectrograph fed by the 4LGSF and its laser tomography adaptive optics module GALACSI. While this so-called narrow-field mode of MUSE was not designed to perform directly imaging of exoplanets and outflows, we show that it can be a game changer to detect and characterize young exoplanets with a prominent emission lines (i.e H{\alpha}, tracer of accretion), at moderate contrasts. These performances are achieved thanks to the combo of a near-diffraction limited PSF and a medium resolution spectrograph and a cross-correlation approach in post-processing . We discuss this in the context of ground and space, infrared and visible wavelengths, preparing for missions like JWST and WFIRST in great synergy and as pathfinder for future ELT/GSMT (Extremely Large and/or Giant Segmented Mirror Telescopes) instruments.

The PDS 70 system has been subject to many studies in the past year following the discovery of two accreting planets in the gap of its circumstellar disk. Nevertheless, the mass accretion rate onto the star is still not well known. Here we determined the stellar mass accretion rate and its variability based on TESS and HARPS observations. The stellar light curve shows a strong signal with a $3.03\pm0.06$ days period, which we attribute to stellar rotation. Our analysis of the HARPS spectra shows a rotational velocity of $v\sin\,i=16.0\pm0.5\,{\rm km\,s^{-1}}$, indicating that the inclination of the rotation axis is $50\pm8$ degrees. This implies that the rotation axes of the star and its circumstellar disk are parallel within the measurement error. We apply magnetospheric accretion models to fit the profiles of the H$\alpha$ line and derive mass accretion rates onto the star in the range of $0.6-2.2\times10^{-10}\,{\rm M_{\odot}yr^{-1}}$, varying over the rotation phase. The measured accretion rates are in agreement with those estimated from NUV fluxes using accretion shock models. The derived accretion rates are higher than expected from the disk mass and planets properties for the low values of the viscous parameter $\alpha$ suggested by recent studies, potentially pointing to an additional mass reservoir in the inner disk to feed the accretion, such as a dead zone. We find that the He I $\lambda$10830 line shows a blueshifted absorption feature, indicative of a wind. The mass-loss rate estimated from the line depth is consistent with an accretion-driven inner disk MHD wind.

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.

We are undertaking a large survey of over thirty disks using the Gemini Planet Imager (GPI) to see whether the observed dust structures match spectral energy distribution (SED) predictions and have any correlation with stellar properties. GPI can observe near-infrared light scattered from dust in circumstellar environments using high-resolution Polarimetric Differential Imaging (PDI) with coronagraphy and adaptive optics. The data have been taken in J and H bands over two years, with inner working angles of 0.08'' and 0.11'' respectively. Ahead of the release of the complete survey results, here we present five objects with extended and irregular dust structures within 2'' of the central star. These objects are: FU Ori; MWC 789; HD 45677; Hen 3-365; and HD 139614. The observed structures are consistent with each object being a pre-main-sequence star with protoplanetary dust. The five objects' circumstellar environments could result from extreme youth and complex initial conditions, from asymmetric scattering patterns due to shadows cast by misaligned disks, or in some cases from interactions with companions. We see complex U_phi structures in most objects that could indicate multiple scattering or result from the illumination of companions. Specific key findings include the first high-contrast observation of MWC 789 revealing a newly-discovered companion candidate and arc, and two faint companion candidates around Hen 3-365. These two objects should be observed further to confirm whether the companion candidates are co-moving. Further observations and modeling are required to determine the causes of the structures.

Context. PDS 70 is a young (5.4 Myr), nearby (~113 pc) star hosting a known transition disk with a large gap. Recent observations with SPHERE and NACO in the near-infrared (NIR) allowed us to detect a planetary mass companion, PDS70b, within the disk cavity. Moreover, observations in H_alpha with MagAO and MUSE revealed emission associated to PDS70b and to another new companion candidate, PDS70c, at a larger separation from the star. Aims. Our aim is to confirm the discovery of the second planet PDS70c using SPHERE at VLT, to further characterize its physical properties, and search for additional point sources in this young planetary system. Methods. We re-analyzed archival SPHERE NIR observations and obtained new data in Y, J, H and K spectral bands for a total of four different epochs. The data were reduced using the data reduction and handling pipeline and the SPHERE data center. We then applied custom routines (e.g. ANDROMEDA and PACO) to subtract the starlight. Results. We re-detect both PDS 70 b and c and confirm that PDS70c is gravitationally bound to the star. We estimate this second planet to be less massive than 5 M Jup and with a T_eff around 900 K. Also, it has a low gravity with log g between 3.0 and 3.5 dex. In addition, a third object has been identified at short separation (~0.12") from the star and gravitationally bound to the star. Its spectrum is however very blue, so that we are probably seeing stellar light reflected by dust and our analysis seems to demonstrate that it is a feature of the inner disk. We, however, cannot completely exclude the possibility that it is a planetary mass object enshrouded by a dust envelope. In this latter case, its mass should be of the order of few tens of M_Earth. Moreover, we propose a possible structure for the planetary system based on our data that, however, cannot be stable on a long timescale.

Protoplanetary disks are known to posses a stunning variety of substructure in the distribution of their mm~sized grains, predominantly seen as rings and gaps (Andrews et al. 2018), which are frequently interpreted as due to the shepherding of large grains by either hidden, still-forming planets within the disk (Zhang et al. 2018) or (magneto-)hydrodynamic instabilities (Flock et al. 2015). The velocity structure of the gas offers a unique probe of both the underlying mechanisms driving the evolution of the disk, the presence of embedded planets and characterising the transportation of material within the disk, such as following planet-building material from volatile-rich regions to the chemically-inert midplane, or detailing the required removal of angular momentum. Here we present the radial profiles of the three velocity components of gas in upper disk layers in the disk of HD 163296 as traced by 12CO molecular emission. These velocities reveal significant flows from the disk surface towards the midplane of disk at the radial locations of gaps argued to be opened by embedded planets (Isella et al. 2016, 2018, Teague et al. 2018, Pinte et al. 2018), bearing striking resemblance to meridional flows, long predicted to occur during the early stages of planet formation (Szulagyi et al. 2014, Morbidelli et al. 2014, Fung & Chiang 2016, Dong et al. 2019). In addition, a persistent radial outflow is seen at the outer edge of the disk, potentially the base of a wind associated with previously detected extended emission (Klaassen et al. 2013).

We report the detection of spiral substructure in both the gas velocity and temperature structure of the disk around TW~Hya, suggestive of planet-disk interactions with an unseen planet. Perturbations from Keplerian rotation tracing out a spiral pattern are observed in the SE of the disk, while significant azimuthal perturbations in the gas temperature are seen in the outer disk, outside 90~au, extending the full azimuth of the disk. The deviation in velocity is either $\Delta v_{\phi} \, / \, v_{\rm kep} \sim 0.1$ or $\Delta v_{z} \, / \, v_{\rm kep} \sim 0.01$ depending on whether the perturbation is in the rotational or vertical direction, while radial perturbations can be ruled out. Deviations in the gas temperature are $\pm 4$ K about the azimuthally averaged profile, equivalent to deviations of $\Delta T_{\rm gas} \, / \, T_{\rm gas} \sim 0.05$. Assuming all three structures can be described by an Archimedean spiral, measurements of the pitch angles of both velocity and temperature spirals show a radially decreasing trend for all three, ranging from 9$^{\circ}$ at 70 au, dropping to 3$^{\circ}$ at 200 au. Such low pitch-angled spirals are not readily explained through the wake of an embedded planet in the location of previously reported at 94 au, but rather require a launching mechanism which results in much more tightly wound spirals. Molecular emission tracing distinct heights in the disk is required to accurately distinguish between spiral launching mechanisms.

While numerical simulations have been playing a key role in the studies of planet-disk interaction, testing numerical results against observations has been limited so far. With the two directly imaged protoplanets embedded in its circumstellar disk, PDS 70 offers an ideal testbed for planet-disk interaction studies. Using two-dimensional hydrodynamic simulations we show that the observed features can be well explained with the two planets in formation, providing strong evidence that previously proposed theories of planet-disk interaction are in action, including resonant migration, particle trapping, size segregation, and filtration. Our simulations suggest that the two planets are likely in 2:1 mean motion resonance and can remain dynamically stable over million-year timescales. The growth of the planets at $10^{-8}-10^{-7}~M_{\rm Jup}~{\rm yr}^{-1}$, rates comparable to the estimates from H$\alpha$ observations, does not destabilize the resonant configuration. Large grains are filtered at the gap edge and only small, (sub-)$\mu$m grains can flow to the circumplanetary disks and the inner circumstellar disk. With the sub-millimeter continuum ring observed outward of the two directly imaged planets, PDS 70 provides the first observational evidence of particle filtration by gap-opening planets. The observed sub-millimeter continuum emission at the vicinity of the planets can be reproduced when (sub-)$\mu$m grains survive over multiple circumplanetary disk gas viscous timescales and accumulate therein. One such possibility is if (sub-)$\mu$m grains grow in size and remain trapped in pressure bumps, similar to what we find happening in circumstellar disks. We discuss potential implications to planet formation in the solar system and mature extrasolar planetary systems.

Advances in high-resolution imaging have revealed H$\alpha$ emission separated from the host star. It is generally believed that the emission is associated with forming planets in protoplanetary disks. However, the nature of this emission is still not fully understood. Here we report a modeling effort of H$\alpha$ emission from the planets around the young star PDS 70. Using standard magnetospheric accretion models previously applied to accreting young stars, we find that the observed line fluxes can be reproduced using a range of parameters relevant to PDS 70b and c, with the mean mass accretion rate of log(${\rm \dot{M}}$) = $-8.0\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ and $-8.1\pm0.6$ M$_{\rm Jup}$ yr$^{-1}$ for PDS 70b and PDS 70c, respectively. Our results suggest that H$\alpha$ emission from young planets can originate in the magnetospheric accretion of mass from the circumplanetary disk. We find that empirical relationships between mass accretion rate and H$\alpha$ line properties frequently used in T Tauri stars are not applicable in the planetary mass regime. In particular, the correlations between line flux and mass accretion rate underpredict the accretion rate by about an order of magnitude, and the width at the 10% height of the line is insensitive to the accretion rate at ${\rm \dot{M}}$ $< 10^{-8}$ M$_{\rm Jup}$ yr$^{-1}$.

Gravitational interactions between a protoplanetary disk and its embedded planet is one of the formation mechanisms of gaps and rings found in recent ALMA observations. To quantify the gap properties measured in not only surface density but also rotational velocity profiles, we run two-dimensional hydrodynamic simulations of protoplanetary disks by varying three parameters: the mass ratio $q$ of a planet to a central star, the ratio of the disk scale height $h_p$ to the orbital radius $r_p$ of the planet, and the viscosity parameter $\alpha$. We find the gap depth $\delta_\Sigma$ in the gas surface density depends on a single dimensionless parameter $K\equiv q^2(h_p/r_p)^{-5}\alpha^{-1}$ as $\delta_\Sigma = (1 + 0.046K)^{-1}$, consistent with the previous results of Kanagawa et al. (2015a). The gap depth $\delta_V$ in the rotational velocity is given by $\delta_V= 0.007 (h_p/r_p) K^{1.38}/(1 +0.06K^{1.03})$. The gap width, in both surface density and rotational velocity, has a minimum of about $4.7 h_p$ when the planet mass $M_p$ is around the disk thermal mass $M_{\text{th}}$, while it increases in a power-law fashion as $M_p/M_{\text{th}}$ increases or decrease from unity. Such a minimum in the gap width arises because spirals from sub-thermal planets have to propagate before they shock the disk gas and open a gap. We compare our relations for the gap depth and width with the previous results, and discuss their applicability to observations.

We present the discovery of a spatially unresolved source of sub-millimeter continuum emission ($\lambda=855$ $\mu$m) associated with a young planet, PDS 70 c, recently detected in H$\alpha$ emission around the 5 Myr old T Tauri star PDS 70. We interpret the emission as originating from a dusty circumplanetary disk with a dust mass between $2\times10^{-3}$ and $4.2 \times 10^{-3}$ Earth masses. Assuming a standard gas-to-dust ratio of 100, the ratio between the total mass of the circumplanetary disk and the mass of the central planet would be between $10^{-4}-10^{-5}$. Furthermore, we report the discovery of another compact continuum source located $0.074''\pm0.013''$ South-West of a second known planet in this system, PDS 70 b, that was previously detected in near-infrared images. We speculate that the latter source might trace dust orbiting in proximity of the planet, but more sensitive observations are required to unveil its nature.

PDS 70b is the most robustly detected young planet imaged in the gap of a transition disk so far, found at a separation of ~195 mas (~22 au) from its host star and at a position angle of about 155 deg. This system is therefore a unique laboratory to characterize the properties of young planetary systems at the stage of their formation. We aim to trace direct and indirect imprints of PDS 70b on the gas and dust emission of the circumstellar disk in order to study the properties of this ~5 Myr young planetary system. We obtained ALMA band 7 observations of PDS 70 in dust continuum and $^{12}CO$ (3-2) and combined them with archival data resulting in an unprecedented angular resolution of about 70 mas (~8 au). We derive an upper limit on circumplanetary material at the location of PDS 70b of ~0.01 $M_{\oplus}$ and find a highly structured circumstellar disk in both dust and gas. The outer dust ring peaks at ~0.65" (74 au) and reveals a possible second unresolved peak at ~0.53" (60 au). The integrated intensity of CO also shows evidence of a depletion of emission at ~0.2" (23 au) with a width of ~0.1" (11 au). The gas kinematics show evidence of a deviation from Keplerian rotation inside $\lesssim$0.8" (91 au). This implies a pressure gradient that can account for the location of the dust ring well beyond the location of PDS 70b. Further in, we detect an inner disk which appears to be connected to the outer disk by a possible bridge feature in the North-West region in both gas and dust. We compare the observations to hydrodynamical simulations including a planet with different masses covering the estimated mass range previously derived from NIR photometry (~5-9 M$_{Jup}$). We find that even a planet with a mass of 10 M$_{Jup}$ may not be sufficient to explain the extent of the wide gap and an additional low-mass companion may be needed to account for the observed disk morphology.

Understanding the origin of the astonishing diversity of exoplanets is a key question for the coming decades. ALMA has revolutionized our view of the dust emission from protoplanetary disks, demonstrating the prevalence of ring and spiral structures that are likely sculpted by young planets in formation. To detect kinematic signatures of these protoplanets and to probe the chemistry of their gas accretion reservoir will require the imaging of molecular spectral line emission at high angular and spectral resolution. However, the current sensitivity of ALMA limits these important spectral studies to only the nearest protoplanetary disks. Although some promising results are emerging, including the identification of the snowlines of a few key molecules and the first attempt at detecting a protoplanet's spiral wake, it is not yet possible to search for these important signatures in a population of disks in diverse environments and ages. Harnessing the tremendous power of (sub)mm observations to pinpoint and characterize the chemistry of planets in formation will require a major increase of ALMA's spectral sensitivity (5-10x), increase in instantaneous bandwidth (2x) at high spectral resolution, and improved angular resolution (2x) in the 2030 era.

The processes that transform gas and dust in circumstellar disks into diverse exoplanets remain poorly understood. One key pathway is to study exoplanets as they form in their young ($\sim$few~Myr) natal disks. Extremely Large Telescopes (ELTs) such as GMT, TMT, or ELT, can be used to establish the initial chemical conditions, locations, and timescales of planet formation, via (1)~measuring the physical and chemical conditions in protoplanetary disks using infrared spectroscopy and (2)~studying planet-disk interactions using imaging and spectro-astrometry. Our current knowledge is based on a limited sample of targets, representing the brightest, most extreme cases, and thus almost certainly represents an incomplete understanding. ELTs will play a transformational role in this arena, thanks to the high spatial and spectral resolution data they will deliver. We recommend a key science program to conduct a volume-limited survey of high-resolution spectroscopy and high-contrast imaging of the nearest protoplanetary disks that would result in an unbiased, holistic picture of planet formation as it occurs.

We present the first images of the transition disk around the close binary system HD 34700A in polarized scattered light using the Gemini Planet Imager instrument on Gemini South. The J and H band images reveal multiple spiral-arm structures outside a large (R = 0.49" = 175 au) cavity along with a bluish spiral structure inside the cavity. The cavity wall shows a strong discontinuity and we clearly see significant non-azimuthal polarization Uphi consistent with multiple scattering within a disk at an inferred inclination ~42deg. Radiative transfer modeling along with a new Gaia distance suggest HD 37400A is a young (~5 Myr) system consisting of two intermediate-mass (~2Msun) stars surrounded by a transitional disk and not a solar-mass binary with a debris disk as previously classified. Conventional assumptions of the dust-to-gas ratio would rule out a gravitational instability origin to the spirals while hydrodynamical models using the known external companion or a hypothetical massive protoplanet in the cavity both have trouble reproducing the relatively large spiral arm pitch angles (~30deg) without fine tuning of gas temperature. We explore the possibility that material surrounding a massive protoplanet could explain the rim discontinuity after also considering effects of shadowing by an inner disk. Analysis of archival Hubble Space Telescope data suggests the disk is rotating counterclockwise as expected from the spiral arm structure and revealed a new low-mass companion at 6.45" separation. We include an appendix which sets out clear definitions of Q, U, Qphi, Uphi, correcting some confusion and errors in the literature.

We present an improved method to measure the rotation curves for disks with non-axisymmetric brightness profiles initially published in Teague et al. (2018a). Application of this method to the well studied AS$~$209 system shows substantial deviations from Keplerian rotation of up to $\pm 5\%$. These deviations are most likely due to perturbations in the gas pressure profile, including a perturbation located at $\approx 250~$au and spanning up to $\approx 50~$au which is only detected kinematically. Modelling the required temperature and density profiles required to recover the observed rotation curve we demonstrate that the rings observed in $\mu$m scattered light are coincident with the pressure maxima, and are radially offset from the rings observed in mm continuum emission. This suggests that if rings in the NIR are due to sub-$\mu$m grains trapped in pressure maxima that there is a vertical dependence on the radius of the pressure minima.

Combining hydrodynamic planet-disk interaction simulations with dust evolution models, we show that protoplanetary disks having a giant planet can reveal diverse morphology in (sub-)millimeter continuum, including a full disk without significant radial structure, a transition disk with an inner cavity, a disk with a single gap and a central continuum peak, and a disk with multiple rings and gaps. Such a diversity originates from (1) the level of viscous transport in the disk which determines the number of gaps a planet can open; (2) the size and spatial distributions of grains determined by the coagulation, fragmentation, and radial drift, which in turn affects the emmisivity of the disk at (sub-)millimeter wavelengths; and (3) the angular resolution used to observe the disk. In particular, our results show that disks having the same underlying gas distribution can have very different grain size/spatial distributions and thus appearance in continuum, depending on the interplay among coagulation, fragmentation, and radial drift. This suggests that proper treatments for the grain growth have to be included in models of protoplanetary disks concerning continuum properties and that complementary molecular line observations are highly desired in addition to continuum observations to reveal the true nature of disks. The fact that a single planet can produce diverse disk morphology emphasizes the need to search for more direct, localized signatures of planets in order to confirm (or dispute) the planetary origin of observed ringed substructures.

We present the first results from a reverberation-mapping campaign undertaken during the first half of 2012, with additional data on one AGN (NGC 3227) from a 2014 campaign. Our main goals are (1) to determine the black hole masses from continuum-Hbeta reverberation signatures, and (2) to look for velocity-dependent time delays that might be indicators of the gross kinematics of the broad-line region. We successfully measure Hbeta time delays and black hole masses for five AGNs, four of which have previous reverberation mass measurements. The values measured here are in agreement with earlier estimates, though there is some intrinsic scatter beyond the formal measurement errors. We observe velocity dependent Hbeta lags in each case, and find that the patterns have changed in the intervening five years for three AGNs that were also observed in 2007.

We present the first kinematical detection of embedded protoplanets within a protoplanetary disk. Using archival ALMA observations of HD 163296, we demonstrate a new technique to measure the rotation curves of CO isotopologue emission to sub-percent precision relative to the Keplerian rotation. These rotation curves betray substantial deviations caused by local perturbations in the radial pressure gradient, likely driven by gaps carved in the gas surface density by Jupiter-mass planets. Comparison with hydrodynamic simulations shows excellent agreement with the gas rotation profile when the disk surface density is perturbed by two Jupiter mass planets at 83 au and 137 au. As the rotation of the gas is dependent on the pressure of the total gas component, this method provides a unique probe of the gas surface density profile without incurring significant uncertainties due to gas-to-dust ratios or local chemical abundances which plague other methods. Future analyses combining both methods promise to provide the most accurate and robust measures of embedded planetary mass. Furthermore, this method provides a unique opportunity to explore wide-separation planets beyond the mm continuum edge and to trace the gas pressure profile essential in modelling grain evolution in disks.

We examine whether various characteristics of planet-driven spiral arms can be used to constrain the masses of unseen planets and their positions within their disks. By carrying out two-dimensional hydrodynamic simulations varying planet mass and disk gas temperature, we find that a larger number of spiral arms form with a smaller planet mass and a lower disk temperature. A planet excites two or more spiral arms interior to its orbit for a range of disk temperature characterized by the disk aspect ratio $0.04\leq(h/r)_p\leq0.15$, whereas exterior to a planet's orbit multiple spiral arms can form only in cold disks with $(h/r)_p \lesssim 0.06$. Constraining the planet mass with the pitch angle of spiral arms requires accurate disk temperature measurements that might be challenging even with ALMA. However, the property that the pitch angle of planet-driven spiral arms decreases away from the planet can be a powerful diagnostic to determine whether the planet is located interior or exterior to the observed spirals. The arm-to-arm separations increase as a function of planet mass, consistent with previous studies; however, the exact slope depends on disk temperature as well as the radial location where the arm-to-arm separations are measured. We apply these diagnostics to the spiral arms seen in MWC 758 and Elias 2-27. As shown in Bae et al. (2017), planet-driven spiral arms can create concentric rings and gaps, which can produce more dominant observable signature than spiral arms under certain circumstances. We discuss the observability of planet-driven spiral arms versus rings and gaps.

Protoplanetary disk simulations show that a single planet can excite more than one spiral arm, possibly explaining recent observations of multiple spiral arms in some systems. In this paper, we explain the mechanism by which a planet excites multiple spiral arms in a protoplanetary disk. Contrary to previous speculations, the formation of both primary and additional arms can be understood as a linear process when the planet mass is sufficiently small. A planet resonantly interacts with epicyclic oscillations in the disk, launching spiral wave modes around the Lindblad resonances. When a set of wave modes is in phase, they can constructively interfere with each other and create a spiral arm. More than one spiral arm can form because such constructive interference can occur for different sets of wave modes, with the exact number and launching position of spiral arms dependent on the planet mass as well as the disk temperature profile. Non-linear effects become increasingly important as the planet mass increases, resulting in spiral arms with stronger shocks and thus larger pitch angles. This is found in common for both primary and additional arms. When a planet has a sufficiently large mass ($\gtrsim$ 3 thermal masses for $(h/r)_p=0.1$), only two spiral arms form interior to its orbit. The wave modes that would form a tertiary arm for smaller mass planets merge with the primary arm. Improvements in our understanding of the formation of spiral arms can provide crucial insights into the origin of observed spiral arms in protoplanetary disks.

As spiral waves driven by a planet in a gaseous disk steepen into a shock, they deposit angular momentum, opening a gap in the disk. This has been well studied using both linear theory and numerical simulations, but so far, only for the primary spiral arm -- the one directly attached to the planet. Using two-dimensional hydrodynamic simulations, we show that the secondary and tertiary arms driven by a planet can also open gaps as they steepen into shocks. The depths of the secondary/tertiary gaps in surface density grow with time in a low viscosity disk ($\alpha = 5\times10^{-5}$), so even low-mass planets (e.g., super-Earth or mini-Neptune-mass) embedded in the disk can open multiple observable gaps, provided that sufficient time has passed. Applying our results to the HL Tau disk, we show that a single 30 Earth-mass planet embedded in the ring at 68.8~au (B5) can reasonably well reproduce the positions of the two major gaps at 13.2 and 32.3 au (D1 and D2), and roughly reproduce two other major gaps at 64.2 and 74.7 au (D5 and D6) seen in the mm continuum. The positions of secondary/tertiary gaps are found to be sensitive to the planetary mass and the disk temperature profile, so with accurate observational measurements of the temperature structure the positions of multiple gaps can be used to constrain the mass of the planet. We also comment on the gaps seen in the TW Hya and HD~163296 disk.

We conjecture that observed protoplanetary disc accretion rates may be explained with low viscosities which could be the result of hydrodynamic turbulence. We show that viscosities parameterized in the usual way with $\alpha \gtrsim 10^{-4}$, comparable to values suggested for hydrodynamic turbulence, can explain the observed accretion rates and lifetimes with plausible inner disc surface densities. Our models are also in better agreement with surface density estimates of the minimum mass solar nebula than models with rapid transport for a given mass accretion rate, such as recent models of accretion driven by magnetic winds. The required surface densities are a natural result of the protostellar infall phase, as long as non-gravitational transport is limited. We argue that, in addition to possible non-ideal magnetic transport due to disc winds possibly modified by the Hall effect, the effects of low-viscosity hydrodynamic accretion deserve more consideration.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

We have recently shown that spiral density waves propagating in accretion disks can undergo a parametric instability by resonantly coupling with and transferring energy into pairs of inertial waves (or inertial-gravity waves when buoyancy is important). In this paper, we perform inviscid three-dimensional global hydrodynamic simulations to examine the growth and consequence of this instability operating on the spiral waves driven by a Jupiter-mass planet in a protoplanetary disk. We find that the spiral waves are destabilized via the spiral wave instability (SWI), generating hydrodynamic turbulence and sustained radially-alternating vertical flows that appear to be associated with long wavelength inertial modes. In the interval $0.3~R_p \leq R \leq 0.7~R_p$, where $R_p$ denotes the semi-major axis of the planetary orbit (assumed to be 5~au), the estimated vertical diffusion rate associated with the turbulence is characterized by $\alpha_{\rm diff} \sim (0.2-1.2) \times 10^{-2}$. For the disk model considered here, the diffusion rate is such that particles with sizes up to several centimeters are vertically mixed within the first pressure scale height. This suggests that the instability of spiral waves launched by a giant planet can significantly disperse solid particles and trace chemical species from the midplane. In planet formation models where the continuous local production of chondrules/pebbles occurs over Myr time scales to provide a feedstock for pebble accretion onto these bodies, this stirring of solid particles may add a time constraint: planetary embryos and large asteroids have to form before a gas giant forms in the outer disk, otherwise the SWI will significantly decrease the chondrule/pebble accretion efficiency.

We present results from a suite of three-dimensional global hydrodynamic simulations which show that spiral density waves propagating in circumstellar disks are unstable to the growth of a parametric instability that leads to break-down of the flow into turbulence. This spiral wave instability (SWI) arises from a resonant interaction between pairs of inertial waves, or inertial-gravity waves, and the background spiral wave. The development of the instability in the linear regime involves the growth of a broad spectrum of inertial modes, with growth rates on the order of the orbital time, and results in a nonlinear saturated state in which turbulent velocity perturbations are of a similar magnitude to those induced by the spiral wave. The turbulence induces angular momentum transport, and vertical mixing, at a rate that depends locally on the amplitude of the spiral wave (we obtain a stress parameter $\alpha \sim 5 \times 10^{-4}$ in our reference model). The instability is found to operate in a wide-range of disk models, including those with isothermal or adiabatic equations of state, and in viscous disks where the dimensionless kinematic viscosity $\nu\le 10^{-5}$. This robustness suggests that the instability will have applications to a broad range of astrophysical disk-related phenomena, including those in close binary systems, planets embedded in protoplanetary disks (including Jupiter in our own Solar System) and FU Orionis outburst models. Further work is required to determine the nature of the instability, and to evaluate its observational consequences, in physically more complete disk models than we have considered in this paper.

The disk surrounding SAO 206462, an 8 Myr-old Herbig Ae star, has recently been reported to exhibit spiral arms, an asymmetric dust continuum, and a dust-depleted inner cavity. By carrying out two-dimensional, two-fluid hydrodynamic calculations, we find that a planetary-mass companion located at the outer disk could be responsible for these observed structures. In this model, the planet excites primary and secondary arms interior to its orbit. It also carves a gap and generates a local pressure bump at the inner gap edge where a vortex forms through Rossby wave instability. The vortex traps radially drifting dust particles, forming a dust-depleted cavity in the inner disk. We propose that the vortex is responsible for the brightest southwestern peak seen in infrared scattered light and sub-millimeter dust continuum emission. In particular, it is possible that the scattered light is boosted as one of the spiral arms passes through the high density vortex region, although the vortex alone may be able to explain the peak. We suggest that a planetary companion with a mass of 10-15 M_J is orbiting SAO 206462 at 100-120 au. Monitoring of the brightest peak over the next few years will help reveal its origin because the spiral arms and vortex will show distinguishable displacement.

We carry out two-fluid, two-dimensional global hydrodynamic simulations to test whether protostellar infall can trigger Rossby wave instability (RWI) in protoplanetry disks. Our results show that infall can trigger the RWI and generate vortices near the outer edge of the mass landing on the disk (i.e. centrifugal radius). We find that the RWI is triggered under a variety of conditions, although the details depend on the disk parameters and the infall pattern. The common key feature of triggering the RWI is the steep radial gradient of the azimuthal velocity induced by the local increase in density at the outer edge of the infall region. Vortices form when the instability enters the nonlinear regime. In our standard model where self-gravity is neglected, vortices merge together to a single vortex within $\sim 20$ local orbital times, and the merged vortex survives for the remaining duration of the calculation ($> 170$ local orbital times). The vortex takes part in outward angular momentum transport, with a Reynolds stress of $\lesssim10^{-2}$. Our two-fluid calculations show that vortices efficiently trap dust particles with stopping times of the order of the orbital time, locally enhancing the dust to gas ratio for particles of the appropriate size by a factor of $\sim 40$ in our standard model. When self-gravity is considered, however, vortices tend to be impeded from merging and may eventually dissipate. We conclude it may well have that protoplanetary disks have favorable conditions for vortex formation during the protostellar infall phase, which might enhance early planetary core formation.

We improve on our previous treatments of long-term evolution of protostellar disks by explicitly solving disk self-gravity in two dimensions. The current model is an extension of the one-dimensional layered accretion disk model of Bae et al. We find that gravitational instability (GI)-induced spiral density waves heat disks via compressional heating (i.e. $P\rm{d}V$ work), and can trigger accretion outbursts by activating the magnetorotational instability (MRI) in the magnetically inert disk dead-zone. The GI-induced spiral waves propagate well inside of gravitationally unstable region before they trigger outbursts at $R \lesssim 1$ AU where GI cannot be sustained. This long-range propagation of waves cannot be reproduced with the previously used local $\alpha$ treatments for GI. In our standard model where zero dead-zone residual viscosity ($\alpha_{\rm rd}$) is assumed, the GI-induced stress measured at the onset of outbursts is locally as large as $0.01$ in terms of the generic $\alpha$ parameter. However, as suggested in our previous one-dimensional calculations, we confirm that the presence of a small but finite $\alpha_{\rm rd}$ triggers thermally-driven bursts of accretion instead of the GI + MRI-driven outbursts that are observed when $\alpha_{\rm rd}=0$. The inclusion of non-zero residual viscosity in the dead-zone decreases the importance of GI soon after mass feeding from the envelope cloud ceases. During the infall phase while the central protostar is still embedded, our models stay in a quiescent accretion phase with $\dot{M}_{acc}\sim10^{-8}-10^{-7}~M_{\odot}~{\rm yr^{-1}}$ over $60~\%$ of the time and spend less than $15~\%$ of the infall phase in accretion outbursts.

The Korean very-long-baseline interferometry (VLBI) network (KVN) and VLBI Exploration of Radio Astrometry (VERA) Array (KaVA) is the first international VLBI array dedicated to high-frequency (23 and 43 GHz bands) observations in East Asia. Here, we report the first imaging observations of three bright active galactic nuclei (AGNs) known for their complex morphologies: 4C 39.25, 3C 273, and M 87. This is one of the initial result of KaVA early science. Our KaVA images reveal extended outflows with complex substructure such as knots and limb brightening, in agreement with previous Very Long Baseline Array (VLBA) observations. Angular resolutions are better than 1.4 and 0.8 milliarcsecond at 23 GHz and 43 GHz, respectively. KaVA achieves a high dynamic range of ~1000, more than three times the value achieved by VERA. We conclude that KaVA is a powerful array with a great potential for the study of AGN outflows, at least comparable to the best existing radio interferometric arrays.

We have carried out the first very long baseline interferometry (VLBI) imaging of 44 GHz class I methanol maser (7_{0}-6_{1}A^{+}) associated with a millimeter core MM2 in a massive star-forming region IRAS 18151-1208 with KaVA (KVN and VERA Array), which is a newly combined array of KVN (Korean VLBI Network) and VERA (VLBI Exploration of Radio Astrometry). We have succeeded in imaging compact maser features with a synthesized beam size of 2.7 milliarcseconds x 1.5 milliarcseconds (mas). These features are detected at a limited number of baselines within the length of shorter than approximately 650 km corresponding to 100 Mlambda in the uv-coverage. The central velocity and the velocity width of the 44 GHz methanol maser are consistent with those of the quiescent gas rather than the outflow traced by the SiO thermal line. The minimum component size among the maser features is ~ 5 mas x 2 mas, which corresponds to the linear size of ~ 15 AU x 6 AU assuming a distance of 3 kpc. The brightness temperatures of these features range from ~ 3.5 x 10^{8} to 1.0 x 10^{10} K, which are higher than estimated lower limit from a previous Very Large Array observation with the highest spatial resolution of ~ 50 mas. The 44 GHz class I methanol maser in IRAS 18151-1208 is found to be associated with the MM2 core, which is thought to be less evolved than another millimeter core MM1 associated with the 6.7 GHz class II methanol maser.

We perform calculations of our one-dimensional, two-zone disk model to study the long-term evolution of the circumstellar disk. In particular, we adopt published photoevaporation prescriptions and examine whether the photoevaporative loss alone, coupled with a range of initial angular momenta of the protostellar cloud, can explain the observed decline of the frequency of optically-thick dusty disks with increasing age. In the parameter space we explore, disks have accreting and/or non-accreting transitional phases lasting of $\lesssim20 %$ of their lifetime, which is in reasonable agreement with observed statistics. Assuming that photoevaporation controls disk clearing, we find that initial angular momentum distribution of clouds needs to be weighted in favor of slowly rotating protostellar cloud cores. Again, assuming inner disk dispersal by photoevaporation, we conjecture that this skewed angular momentum distribution is a result of fragmentation into binary or multiple stellar systems in rapidly-rotating cores. Accreting and non-accreting transitional disks show different evolutionary paths on the $\dot{M}-R_{\rm wall}$ plane, which possibly explains the different observed properties between the two populations. However, we further find that scaling the photoevaporation rates downward by a factor of 10 makes it difficult to clear the disks on the observed timescales, showing that the precise value of the photoevaporative loss is crucial to setting the clearing times. While our results apply only to pure photoevaporative loss (plus disk accretion), there may be implications for models in which planets clear disks preferentially at radii of order 10 AU.

We extend the one-dimensional, two-zone models of long-term protostellar disk evolution with infall of Zhu et al. to consider the potential effects of a finite viscosity in regions where the ionization is too low for the magnetorotational instability (MRI) to operate (the "dead zone"). We find that the presence of a small but finite dead zone viscosity, as suggested by simulations of stratified disks with MRI-active outer layers, can trigger inside-out bursts of accretion, starting at or near the inner edge of the disk, instead of the previously-found outside-in bursts with zero dead zone viscosity, which originate at a few AU in radius. These inside-out bursts of accretion bear a qualitative resemblance to the outburst behavior of one FU Ori object, V1515 Cyg, in contrast to the outside-in burst models which more closely resemble the accretion events in FU Ori and V1057 Cyg. Our results suggest that the type and frequency of outbursts are potentially a probe of transport efficiency in the dead zone. Simulations must treat the inner disk regions, $R \lesssim 0.5$ AU, to show the detailed time evolution of accretion outbursts in general and to observe the inside-out bursts in particular.

We report simultaneous multi-frequency observing performance at 22 and 43 GHz of the 21-m shaped-Cassegrain radio telescopes of the Korean VLBI Network (KVN). KVN is the first millimeter-dedicated VLBI network in Korea having a maximum baseline length of 480 km. It currently operates at 22 and 43 GHz and planed to operate in four frequency bands, 22, 43, 86, and 129 GHz. The unique quasioptics of KVN enable simultaneous multi-frequency observations based on efficient beam filtering and accuarate antenna-beam alignment at 22 and 43 GHz. We found that the offset of the beams is within <5 arcseconds over all pointing directions of antenna. The dual polarization, cooled HEMT receivers at 22 and 43 GHz result in receiver noise temperatures less than 40 K at 21.25-23.25 GHz and 80 K at 42.11-44.11 GHz. The pointing accuracies have been measured to be 3 arcseconds in azimuth and elevation for all antennas. The measured aperture efficiencies are 65%(K)/67%(Q), 62%(K)/59%(Q), and 66%(K)/60%(Q) for the three KVN antennas, KVNYS, KVNUS, and KVNTN, respectively. The main-beam efficiencies are measured to be 50%(K)/52%(Q), 48%(K)/50%(Q), and 50%(K)/47%(Q) for KVNYS, KVNUS, and KVNTN, respectively. The estimated Moon efficiencies are 77%(K)/90%(Q), 74%(K)/79%(Q), and 80%(K)/86%(Q) for KVNYS, KVNUS, KVNTN, respectively. The elevation dependence of the aperture efficiencies is quite flat for elevations > 20 degrees.

We report a multi-epoch, simultaneous 22 GHz H2O and 44 GHz class I CH3OH maser line survey towards 180 intermediate-mass young stellar objects, including 14 Class 0, 19 Class I objects, and 147 Herbig Ae/Be stars. We detected H2O and CH3OH maser emission towards 16 (9%) and 10 (6%) sources with one new H2O and six new CH3OH maser sources. The detection rates of both masers rapidly decrease as the central (proto)stars evolve, which is contrary to the trends in high-mass star-forming regions. This suggests that the excitations of the two masers are closely related to the evolutionary stage of the central (proto)stars and the circumstellar environments. H2O maser velocities deviate on average 9 km s^-1 from the ambient gas velocities whereas CH3OH maser velocities match quite well with the ambient gas velocities. For both maser emissions, large velocity differences (|v_{H2O} - v_{sys} | > 10 km s^-1 and |v_{CH3OH} - v_{sys}| > 1 km s^-1) are mostly confined to Class 0 objects. The formation and disappearance of H2O masers is frequent and their integrated intensities change by up to two orders of magnitude. In contrast, CH3OH maser lines usually show no significant change in intensity, shape, or velocity. This is consistent with the previous suggestion that H2O maser emission originates from the base of an outflow while 44 GHz class I CH3OH maser emission arises from the interaction region of the outflow with the ambient gas. The isotropic maser luminosities are well correlated with the bolometric luminosities of the central objects. The fitted relations are L_{H2O} = 1.71 * 10^{-9} (L_{bol})^{0.97} and L_{CH3OH} = 1.71 * 10^{-10} (L_{bol})^{1.22}.