No Arabic abstract
We report observations of the cyanide anion, CN, in the disk around TW~Hya covering the $N=1-0$, $N=2-1$ and $N=3-2$ transitions. Using line stacking techniques, 24 hyperfine transitions are detected out of the 30 within the observed frequency ranges. Exploiting the super-spectral resolution from the line stacking method reveals the splitting of hyperfine components previously unresolved by laboratory spectroscopy. All transitions display a similar emission morphology, characterized by an azimuthally symmetric ring, peaking at $approx 45$~au (0.75), and a diffuse outer tail extending out to the disk edge at $approx 200$~au. Excitation analyses assuming local thermodynamic equilibrium (LTE) yield excitation temperatures in excess of the derived kinetic temperatures based on the local line widths for all fine structure groups, suggesting assumptions of LTE are invalid. Using the 0D radiative transfer code RADEX, we demonstrate that such non-LTE effects may be present when the local H$_2$ density drops to $10^{7}~{rm cm^{-3}}$ and below. Comparison with models of TW~Hya find similar densities at elevated regions in the disk, typically $z , / , r gtrsim 0.2$, consistent with model predictions where CN is formed via vibrationally excited H$_2$ in the disk atmospheric layers where UV irradiation is less attenuated.
The recent detection of gas-phase methanol (CH$_3$OH) lines in the disc of TW Hya by Walsh et al. provided the first observational constraints on the complex O-bearing organic content in protoplanetary discs. The emission has a ring-like morphology, with a peak at $sim 30-50$ au and an inferred column density of $sim 3-6times10^{12}$ cm$^{-2}$. A low CH$_3$OH fractional abundance of $sim 0.3-4times 10^{-11}$ (with respect to H$_2$) is derived, depending on the assumed vertical location of the CH$_3$OH molecular layer. In this study, we use a thermo-chemical model of the TW Hya disc, coupled with the ALCHEMIC gas-grain chemical model, assuming laboratory-motivated, fast diffusivities of the surface molecules to interpret the CH$_3$OH detection. Based on this disc model, we performed radiative transfer calculations with the LIME code and simulations of the observations with the CASA simulator. We found that our model allows to reproduce the observations well. The CH$_3$OH emission in our model appears as a ring with radius of $sim60$ au. Synthetic and observed line flux densities are equal within the rms noise level of observations. The synthetic CH$_3$OH spectra calculated assuming local thermodynamic equilibrium (LTE) can differ by up to a factor of 3.5 from the non-LTE spectra. For the strongest lines, the differences between LTE and non-LTE flux densities are very small and practically negligible. Variations in the diffusivity of the surface molecules can lead to variations of the CH$_3$OH abundance and, therefore, line flux densities by an order of magnitude.
We present a detailed analysis of the spatially and spectrally resolved 12CO J=2-1 and J=3-2 emission lines from the TW Hya circumstellar disk, based on science verification data from the Atacama Large Millimeter/Submillimeter Array (ALMA). These lines exhibit substantial emission in their high-velocity wings (with projected velocities out to 2.1 km/s, corresponding to intrinsic orbital velocities >20 km/s) that trace molecular gas as close as 2 AU from the central star. However, we are not able to reproduce the intensity of these wings and the general spatio-kinematic pattern of the lines with simple models for the disk structure and kinematics. Using three-dimensional non-local thermodynamic equilibrium molecular excitation and radiative transfer calculations, we construct some alternative models that successfully account for these features by modifying either (1) the temperature structure of the inner disk (inside the dust-depleted disk cavity; r < 4 AU); (2) the intrinsic (Keplerian) disk velocity field; or (3) the distribution of disk inclination angles (a warp). The latter approach is particularly compelling because a representative warped disk model qualitatively reproduces the observed azimuthal modulation of optical light scattered off the disk surface. In any model scenario, the ALMA data clearly require a substantial molecular gas reservoir located inside the region where dust optical depths are known to be substantially diminished in the TW Hya disk, in agreement with previous studies based on infrared spectroscopy. The results from these updated model prescriptions are discussed in terms of their potential physical origins, which might include dynamical perturbations from a low-mass companion with an orbital separation of a few AU.
We obtain high spatial and spectral resolution images of the CO J=2-1, CN N=2-1 and CS J=5-4 emission with ALMA in Cycle~2. The radial distribution of the turbulent broadening is derived with three approaches: two `direct and one modelling. The first requires a single transition and derives Tex{} directly from the line profile, yielding a vturb{}. The second assumes two different molecules are co-spatial thus their relative linewidths allow for a calculation of Tkin{} and vturb{}. Finally we fit a parametric disk model where physical properties of the disk are described by power laws, to compare our `direct methods with previous values. The two direct methods were limited to the outer $r > 40$~au disk due to beam smear. The direct method found vturb{} ranging from $approx$~vel{130} at 40~au, dropping to $approx$~vel{50} in the outer disk, qualitatively recovered with the parametric model fitting. This corresponds to roughly $0.2 - 0.4~c_s$. CN was found to exhibit strong non-LTE effects outside $r approx 140$~au, so vturb{} was limited to within this radius. The assumption that CN and CS are co-spatial is consistent with observed linewidths only within $r lesssim 100$~au, within which vturb{} was found to drop from vel{100} ($approx~0.4~c_s$) to nothing at 100~au. The parametric model yielded a near constant vel{50} for CS ($0.2 - 0.4~c_s$). We demonstrate that absolute flux calibration is and will be the limiting factor in all studies of turbulence using a single molecule. The magnitude of the dispersion is comparable with or below that predicted by the magneto-rotational instability theory. A more precise comparison would require to reach an absolute calibration precision of order 3%, or to find a suitable combination of light and heavy molecules which are co-located in the disk.
Despite their importance in the star formation process, measurements of magnetic field strength in proto-planetary discs remain rare. While linear polarisation of dust and molecular lines can give insight into the magnetic field structure, only observations of the circular polarisation produced by Zeeman splitting provide a direct measurement of magnetic field strenghts. One of the most promising probes of magnetic field strengths is the paramagnetic radical CN. Here we present the first Atacama Large Millimeter/submillimeter Array (ALMA) observations of the Zeeman splitting of CN in the disc of TW~Hya. The observations indicate an excellent polarisation performance of ALMA, but fail to detect significant polarisation. An analysis of eight individual CN hyperfine components as well as a stacking analysis of the strongest (non-blended) hyperfine components yields the most stringent limits obtained so far on the magnetic field strength in a proto-planetary disc. We find that the vertical component of the magnetic field $|B_z|<0.8$~mG ($1sigma$ limit). We also provide a $1sigma$ toroidal field strength limit of $<30$~mG. These limits rule out some of the earlier accretion disc models, but remain consistent with the most recent detailed models with efficient advection. We detect marginal linear polarisation from the dust continuum, but the almost purely toroidal geometry of the polarisation vectors implies that his is due to radiatively aligned grains.
We analyze high angular resolution ALMA observations of the TW Hya disk to place constraints on the CO and dust properties. We present new, sensitive observations of the $^{12}$CO $J = 3-2$ line at a spatial resolution of 8 AU (0farcs14). The CO emission exhibits a bright inner core, a shoulder at $rapprox70$ AU, and a prominent break in slope at $rapprox90$ AU. Radiative transfer modeling is used to demonstrate that the emission morphology can be reasonably reproduced with a $^{12}$CO column density profile featuring a steep decrease at $rapprox15$ AU and a secondary bump peaking at $rapprox70$ AU. Similar features have been identified in observations of rarer CO isotopologues, which trace heights closer to the midplane. Substructure in the underlying gas distribution or radially varying CO depletion that affects much of the disks vertical extent may explain the shared emission features of the main CO isotopologues. We also combine archival 1.3 mm and 870 $mu$m continuum observations to produce a spectral index map at a spatial resolution of 2 AU. The spectral index rises sharply at the continuum emission gaps at radii of 25, 41, and 47 AU. This behavior suggests that the grains within the gaps are no larger than a few millimeters. Outside the continuum gaps, the low spectral index values of $alphaapprox 2$ indicate either that grains up to centimeter size are present, or that the bright continuum rings are marginally optically thick at millimeter wavelengths.