Inferring giant planets from ALMA millimeter continuum and line observations in (transition) disks

• Published in: Astronomy and astrophysics (Berlin) 2018-04, Vol.612, p.A104
• Main Authors: Facchini, S., Pinilla, P., van Dishoeck, E. F., de Juan Ovelar, M.
• Format: Article
• Language: English
• Subjects: astrochemistry; Chemical speciation; Computer simulation; Density; Dependence; Dust; Emissions; Evolution; Gas temperature; Holes; Infrared imagery; Infrared signatures; Mathematical models; Organic chemistry; Parameters; Planet formation; Planetary orbits; Planets; planet–disk interactions; protoplanetary disks; Protoplanets; Radiative transfer; Radio telescopes; Rotation; Simulation; submillimeter: planetary systems; Two dimensional models; Viscosity
• ISSN: 0004-6361; 1432-0746
• DOI: 10.1051/0004-6361/201731390

Context. Radial gaps or cavities in the continuum emission in the IR-mm wavelength range are potential signatures of protoplanets embedded in their natal protoplanetary disk are. Hitherto, models have relied on the combination of mm continuum observations and near-infrared scattered light images to put constraints on the properties of embedded planets. Atacama Large Millimeter/submillimeter Array (ALMA) observations are now probing spatially resolved rotational line emission of CO and other chemical species. These observations can provide complementary information on the mechanism carving the gaps in dust and additional constraints on the purported planet mass. Aims. We investigate whether the combination of ALMA continuum and CO line observations can constrain the presence and mass of planets embedded in protoplanetary disks. Methods. We post-processed azimuthally averaged 2D hydrodynamical simulations of planet-disk models, in which the dust densities and grain size distributions are computed with a dust evolution code that considers radial drift, fragmentation, and growth. The simulations explored various planet masses (1 MJ ≤ Mp ≤ 15 MJ) and turbulent parameters (10−4 ≤ α ≤ 10−3). The outputs were then post-processed with the thermochemical code DALI, accounting for the radially and vertically varying dust properties. We obtained the gas and dust temperature structures, chemical abundances, and synthetic emission maps of both thermal continuum and CO rotational lines. This is the first study combining hydrodynamical simulations, dust evolution, full radiative transfer, and chemistry to predict gas emission of disks hosting massive planets. Results. All radial intensity profiles of 12CO, 13CO, and C18O show a gap at the planet location. The ratio between the location of the gap as seen in CO and the peak in the mm continuum at the pressure maximum outside the orbit of the planet shows a clear dependence on planet mass and is independent of disk viscosity for the parameters explored in this paper. Because of the low dust density in the gaps, the dust and gas components can become thermally decoupled and the gas becomes colder than the dust. The gaps seen in CO are due to a combination of gas temperature dropping at the location of the planet and of the underlying surface density profile. Both effects need to be taken into account and disentangled when inferring gas surface densities from observed CO intensity profiles; otherwise, the gas surface density