No Arabic abstract
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 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.
We report the detection of spiral substructure in both the gas velocity and temperature structure of the disk around TW~Hya, suggestive of planet-disk interactions with an unseen planet. Perturbations from Keplerian rotation tracing out a spiral pattern are observed in the SE of the disk, while significant azimuthal perturbations in the gas temperature are seen in the outer disk, outside 90~au, extending the full azimuth of the disk. The deviation in velocity is either $Delta v_{phi} , / , v_{rm kep} sim 0.1$ or $Delta v_{z} , / , v_{rm kep} sim 0.01$ depending on whether the perturbation is in the rotational or vertical direction, while radial perturbations can be ruled out. Deviations in the gas temperature are $pm 4$ K about the azimuthally averaged profile, equivalent to deviations of $Delta T_{rm gas} , / , T_{rm gas} sim 0.05$. Assuming all three structures can be described by an Archimedean spiral, measurements of the pitch angles of both velocity and temperature spirals show a radially decreasing trend for all three, ranging from 9$^{circ}$ at 70 au, dropping to 3$^{circ}$ at 200 au. Such low pitch-angled spirals are not readily explained through the wake of an embedded planet in the location of previously reported at 94 au, but rather require a launching mechanism which results in much more tightly wound spirals. Molecular emission tracing distinct heights in the disk is required to accurately distinguish between spiral launching mechanisms.
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.
We report two-dimensional spectroastrometry of Br gamma emission of TW Hya to study the kinematics of the ionized gas in the star-disk interface region. The spectroastrometry with the integral field spectrograph SINFONI at the Very Large Telescope is sensitive to the positional offset of the line emission down to the physical scale of the stellar diameter (~0.01 AU). The centroid of Br gamma emission is displaced to the north with respect to the central star at the blue side of the emission line, and to the south at the red side. The major axis of the centroid motion is P.A.= -20 degrees, which is nearly equal to the major axis of the protoplanetary disk projected on the sky, previously reported by CO sub millimeter spectroscopy (P.A.= -27 degrees) The line-of-sight motion of the Br gamma emission, in which the northern side of the disk is approaching toward us, is also consistent with the direction of the disk rotation known from the CO observation. The agreement implies that the kinematics of Br gamma emission is accounted for by the ionized gas in the inner edge of the disk. A simple modeling of the astrometry, however, indicates that the accretion inflow similarly well reproduces the centroid displacements of Br gamma, but only if the position angles of the centroid motion and the projected disk ellipse is a chance coincidence. No clear evidence of disk wind is found.
For over a decade, the structure of the inner cavity in the transition disk of TW Hydrae has been a subject of debate. Modeling the disk with data obtained at different wavelengths has led to a variety of proposed disk structures. Rather than being inconsistent, the individual models might point to the different faces of physical processes going on in disks, such as dust growth and planet formation. Our aim is to investigate the structure of the transition disk again and to find to what extent we can reconcile apparent model differences. A large set of high-angular-resolution data was collected from near-infrared to centimeter wavelengths. We investigated the existing disk models and established a new self-consistent radiative-transfer model. A genetic fitting algorithm was used to automatize the parameter fitting. Simple disk models with a vertical inner rim and a radially homogeneous dust composition from small to large grains cannot reproduce the combined data set. Two modifications are applied to this simple disk model: (1) the inner rim is smoothed by exponentially decreasing the surface density in the inner ~3 AU, and (2) the largest grains (>100 um) are concentrated towards the inner disk region. Both properties can be linked to fundamental processes that determine the evolution of protoplanetary disks: the shaping by a possible companion and the different regimes of dust-grain growth, respectively. The full interferometric data set from near-infrared to centimeter wavelengths requires a revision of existing models for the TW Hya disk. We present a new model that incorporates the characteristic structures of previous models but deviates in two key aspects: it does not have a sharp edge at 4 AU, and the surface density of large grains differs from that of smaller grains. This is the first successful radiative-transfer-based model for a full set of interferometric data.