أهداف: تم دراسة استجابة دكس البروتوبلانيتي لانفجارات مناخية من المدة المختلفة مع الغرض من تحديد تأثير الانفجارات على قوة ومستدامة الاضطراب الجاذبي في الدكس. تم التركيز بشكل خاص على توزيع الغاز والغبار الناضج (من 1 ملم إلى بضعة سنتيمترات) أثناء وبعد الانفجار.
طرق: تم استخدام تحليلات الهيدروديناميكا الرقمية باستخدام رمز FEOSAD لدراسة ديناميكية الغاز والغبار في حد الدكس الخفيف. كما تم اعتماد الارتجاع الضوئي للانفجارات المناخية مع درجات الشدة المستمرة والنازولة لتمثيل حالتي الحدود الطبيعية لانفجارات نوع FU-Orionis.
نتائج: التأثير القصير الأمد للانفجار هو تقليص قوة الاضطراب الجاذبي بالحرارة وتوسيع الدكس. أطول الانفجارات مع مدة مقاربة للدوران الحلزوني للحلزون الذي يشكل الدكس الأصلي يمكن أن يذهب إلى تحلل النمط الذي كان موجودا في الغاز في نهاية الانفجار، بينما لا يمكن لأقصر الانفجارات سوى تضعيف النمط الحلزوني. يختلف رد الغبار الناضج على الانفجار. يتحول التوزيع الحلزوني الأصلي مع الأفق المتضخم في المناطق بين الأذرع إلى توزيع حلقي مع الأفق المتضخم في الفجوات العميقة. يتم تحقيق هذا التحول بشكل أكثر تأكيدا لأطول الانفجارات المدتية. التأثير الطويل الأمد للانفجار يعتمد على الظروف الأصلية للدكس في بداية الانفجار. في بعض الحالات، يبدأ التفريع النشط في الدكس بعد آلاف السنين من بداية الانفجار، الذي لم يكن موجودا في النموذج دون الانفجارات (مختصر).
Aims: Response of a protoplanetary disk to luminosity bursts of various duration is studied with the purpose to determine the effect of the bursts on the strength and sustainability of gravitational instability in the disk. A special emphasis is paid to the spatial distribution of gas and grown dust (from 1 mm to a few cm) during and after the burst. Methods: Numerical hydrodynamics simulations using the FEOSAD code were employed to study the dynamics of gas and dust in the thin-disk limit. Dust-to-gas friction including back reaction and dust growth were also considered. Bursts of various duration (from 100 to 500 yr) were initiated in accordance with a thermally ignited magnetorotational instability. Luminosity curves for constant- and declining-magnitude bursts were adopted to represent two typical limiting cases for FU-Orionis-type eruptions. Results: The short-term effect of the burst is to reduce the strength of gravitational instability by heating and expanding the disk. The longest bursts with duration comparable to the revolution period of the spiral can completely dissolve the original two-armed spiral pattern in the gas disk by the end of the burst, while the shortest bursts only weaken the spiral pattern. The reaction of grown dust to the burst is somewhat different. The spiral-like initial distribution with deep cavities in the inter-armed regions transforms into a ring-like distribution with deep gaps. This transformation is most expressed for the longest-duration bursts. The long-term effect of the burst depends on the initial disk conditions at the onset of the burst. In some cases, vigorous disk fragmentation sets in several thousand years after the burst, which was absent in the model without the bursts (abridged).
High-resolution observations of edge-on proto-planetary disks in emission from molecular species sampling different critical densities and formation pathways offer the opportunity to trace the vertical chemical and physical structures of protoplanetary disks. Based on analysis of sub-arcsecond resolution Atacama Large Millimeter Array (ALMA) archival data for the edge-on Flying Saucer disk (2MASS J16281370-2431391), we establish the vertical and radial differentiation of the disk CN emitting regions with respect to those of $^{12}$CO and CS, and we model the disk physical conditions from which the CN emission arises. We demonstrate that the disk $^{12}$CO (2-1), CN (2-1), and CS J=5-4 emitting regions decrease in scale height above the midplane, such that 12CO, CN, and CS trace layers of increasing density and decreasing temperature. We find that at radii > 100 au from the central star, CN emission arises predominantly from intermediate layers, while in the inner region of the disk, CN appears to arise from layers closer to the midplane. We investigate disk physical conditions within the CN emitting regions, as well as the ranges of CN excitation temperature and column density, via RADEX non-LTE modeling of the three brightest CN hyperfine lines. Near the disk midplane, where we derive densities nH2 ~10$^{7}$ cm$^{-3}$ at relatively low T$_{kin}$ (~12 K), we find that CN is thermalized, while sub-thermal, non-LTE conditions appear to obtain for CN emission from higher (intermediate) disk layers. We consider whether and how the particular spatial location and excitation conditions of CN emission from the Flying Saucer can be related to CN production that is governed, radially and vertically, by the degree of irradiation of the flared disk by X-rays and UV photons from the central star.
Determining the gas density and temperature structures of protoplanetary disks is a fundamental task to constrain planet formation theories. This is a challenging procedure and most determinations are based on model-dependent assumptions. We attempt a direct determination of the radial and vertical temperature structure of the Flying Saucer disk, thanks to its favorable inclination of 90 degrees. We present a method based on the tomographic study of an edge-on disk. Using ALMA, we observe at 0.5$$ resolution the Flying Saucer in CO J=2-1 and CS J=5-4. This edge-on disk appears in silhouette against the CO J=2-1 emission from background molecular clouds in $rho$ Oph. The combination of velocity gradients due to the Keplerian rotation of the disk and intensity variations in the CO background as a function of velocity provide a direct measure of the gas temperature as a function of radius and height above the disk mid-plane. The overall thermal structure is consistent with model predictions, with a cold ($< 15-12 $~K), CO-depleted mid-plane, and a warmer disk atmosphere. However, we find evidence for CO gas along the mid-plane beyond a radius of about 200,au, coincident with a change of grain properties. Such a behavior is expected in case of efficient rise of UV penetration re-heating the disk and thus allowing CO thermal desorption or favoring direct CO photo-desorption. CO is also detected up to 3-4 scale heights while CS is confined around 1 scale height above the mid-plane. The limits of the method due to finite spatial and spectral resolutions are also discussed. This method appears to be very promising to determine the gas structure of planet-forming disks, provided that the molecular data have an angular resolution which is high enough, of the order of $0.3 - 0.1$ at the distance of the nearest star forming regions.
We report FUV, optical, and NIR observations of three T Tauri stars in the Orion OB1b subassociation with H$alpha$ equivalent widths consistent with low or absent accretion and various degrees of excess flux in the mid-infrared. We aim to search for evidence of gas in the inner disk in HST ACS/SBC spectra, and to probe the accretion flows onto the star using H$alpha$ and He I $lambda$10830 in spectra obtained at the Magellan and SOAR telescopes. At the critical age of 5 Myr, the targets are at different stages of disk evolution. One of our targets is clearly accreting, as shown by redshifted absorption at free-fall velocities in the He I line and wide wings in H$alpha$; however, a marginal detection of FUV H$_2$ suggests that little gas is present in the inner disk, although the spectral energy distribution indicates that small dust still remains close to the star. Another target is surrounded by a transitional disk, with an inner cavity in which little sub-micron dust remains. Still, the inner disk shows substantial amounts of gas, accreting onto the star at a probably low, but uncertain rate. The third target lacks both a He I line or FUV emission, consistent with no accretion or inner gas disk; its very weak IR excess is consistent with a debris disk. Different processes occurring in targets with ages close to the disk dispersal time suggest that the end of accretion phase is reached in diverse ways.
We present 870 $mu$m ALMA observations of polarized dust emission toward the Class II protoplanetary disk IM Lup. We find that the orientation of the polarized emission is along the minor axis of the disk, and that the value of the polarization fraction increases steadily toward the center of the disk, reaching a peak value of ~1.1%. All of these characteristics are consistent with models of self-scattering of submillimeter-wave emission from an optically thin inclined disk. The distribution of the polarization position angles across the disk reveals that while the average orientation is along the minor axis, the polarization orientations show a significant spread in angles; this can also be explained by models of pure scattering. We compare the polarization with that of the Class I/II source HL Tau. A comparison of cuts of the polarization fraction across the major and minor axes of both sources reveals that IM Lup has a substantially higher polarization fraction than HL Tau toward the center of the disk. This enhanced polarization fraction could be due a number of factors, including higher optical depth in HL Tau, or scattering by larger dust grains in the more evolved IM Lup disk. However, models yield similar maximum grain sizes for both HL Tau (72 $mu$m) and IM Lup (61 $mu$m, this work). This reveals continued tension between grain-size estimates from scattering models and from models of the dust emission spectrum, which find that the bulk of the (unpolarized) emission in disks is most likely due to millimeter (or even centimeter) sized grains.
The protoplanetary disk around HL Tau is so far the youngest candidate of planet formation, and it is still embedded in a protostellar envelope with a size of thousands of au. In this work, we study the gas kinematics in the envelope and its possible influence on the embedded disk. We present our new ALMA cycle 3 observational results of HL Tau in the 13CO (2-1) and C18O (2-1) emission at resolutions of 0.8 (110 au), and we compare the observed velocity pattern with models of different kinds of gas motions. Both the 13CO and C18O emission lines show a central compact component with a size of 2 (280 au), which traces the protoplanetary disk. The disk is clearly resolved and shows a Keplerian motion, from which the protostellar mass of HL Tau is estimated to be 1.8+/-0.3 M$_odot$, assuming the inclination angle of the disk to be 47 deg from the plane of the sky. The 13CO emission shows two arc structures with sizes of 1000-2000 au and masses of 3E-3 M$_odot$ connected to the central disk. One is blueshifted and stretches from the northeast to the northwest, and the other is redshifted and stretches from the southwest to the southeast. We find that simple kinematical models of infalling and (counter-)rotating flattened envelopes cannot fully explain the observed velocity patterns in the arc structures. The gas kinematics of the arc structures can be better explained with three-dimensional infalling or outflowing motions. Nevertheless, the observed velocity in the northwestern part of the blueshifted arc structure is ~60-70% higher than the expected free-fall velocity. We discuss two possible origins of the arc structures: (1) infalling flows externally compressed by an expanding shell driven by XZ Tau and (2) outflowing gas clumps caused by gravitational instabilities in the protoplanetary disk around HL Tau.