<figure><img src="image_1.png" /></figure>

HAL
open science

# Kinematic Evidence for an Embedded Protoplanet in a Circumstellar Disk

C. Pinte, D. J. Price, F. Ménard, G. Duchêne, W. R. F. Dent, T. Hill, I. de Gregorio-Monsalvo, A. Hales, D. Mentiplay

► To cite this version:

C. Pinte, D. J. Price, F. Ménard, G. Duchêne, W. R. F. Dent, et al. Kinematic Evidence for an Embedded Protoplanet in a Circumstellar Disk. The Astrophysical Journal Letters, 2018, 860, 10.3847/2041-8213/aac6dc. insu-03693554

HAL Id: insu-03693554

https://insu.hal.science/insu-03693554v1

Submitted on 1 Sep 2025

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_2.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_3.png" /></figure>

# Kinematic Evidence for an Embedded Protoplanet in a Circumstellar Disk

C. Pinte¹,², D. J. Price¹, F. Ménard², G. Duchêne²,³, W. R. F. Dent⁴,

T. Hill⁴, I. de Gregorio-Monsalvo⁴, A. Hales⁴,⁵, and D. Ment play¹,⁰

¹ Monash Centre for Astrophysics (MoCA) and School of Physics and Astronomy, Monash University, Clayton Vic 3800, Australia; christophe.pinte@monash.edu

² Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France

³ Astronomy Department, University of California, Berkeley, CA 94720-3411, USA

⁴ Atacama Large Millimeter/Submillimeter Array, Joint ALMA Observatory, Alonso de Córdova 3107, Vitacura 763-0355, Santiago, Chile

⁵ National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903-2475, USA

Received 2018 April 22; revised 2018 May 18; accepted 2018 May 19; published 2018 June 13

## Abstract

Disks of gas and dust surrounding young stars are the birthplace of planets. However, the direct detection of protoplanets forming within disks has proved elusive to date. We present the detection of a large, localized deviation from Keplerian velocity in the protoplanetary disk surrounding the young star HD 163296. The observed velocity pattern is consistent with the dynamical effect of a two-Jupiter-mass planet orbiting at a radius ≈260 au from the star.

**Key words:** protoplanetary disks – radiative transfer – stars: formation – stars: individual (HD 163296) – submillimeter: planetary systems – techniques: interferometric

## 1. Introduction

Direct observations of forming planets in protoplanetary disks is the ultimate goal of disk studies. The disk usually outshines the planet, requiring observations at high contrast and angular resolution. Detections by direct imaging have been reported in several disks: HD 100546 (Quanz et al. 2013a, 2015; Brittain et al. 2014; Currie et al. 2015), LkCa 15 (Kraus & Ireland 2012; Sallum et al. 2015), HD 169142 (Quanz et al. 2013b; Biller et al. 2014; Reggiani et al. 2014), and MWC 758 (Reggiani et al. 2018). Yet, most of the detections to date have been subsequently challenged (e.g., Thalmann et al. 2015, 2016; Rameau et al. 2017; Ligi et al. 2018). The quest continues.

An alternative approach is to search for indirect signatures imprinted by planets on their host disk. The Atacama Large Millimeter Array (ALMA) and adaptive optics systems have revealed a variety of structures: gaps and rings (ALMA Partnership et al. 2015; Andrews et al. 2016; Isella et al. 2016) and spirals (e.g., Benisty et al. 2015; Stolker et al. 2016) that could be signposts of planets, but numerous other explanations also exist that do not require planets (e.g., Takahashi & Inutsuka 2014; Flock et al. 2015; Gonzalez et al. 2015; Lorén-Aguilar & Bate 2015; Zhang et al. 2015; Béthune et al. 2016). Embedded planets in circumstellar disks will launch spiral waves at Lindblad resonances both inside and outside of their orbit (e.g., Ogilvie & Lubow 2002), disturbing the local Keplerian velocity pattern. Hydrodynamic simulations show that the impact on the velocity pattern should be detectable by high spectral resolution ALMA line observations (Perez et al. 2015). Deviations from Keplerian rotation have been detected around circumbinary disks, with streamers at near free-fall velocities (Casassus et al. 2015; Price et al. 2018) and radial flows or warps (Walsh et al. 2017).

HD 163296 is a ~4.4 Myr old Herbig Ae star located at a distance of 101.5 ± 1.2 pc from the Sun (Gaia Collaboration et al. 2018). We rescaled all relevant quantities from previous papers based on the new Gaia distance. HD 163296 has a mass of 1.9 $M_\odot$ (e.g., Flaherty et al. 2015), a luminosity of 25 $L_\odot$ (Natta et al. 2004), and an A1Ve spectral type, with effective temperature 9300 K. Observations with the Hubble Space Telescope (HST) revealed a disk in scattered light that extends as far out as 375 au (Grady et al. 2000). Interestingly, Grady et al. (2000) inferred the

presence of a giant planet at ≈270 au based on the gap observed in scattered light at that radius. de Gregorio-Monsalvo et al. (2013) presented ALMA data and showed that the gaseous component of the disk extends to distances of at least $R_{\text{out-CO}} = 415$ au in CO while the continuum is detected only to $R_{\text{out-Dust}} = 200$ au. Higher resolution ALMA imaging revealed a bright inner disk component within the inner 0.′′5, and a spectacular series of three rings at ≈65 au, 100 au, with a fainter ring at 160 au (Isella et al. 2016).

In this Letter, we present the detection of a local deviation from the Keplerian velocity pattern found in high spectral resolution ALMA imaging. By comparing with models we find this to be consistent with the presence of a few-Jupiter-mass protoplanet in the disk.

## 2. Observations and Data Reduction

We use archival ALMA data. Observations were performed on 2012 June 9, 11, 22, and July 6 at Band 7 (2011.0.000010. SV), and on 2015 August 5, 8, and 9 at Band 6 (2013.1.00601. S). A complete description of the data was presented in de Gregorio-Monsalvo et al. (2013) and Isella et al. (2016). For the Band 7 data, we re-used the maps produced by de Gregorio-Monsalvo et al. (2013), with a 0.′′52 × 0.′′38 beam at PA = 82°, and a channel width of 0.11 km s⁻¹.

We used CASA scripts provided by ALMA to calibrate the Band 6 data. As the data from the night of August 9 showed significantly higher noise and flux levels, we selected only the data from August 5 and 8 for the analysis. We performed three successive rounds of phase self-calibration, the last with solutions calculated for each individual integration (6 s), followed by a phase and amplitude self-calibration. The continuum self-calibration solutions were applied to the CO lines. Imaging was performed at 0.1 km s⁻¹ resolution, using Briggs weighting with a robust parameter of −0.5 to obtain a synthesized beam of 0.′′28 × 0.′′18 at PA = −88°. We did not subtract the continuum emission in order to avoid underestimating the gas temperature and affecting the apparent morphology of the emission (e.g., Weaver et al. 2018). At the location of the detected velocity deviation, continuum emission is negligible, and an analysis on continuum-subtracted data would lead to the same results.
<figure><img src="image_4.png" /><figcaption>**Figure 1.** Kinematic asymmetry in HD 163296. Band 6 continuum emission (top left) and channel map of $^{12}$CO line emission at +1 km s$^{-1}$ from the systemic velocity (top right, with a close-up shown in the bottom right) shows a distinct “kink” in the emission (highlighted by the dotted circle). Comparison with the continuum emission (top left) locates this outside of the outermost dust ring. The corresponding emission on the opposite side of the disk (bottom left; showing −1 km s$^{-1}$ channel) shows no corresponding feature, indicating the disturbance to the flow is localized in both radius and azimuth. The channel width is $\Delta v = 0.1$ km s$^{-1}$. The white contour shows the 5-$\sigma$ ($\sigma = 0.1$ mJy beam$^{-1}$) level of the continuum map. The dashed line is the expected location of the isovelocity curve on the upper surface of a disk with an opening angle of 15$^{\circ}$ and an inclination of 45$^{\circ}$. Dotted lines in the bottom-right figure indicate 15% deviations (≈0.4 km s$^{-1}$) from Keplerian flow around the star. The potential planet location is marked by a cyan dot, assuming it is located in the midplane.</figcaption></figure>

### 3. Results and Analysis

The disk shows the typical butterfly pattern of disks in Keplerian rotation (de Gregorio-Monsalvo et al. 2013; Rosenfeld et al. 2013). In a given channel, the emission is concentrated along an isovelocity curve, corresponding to the region of the disk where the projected velocity is equal to the channel velocity. The emission from the upper and lower surfaces—above and

below the midplane as seen by the observer—and from the near and far sides of these surfaces, is well separated (Figure 1, and schematic view in Figure 3).

In a recent paper (Pinte et al. 2018) we showed how to reconstruct the position and velocity of each of the CO layers, for disks at intermediate inclination, by simple geometrical arguments based on the emission in each channel map. HD 163296 displays
<figure><img src="image_5.png" /><figcaption>**Figure 2.** Channel maps around the detected deviation from Keplerian velocity. The “kink” is most visible in channels at velocities between 0.8 and 1.2 km s⁻¹ (top row) and is also seen in the *J* = 3–2 transition in similar velocity channels (bottom row) indicating it is localized in both space and velocity.</figcaption></figure>

a similar scale height and velocity profile to the T Tauri star IM Lupi (Pinte et al. 2018), with a flared CO emitting surface and decreasing velocities and temperature with radius (C. Pinte et al. 2018, in preparation).

Significantly, HD 163296 shows an asymmetry between the southeast and northwest sides of the disk at a cylindrical radius of ≈260 au, outside of the third dust ring seen in continuum emission. This asymmetry is most evident in channels at a projected velocity of ≈6.8 ± 0.2 km s⁻¹ (≈1 km s⁻¹ from the systemic velocity). Figure 1 shows the corresponding individual velocity channels. The emission feature—highlighted by the dotted circle—corresponds to a kink in the upper-surface isovelocity curve northwest of the central object at velocities close to *dv* = +1 km s⁻¹. The symmetric channel (*dv* = −1 km s⁻¹) shows a smooth Keplerian profile to the southeast. We detect a similar deformation of the isovelocity curves at the same location in both ¹²CO *J* = 2–1 and 3–2 transitions (Figure 2). While it is not as obvious in the Band 7 Early-Science data due to the limited spatial resolution, the deformation of the isovelocity curve is present and could already be seen, with the benefit of hindsight from our Band 6 detection, in de Gregorio-Monsalvo et al. (2013) and Rosenfeld et al. (2013; their Figures 3 and 2, respectively). The asymmetry is not detectable in the less abundant isotologues ¹³CO and C¹⁸O, where the emission is more diffuse and fainter because of the lower optical depth.

The deformation of the emission is localized to an area approximately 0.5 in size (indicated by the dotted circle in Figures 1 and 2) and to channel maps at velocities between 0.8 and 1.2 km s⁻¹ from the systemic velocity (top row of Figure 2). This argues for a localized perturbation and excludes an origin from any large-scale structure in the disk.

## 4. Models and Discussion

The detected asymmetry matches our expectations for a local deviation from Keplerian velocity caused by a massive body embedded in the disk. A local deviation of ≈0.4 km s⁻¹ is enough to reproduce the observed spatial shift (Figure 1, bottom-right panel). The dotted lines shown in the bottom-right

<figure><img src="image_6.png" /><figcaption>**Figure 3.** Geometry of the inclined and flared disk, showing a schematic of the expected emission from two infinitely thin emitting surfaces. Green shows the emission from the lower surface of the disk, and red shows the upper surface. We added a 10% deviation in azimuthal velocity north of the star, which appears as a “kink” in the line emission. Emission is only seen when the projected velocity matches the channel velocity, producing the characteristic “butterfly” shape. Emission is preferentially seen on the upper surface of the disk due to the higher inclination with respect to the line of sight.</figcaption></figure>

panel delineate what would be ≈15% deviations in the local velocity field, which is the approximate extent of the deviation from Keplerian rotation. Most significantly, the shape of the deviation in the emission maps is similar to the prediction by Perez et al. (2015) for the kinematic signatures of an embedded planet, where the wake of the spiral generated by the planet was shown to produce a kink in the emission due to the deviation from the Keplerian rotation around the central star.

The basic feature of the channel maps can be explained with a simple model assuming emission from two infinitely thin emitting surfaces. Figure 3 shows the expected emission arising from such a model, showing the butterfly signature from the
<figure><img src="image_7.png" /><figcaption>**Figure 4.** Left panel: surface density in 3D hydrodynamics simulations of the HD 163296 disk, shown after 35 orbits of a 2 $M_{\text{jup}}$ planet and viewed at a face-on inclination. Dots mark the star and planet. Right panel: deviation of the azimuthal velocity from Keplerian velocity.</figcaption></figure>

<figure><img src="image_8.png" /></figure>

disk. Asymmetries of the velocity field, added in an ad hoc manner in the model for illustrative purposes, are evident as small bumps on the line emissions.

To go beyond this simple model and infer the mass of the putative planet, we performed a series of 3D global simulations using the PHANTOM Smoothed Particle Hydrodynamics (SPH) code (Price et al. 2017). We adopted the gas disk parameters from de Gregorio-Monsalvo et al. (2013). We employed gas-only simulations, ignoring the effect of dust, using 1 million SPH particles and a central mass of $1.9 M_{\odot}$. The inner radius of the disk in our model was set to 50 au (mainly to speed up the calculations as the inner disk is irrelevant for our present purpose), with an initial outer radius set to 500 au. We set the gas mass between those radii to $10^{-2} M_{\odot}$, and use an exponentially tapered power-law surface density profile with a critical radius of 100 au, power-law index of $p = -1.0$, and an exponent $\gamma = 0.8$. The disk aspect ratio was set to 0.08 at 50 au, with a vertically isothermal profile. We set the artificial viscosity in the code in order to obtain an average Shakura & Sunyaev (1973) viscosity of $10^{-3}$ (Lodato & Price 2010), in agreement with the upper limits found by Flaherty et al. (2015, 2017).

We embedded a single planet in the disk orbiting at 260 au with a mass of either 1, 2, 3, or 5 $M_{\text{jup}}$. We used sink particles (Bate et al. 1995) to represent the star and planet. We set the accretion radius of the planet to half of the Hill radius (7.05, 8.85, 10.15, and 12 au, respectively), with an accretion radius of 10 au for the central star. The model surface density is plotted in Figure 4 for the 2 $M_{\text{jup}}$ planet. We evolved the models for 35 orbits of the planet (≈100,000 years), which is sufficient to establish the flow pattern around the planet.

To compute the temperature and synthetic line maps, we used the MCFOST Monte Carlo radiative transfer code (Pinte et al. 2006, 2009), assuming $T_{\text{gas}} = T_{\text{dust}}$, and local thermodynamic equilibrium as we are looking at low-$J$ CO lines. The central star was represented by a sphere of radius $2.1 R_{\odot}$, radiating isotropically with a Kurucz spectrum at 9250 K. We used a Voronoi tesselation where each cell corresponds to an SPH particle, avoiding the need to interpolate between the SPH and radiative transfer codes. We set the CO abundance following the prescription in Appendix B of Pinte et al. (2018) to account for freeze-out where $T < 20$ K, as well as photo-dissociation and photo-desorption in locations where the UV radiation is high. We adopted a turbulent velocity of 50 m s$^{-1}$,

consistent with the upper limits found by Flaherty et al. (2015, 2017). We assumed a population of astrosilicate (Draine 2003) grains with sizes ranging from 0.03 to 1000 µm and following a power-law $dn(a) \propto a^{-3.5}da$, a gas-to-dust ratio of 100, and computed the dust optical properties using Mie theory.

Figure 5 presents the predicted emission in $^{12}$CO $J = 2-1$ of our theoretical models for four different planet masses. A 2 $M_{\text{jup}}$ planet appears to reproduce a deformation of the $^{12}$CO isovelocity curve that is consistent with the observations. At 1 $M_{\text{jup}}$, the planet only produces a small deformation that is barely visible in the channel maps, while a more massive planet triggers a strong spiral arm that would have been detected in channel maps at least up to $0.5 \text{ km s}^{-1}$ from the nominal velocity of 1 km s$^{-1}$. The twisted emission in the channel maps is a direct consequence of deviation from Keplerian velocity generated by the planet along the wake of the spiral arms (Figure 4, right panel). Perez et al. (2015) also predicted that the circumplanetary disk can be detected as a compact emission separated in velocity from the circumstellar disk emission. The circumplanetary disk radius is about one-third of the Hill radius (e.g., Ayliffe & Bate 2009). A 2–5 Jupiter mass planet would produce a circumplanetary disk with a diameter smaller than 6–8 au, respectively. At the current spatial resolution of the ALMA observations, its flux will be diluted in the beam (≈20 au).

Note that for the adopted disk parameters, the planet migrates by about 30 au during the simulation, and the synthetic maps display the velocity deviation slightly closer to the star than in the data. At this rate, the planet would reach the star in about 1 Myr (though we overestimate the migration rate by a factor of 2–3 due to the relatively large sink particle radius that we adopted; see Ayliffe & Bate 2010). If the detection is confirmed, the survival of such an embedded planet could put additional constraints on the disk surface density profile and viscosity.

Grady et al. (2000) detected a gap in the scattered light images with HST/STIS at 260 au, and estimated the mass of a potential planet to be 0.4 $M_{\text{jup}}$ (based on some simple analytical derivation). Isella et al. (2016) also detected a small dip in the integrated CO brightness profile at $\approx 2''/2$ (see their Figure 1, or Figure 5 in Liu et al. 2018). In our model, the gap appears in scattered light for a planet mass larger than 2 $M_{\text{jup}}$, but remains undetected in the synthetic CO maps. The final profile of a planetary gap establishes itself on a viscous timescale
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Comparison of $^{12}$CO J = 2–1 ALMA observations (top row) with synthetic channel maps from our 3D hydrodynamics calculations with embedded planets of 1, 2, 3, and 5 $M_{Jup}$ (from top to bottom). Channel width is 0.1 km s$^{-1}$. Synthetic maps were convolved to a Gaussian beam to match the spatial resolution of the observations.</figcaption></figure>

(thousands of orbits with a viscosity of $10^{-3}$), however, and the gap width and depth in our models are only lower limits.

The effect of the planet appears fainter in the $^{13}$CO channels maps than in the $^{12}$CO maps, even if the planet is located in the midplane. This is due to optical depth and vertical temperature gradient effects: the $^{12}$CO is coming from a vertically narrow and warm layer above the midplane, while the $^{13}$CO is originating from a deeper, thicker layer, where the disk is almost vertically isothermal, resulting in a uniform emissivity that washes out some of the kinematics signal.

Are we seeing the signature of an embedded planet? Can we exclude wishful thinking? The strongest evidence is that the perturbation to the disk kinematics is highly localized in both space and velocity. This excludes any mechanism that merely

produces axisymmetric rings in disks. This excludes, for example, ice lines (Zhang et al. 2015), self-induced dust traps (Gonzalez et al. 2015), instabilities (Takahashi & Inutsuka 2014; Lorén-Aguilar & Bate 2015), and zonal flows (Flock et al. 2015). A spiral wave could in principle result from the disk self-gravity, but multiple, large-scale spirals would be expected in that case (e.g., Dipierro et al. 2015), which the localized deviation seen in HD 163296 would seem to exclude.

The localized nature of the kinematic perturbation, that it occurs close to the gap found by Grady et al. (2000), and the similarity to the signatures predicted by our hydrodynamic models is strong evidence for a young protoplanet in a gas-rich disk. However, confirmation with direct imaging is the only way to be sure. The relatively large planet mass and its known
location in the disk means direct imaging follow-up might be able to detect it, depending on how embedded it is in the disk. So far, no point source has been detected at the location of the potential planet with near-IR adaptive optics systems. A 2 *M*<sub>Jup</sub> planet is consistent with the upper limits (for an unobscured planet) obtained by adaptive optics systems, such as SPHERE (Muro-Arena et al. 2018) and Keck L' (Guidi et al. submitted). Using the formalism developed in Pinte et al. (2018), we find that the velocity kink is located at a distance of ≈260 au, and an elevation above the midplane of ≈70 au. Assuming that the potential planet is located in the midplane, it would be at a projected distance of 2.3 ± 0.2 au and PA = -3° ± 5° from the star, where we estimated the uncertainties by locating the velocity deviation with half of a beam accuracy. If the planet orbit is slightly inclined compared to the disk's plane, the position on the sky will be shifted along a line going from the northeast to the southwest directions.

Can massive planets form at a distance of 250 au from the star? The location of giant planets in the outer regions of disks would be broadly consistent with gravitational instability. On the other hand, the timescale for core accretion may also be reasonable given that HD 163296 is a relatively old disk (≈5 Myr). The planet may also have undergone outward migration, depending upon on the initial profile of the disk. It is beyond the scope of this Letter to speculate further.

## 5. Summary

We detected a 15% deviation from Keplerian flow around the star in the disk around HD 163296. The deviation was detected in both Band 6 and Band 7 in two different transitions of ¹²CO and matches the kinematic signature predicted for an embedded protoplanet.

Comparing the observations to a series of 3D hydrodynamic and radiative transfer models of embedded planets suggests that the kinematic feature is caused by a planet of ≈2 *M*<sub>Jup</sub> in the midplane. Such a planet would be located at a distance of ≈260 au.

This Letter makes use of the following ALMA data: ADS/JAO.ALMA#2011.0.00010.SV and ADS/JAO.ALMA#2013.1.00601.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. This work was performed on the ozSTAR national facility at Swinburne University of Technology. ozSTAR is funded by Swinburne and the Australian Government's Education Investment Fund. We thank the anonymous referee for insightful comments and suggestions. C.P. and D.J.P. acknowledge funding from the Australian Research Council via FT170100040, FT130100034, and DP180104235. F.M. and C.P. acknowledge funding from ANR of France (ANR-16-CE31-0013).

Facility: ALMA.

Software: CASA (McMullin et al. 2007), phantom (Price et al. 2017), splash (Price 2007), mcfost (Pinte et al. 2006, 2009).

## ORCID iDs

C. Pinte https://orcid.org/0000-0001-5907-5179

D. J. Price https://orcid.org/0000-0002-4716-4235

F. Ménard https://orcid.org/0000-0002-1637-7393

G. Duchêne https://orcid.org/0000-0002-5092-6464

T. Hill https://orcid.org/0000-0003-3253-1255

I. de Gregorio-Monsalvo https://orcid.org/0000-0003-4518-407X

A. Hales https://orcid.org/0000-0001-5073-2849

D. Mentiplay https://orcid.org/0000-0002-5526-8798

## References

ALMA Partnership, Brogan, C. L., & Pérez, L. M. 2015, ApJL, 808, L3

Andrews, S. M., Wilner, D. J., Zhu, Z., et al. 2016, ApJL, 820, L40

Ayliffe, B. A., & Bate, M. R. 2009, MNRAS, 397, 657

Ayliffe, B. A., & Bate, M. R. 2010, MNRAS, 408, 876

Bate, M. R., Bonnell, I. A., & Price, N. M. 1995, MNRAS, 277, 362

Benisty, M., Juhasz, A., Boccaletti, A., et al. 2015, A&A, 578, L6

Béthune, W., Lesur, G., & Ferreira, J. 2016, A&A, 589, A87

Biller, B. A., Males, J., Rodigas, T., et al. 2014, ApJL, 792, L22

Brittain, S. D., Carr, J. S., Najita, J. R., Quanz, S. P., & Meyer, M. R. 2014, ApJ, 791, 136

Casassus, S., Marino, S., Pérez, S., et al. 2015, ApJ, 811, 92

Currie, T., Cloutier, R., Brittain, S., et al. 2015, ApJL, 814, L27

de Gregorio-Monsalvo, I., Ménard, F., Dent, W., et al. 2013, A&A, 557, A133

Dipierro, G., Pinilla, P., Lodato, G., & Testi, L. 2015, MNRAS, 451, 974

Draine, B. T. 2003, ApJ, 598, 1017

Flaherty, K. M., Hughes, A. M., Rose, S. C., et al. 2017, ApJ, 843, 150

Flaherty, K. M., Hughes, A. M., Rosenfeld, K. A., et al. 2015, ApJ, 813, 99

Flock, M., Ruge, J. P., Dzyurkevich, N., et al. 2015, A&A, 574, A68

Gaia Collaboration, Brown, A. G. A., Vallenari, A., et al. 2018, arXiv:1804.09365

Gonzalez, J.-F., Laibe, G., Maddison, S. T., Pinte, C., & Ménard, F. 2015, MNRAS, 454, L36

Grady, C. A., Devine, D., Woodgate, B., et al. 2000, ApJ, 544, 895

Isella, A., Guidi, G., Testi, L., et al. 2016, PhRvL, 117, 251101

Kraus, A. L., & Ireland, M. J. 2012, ApJ, 745, 5

Ligi, R., Vigan, A., Gratton, R., et al. 2018, MNRAS, 473, 1774

Liu, S.-F., Jin, S., Li, S., Isella, A., & Li, H. 2018, arXiv:1803.05437

Lodato, G., & Price, D. J. 2010, MNRAS, 405, 1212

Lorén-Aguilar, P., & Bate, M. R. 2015, MNRAS, 453, L78

McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007, in ASP Conf. Ser. 376, Astronomical Data Analysis Software and Systems XVI, ed. R. A. Shaw, F. Hill, & D. J. Bell (San Francisco, CA: ASP), 127

Muro-Arena, G. A., Dominik, C., Waters, L. B. F. M., et al. 2018, arXiv:1802.03328

Natta, A., Testi, L., Neri, R., Shepherd, D. S., & Wilner, D. J. 2004, A&A, 416, 179

Ogilvie, G. I., & Lubow, S. H. 2002, MNRAS, 330, 950

Perez, S., Dunhill, A., Casassus, S., et al. 2015, ApJL, 811, L5

Pinte, C., Harries, T. J., Min, M., et al. 2009, A&A, 498, 967

Pinte, C., Ménard, F., Duchêne, G., et al. 2018, A&A, 609, A47

Pinte, C., Ménard, F., Duchêne, G., & Bastien, P. 2006, A&A, 459, 797

Price, D. J. 2007, PASA, 24, 159

Price, D. J., Cuello, N., Pinte, C., et al. 2018, MNRAS, 477, 1270

Price, D. J., Wurster, J., Nixon, C., et al. 2017, arXiv:1702.03930

Quanz, S. P., Amara, A., Meyer, M. R., et al. 2013a, ApJL, 766, L1

Quanz, S. P., Amara, A., Meyer, M. R., et al. 2015, ApJ, 807, 64

Quanz, S. P., Avenhaus, H., Buenzli, E., et al. 2013b, ApJL, 766, L2

Rameau, J., Follette, K. B., Pueyo, L., et al. 2017, AJ, 153, 244

Reggiani, M., Christiaens, V., Absil, O., et al. 2018, A&A, 611, A74

Reggiani, M., Quanz, S. P., Meyer, M. R., et al. 2014, ApJL, 792, L23

Rosenfeld, K. A., Andrews, S. M., Hughes, A. M., Wilner, D. J., & Qi, C. 2013, ApJ, 774, 16

Sallum, S., Follette, K. B., Eisner, J. A., et al. 2015, Natur, 527, 342

Shakura, N. I., & Sunyaev, R. A. 1973, A&A, 24, 337

Stolker, T., Dominik, C., Avenhaus, H., et al. 2016, arXiv

Takahashi, S. Z., & Inutsuka, S.-i 2014, ApJ, 794, 55

Thalmann, C., Janson, M., Garufi, A., et al. 2016, ApJL, 828, L17

Thalmann, C., Mulders, G. D., Janson, M., et al. 2015, ApJL, 808, L41

Walsh, C., Daley, C., Facchini, S., & Juhasz, A. 2017, A&A, 607, A114

Weaver, E., Isella, A., & Boehler, Y. 2018, ApJ, 853, 113

Zhang, K., Blake, G. A., & Bergin, E. A. 2015, ApJL, 806, L7