Issue |
A&A
Volume 632, December 2019
|
|
---|---|---|
Article Number | L12 | |
Number of page(s) | 9 | |
Section | Letters to the Editor | |
DOI | https://doi.org/10.1051/0004-6361/201936750 | |
Published online | 11 December 2019 |
Letter to the Editor
Multiple nitrogen reservoirs in a protoplanetary disk at the epoch of comet and giant planet formation
1
Univ. Grenoble Alpes, CNRS, IPAG, 38000 Grenoble, France
e-mail: pierre.hily-blant@univ-grenoble-alpes.fr
2
IRAM, 300 Rue de la Piscine, Grenoble, France
3
Chester F. Carlson Center for Imaging Science, School of Physics & Astronomy, and Laboratory for Multiwavelength Astrophysics, Rochester Institute of Technology, 54 Lomb Memorial Drive, Rochester, NY 14623, USA
Received:
20
September
2019
Accepted:
4
November
2019
The isotopic ratio of nitrogen measured in primitive Solar System bodies shows a broad range of values, the origin of which remains unknown. One key question is whether these isotopic reservoirs of nitrogen predate the comet formation stage or are posterior to it. Another central question is elucidating the processes that can produce the observed variations in the 14N/15N isotopic ratio. Disks that orbit pre-main-sequence (T Tauri) stars provide unique opportunities for observing the chemical content of analogs of the protosolar nebula and therefore for building a comprehensive scenario that can explain the origin of nitrogen in the Solar System and in planet-forming disks. With ALMA, it has become possible to measure isotopic ratios of nitrogen-bearing species in such environments. We present spectrally and spatially resolved observations of the hyperfine structure of the 4−3 rotational transition of HCN and its main isotopologs H13CN and HC15N in the disk orbiting the 8 Myr old T Tauri star TW Hya. The sensitivity allows directly measuring the HCN/H13CN and HCN/HC15N abundance ratios with minimal assumptions. Averaged spatially over the disks, the ratios are 86 ± 4 and 223 ± 21, respectively. The latter value is significantly lower than the CN/C15N ratio of 323 ± 30 in this disk and thus provides the first evidence that two isotopic reservoirs of nitrogen are present in a disk at the stage of giant planet and comet formation. Furthermore, we find clear evidence for an increase in the ratio of HCN to HC15N with radius. The ratio in the outer disk, at 45 au, is 339 ± 28, in excellent agreement with direct measurements in the local interstellar medium, and with the bulk nitrogen isotopic ratio predicted from galactic evolution calculations. In the comet formation region at r = 20 au, the ratio is a factor ≈3 lower, 121 ± 11. This radial increase qualitatively agrees with the scenario in which selective photodissociation of N2 is the dominant fractionation process. However, our isotopic ratios and kinetic temperature of the HCN-emitting layers quantitatively disagree with models of nitrogen chemistry in disks.
Key words: comets: general / protoplanetary disks
© P. Hily-Blant et al. 2019
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1. Introduction
The isotopic ratio of nitrogen, 14N/15N, in various bodies of the Solar System (meteorites, comets, planets, etc.) shows the largest variations among the most abundant constituents (carbon and oxygen), with values ranging from ∼50 in submicron grains that are immersed in a chondrite matrix (so-called hotspots) to 441 in the Sun and Jupiter (Bonal et al. 2010; Marty et al. 2011; Fouchet et al. 2000; Hily-Blant et al. 2013). When, where, and how these isotopic reservoirs of nitrogen formed during the evolution from a molecular cloud to planetary systems remains to be determined. Comets provide a special case among primitive cosmomaterials: regardless of the cometary type (short or long period) and of the carrier of nitrogen that is observed in their coma (CN, HCN, NH2, NH3, N2, or NO), the isotopic ratio is consistently found to be ≈140 (Wampfler et al. 2018; Hily-Blant et al. 2017, hereafter Paper I, and references therein). This is lower by a factor three than the ratio of the bulk proto-Sun, which is 441. This well-established observational fact remains unexplained.
Evolved protoplanetary disks at the stage of giant planet and comet formation provide unique opportunities for investigating how and when this 15N-rich reservoir of nitrogen was built up. The nitrogen isotopic ratio was measured via rare isotopologs1 of HCN, specifically, H13CN and HC15N, for a small sample of disks with the ALMA interferometer (Guzmán et al. 2017), yielding an HCN/HC15N between 83 ± 32 and 156 ± 71, and a mean value of 111 ± 19, under the assumption that HCN/H13CN = 70. On the other hand, a direct measurement, one that does not involve assuming a 12C/13C ratio, of CN/C15N in the TW Hya system yielded a significantly higher ratio, 323 ± 30 (Paper I). Based on this value and on directly measured values in the dense local interstellar medium (ISM), a value of 14N/15N = 330 ± 30 was proposed for the present-day bulk solar neighborhood. This also agrees excellently well with predictions from the galactic chemical evolution model of Romano et al. (2017) for the local ISM. A new picture thus emerged in which HCN and CN trace two isotopic reservoirs of nitrogen that are present in analogs of the protosolar nebula (PSN) at the time of comet formation, with HCN representing a 15N-enriched reservoir of nitrogen similar to that recorded in Solar System comets. However, HCN/HC15N was only measured indirectly, and moreover, this was not done in TW Hya, therefore clear-cut evidence for this scenario was lacking.
To elucidate how the various nitrogen reservoirs in the Solar System were formed requires identifying the processes that lead to variations of 14N/15N (Hily-Blant et al. 2013, 2017; Füri & Marty 2015). The efficiency of chemical mass fractionation, which is the only process expected in cold and shielded gas (Terzieva & Herbst 2000), has recently been called into question both theoretically (Roueff et al. 2015; Wirström & Charnley 2018) and observationally by measuring HCN/HC15N in the L1498 pre-stellar core (Magalhães et al. 2018). The obtained value, 338 ± 28, agrees very well with the present-day elemental ratio of ∼330 proposed in Paper I and thus supports the prediction that chemical fractionation is not an efficient process even in cold pre-stellar cores. Further credence was provided by subsequent direct measurements in HC3N and HC5N toward similar sources (Hily-Blant et al. 2018; Taniguchi & Saito 2017). The concordance between observations and models may be only an apparent one, however, because sharp deviations are recurrently obtained for N2H+ (Bizzocchi et al. 2013; Redaelli et al. 2018).
The other fractionation process, selective photodissociation of N2, which is expected in unshielded gas, has been investigated in the context of protoplanetary disks and molecular clouds (Heays et al. 2014; Visser et al. 2018; Furuya & Aikawa 2018). The opacity of N15N is far lower than that of N2, therefore more 15N atoms are released upon dissociation than 14N. This leads to an atomic isotopic ratio that is lower than the bulk. Gas-phase and/or surface chemistry then propagates this deviation into other species. A correlation (positive or negative) between the UV flux and the 14N/15N is therefore expected in some species. For HCN in disks, the HCN/HC15N is predicted to increase with decreasing UV flux.
In this Letter, we aim at providing the requisite evidence for two distinct N reservoirs within the disk that orbits the nearby T Tauri star TW Hya through an analysis of the emission of HCN and isotopologs, based on archival ALMA data (2016.1.00629.S, Cleeves et al.). We also explore the possibility of a radial gradient in the HCN/HC15N to determine whether selective photodissociation is the driving fractionation process in PSN analogs.
2. Observations and data reduction
The (4−3) rotational transitions of HCN, H13CN, and HC15N at 354.5, 345.4, and 344.2 GHz, respectively, have been mapped with the ALMA interferometer with a spectral resolution of 61 KHz or 0.053 km s−1 (see Appendix A). All the lines were observed simultaneously, thus mitigating cross-calibration biases. The standard pipeline of CASA 5.4.0 was used to generate the complex visibilities, which were then self-calibrated in Gildas2. The final sensitivities are 2.3, 1.9, and 2.0 mJy beam−1 per 61 kHz channel for HCN, H13CN, and HC15N respectively. The resulting integrated intensity maps are shown in Fig. 1.
Fig. 1. Integrated intensity maps of HCN, H13CN, and HC15N (from left to right) in Jy beam km s−1. The maps are centered at α = 11:01:51.81150, δ = −34:42:17.2636 (J2000.0). The visibilities were projected on the same spatial grid and convolved to a synthesized beam of HPBW 0.26″ × 0.21″. The emission is centered on the systemic velocity of 2.83 km s−1 and is integrated over 5 km s−1 for HCN and 2 km s−1 for the less abundant isotopologs. |
The HCN(4−3) spectra (see Fig. 2) contain two weak hyperfine (hf) transitions of equal intensity, each carrying 2.1% of the total flux, and shifted by −1.672 and +1.360 km s−1 with respect to our reference frequency 354505.477 MHz. They surround a strong, non-Gaussian, central feature made of three overlapping hf transitions representing 95.8% of the flux. The intrinsically weaker intensity of H13CN and HC15N compared to HCN leads to one single feature only, which in the case of H13CN carries 95.8% of the total flux. For HC15N, the splitting is entirely unresolved and the single feature thus carries 100% of the rotational transition flux (see Table A.1).
Fig. 2. Spectra of HCN, H13CN, and HC15N (from top to bottom) from the Keplerian mask at 2.83 km s−1 (systemic velocity) and averaged within annuli centered at radii r = 5 to 55 au. The specific intensity scale is K using K/(Jy beam−1) conversion factors of 178.7, 188.1, and 189.0 for HCN, H13CN, and HC15N, respectively. |
3. Results
The conversion of line flux into abundance ratios is a simple procedure provided the lines are optically thin and the excitation temperatures are known. Even for optically thin transitions, however, the overlap in velocity space of lines that are emitted at different locations in the disk, which is a result of the Keplerian rotation of the disk, complicates the conversion from flux into abundance ratios because this mixes different regions of the disk, broadens the lines, and may also weaken the assumption that the emitting species are cospatial. The latter is especially true in the case of species with very different abundances, such as nitrogen isotopologs. Mitigating both the optical depth and Keplerian blurring effects thus offers clear advantages.
The nearly face-on configuration of the TW Hya disk naturally reduces the Keplerian blurring. In addition, the contribution from different locations along the line of sight to the emerging lines is limited to the vertical direction. Moreover, the spectrally resolved hf multiplet of HCN provides two weak lines that are most likely optically thin. Their flux is thus directly proportional to their column density.
Because the emission of all three isotopologs is well detected and spatially resolved, one goal might be measuring the HCN isotopic ratios as a function of radius, provided that line broadening can be minimized. This implies that the spatial area over which spectra are averaged needs to be limited. Another complementary approach consists of deriving a spatially averaged isotopic ratio by measuring line flux ratios over the whole disk. This method provides higher sensitivity, but Keplerian line broadening must still be controlled. We applied both approaches and detail them in the following.
3.1. Spatially averaged line fluxes
In order to limit the Keplerian line-broadening, the spectra were averaged over several groups of pixels that were selected for their closely similar projected Keplerian velocity. In practice, a 2D mask was thus determined by selecting the disk area where the 354 507.455 MHz hyperfine emission of HCN is higher than 40 mJy beam−1 in each velocity channel. Although the overlap of different Keplerian velocities is well controlled, it is still noticeable close to the disk center, and we further masked out the spectra located at radii smaller than 0.3″. The resulting mask recovered the usual dipole-like pattern, but with a central hole due to the radial filter (see Fig. B.1). As a result, the spectra averaged within each Keplerian velocity mask have a high signal-to-noise ratio (see Fig. B.2), allowing the two weak hf transitions of HCN to be properly fitted. The line fluxes of the two weak hf lines at 354 503.869 and 354 507.455 MHz were measured through Gaussian fitting and their averaged value, corrected for their summed relative intensity (4.2%), was then used as an estimator of the HCN flux (see Table B.1).
3.2. Radially dependent fluxes
In order to measure the flux ratios as a function of distance from the central T Tauri star and still avoid Keplerian line-broadening, a radial mask was superimposed on the previously described Keplerian masks. To calculate disk radial offsets in au, we adopted a distance of 59.5 ± 0.9 pc to TW Hya (Gaia Collaboration 2018). Spectra contained within concentric annuli, centered at radii r = 5 − 55 au in steps of 5 au (one-third of the half-power beam width, HPBW), and of constant thickness δr = 5 au, were averaged. The outer radius was imposed by the HC15N detection limit (see Fig. 1). The number of spectra in each annulus increased with r, which compensated for the natural decrease of the signal-to-noise ratio with radius. In the process, spectra were weighted according to their flux rms. As in the spatially averaged approach, the line fluxes of HCN, H13CN, and HC15N were measured through Gaussian fitting.
Radially averaged spectra corresponding to the Keplerian mask at 2.83 km s−1 are shown in Fig. 2. The velocity confusion at small radii is evident, while at radii above 20 au, it is efficiently removed by our masking procedure. Average spectra in all Keplerian channels are displayed in Fig. B.3. While the peak intensity of HCN hardly changes beyond 20 au, the peak intensities of H13CN and HC15N decrease when moving outward. Also, all lines become narrower as the radius increases (see Table B.3). For HCN, the FWHM is 1.5 ± 0.7 km s−1 at r ≤ 0.2″, decreasing to 0.34 ± 0.02 km s−1 at 0.9″.
4. Isotopic ratios
Because the three molecules have very similar spectroscopic constants, their rotational flux ratios are equal to their column density ratios to within 1%, provided the lines are optically thin and have the same excitation temperature. Under these assumptions, which we verify below, we obtain the isotopic ratios directly from the fluxes. The results are shown in Fig. 3 and are listed in Tables 1 and B.2.
Fig. 3. Left: HCN/HC15N and HCN/H13CN isotopic ratios in each Keplerian channel derived from the spatially averaged procedure (see Sect. 3.1). The colored bands indicate the average value ±1σ. The spatially averaged CN/C15N ratio from Paper I is shown as a hatched rectangle. Right: radial dependence (see Sect. 3.2) of the HCN/H13CN and HCN/HC15N ratios. Each box is 5 au wide (one-third of the HPBW) and ±1σ in height. The full line shows the best fit of a power law to the 15N ratio, with a slope of 1.3. The gray shaded areas are the present-day 14N/15N and 12C/13C elemental ratios in the local ISM (330 ± 30 (Paper I) and 65 ± 10 Halfen et al. 2017, respectively). |
For HCN/HC15N, we find ratios between 185 and 252 across the Keplerian channels (typical uncertainty of 5−10%) with an average isotopic ratio 223 ± 21. Here and elsewhere, uncertainties on the absolute flux include a conservative 10% error (statistical) from amplitude calibration. The quoted uncertainty on average values is the largest of the weighted uncertainties and the dispersion of the averaged values. The lowest values (185, 186, and 198) are obtained in the three central Keplerian channels, and the average ratio in the eight remaining channels is 233 ± 12. HCN/H13CN is more uniform, with values ranging from 82 to 92, which averages to 86 ± 4. The corresponding average of H13CN/HC15N is 2.6 ± 0.2, with values ranging from 2.22 to 2.75. The low values of HCN/HC15N might arise because the opacity of HCN is substantial even in its weak hf lines. However, the almost constant HCN/H13CN suggests that these lines are instead optically thin. Still, assuming a constant isotopic ratio of 233 and adopting the usual expression τ/(1 − exp(−τ)) for the opacity correction for the integrated intensity, with τ the opacity (Goldsmith & Langer 1999), an opacity of 0.4 for the weakest hf lines, corresponding to a rotational opacity τ12 ≈ 20, would be needed to explain the flux ratios in these three channels. Alternatively, assuming a constant H13CN/HC15N ratio of 2.7 based on the outer channels (see Table B.2), the total H13CN opacity would be τ13 = 0.3 in the three middle channels. It is thus likely that the opacity of the weakest hf line of HCN(4−3) and of H13CN(4−3) is ≈0.3, which is not entirely negligible, but still fulfills our optically thin assumption. We note that in these three channels, the HCN/H13CN would change by less than 4% if these corrections were applied.
The radially dependent approach shows that the HCN-to-HC15N isotopic ratio varies with radius while HCN/H13CN remains constant within 1σ. More specifically, the former ratio increases from 121 at 20 au to 339 at 45 au, and may decrease at larger distances. The nitrogen isotopic ratio in HCN close to the star is thus a factor 2.7 lower than the spatially averaged CN/C15N, 323 ± 30, whereas between 45 and 50 au, the two ratios become consistent within 1σ. The ratios change by less than 4% when the three central Keplerian channels are excluded. The best fit of a power law to the HCN/HC15N gives 117 ± 10×(r/20 au)1.3 ± 0.13, indicating a shallower dependence with radius than a pure geometrical dilution of the UV flux from the central star (∝r−2).
A key assumption in our analysis is that a single excitation temperature characterizes the weak hf lines of HCN and the overlapping hf lines of H13CN and of HC15N. However, line overlap may lead to unequal excitation temperatures among the hf manifold of HCN (Magalhães et al. 2018). We have tested our assumption (see Appendix C) by computing the excitation temperatures of HCN, H13CN, and HC15N lines while taking line overlap in the excitation into account. Our calculations assume uniform physical conditions and cover the density and gas temperature regimes from 106 to 1011 cm−3 and 5 to 30 K, respectively. We found that deviations from a single excitation temperature are only obtained at densities below 3 × 107 cm−3 (see Fig. C.1). On the other hand, the full width at half-maximum (FWHM) of the HCN, H13CN, and HC15N lines at a radius of 0.6″ (≈35 au) and beyond is ≈0.35 km s−1, which places an upper limit on the kinetic temperature of (FWHM/2.35)2 × μmH/(2kB)≈35 K (Teague et al. 2016). Referring to the TW Hya model of van’t Hoff et al. (2017), this upper limit indicates that the density is at these radii higher than 109 cm−3, which means that it is in a regime where lines are expected to be thermalized. This validates our basic assumption. Our upper limit on the kinetic temperature is also well below the gas temperature of ∼200 K where the HCN abundance is the highest in the generic models of Visser et al. (2018).
5. Discussion
The spatially averaged directly measured CN/C15N and HCN/HC15N in the TW Hya circumstellar disk provide clear-cut confirmation of at least two isotopic reservoirs of nitrogen, as anticipated by Paper I. In CN, the ratio is 1.4 ± 0.2 times higher than in HCN, and their difference is 90 ± 36. Based on the ∼8 Myr age of TW Hya, these results indicate that multiple isotopic reservoirs are present at the key protoplanetary disk evolutionary stages when comets and giant planets are likely forming.
In the comet formation zone between 20 and 30 au, the HCN/HC15N is 121 ± 11, which is a factor of 2.8 compared to the outer disk at 55 au, where the isotopic ratio of HCN matches the revised elemental ratio of 330 well. The 15N-enrichment in the region inside of 30 au is reminiscent of the threefold 15N-enrichment of comets compared to the bulk proto-Sun, providing a strong indication that HCN traces the fractionated reservoir that is recorded by comets. Guzmán et al. (2017) also found marginal evidence for an increase of the HCN-based 14N/15N ratio with radius in the inner ∼50 au of the V4046 Sgr disk. We also note that the inner disk ratio confirms the mean HCN/HC15N ratio of 111 ± 19 (Paper I) that was obtained indirectly by Guzmán et al. (2017) assuming HCN/H13CN = 70 ± 20. Adopting our measured averaged of 86 for the latter would bring the former ratio to 136 ± 23. We also very tentatively observe a decrease of the ratio beyond 55 au, which could indicate chemical processes in a gas not exposed to strong UV radiation and where depletion is limited. Measurements with a high signal-to-noise ratio beyond 55 au are needed to confirm this behavior, however.
The average HCN/H13CN, 86 ± 4, is significantly higher than the 12C/13C ratio of 65 ± 10 in the local ISM (Halfen et al. 2017), but agrees surprisingly well with the value of 89 ± 1 in the PSN (Clayton & Nittler 2004). This indicates that HCN is slightly depleted in 13C in TW Hya, opposite to what was found in the L1498 pre-stellar core (Magalhães et al. 2018), suggesting that the underlying process may not be the same in these two environments. The depletion of HCN in 13C is at variance with the models of Visser et al. (2018), however, which instead predict that HCN is enriched in 13C at radii below ∼100 au. Optically thick H13CN emission could artificially increase the HCN/H13CN provided that the weak HCN hf transition remains thin. However, the H13CN opacity required to decrease the ratio to the elemental value is 0.55, which is higher than the value of 0.3 derived previously for the central Keplerian channels. This also implies that the total opacity of HCN would be 35, or 0.75 for the weakest HCN hf, which is twice higher than the value of 0.3 that is required to bring the spatially averaged HCN/HC15N to a constant value (see Sect. 4). We note that high 12C/13C ratios have already been measured in carbon monoxide from ices toward young stellar objects (Smith et al. 2015).
The increasing HCN/HC15N with radius may indicate that selective photodissociation represents the dominant fractionation process in disks (Visser et al. 2018). In these models, the HCN-to-HC15N column density ratio first decreases with radius before it increases back to the elemental ratio as a consequence of UV flux attenuation that is due to a combination of geometrical dilution and dust extinction. We note that the minimum ratio in the grid of models of Visser et al. (2018) is ≈240−280 for disks around T Tauri stars, which would translate into 180−210 when 330 (rather than 441) is assumed for the elemental ratio. This is still higher than the lowest ratio reported here and might indicate that the UV flux is underestimated in these models. However, dedicated models of TW Hya are required because the precise radial dependence of the isotopic ratio strongly depends on the physical parameters of the model (disk mass, grain size, and dust-to-gas mass ratio).
Overall, the increasing HCN/HC15N from 121 up to the bulk value in the outer disk supports a scenario in which a single reservoir of interstellar, not fractionated, nitrogen dominates at large radii in PSN analogs. This nitrogen reservoir is enriched in 15N by the selective photodissociation of N2 at smaller distance from the central star.
IRAM memo available at http://www.iram.fr/IRAMFR/ARC/documents/filler/casa-gildas.pdf
Acknowledgments
PHB thanks D. Petry for his help in processing the ALMA data, and A. Faure for useful comments and discussions. We also thank the referee for valuable comments that greatly clarified the paper. This study makes use of the ADS/JAO.ALMA#2016.1.00629.S data. 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. This work has received financial support from the Institut Universitaire de France and from Labex OSUG@2020. JHK’s research on protoplanetary disks orbiting young stars near Earth is supported by NASA Exoplanets Research Program grants NNX16AB43G and 80NSSC19K0292 to RIT. This study made use of the CDMS database (Müller et al. 2005).
References
- Beckwith, S. V. W., & Sargent, A. I. 1993, ApJ, 402, 280 [NASA ADS] [CrossRef] [Google Scholar]
- Bizzocchi, L., Caselli, P., Leonardo, E., & Dore, L. 2013, A&A, 555, A109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bonal, L., Huss, G. R., Krot, A. N., et al. 2010, Geochim. Cosmochim. Acta, 74, 6590 [NASA ADS] [CrossRef] [Google Scholar]
- Clayton, D. D., & Nittler, L. R. 2004, ARA&A, 42, 39 [NASA ADS] [CrossRef] [Google Scholar]
- Daniel, F., & Cernicharo, J. 2008, A&A, 488, 1237 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Fouchet, T., Lellouch, E., Bézard, B., et al. 2000, Icarus, 143, 223 [NASA ADS] [CrossRef] [Google Scholar]
- Füri, E., & Marty, B. 2015, Nat. Geosc., 8, 515 [CrossRef] [Google Scholar]
- Furuya, K., & Aikawa, Y. 2018, ApJ, 857, 105 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Goldsmith, P. F., & Langer, W. D. 1999, ApJ, 517, 209 [NASA ADS] [CrossRef] [Google Scholar]
- Guzmán, V. V., Öberg, K. I., Huang, J., Loomis, R., & Qi, C. 2017, ApJ, 836, 30 [NASA ADS] [CrossRef] [Google Scholar]
- Halfen, D. T., Woolf, N. J., & Ziurys, L. M. 2017, ApJ, 845, 158 [NASA ADS] [CrossRef] [Google Scholar]
- Heays, A. N., Visser, R., Gredel, R., et al. 2014, A&A, 562, A61 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hily-Blant, P., Bonal, L., Faure, A., & Quirico, E. 2013, Icarus, 223, 582 [Google Scholar]
- Hily-Blant, P., Magalhaes, V., Kastner, J., et al. 2017, A&A, 603, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hily-Blant, P., Faure, A., Vastel, C., et al. 2018, MNRAS, 480, 1174 [NASA ADS] [CrossRef] [Google Scholar]
- Magalhães, V. S., Hily-Blant, P., Faure, A., Hernandez-Vera, M., & Lique, F. 2018, A&A, 615, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marty, B., Chaussidon, M., Wiens, R. C., Jurewicz, A. J. G., & Burnett, D. S. 2011, Science, 332, 1533 [NASA ADS] [CrossRef] [Google Scholar]
- Müller, H. S. P., Schlöder, F., Stutzki, J., & Winnewisser, G. 2005, J. Mol. Struct., 742, 215 [NASA ADS] [CrossRef] [Google Scholar]
- Redaelli, E., Bizzocchi, L., Caselli, P., et al. 2018, A&A, 617, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Romano, D., Matteucci, F., Zhang, Z.-Y., Papadopoulos, P. P., & Ivison, R. J. 2017, MNRAS, 470, 401 [NASA ADS] [CrossRef] [Google Scholar]
- Roueff, E., Loison, J. C., & Hickson, K. M. 2015, A&A, 576, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Smith, R. L., Pontoppidan, K. M., Young, E. D., & Morris, M. R. 2015, ApJ, 813, 120 [NASA ADS] [CrossRef] [Google Scholar]
- Taniguchi, K., & Saito, M. 2017, PASJ, 69, L7 [NASA ADS] [CrossRef] [Google Scholar]
- Teague, R., Guilloteau, S., Semenov, D., et al. 2016, A&A, 592, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Terzieva, R., & Herbst, E. 2000, MNRAS, 317, 563 [NASA ADS] [CrossRef] [Google Scholar]
- van’t Hoff, M. L. R., Walsh, C., Kama, M., Facchini, S., & van Dishoeck, E. F. 2017, A&A, 599, A101 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Visser, R., Bruderer, S., Cazzoletti, P., et al. 2018, A&A, 615, A75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wampfler, S., Rubin, M., Altwegg, K., et al. 2018, in COSPAR Meeting, 42nd COSPAR Scientific Assembly, 42, B1.3–21–18 [Google Scholar]
- Wirström, E. S., & Charnley, S. B. 2018, MNRAS, 474, 3720 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Observations
The observations were performed with the ALMA interferometer (2016.1.00629.S, PI:I. Cleeves) and were spread over six epochs (2016 December 30, 2017 July 5, 2017 July 10, 2017 July 15, 2017 July 21, and 2017 July 22). J1037−2934 was used as a bandpass, flux, and phase calibrator. The phase was also calibrated with J1103−3251, and J1058+0133 was also used for bandpass calibration.
H13CN and HC15N were tuned in the lower side-band and HCN in the upper side-band, allowing simultaneous observations for the three species. We used the standard pipeline of CASA 5.4.0 to generate the first sets of complex visibilities, which were then exported into the Gildas format using the Gildas-Casa filer of Guilloteau, Chapillon et al. (see the IRAM memo3). Self-calibration was then applied to generate the final set of visibilities.
Spectroscopy shows that the three hyperfine (hf) transitions of HCN(4−3) are separated by less than 0.13 km s−1 and are not resolved. Their total relative intensity is 95.8%. The hf transition at 354505.846 MHz is not detected given its 0.03% relative strength. The properties of the hf lines of HCN and H13CN are summarized in Table A.1, while for HC15N, the unresolved hyperfine splitting is not given. We also show in Fig. B.1 channel maps of HCN strong central feature, H13CN, and HC15N. The velocity range covers the main and overlapping hf lines of HCN.
Spectroscopic information on the observed transitions.
Appendix B: Keplerian masks
B.1. Results from the spatially averaged procedure
In Fig. B.1 (top row), we show the Keplerian mask defined as explained in the main text. Each mask is overlaid in the channel maps of the three isotopologs, showing the usual dipole-like pattern (Beckwith & Sargent 1993). Comparing the rows showing the main feature and the weak hf line of HCN, we note that the extended HCN emission is filtered out in the present study, as is the case of H13CN, although to a lesser extent. The spectra averaged within each mask, shown in Fig. B.2, were fit with Gaussian profiles, the results of which are summarized in Table B.1. The obtained fluxes are then used to derive the spatially averaged isotopic ratios listed in Table B.2.
Fig. B.1. Keplerian masks (top row) overlaid on (from top to bottom) the HCN (main feature and weak hf line at 345507.455 MHz), H13CN, and HC15N, channel maps. The masks (top row) were obtained for a threshold of 40 mJy beam−1 applied to the weak hf transition at 354507.455 MHz. |
Results of the Gaussian fits to the spatially averaged spectra of Fig. B.2.
Isotopic ratios from the spatially averaged spectra.
Fig. B.2. Spatially averaged spectra (Jy beam−1) within each Keplerian velocity channel. HCN, H13CN, and HC15N are shown from top to bottom. The Keplerian channels are those of Fig. B.1. |
B.2. Radial dependence procedure
By combining two masks, one in velocity space (Keplerian mask) and one in the radial direction, we are able to study the radial dependence of the 14N/15N isotopic ratio in HCN toward the TW Hya system. In Fig. 2, we show the spectra averaged within each annulus, but we concentrate on the central Keplerian channel at 2.83 km s−1 (systemic velocity). Figure B.3 is complementary in that it shows the spectra averaged at three radii, but within all Keplerian masks.
Fig. B.3. Spectra of HCN (black), H13CN (red, ×5), and HC15N (green, ×10) averaged within each Keplerian velocity channel (see, e.g., Table B.1) and within annuli of constant thickness 5 au, centered at 30, 40, and 50 au (from top to bottom). The specific intensity scale is K using K/(Jy beam−1) conversion factors of 178.7, 188.1, and 189.0 for HCN, H13CN, and HC15N, respectively. The spectra are shifted vertically for the sake of legibility. |
As in the radially averaged approach, the spectral line intensity and shape at any radius do not change appreciably in the various Keplerian velocity channels. In contrast, the dependence on radius is significant, especially for H13CN and most importantly for HC15N. As can be seen in Table B.3 and as mentioned in the main text, all lines become narrower as the radius increases (see Table B.3).
Radial dependence of the line width (FWHM, km s−1) of HCN, H13CN, and HC15N.
B.3. Radially dependent line width
Gaussian fits to Keplerian-filtered spectra averaged within concentric annuli of thickness 0.1″ were performed in order to measure the dependence of the FWHM of the HCN, H13CN, and HC15N lines. The results are summarized in Table B.3. They are used to infer an upper limit on the kinetic temperature following Paper I and Teague et al. (2016).
Appendix C: Single excitation temperature assumption
The main assumption in our analysis is that of a single excitation temperature characterizing the weak hf lines of HCN and the overlapping hf lines of H13CN and of HC15N. A thorough investigation of this assumption can only be made with two-dimensional radiative transfer taking into account the overlap of hf transitions in the level population evaluation, which in turn requires a complete physical and kinematic structure of the disk. This type of analysis has been made in the context of pre-stellar starless cores by Magalhães et al. (2018), in spherical geometry, with the 1D ALICO model (Daniel & Cernicharo 2008).
In the present context, we were unable to make detailed calculations like this. Instead, we used our 1D ALICO code to investigate the density and kinetic gas temperature at which departures from the single excitation temperature assumption may be obtained. To do so, the level populations were computed while treating the effect of hf line overlap in the excitation of HCN, H13CN, and HC15N. Our grid of models encompasses a range of H2 density (106 − 1011 cm−3) and kinetic temperature (from 5 to 30 K) relevant to the TW Hya orbiting disk. The results of these calculations in terms of Tex are shown in Fig. C.1. There are significant departures at densities between 106 cm−3 and 108 cm−3. At higher density, all lines are thermalized (models at nH2 higher than 109 cm−3 are not shown). From the generic chemical models of Visser et al. (2018), it appears that HCN is concentrated at radii r = 5 − 50 au and scale heights z/r = 0.2 to 0.3, where the gas kinetic temperature is ≈200 K and the density in the range 106 to 108 cm−3, which would indicate that the single excitation temperature assumption fails.
Fig. C.1. Excitation temperature of the six hf lines of HCN (see Table A.1) and the strongest hf lines of H13CN (red) and HC15N (green). The four overlapping hf lines of HCN are shown as open circles and the two weakest resolved lines with filled dots. |
On the other hand, the FWHM of the HCN, H13CN, and HC15N lines at a radius of 0.6″ (≈35 au) and beyond is ≈0.35 km s−1 (see Table B.3), placing a conservative upper limit on the kinetic temperature of (Teague et al. 2016)
Referring to the tailored TW Hya model of van’t Hoff et al. (2017), this upper limit, applied to the relevant range of radii, indicates that the density is higher than 109 cm−3, in a regime where lines should be thermalized. This validates our single excitation temperature assumption. Our upper limit on the kinetic temperature is also well below the gas temperature of ∼200 K where the HCN abundance is the highest in the models of Visser et al. (2018), although we note that these models are generic and are not tailored to TW Hya.
All Tables
All Figures
Fig. 1. Integrated intensity maps of HCN, H13CN, and HC15N (from left to right) in Jy beam km s−1. The maps are centered at α = 11:01:51.81150, δ = −34:42:17.2636 (J2000.0). The visibilities were projected on the same spatial grid and convolved to a synthesized beam of HPBW 0.26″ × 0.21″. The emission is centered on the systemic velocity of 2.83 km s−1 and is integrated over 5 km s−1 for HCN and 2 km s−1 for the less abundant isotopologs. |
|
In the text |
Fig. 2. Spectra of HCN, H13CN, and HC15N (from top to bottom) from the Keplerian mask at 2.83 km s−1 (systemic velocity) and averaged within annuli centered at radii r = 5 to 55 au. The specific intensity scale is K using K/(Jy beam−1) conversion factors of 178.7, 188.1, and 189.0 for HCN, H13CN, and HC15N, respectively. |
|
In the text |
Fig. 3. Left: HCN/HC15N and HCN/H13CN isotopic ratios in each Keplerian channel derived from the spatially averaged procedure (see Sect. 3.1). The colored bands indicate the average value ±1σ. The spatially averaged CN/C15N ratio from Paper I is shown as a hatched rectangle. Right: radial dependence (see Sect. 3.2) of the HCN/H13CN and HCN/HC15N ratios. Each box is 5 au wide (one-third of the HPBW) and ±1σ in height. The full line shows the best fit of a power law to the 15N ratio, with a slope of 1.3. The gray shaded areas are the present-day 14N/15N and 12C/13C elemental ratios in the local ISM (330 ± 30 (Paper I) and 65 ± 10 Halfen et al. 2017, respectively). |
|
In the text |
Fig. B.1. Keplerian masks (top row) overlaid on (from top to bottom) the HCN (main feature and weak hf line at 345507.455 MHz), H13CN, and HC15N, channel maps. The masks (top row) were obtained for a threshold of 40 mJy beam−1 applied to the weak hf transition at 354507.455 MHz. |
|
In the text |
Fig. B.2. Spatially averaged spectra (Jy beam−1) within each Keplerian velocity channel. HCN, H13CN, and HC15N are shown from top to bottom. The Keplerian channels are those of Fig. B.1. |
|
In the text |
Fig. B.3. Spectra of HCN (black), H13CN (red, ×5), and HC15N (green, ×10) averaged within each Keplerian velocity channel (see, e.g., Table B.1) and within annuli of constant thickness 5 au, centered at 30, 40, and 50 au (from top to bottom). The specific intensity scale is K using K/(Jy beam−1) conversion factors of 178.7, 188.1, and 189.0 for HCN, H13CN, and HC15N, respectively. The spectra are shifted vertically for the sake of legibility. |
|
In the text |
Fig. C.1. Excitation temperature of the six hf lines of HCN (see Table A.1) and the strongest hf lines of H13CN (red) and HC15N (green). The four overlapping hf lines of HCN are shown as open circles and the two weakest resolved lines with filled dots. |
|
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.