2017A&A...598A.118C

Context. Quantifying the gas surface density inside the dust cavities and gaps of transition disks is important to establish their origin.
Aims. We seek to constrain the surface density of warm gas in the inner disk of HD 139614, an accreting 9Myr Herbig Ae star with a (pre-)transition disk exhibiting a dust gap from 2.3±0.1 to 5.3±0.3AU.
Methods. We observed HD 139614 with ESO/VLT CRIRES and obtained high-resolution (R∼90000) spectra of CO ro-vibrational emission at 4.7µm. We derived constraints on the disk's structure by modeling the CO isotopolog line-profiles, the spectroastrometric signal, and the rotational diagrams using grids of flat Keplerian disk models.
Results. We detected υ = 1-0 ¹²CO, 2-1 ¹²CO, 1-0 ¹³CO, 1-0 C¹⁸O, and 1-0 C¹⁷O ro-vibrational lines. Lines are consistent with disk emission and thermal excitation. ¹²CO υ=1-0 lines have an average width of 14km/s, T_gas of 450K and an emitting region from 1 to 15AU. ¹³CO and C¹⁸O lines are on average 70 and 100K colder, 1 and 4km/s narrower than ¹²CO υ=1-0, and are dominated by emission at R≥6AU. The ¹²CO υ=1-0 composite line-profile indicates that if there is a gap devoid of gas it must have a width narrower than 2AU. We find that a drop in the gas surface density (δ_gas) at R<5-6AU is required to be able to simultaneously reproduce the line-profiles and rotational diagrams of the three CO isotopologs. Models without a gas density drop generate ¹³CO and C¹⁸O emission lines that are too broad and warm. The value of δ_gas can range from 10^–2 to 10^–4 depending on the gas-to-dust ratio of the outer disk. We find that the gas surface density profile at 1<R<6AU is flat or increases with radius. We derive a gas column density at 1<R<6AU of N_H=3x10¹⁹-10²¹cm^–2 (7x10^–5-2.4x10^–3g/cm²) assuming N_CO=10^–4N_H. We find a 5σ upper limit on the CO column density N_CO at R≤1AU of 5x10¹⁵cm^–2 (N_H≤5x10¹⁹cm^–2).
Conclusions. The dust gap in the disk of HD 139614 has molecular gas. The distribution and amount of gas at R≤6AU in HD 139614 is very different from that of a primordial disk. The gas surface density in the disk at R≤1AU and at 1<R<6AU is significantly lower than the surface density that would be expected from the accretion rate of HD 139614 (10^–8M_☉/yr) assuming a standard viscous α-disk model. The gas density drop, the non-negative density gradient in the gas inside 6AU, and the absence of a wide (>2AU) gas gap, suggest the presence of an embedded <2M_J planet at around 4AU.