CO desorption in the AS 209 disk

Evidence for a CO desorption front in the outer AS 209 disk

Abstract

Millimeter observations of CO isotopologues are often used to make inferences about protoplanetary disk gas density and temperature structures. The accuracy of these estimates depends on our understanding of CO freezeout and desorption from dust grains. Most models of these processes indicate that CO column density decreases monotonically with distance from the central star due to a decrease in gas density and freezeout beyond the CO snowline. We present ALMA Cycle 2 observations of CO, CO, and CO emission that instead suggest CO enhancement in the outer disk of T Tauri star AS 209. Most notably, the CO emission consists of a central peak and a ring at a radius of (120 AU), well outside the expected CO snowline. We propose that the ring arises from the onset of CO desorption near the edge of the millimeter dust disk. CO desorption exterior to a CO snowline may occur via non-thermal processes involving cosmic rays or high-energy photons, or via a radial thermal inversion arising from dust migration.

astrochemistry—ISM: molecules—protoplanetary disks—radio lines: ISM

1 Introduction

Characterizing protoplanetary disk structures is essential for understanding planet formation. Disk properties are often inferred by modeling millimeter observations of CO isotopologues (e.g. Dartois et al. 2003; Rosenfeld et al. 2013a; Williams & Best 2014). The vertical structure of a disk is typically described in terms of a surface layer in which strong UV radiation photodissociates molecules, an intermediate warm molecular layer, and a cold midplane (Aikawa et al. 2002). The relationship between CO abundances and overall gas density structure is tied to the disk temperature profile. Near the star, CO is in the gas phase, but in the cold outer disk midplane, CO is expected to be frozen out onto dust grains. The boundary separating these two regions is the CO snowline, thought to influence the formation locations and compositions of planets and planetesimals (e.g., Öberg et al. 2011a; Qi et al. 2013; Ali-Dib et al. 2014). While the fractional abundance of gas-phase CO is generally assumed to be constant in the warm molecular layer, depletion due to photodissociation and freezeout in other layers must be accounted for in order to avoid underestimating the disk mass or surface density (Qi et al. 2011). Because CO is optically thick in disks, CO and CO observations have been used to derive disk masses, but these may be substantially underestimated if selective photodissociation is neglected (Miotello et al. 2014). Comparisons of gas mass estimates from HD and CO observations for the TW Hya disk also indicate that the CO:H ratio in the warm molecular layer may not follow typical ISM abundances (Bergin et al. 2013; Favre et al. 2013).

Desorption in the outer disk midplane, preventing the complete freezeout of CO exterior to the snowline, may be an additional complicating factor in characterizing disk structure with CO observations. Two mechanisms may return CO to the gas phase in the outer disk. First, non-thermal desorption may occur when disk densities decrease enough for cosmic rays, X-rays, and UV radiation to penetrate to the midplane (e.g. Willacy 2007; Walsh et al. 2010). Second, the inward drift of large grains or external photoevaporation may invert the radial temperature profile, leading to thermal CO desorption (Cleeves 2016; Facchini et al. 2016). Öberg et al. (2015) presented indirect evidence for CO desorption in the outer disk of T Tauri star IM Lup with a detection of double DCO rings, attributing the outer ring to increased production of DCO due to UV photodesorption of CO ice.

In this letter, we present evidence for a CO desorption front in the outer disk of AS 209, based primarily on the emission pattern of CO. In sections 2 and 3, we describe ALMA observations of CO, CO, and CO transitions in the disk around AS 209, a 1.6 Myr T Tauri star with a mass of 0.9 (Andrews et al. 2009). The system is thought to be part of the Ophiuchus star-forming region about 120 pc away (Loinard et al. 2008). In section 4, we simulate observations with a toy model to demonstrate that a CO abundance enhancement in the outer disk is consistent with the different isotopologue emission morphologies. Section 5 discusses approaches for further characterizing CO desorption in disks and implications for disk structure inference.

2 Observations and Data Reduction

AS 209 (J2000.0 R.A. , decl. ) was observed with the Atacama Large Millimeter/Submillimeter Array on 2014 July 2 (project code ADS/JAO.ALMA#2013.1.00226), with 21 minutes on source. The configuration consisted of thirty-four 12 m antennae, with baselines between 20 and 650 m. Observations were set up with thirteen Band 6 spectral windows (SPWs). Twelve (including those with CO isotopologue lines) had spectral resolutions of 61 kHz and bandwidths of 59 MHz, while the thirteenth had a resolution of 122 kHz and bandwidth of 469 MHz. An earlier reduction of the 1.4 mm continuum and CO data was published in Huang & Öberg (2015) as part of an ND analysis.

ALMA/NAASC staff provided a calibration script using the quasar J1733-1304 for bandpass and phase calibration and Titan for flux calibration. It was modified to scale visibility weights properly1 and executed in CASA 4.4.0 to calibrate visibilities. The CO isotopologue SPWs were phase self-calibrated with solutions obtained from averaging six SPWs free of strong line emission. After subtracting the continuum in the uv-plane, each line was imaged and CLEANed with Keplerian rotation masks. The 1.4 mm dust continuum noise level is mJy beam, and the ALMA systematic flux uncertainty is 10. The continuum flux density, obtained by integrating interior to the 3 contour, is 252 25 mJy, which is consistent with Submillimeter Array observations of AS 209 presented in Öberg et al. (2011b). Channel rms for 0.1 km s bins and integrated flux values are listed in Table 1.

Transition Rest Frequency E Beam (P. A.) Channel rms Integrated Flux
(GHz) (K) (mJy beam) (Jy km s)
CO 230.53800 16.6 10 7.50.8
CO 220.39868 15.9 15 2.10.2
CO 219.56035 15.8 8 0.540.05
Table 1: Summary of Line Observations

3 Observational Results

Figure 1: AS 209 disk emission maps and radial intensity profiles. Left column: Top left panel shows the 1.4 mm continuum intensity, with [4, 8, 16, 32…] contours, where mJy beam. The next three panels are CO, CO, and CO integrated intensity maps. Line emission color bars start at 2. Synthesized beams are shown in each panel’s lower left. Red crosses mark the continuum peak position. Right column: Deprojected and normalized continuum, CO, CO and CO radial intensity profiles. Purple ribbons show the standard deviation at each radius.
Figure 2: Channel maps for the AS 209 disk CO transition. Synthesized beams are shown in the lower left and channel velocities (km s) in the lower right of each panel. Red crosses mark the phase center. The red ellipse traces a projected radius, assuming P.A. = 86 and inclination = 38. Contours are [3, 5, 7], with = 5 mJy beam.

Column one of Figure 1 shows the 1.4 mm dust continuum intensity and CO, CO, and CO integrated intensity maps (summed from -2.0 to 11.0 km s). Column two of Figure 1 shows corresponding deprojected and azimuthally averaged intensity profiles. The position angle and inclination of 86 and 38, respectively, are adopted from Andrews et al. 2009. The continuum intensity is normalized to its peak value, 71.7 mJy beam (synthesized beam: ). Line profiles are normalized to the peak integrated intensities of 881, 266, and 91 mJy beam km s for CO, CO, and CO, respectively. The rms in the integrated intensity maps are 23, 27, and 14 mJy beam km s for CO, CO, and CO, respectively.

The continuum and CO emission are centrally peaked and decrease monotonically with radius, although there is significant CO emission outside the continuum detection threshold. The CO emission is weaker on the west side, likely due to cloud contamination (Öberg et al. 2011b). In contrast to CO, the CO profile bulges outward at (120 AU). Like the other isotopologues, the CO emission is centrally peaked, but also has a ring at coinciding with the CO “bulge.”

Figure 2 shows CO channel maps. A red ellipse traces a projected radius. Line wings contribute to the central peak observed in the integrated intensity map, but in channels near the systemic velocity of km s, little emission is present within of the disk center. Like the CO integrated intensity map, the channel maps demonstrate that the emission traces out a ring at .

4 A toy model for CO desorption in the outer disk

A CO desorption front, the onset of desorption in the outer disk outside a CO snowline, may produce an emission ring. To explore qualitatively how such a CO distribution can affect the emission profiles of different isotopologues, we simulated observations of a T Tauri disk using a toy model adapted from parametric gas density and temperature models that have been applied to CO emission in many disks (e.g. Rosenfeld et al. 2013a; Williams & Best 2014; Czekala et al. 2015). Of the isotopologues observed in the AS 209 disk, we assume that CO, the most optically thin, best traces the CO column density. Because the millimeter continuum emission from the AS 209 disk is optically thin at all radii (Pérez et al. 2012; Tazzari et al. 2016), dust opacity is expected to have negligible effects on the CO emission. Our CO surface density model is therefore motivated by the continuum-subtracted CO radial intensity profile.

The gas surface density of protoplanetary disks is often modeled with the Lynden-Bell and Pringle (1974) similarity solution for a viscous accretion disk. We adapt this for our model CO surface density profile (in cylindrical coordinates) by adding a Gaussian ring in the outer disk to simulate a large-scale return of CO into the gas phase:

(1)

Following standard prescriptions for gas abundance distributions (e.g. Rosenfeld et al. 2013b), we scale CO abundances vertically with the midplane pressure scale height:

(2)

with

(3)

where is the mass of atomic hydrogen, and is the mean gas particle mass.

Similarly to Rosenfeld et al. (2013a) and Dartois et al. (2003), we model the vertical temperature gradient in cylindrical coordinates as

(4)

where

(5)
(6)

and are the atmosphere and midplane temperatures, respectively, at AU. The temperature gradient scale height is

(7)

From a power-law fit to the midplane dust temperature that Andrews et al. (2009) calculated for the AS 209 disk, K and = 0.48. Since the atmosphere gas temperature profile of the AS 209 disk has not been constrained, we set and K, similar to values Dartois et al. (2003) found for the disk around T Tauri star DM Tau. After some experimentation to determine parameter values that could reproduce the main features of the observations, we set , , AU, , AU, and AU.

We use this model and the LAMDA database molecular data (Schöier et al. 2005; Yang et al. 2010) to calculate CO, CO, and CO intensities with the radiative transfer code RADMC-3D (Dullemond 2012). We fix the CO/CO and CO/CO ratios to the ISM values of 69 and 557, respectively (Wilson 1999). Local thermal equilibrium was assumed because CO and its isotopologues have low critical densities relative to disk gas densities (Pavlyuchenkov et al. 2007). Keplerian gas velocities were assumed, and the turbulent broadening parameter was set to 0.01 km s, on par with the TW Hya disk upper limit Hughes et al. (2011) derived. After computing sky brightness images at the orientation of the AS 209 disk, model visibilities were produced at the same spatial frequencies as the AS 209 data by using Python package vis_sample.2 Model visibilities were imaged and CLEANed in CASA 4.4.0, then integrated across the same velocity ranges as the observations to produce integrated intensity maps. Observations were then simulated by adding Gaussian noise to the model visibilities so that the rms in the model image cubes was comparable to that of the data image cubes. Temperature and column density profiles and integrated intensity maps are shown for the model with the Gaussian CO ring in Figure 3. For comparison, we also model CO abundances without an outer ring by setting in Eq. 2, holding other parameters the same. Corresponding structure plots and integrated intensity maps are shown in Figure 4.

Figure 3: Model structures and resulting integrated intensity maps for the model with an outer CO ring ( in eq. 2). Top left: Temperature gradient profile up to four pressure scale heights (Eq. 3) above the midplane. Red line denotes 20 K isotherm. Bottom left: CO column density (bottom). Middle: Integrated intensity maps for CO, CO, and CO . Right column: Corresponding maps with gaussian noise added to model visibilities. Synthesized beams are drawn in the lower left. Red crosses mark phase centers.
Figure 4: Same as Figure 3, but for a model with monotonically decreasing CO column densities ( in eq. 2).

For either model, the CO integrated emission is centrally peaked, in line with observations. However, only the model with an outer CO ring, shown in Figure 3, reproduces the central peak and ring in the CO emission, and the shoulder in CO emission. The morphological similarities between the data and model with an outer CO ring suggest that excess CO in the outer disk in conjunction with line opacity can account for the variation in emission patterns for the AS 209 disk’s CO isotopologues. The model results highlight that a non-monotonic CO radial column density profile can be obscured in (partially) thick CO and CO observations, and thus direct evidence of CO desorption in the outer disk would often require observations of rarer, optically thin isotopologues such as CO.

5 Discussion

5.1 Origin of the CO isotopologue emission morphology

The outer emission ring in CO observed in the AS 209 disk suggests that its CO abundances are enhanced at large radii. In the radial intensity profile of CO in Figure 1, the ring peaks at 120 AU. While no snowline estimate has been reported in the literature for the AS 209 disk, the CO snowline has been estimated to lie at 30 AU for the TW Hya disk and 90 AU for the HD 163296 disk, corresponding to midplane temperatures of 17 and 25 K, respectively (Qi et al. 2013, 2015). TW Hya is a T Tauri star with –4000 K and = 0.19–0.23 (Webb et al. 1999; Vacca & Sandell 2011). AS 209 is somewhat warmer ( K and ), suggesting that its CO snowline should lie at a slightly larger radius compared to TW Hya. As a Herbig Ae star, HD 163296 has a much warmer disk, so its CO snowline radius can be treated as an upper bound for AS 209. Furthermore, the AS 209 disk midplane dust temperature profile derived by Andrews et al. (2009) drops to 25 K at 35 AU and 17 K at 80 AU, suggesting that the CO snowline lies within this range. The decrement in emission between the central peak and outer ring of CO in the AS 209 disk occurs at (80 AU), consistent with the expected location of initial CO freezeout. Meanwhile, the CO ring in the AS 209 disk occurs well outside the expected CO snowline.

Based on T Tauri disk models indicating that radial drift could increase midplane dust temperatures at large radii, Cleeves (2016) suggested that subsequent thermal desorption of CO in the outer disk could cause CO enhancement followed by a second CO snowline once temperatures began falling again. Pérez et al. (2012) found that larger dust grains emitted from a more compact region than small grains in the AS 209 disk, providing evidence of radial drift of solids. Cleeves’ models and Pérez et al.’s observations suggest that the CO emission morphology in the AS 209 disk can plausibly be explained by thermal CO desorption in the outer disk, with CO snowlines both interior and exterior to the desorption ring. Non-thermal desorption may also substantially enhance outer disk CO. Further observations and modeling are necessary to constrain the relative contributions of thermal versus non-thermal processes.

Simulated observations from our toy model of CO desorption in the outer disk indicate that isotopologue opacity differences in the AS 209 disk can explain why there is a CO ring, but only a shoulder in CO and a monotonically decreasing intensity profile in CO. Selective photodissociation is a less likely alternative. Models from Miotello et al. (2014) indicate that selective photodissociation leads to lower CO intensities in the outer disk compared to models without, whereas we observe an increase in CO intensity in the outer disk of AS 209.

Another possibility to examine is that dust opacity creates the appearance of a ring in optically thin CO emission. The peak brightness temperature of the 1.4 mm AS 209 continuum is only 7 K, which is consistent with the findings from Pérez et al. 2012 and Tazzari et al. 2016 of an optically thin disk. However, Öberg et al. 2015 suggested that such low brightness temperatures could also be consistent with dust concentrated into optically thick rings much narrower than the beam. While existing observations do not point to this scenario, it could be tested with ALMA at higher spatial resolutions.

5.2 Implications for disk structure modeling

Disk structure inferences made from CO observations depend on assumptions about chemistry (e.g. Favre et al. 2013; Miotello et al. 2014; Williams & Best 2014). Common parametric disk gas structure models stipulate that CO emission decreases monotonically with radius (Rosenfeld et al. 2013a, b), in line with most CO observations and chemical models (Aikawa & Herbst 2001; Willacy 2007; Walsh et al. 2012). However, the optically thick lower- transitions of CO in disks could easily obscure CO substructures, as indicated by our toy model. Although CO isotopologue emission rings have been observed in disks, they have been primarily associated with systems with prominent dust cavities, suggesting that gas and dust cleared for similar reasons (Bruderer et al. 2014; van der Marel et al. 2016; Tang et al. 2016). The non-monotonicity of the CO radial intensity profile for the AS 209 disk, as well as the double DCO rings of the IM Lup disk in Öberg et al. (2015), are unexpected because neither disk shows millimeter dust gaps. Since the outer CO emission ring observed in the AS 209 disk seems inconsistent with the standard parameterizations used to model CO emission in disks, the observations presented in this paper indicate that desorption effects should also be accounted for when tracing disk structure with CO.

5.3 Developing further constraints on desorption

Given the evidence of CO desorption in the AS 209 disk in this work, as well as observations of IM Lup (Öberg et al. 2015), we expect that signatures of CO desorption and second snowlines will be apparent in the outer regions of other disks observed at high resolution with ALMA. Recent observations of CO and CO in the TW Hya disk indicate an outer ring, which Schwarz et al. (2016) suggested could be due to desorption. Because optically thick CO and even CO transitions could obscure significant features in CO abundance profiles, especially for more massive, gas-rich disks, sensitive high-resolution observations of CO and CO can place further constraints on desorption mechanisms and better inform CO emission parameterizations. Öberg et al. (2015) suggested that observing multiple transitions would provide better estimates of gas temperature in the outer disk, thereby distinguishing between thermal and non-thermal desorption. Better temperature measurements will also be necessary to derive CO column densities more reliably in the outer disk. Furthermore, NH is often used to locate the CO snowline because its formation is inhibited and its destruction accelerated by gas-phase CO (Qi et al. 2013). If abundant CO returns to the gas phase near the midplane, then NH may trace this desorption front.

5.4 Summary

We present ALMA observations of CO, CO, and CO emission in the AS 209 disk. The CO emission ring is interpreted as the onset of large-scale CO desorption in the outer disk. Based on simulated observations from a toy model, we propose that CO desorption in the outer disk, along with line opacity differences, can explain variations in emission patterns among the isotopologues. The CO isotopologue emission in AS 209 is not described well by common parametric models used to infer disk structure from CO observations, indicating that CO is not a straightforward tracer of total disk gas properties. To constrain the effects of desorption on CO abundances, especially on the creation of outer disk CO desorption fronts or secondary snowlines, high-resolution observations of rare, optically thin CO isotopologues are crucial.

We thank the referee for comments improving this paper, Ian Czekala, Ilse Cleeves, and Meredith MacGregor for helpful discussions, and Ryan Loomis for access to vis_sample. This paper makes use of ALMA data ADS/JAO.ALMA#2013.1.00226.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), 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. KIÖ acknowledges funding from the Alfred P. Sloan Foundation and the Packard Foundation. This material is based upon work supported by the National Science Foundation Graduate Research Fellowship under Grant No. DGE-1144152. \softwareCASA 4.4.0, RADMC-3D, pwkit (https://github.com/pkgw/pwkit), vis_sample (https://github.com/AstroChem/vis_sample), cubehelix (https://github.com/jradavenport/cubehelix)

Footnotes

  1. See
    https://casaguides.nrao.edu/index.php/DataWeightsAndCombination
  2. Available at https://github.com/AstroChem/vis_sample

References

  1. Aikawa, Y., & Herbst, E. 2001, A&A, 371, 1107
  2. Aikawa, Y., van Zadelhoff, G. J., van Dishoeck, E. F., & Herbst, E. 2002, A&A, 386, 622
  3. Ali-Dib, M., Mousis, O., Petit, J.-M., & Lunine, J. I. 2014, ApJ, 793, 9
  4. Andrews, S. M., Wilner, D. J., Hughes, A. M., Qi, C., & Dullemond, C. P. 2009, ApJ, 700, 1502
  5. Bergin, E. A., Cleeves, L. I., Gorti, U., et al. 2013, Nature, 493, 644
  6. Bruderer, S., van der Marel, N., van Dishoeck, E. F., & van Kempen, T. A. 2014, A&A, 562, A26
  7. Cleeves, L. I. 2016, ApJ, 816, L21
  8. Czekala, I., Andrews, S. M., Jensen, E. L. N., et al. 2015, ApJ, 806, 154
  9. Dartois, E., Dutrey, A., & Guilloteau, S. 2003, A&A, 399, 773
  10. Dullemond, C. P. 2012, RADMC-3D: A multi-purpose radiative transfer tool, Astrophysics Source Code Library, ascl:1202.015
  11. Facchini, S., Clarke, C. J., & Bisbas, T. G. 2016, MNRAS, 457, 3593
  12. Favre, C., Cleeves, L. I., Bergin, E. A., Qi, C., & Blake, G. A. 2013, ApJ, 776, L38
  13. Huang, J., & Öberg, K. I. 2015, ApJ, 809, L26
  14. Hughes, A. M., Wilner, D. J., Andrews, S. M., Qi, C., & Hogerheijde, M. R. 2011, ApJ, 727, 85
  15. Loinard, L., Torres, R. M., Mioduszewski, A. J., & Rodríguez, L. F. 2008, ApJ, 675, L29
  16. Lynden-Bell, D., & Pringle, J. E. 1974, MNRAS, 168, 603
  17. Miotello, A., Bruderer, S., & van Dishoeck, E. F. 2014, A&A, 572, A96
  18. Öberg, K. I., Furuya, K., Loomis, R., et al. 2015, ApJ, 810, 112
  19. Öberg, K. I., Murray-Clay, R., & Bergin, E. A. 2011a, ApJ, 743, L16
  20. Öberg, K. I., Qi, C., Fogel, J. K. J., et al. 2011b, ApJ, 734, 98
  21. Pavlyuchenkov, Y., Semenov, D., Henning, T., et al. 2007, ApJ, 669, 1262
  22. Pérez, L. M., Carpenter, J. M., Chandler, C. J., et al. 2012, ApJ, 760, L17
  23. Qi, C., D’Alessio, P., Öberg, K. I., et al. 2011, ApJ, 740, 84
  24. Qi, C., Öberg, K. I., Andrews, S. M., et al. 2015, ApJ, 813, 128
  25. Qi, C., Öberg, K. I., Wilner, D. J., et al. 2013, Science, 341, 630
  26. Rosenfeld, K. A., Andrews, S. M., Hughes, A. M., Wilner, D. J., & Qi, C. 2013a, ApJ, 774, 16
  27. Rosenfeld, K. A., Andrews, S. M., Wilner, D. J., Kastner, J. H., & McClure, M. K. 2013b, ApJ, 775, 136
  28. Schöier, F. L., van der Tak, F. F. S., van Dishoeck, E. F., & Black, J. H. 2005, A&A, 432, 369
  29. Schwarz, K. R., Bergin, E. A., Cleeves, L. I., et al. 2016, ArXiv e-prints, arXiv:1512.05679
  30. Tang, Y.-W., Dutrey, A., Guilloteau, S., et al. 2016, ApJ, 820, 19
  31. Tazzari, M., Testi, L., Ercolano, B., et al. 2016, A&A, 588, A53
  32. Vacca, W. D., & Sandell, G. 2011, ApJ, 732, 8
  33. van der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2016, A&A, 585, A58
  34. Walsh, C., Millar, T. J., & Nomura, H. 2010, ApJ, 722, 1607
  35. Walsh, C., Nomura, H., Millar, T. J., & Aikawa, Y. 2012, ApJ, 747, 114
  36. Webb, R. A., Zuckerman, B., Platais, I., et al. 1999, ApJ, 512, L63
  37. Willacy, K. 2007, ApJ, 660, 441
  38. Williams, J. P., & Best, W. M. J. 2014, ApJ, 788, 59
  39. Wilson, T. L. 1999, Reports on Progress in Physics, 62, 143
  40. Yang, B., Stancil, P. C., Balakrishnan, N., & Forrey, R. C. 2010, ApJ, 718, 1062
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
""
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
   
Add comment
Cancel
Loading ...
278669
This is a comment super asjknd jkasnjk adsnkj
Upvote
Downvote
""
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters
Submit
Cancel

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test
Test description