{'fraction_non_alphanumeric': 0.03862735553855782, 'fraction_numerical': 0.03563729921424101, 'mean_word_length': 4.241071428571429, 'pattern_counts': {'":': 0, '<': 1, '<?xml version=': 0, '>': 0, 'https://': 0, 'lorem ipsum': 0, 'www.': 0, 'xml': 0}, 'pii_count': 0, 'substrings_counts': 0, 'word_list_counts': {'cursed_substrings.json': 0, 'profanity_word_list.json': 0, 'sexual_word_list.json': 0, 'zh_pornsignals.json': 0}}
ALMA observations of the Th 28 protostellar disk - A new example of counter-rotation between disk and optical jet

Recently, differences in Doppler shifts across the base of four close classical T Tauri star jets have been detected with the HST in optical and near-ultraviolet (NUV) emission lines, and these Doppler shifts were interpreted as rotation signatures under the assumption of steady state flow. To support this interpretation, it is necessary that the underlying disks rotate in the same sense. Agreement between disk rotation and jet rotation determined from optical lines has been verified in two cases and rejected in one case. Meanwhile, the NUV lines, which may trace faster and more collimated inner spines of the jet than optical lines, either agree or show no clear indication. We propose to perform this test on the fourth system, Th 28. We present ALMA high angular resolution Band 7 continuum, 12CO(3-2) and 13CO(2-1) observations of the circumstellar disk around the T Tauri star Th 28. We were able to detect, in CO and continuum, clear signatures of a disk in Keplerian rotation around Th28. The 12CO emission is resolved, allowing us to derive estimates of disk position angle and inclination. The large velocity separation of the peaks in 12CO, combined with the resolved extent of the emission, indicate a central stellar mass in the range 1-2 Msun. The rotation sense of the disk is well detected in both 13CO and 12CO emission lines, and this direction is opposite to that implied by the transverse Doppler shifts measured in the optical lines of the jet. The Th 28 system is the second system where counter-rotation between the disk and the optical jet is detected. These findings imply either that optical transverse velocity gradients detected with HST do not trace jet rotation or that modeling the flow with the steady assumption is not valid. In both cases jet rotation studies that rely solely on optical lines are not suitable to derive the launching radius of the jet.


Introduction
The interplay between accretion and ejection of matter is believed to be a crucial element in the formation of stars. Yet, the exact link between jets and accretion disks is still a critical issue in contemporary astrophysics ). The most attractive possibility is a fundamental connection via angular momentum transfer from disk to jet by means of magnetocentrifugal forces, so that circumstellar material may continue to accrete onto the central object (e.g., Blandford & Payne 1982). Exactly where this transfer occurs and how it impacts the disk physics is, however, still hotly debated (Ferreira et al. 2006;Pudritz et al. 2007;Shang et al. 2007;Romanova et al. 2009;Cabrit 2009). Owing to their low extinction and their proximity (≈ 150 pc), accreting classical T Tauri stars (hereafter CTTSs) offer a unique laboratory to elucidate the connection between accretion and ejection.
A recent major breakthrough was the detection with the STIS instrument mounted on the HST of transverse Doppler shifts of 5 -30 km.s −1 in optical and near-ultraviolet (hereafter NUV) emission lines on scales of 50 AU across the base of four close CTTS jets Woitas et al. 2005;Coffey et al. 2004Coffey et al. , 2007. If these shifts are due to jet rotation in a steady state outflow, they would imply that jets are magnetically launched from the disk surface at relatively small-scale radii of 0.1-1.6 AU, thus solving the long-standing problem of the jet origin in CTTS Anderson et al. 2003;Coffey et al. 2004;Woitas et al. 2005). However, effects other than rotation (e.g., precession; see Cerqueira et al. 2006) could cause transverse Doppler shifts similar to those detected in T Tauri jets. It is crucial that the underlying disks rotate in the same sense as the jet to support the jet rotation interpretation. This comparison was conducted in three of the four systems mentioned above (CW Tau, DG Tau, and RW Aur), plus in RY Tau. In RY Tau, the study was inconclusive as no systematic Doppler shifts were found in the jet; see Coffey et al. (2015). The sense of rotation of the jets of CW Tau and DG Tau (Coffey et al. 2007), traced by optical lines, coincides with the sense of rotation of their disks (Testi et al. 2002;Cabrit private comm., respectively). The jet rotation sense traced in near ultra-violet (NUV) emission lines at higher flow velocities is coherent for DG Tau, but these emission lines are below the detection threshold for CW Tau (Coffey et al. 2007). In RW Aur, the optical jet and disk rotations appear opposite (see Coffey et al. 2004 andCabrit et al. 2006). Here, jet rotation traced in the NUV was found to be compatible with the disk rotation sense (i.e., opposite to the optical jet rotation), but could only be measured in one lobe and at one of the two epochs (spaced by six months; see Coffey et al. 2012); because of its complexity and variability the authors declared RW Aur as unsuitable for jet rotation studies. We present the confrontation of jet and disk rotation senses for the fourth system, Th 28.   integrated from -15 to 20 km.s −1 . Contours start at 5 σ with 20 σ steps. Middle: moment zero of the 13 CO(3-2) integrated from -15 to 20 km.s −1 . Contours start at 5 σ with 5 σ steps. Right: continuum map. Contours start at 5 σ with 20 σ steps. Noise level at 1 σ is indicated in the upper right corner of each panel; the unit is mJy/beam.km/s in the two left panels and mJy/beam in the right panel. The black ellipses on the three panels represent the FWHM intensity contours. White hatched ellipses represent the clean beam FWHM.
The T Tauri star ThA 15-28 (hereafter Th 28; The 1962), also known as Sz 102 (Schwartz 1977) or Krautter's star, is a young member of the Lupus 3 association. Although uncertain, as the photospheric spectrum is veiled, a spectral type G8 to K2 was quoted for the central object of Th 28 (Graham & Heyer 1988;Hughes et al. 1994;Mortier et al. 2011). Even for a K star, the stellar emission from Th 28 is significantly underluminous: 0.03 L (Mortier et al. 2011), which suggests that the disk is observed close to edge-on and obscures the central star. The derived disk luminosity, inferred from the spectral energy distribution (SED) integrated excess above the stellar photosphere, is ∼0.13 L (Mortier et al. 2011;Merín et al. 2008). We adopt the recent value of d=185 +11 −10 pc derived by Galli et al. (2013) for the Lupus 3 cloud, with an improved convergent point search method, for the distance to the source.
Th 28 drives a bright and striking bipolar jet whose axis lies at PA=98 • (Krautter 1986;Graham & Heyer 1988). The large proper motions derived for the distant knots are compatible with a close to edge-on geometry (e.g., Comerón & Fernández 2010). Th 28 is considered one of the best case for jet rotation detection in the HST sample, as transverse velocity gradients are well detected in both lobes of its optical atomic jet and are consistent in different optical lines (Coffey et al. 2004). In optical, the jet rotation sense is therefore well determined with a clockwise direction looking down the blueshifted jet lobe toward the source. The situation for NUV lines were inconclusive (Coffey et al. 2007). Only one measurement on the receding jet, that is, however, consistent with the sense of rotation derived in optical, exceeded the incertitude level.
In this paper, we report ALMA band 7 CO and continuum observations at 0.4 -0.5 angular resolution of the Th 28 system, aimed to detect the disk and determine its sense of rotation. We detail our observations and data reduction in Sect. 2, and demonstrate that the CO and continuum emissions clearly trace a disk in Keplerian rotation around Th 28 in Sect. 3. The Sect. 4 offers a detailed analysis of the disk properties. In Sect. 5 we discuss the implications for jet launching models and central source properties. We summarize our conclusions in Sect. 6.

Observations and data reduction
The characteristics of our lines and continuum observations are detailed below. The resulting beam sizes and sensitivities are summarized in Table 1.
The 12 CO(J=3→2) and 13 CO(J=3→2) emission lines plus the continuum emission of Th 28 were observed using ALMA Band 7 (275-373 GHz). The observations were performed on April 28 2014 in two tracks of ∼ 40 minutes each during the ALMA cycle 1 campaign. The Band 7 data contained two spectral windows of 468.75 MHz bandwidth each (3840 channels) that were tuned at 345.818 GHz and 330.609 GHz. This allowed us to simultaneously cover the 12 CO(J=3→2) and 13 CO(J=3→2) at 345.796 GHz and 330.588 GHz, respectively. A third spectral window of 2 GHz bandwidth in 128 channels centered at 332.014 GHz was dedicated to the detection of the continuum emission from the Th 28 dust disk. The phase center of the observations was taken to be the position of Th 28 as seen with HST: α(J2000) = 16 h 08 m 29 s .72 and δ(J2000) = -39 • 03 11 . 01. The data were taken using the cycle 1 semiextended configuration of ALMA with baselines ranging from 20 m to 480 m.
The data were reduced using the common astronomy software application (hereafter CASA; see McMullin et al. 2006). We performed an initial correction for rapid atmospheric variations at each antenna using water vapor radiometer data and corrected for the time and frequency dependence of the system temperatures. Bandpass calibration was performed on the quasars J1733-1304 and J1517-2422. The primary flux calibration was made using Titan and J1517-2422. Quasar J1534-3526 was used in both observations to calibrate the time variation of the complex gains. Based on the dispersion between the fluxes derived for the phase calibrator in each observing session (using the two different amplitude calibrator), we estimate the absolute flux calibration to be accurate within ∼10%. Owing to residual inconsistency in the phase calibration between the two tracks, we used the continuum spectral window of one track to further phase calibrate all spectral windows. Imaging was carried out using the CLEAN algorithm of CASA. We used the Briggs weighting with The mean frequency was calculated assuming a S(ν) ∝ ν −3 emission spectra accurately describes the ISM SED slope in the 2 GHz spectral window.
The rms presented here report 1σ rms of the integrated maps shown in Fig. 1.
a robust parameter of 0.5 1 . The channel spacing of 0.122 MHz resulted in a maximum spectral resolution of ∼0.11 km.s −1 that we degraded down to 0.4 km.s −1 to detect weak 12 CO and 13 CO emission. The resulting root mean square (RMS) noise level is of 4 mJy/beam per 0.4 km.s −1 channel in 12 CO(3-2) and 13 CO(3-2). The continuum displays a RMS noise of 0.08 mJy/beam over its 2 GHz bandwidth. Panels (a), (b), and (c) in Figure 1 show the 12 CO(3-2), 13 CO(3-2), and continuum, respectively, which will be described in Sect 3.

CO emissions probing a disk in Keplerian rotation
Left panels of Figure  We detail each of these items below.
(a) Compact and elongated CO emission: The 12 CO emission shows an elliptical shape with a main axis close to the northsouth direction, roughly perpendicular to the main axis of the beam. It is therefore clearly resolved in the north-south direction. Indeed, the emitting area with fluxes above half maximum emission is an ellipse of 0.68 ×0.48 at PA≈10 • (black ellipse on Fig. 1 The PA of the long axis of the CO emissions (see (a) and (b)) are therefore perpendicular, within the uncertainties, to the PA of the optical blueshifted jet (98 • ).
(c) Double-peaked profile of the CO lines: The integrated CO line profiles are double peaked with symmetric wings (see Fig. 2-Top), exhibiting a strong resemblance to the typical profile shape of rotating Keplerian disks in CTTS (Duvert et al. 2000;Guilloteau & Dutrey 1994). The 12 CO and 13 CO show peak fluxes of ∼0.45 Jy and 0.04 Jy, respectively. The apparent separations of the two velocity peaks toward Th 28 is larger in 13 CO than in 12 CO: ∆v = 8.7±0.1 km.s −1 and 7.1±0.3 km.s −1 , respectively. Such a behavior is expected for a Keplerian disk, since the intensity peaks roughly correspond to the outer radius of the disk, which should be larger in the more optically thick 12 CO(2-1) line.
(d) The position-velocity (PV) diagram of the 12 CO emission line along the PA=7.3 • of the disk is shown in Fig. 3. It further confirms the hypothesis of a disk in rotation with a quasi-perfect point symmetry both in position and velocity. Such a pattern is indeed indicative of a globally axisymmetric rotating structure around Th 28. We have superimposed Keplerian curves for stellar central masses M =1, 2, and 3 M (see Sect. 5.2 for a more detailed analysis), assuming the disk has an inclination of 73 degrees.
With the supporting evidence presented above, we claim to observe a protoplanetary disk in Keplerian rotation around Th 28. Bottom left: first moment of the 12 CO(3-2) emission line. The black arrow indicates the PA of the disk as seen in 12 CO. The latter is used to build the PV diagram shown in Figure 3. The red and blue arrows represent the jet axis of Th 28, while the line-of-sight vectors represent their rotation sense in optical tracers derived by Coffey et al. (2004). Bottom right: first moment of the 13 CO(3-2) emission line. In these two first-moment maps, a clear velocity gradient due to the rotation of the circumstellar disk of Th 28 is seen.  signal-to-noise (S/N) ratio of ∼150, we report the unambiguous detection of continuum emission in the circumstellar disk of Th 28. The comparison of the beam (0.50 ×0.40 with PA of 88.8 • ) to the emission area with fluxes above half of the peak flux (0.50 ×0.41 with PA of 88.0 • ) indicates unresolved continuum emission from Th 28 at 0.90 mm. The mm-dust emission is therefore compact, in accordance with an emission originating from a disk.
As for the line results (see Sect. 3.1), to free ourselves of deconvolution uncertainties, we performed a fit of the disk parameters in the uv plane. Keeping in mind that the continuum appears unresolved in the image plane, it reveals a disk of ∼0.26 ×0.10 at PA=15±4 • (see Table 2), which is roughly compatible with those derived in 12 CO & 13 CO.
Mass of the disk: If we assume optically thin emission from dust grains at T =20 K, a dust opacity coefficient k ν = 0.1cm 2 .g −1 × (ν/1000GHz) β (very similar to that advocated by Hildebrand 1983, including a gas to dust ratio of 100 by mass) with β ∼ 0.5 (Beckwith & Sargent 1991), our total flux density of 13.3 mJy implies a total mass of ∼8.5×10 −4 M at the distance of 185 pc of Th 28. This circumstellar disk mass is relatively low but is not inconsistent with known CTTS disk masses (4×10 −3 M with a 0.72 dex standard deviation; Andrews & Williams 2005). Also, this value is highly sensible to the assumptions made on the parameters listed above. For instance a variation of −20% (resp. +20%) of the dust temperature implies an augmentation (resp. diminution) of the disk total mass to 1.2×10 −3 M (resp. 6.6×10 −4 M ), while varying the opacity index β by minus or plus 20% changes the total disk mass to 7.6×10 −4 M and 9.4×10 −4 M , respectively.

V lsr of the source
To determine the V lsr of Th 28, we fitted a single Gaussian component to the high-velocity wings of the 12 CO and 13 CO profiles (see top of Fig. 2 and Table 2). We derived a central velocity of V lsr ≈ 2.92 km.s −1 in 12 CO and V lsr ≈ 3.14 km.s −1 in 13 CO. In the following, we adopt a V lsr =3±0.1 km.s −1 for Th 28. This value is compatible with the range of V lsr observed toward the Lupus 3 cloud (James et al. 2006).

Disk inclination
The CO emission aspect ratio derived from the 2D-Gaussian fits in the uv plane provide estimates of the inclination to the line of sight of 73±1 • and 73±9 • from the 12 CO and 13 CO emission lines, respectively (see Table 2). This inclination is comparable but somewhat lower than the inclination of 82.3±2 • for the optical jet determined from the proper motions of distant knots (see Appendix). However, the CO emission may be underestimating the true disk inclination. Indeed, CO emissions are expected to arise above the disk plane at typical altitudes of 2-3 disk scale heights (Dartois et al. 2003). For a close to edge-on geometry, the apparent ellipse in the plane of the sky has a larger semiminor axis than a flat disk, mimicking a less inclined flat disk.

Disk rotation sense
The lower panels of Figure 2 show the first moment map of the 12 CO(3-2) and 13 CO(3-2) emission lines. A velocity gradient along the PA of the disk is striking from those plots. They reveal that the northeastern part of the disk is redshifted, while the southwestern part is blueshifted. We therefore report unambiguous detection of the sense of rotation of the Th 28 disk. The derived sense of rotation is opposite that derived for the optical jet with the HST (see Figure 2). We discuss the implications of this result for jet launching models in the following section.

12 CO versus 13 CO
To derive the overall flux ratio of the 13 CO emission line with respect to the 12 CO emission line, we integrated both lines (i) in the velocity range -15 km.s −1 to 20 km.s −1 and (ii) over the area where the 13 CO displays fluxes above 5σ (see Fig. 1-Middle). In so doing, we achieve integrated fluxes of ∼3.75 Jy.km.s −1 and ∼0.60 Jy.km.s −1 for the 12 CO and 13 CO, respectively. The integrated flux ratio, S( 13 CO)/S( 12 CO), is therefore of ∼0.16. The discrepancy between the integrated flux ratio (0.16) and the fractional abundance of CO isotopes ([ 13 CO]/[ 12 CO] ∼0.017 in the ISM) implies the 12 CO to be optically thick. Nevertheless, an optically thick 12 CO line profile supplies a valuable indication on the temperature profile in the disk (Beckwith & Sargent 1993, hereafter BS93); see Section 4.5.

Temperature profile (T(K) versus radius) in the disk
Since the 12 CO emission line is optically thick (see above), its spectrum is a direct probe of the temperature as a function of radius. Following the analysis developed by BS93, if the temperature profile follows a power law in radius (r), as T ∝ r −q , the flux density is roughly proportional to v 3q−5 obs in the high-velocity wings. Performing the fit of the flux density on the high-velocity blue wing [-8.0 km.s −1 ; -2.4 km.s −1 ] of the 12 CO emission line, we find an index q=0.8±0.2. Under these assumptions, the temperature profile of Th 28 is steeper than the typical values derived in CTTS disks (0.5-0.75; see Isella et al. 2009;Beckwith & Sargent 1991).

Discussion
Our ALMA observations of Th 28 reveal, for the first time, the millimeter CO and continuum emission from its accretion disk. We first discuss the implications brought by these observations for jet rotation studies. We then briefly discuss additional first order constraints on the Th 28 system.

Implications for jet rotation studies
Our high S/N ratio and well-behaved ALMA CO observations unambiguously establish that the Th 28 disk rotates in a sense opposite to that derived with the HST for both lobes of the Th 28 atomic optical jet by Coffey et al. (2004). Now, consistency between disk and jet rotation has been investigated in five systems. One was not informative, since the jet rotation sense could not be detected (RY Tau). Two of them (DG Tau, CW Tau) show consistent disk and optical jet rotation sense, while the two others (RW Aur and now Th 28) show inconsistent rotation sense. However, RW Aur-A is peculiar since its star shows periodic radial velocity variations possibly signaling a close companion (Petrov et al. 2001), which might confuse its jet dynamics, for example, through jet precession/wiggling or shock asymmetries (Cerqueira et al. 2006). On the contrary, the Th 28 jet was one of the best cases of jet rotation detection, with different optical lines and both lobes of the jet giving similar transverse velocity gradients. So, the discrepancy observed in that most favorable case between jet and disk rotation suggests broader implications for optical jet rotation studies. We discuss below three possible explanations for the observed discrepancy: 1) The disk rotation is probed on radial distances r ≥ 20 AU from the source, while the jet is expected to originate from much smaller disk radii, r < 5 AU ). One could imagine either a true change of rotation sense in the inner disk regions or a warp between the inner and outer disk that could induce an artificial change of rotation sense in the outer disk due to projection effects. Although we cannot fully exclude these possibilities, they seem unlikely as both the inclination and the PA of the outer CO disk axis agree well with that of the atomic jet, indicating no significant change of disk geometry between AU to tens of AU scales. In addition, the high inclination of the outer disk rotation axis means that line-of-sight velocities are very close to the underlying disk rotation velocities and that projection effects should be minimized. Moreover, the disk of Th 28 does not show any evidence of perturbation in its PV diagram at high velocities that could betray such behavior.
2) The observed transverse velocity gradients in jets do trace rotation but in a nonstationary part of the flow. Indeed, shocks and/or nonstationary ejection may significantly alter the apparent rotation of the flow, either amplifying it (Fendt et al. 2011) or even reversing it (Sauty et al. 2012) owing to the exchange of angular momentum between the flow and the magnetic field. Also, Staff et al. (2015) recently presented 3D magnetohydrodynamic simulations of disk winds where the less opened magnetic field configuration results in a wide, two-component jet. They show that kink mode creates a narrow corkscrew-like jet with regions in which jet rotation at ±2 AU from its axis is opposite to the disk. However such counter-rotating jet signatures would remain unresolved at the current achievable angular resolution of 0.1 (=14 AU) with the HST in the optical domain.
3) Optical transverse velocity gradients in jets do not trace rotation in the jet body but rather a departure from axisymmetry. Indeed, high angular resolution observations show evidence of small-scale jet axis wiggling and/or asymmetric shock fronts (see, e.g., White et al. 2014;Coffey et al. 2015). White et al. (2014) recently demonstrated in the DG Tau case that spurious rotation signatures on the same magnitude as those observed in the optical can be mimicked if the jet axis wiggling is not properly taken into account. Moreover, transverse velocity gradients traced by optical emission lines probe an intermediate velocity component (IVC) of the flow, the nature of which is still unclear. White et al. (2014) proposed that the IVC in DG Tau traces a dragged mixing layer between a fast inner jet and an outer slower wind. Therefore, transverse gradients detected in the NUV domain, which probe a faster and more collimated component of the flow, may trace jet rotation more reliably. The NUV velocity gradients detected in two sources so far (DG Tau and RW Aur) appear to rotate in the same sense as their underlying disks. However, RW Aur shows a time-variable behavior that renders it unsuitable for jet rotation studies. In the case of Th 28 no sense of rotation could be reliably determined in the NUV. Difficulties were encountered in this plane-of-sky jet because of a large, low-velocity absorption feature that closely coincides with the jet radial velocity. Further studies are clearly needed to establish whether observed NUV velocity gradients do indeed trace jet rotation.
In any case, the implications are that the velocity gradients derived from optical lines in T Tauri jets cannot be used to infer jet launching radius using the method proposed by Bacciotti et al. (2002); Anderson et al. (2003); Ferreira et al. (2006). Indeed, application of this method requires that the flow is both stationary and axisymmetric. Both assumptions seem inconsistent with the discrepancy found between the optical jet and disk rotation sense in Th 28.

Implications for the nature of Th 28
Comerón & Fernández (2010) derived a central stellar mass of 0.6-0.9 M from the maximum velocity of 450 km.s −1 observed in the redshifted wings of the Hα emission line profile, assuming that this velocity traces the free-fall velocity of infalling gas at the base of the accretion columns. In the following, we argue that the combination of the spatial extension and peak velocity separation observed in 12 CO indicates a stellar mass more in the range 1-2 M .
An accurate derivation of the central stellar mass would require a detailed disk modeling, with radiative transfer effects properly taken into account, and fits of the spectral energy distribution. Such a detailed analysis goes beyond the scope of the present paper. Instead, we propose a simple derivation of the central stellar mass standing on the 12 CO emission line. To meaningfully constrain the central stellar mass from the position-velocity diagram, it is critical to take projection and beam convolution effects into account. To do so, we build a simple geometrical model of an optically thick flat disk with a radial power law temperature distribution of index q and we convolve this model with the effective ALMA beam of the 12 CO observations. We fix the disk inclination to i=80 • and the power law index q of the temperature profile to q=0.8, which is most compatible with the shape of the wings in the 12 CO profile (see Sect. 4.5). In order to reproduce the observed 12 CO brightness profile along the disk major axis, a disk outer radius of R out =0.5 is required, corresponding to 92.5 AU at the adopted distance of d=185 pc for Th 28. The observed velocity separation between the peaks is then best fitted with a stellar mass of 1.6 M . The model described above depends on the index of the temperature profile, q, whose value is uncertain by 25 % (see Sect. 4.5). Spreading this uncertainty in the model, we derive an uncertainty on the mass of the central star of ±0.2 M . Moreover, a central stellar mass in the range 1-2 M appears more compatible with a spectral type of early K as quoted in the literature for this object. Using the pre-main-sequence tracks of Siess et al. (2000), a spectral type of K2 (T e f f = 4900-5000 K) corresponds to M star ≥ 2 M for ages less than 3.5 Myrs.
Concerning the low luminosity of the central object, our best fit to the CO emission profiles gives V lsr =3±0.1 km.s −1 , which is consistent with the range of values observed in the Lupus 3 cloud (James et al. 2006). It is then unlikely that the low luminosity of the central Th 28 star could be related to a larger distance. One would expect before extinction a ∼1.9 L luminosity from a 1.5 M protostar of age ≤ 3.5 Myr (Siess et al. 2000) 3 . Therefore an extinction of 4.5 mag would be sufficient to attenuate the star down to the luminosity observed for Th 28 (0.03 L , see Mortier et al. 2011). Such an extinction appears easily achievable from a close to edge-on system and our ALMA observations confirm a large inclination to the line of sight (i ≥ 73 • ). This is also consistent with the fast proper motion of the jet knots at the distance of Lupus 3 when compared to their radial velocities (see Sect. 4.2).
• The 12 CO emission is clearly resolved and elongated along PA=7.3 • . The morphology of the disk seen in CO matches very well the morphology derived for the large-scale atomic jet. The PA of the disk is perpendicular to that of the jet, and both inclinations are comparable. • The 12 CO shows kinematics that are consistent with a disk in Keplerian rotation. Indeed its double-peak integrated profile, plus its PV diagram along the PA of the disk, are characteristic of Keplerian rotation. • We derive a V lsr =3±0.1 km.s −1 for the central source in agreement with the range of values observed in the Lupus 3 cloud.

see also Webpage Siess
• We constrain the power law index, q, of the temperature distribution to q 0.8 and the R out to 95 AU. • The combination of large R out and peak velocity separation in 12 CO of ∆v=7 km.s −1 is best reproduced with a central stellar mass of 1.4-1.8 M . This is also consistent with early K spectral type estimates for this source. • The rotation sense of the disk is well detected with our ALMA CO observations and this direction is opposite that of the transverse velocity shifts previously detected with HST in the optical jet. This discrepancy in rotation senses implies that velocity gradients measured in optical lines cannot be used to infer launching radii. The NUV domain, which probes a faster and more collimated inner part of the jet, is likely more reliable to measure jet rotation from Doppler gradients.