---
author:
- 'Á. Sánchez-Monge'
- 'R. Cesaroni'
- 'M. T. Beltrán'
- 'M. S. N. Kumar'
- 'T. Stanke'
- 'H. Zinnecker'
- 'S. Etoka'
- 'D. Galli'
- 'C. A. Hummel'
- 'L. Moscadelli'
- 'T. Preibisch'
- 'T. Ratzka'
- 'F. F. S. van der Tak'
- 'S. Vig'
- 'C. M. Walmsley'
- 'K.-S. Wang'
date: 'Received date; accepted date'
title: |
     A candidate circumbinary Keplerian disk in 35:\
    A study with ALMA 
---

Introduction {#sint}
============

While different scenarios have been proposed to explain the formation of high-mass (i.e. OB-type) stars (monolithic collapse in a turbulence-dominated core – Krumholz et al. [@krumholz2009]; competitive accretion driven by a stellar cluster – Bonnell & Bate [@bonnellbate2006]; Bondi-Hoyle accretion – Keto [@keto2007]; see also the review by Zinnecker & Yorke [@ziyo]), all of them predict the formation of circumstellar disks. It is thus surprising that only a handful of disk candidates have been observed in association with massive (proto)stars. As a matter of fact, despite many observational efforts, convincing evidence of disks has been found only around early B-type (proto)stars, while circumstellar disks around O-type stars remain elusive (Wang et al. [@wang2012]; Cesaroni et al. [@cesaroni2007] and references therein). Moreover, a detailed investigation of the disk properties, comparable to that performed in disks around low-mass stars (e.g. Dutrey et al. [@dutrey2007]) is still missing due to the large distances of OB-type (proto)stars and the limited angular resolution at (sub)millimeter wavelengths. With the advent of the Atacama Large Millimeter Array (ALMA) the situation is bound to improve dramatically, as resolutions $\ll$1 will be easily obtained.

With this in mind, we performed ALMA Cycle 0 observations of two IR sources containing B-type (proto)stars. These were chosen on the basis of their luminosities (on the order of $10^4~L_\odot$), presence of bipolar nebulosities/outflows, detection of broad line wings in typical jet/outflow tracers (SiO), and strong emission in hot molecular core (HMC) tracers (such as methyl cyanide, ). Here we present the most important results obtained for one of the two sources, 35.

35 is a well known star forming region located at a distance of $2.19_{-0.20}^{+0.24}$ kpc (Zhang et al. [@zhang2009]), with a luminosity of $\sim$$3\times10^4~L_\odot$[^1]. The region is characterized by the presence of a butterfly shaped reflection nebula oriented NE–SW (see Fig. \[flarge\]), as well as a bipolar molecular outflow in the same direction, observed in  by Dent et al. ([@dent1985a]), Gibb et al. ([@gibb2003]; hereafter GHLW), Birks et al. ([@birks2006]; hereafter BFG), and López-Sepulcre et al. ([@lopezsepulcre2009]). The (1–0) line emission appears to trace also a N–S collimated flow (see BFG), coinciding with a thermal radio jet (Heaton & Little [@heatonlittle1988]; GHLW) seen also at IR wavelengths (Dent et al. [@dent1985b]; Walther et al. [@walther1990]; Fuller et al. [@fuller2001]; De Buizer [@debuizer2006]; Zhang et al. [@zhang2013]). It has been proposed that the poorly collimated NE–SW outflow and the N–S jet could be manifestations of the same bipolar flow undergoing precession (Little et al. [@little1998]). However, evidence for multiple outflows in this region is provided by SiO, , , and  line observations (GHLW; Lee et al. [@lee2012]).

A molecular clump elongated perpendicular to the NE–SW outflow has been mapped in dense gas tracers (, CS), whose emission exhibits a velocity gradient from NW to SE (Little et al. [@little1985]; Brebner et al. [@brebner1987]). This was first interpreted as a large ($\sim$1 or 0.6 pc) disk/toroid rotating about the NE–SW outflow axis, but GHLW, on the basis of their  and  observations, propose that this is actually a fragmented rotating envelope containing multiple young stellar objects (YSOs). Indeed, GHLW identify a core at the center of the outflow and another core, named G35MM2, offset to the SE.

Observations and results {#sobs}
========================

35 was observed with ALMA Cycle 0 at 350 GHz in May and June 2012, with baselines in the range 36–400 m, providing sensitivity to structures $\le$2 . The digital correlator was configured in 4 spectral windows (with dual polarization) of 1875 MHz and 3840 channels each (covering the ranges 334.85–338.85 GHz and 346.85–350.85 GHz), providing a resolution of $\sim$0.4 km s$^{-1}$. Flux, gain, and bandpass calibrations were obtained through observations of Neptune and J1751$+$096. The data were calibrated and imaged using CASA. A continuum map was obtained from line-free channels and subtracted from the data. The synthesized beam is $0\farcs51\times0\farcs46$, P.A.=48. The rms noise is $\sim$6 mJy beam$^{-1}$ for individual line channels, while in the continuum image it is $\sim$1.8 mJy beam$^{-1}$, implying a S/N of only $\sim$100. The latter indicates a reduced dynamic range.

In Fig. \[flarge\], we present the map of the 350 GHz continuum emission overlaid on an enhanced resolution Spitzer/IRAC image at 4.5  extracted from the GLIMPSE survey (Benjamin et al. [@benj]). The sub-millimeter continuum emission is clearly tracing an elongated structure across the waist of the butterfly shaped nebula. In all likelihood, we are detecting the densest part of the flattened molecular structure observed on a larger scale by Little et al. ([@little1985]), Brebner et al. ([@brebner1987]), and GHLW.

Along the elongated structure a chain of at least 5 cores is seen (Fig. \[fcont\]), lending support to GHLW’s idea that one is dealing with a fragmented structure instead of the smooth disk/toroid hypothesized by Little et al. ([@little1985]). We stress that the angular resolution of our maps ($\sim$7 times better than previous (sub)millimeter observations) reveals that the YSOs powering the outflow(s) lie inside cores A and/or B (see Fig. \[fcont\]), because these two are the only cores located close to the geometrical center of the bipolar nebula. In particular, core B lies along the N–S jet traced by the IR and radio emission, and coincides with one of the free-free sources detected by GHLW. The hypothesis by BFG that core G35MM2 could be driving the NE–SW outflow is not convincing, as this core is off-center, lying right at the border of the waist outlined by the IR emission.

Methyl cyanide emission (as well as other typical HMC tracers) is clearly detected only towards cores A and B, and marginally towards G35MM2. Emission by vibrationally excited lines of  also indicates that cores A and B could be hosting massive stars. Core B coincides also with a compact free-free continuum source detected by GHLW at 6 and 3.6 cm, and by Codella et al. ([@codella2010]) at 1.3 cm. This emission could be part of the N–S thermal radio jet or might be coming from an  region ionized by an embedded early-type star. We will discuss this possibility in Sect. \[sdis\]. Faint radio emission at a $5.6\sigma$ level is detected also towards core A, consistent with the presence of embedded star(s).

We now investigate the gas velocity field in the two cores by computing the first moment of a  line, a dense gas tracer. Figure \[fvelo\] plots the result for the (19–18) $K$=2 line, with overlaid the line emission averaged over the same velocity interval used to calculate the first moment. Note that the mean velocity of core B ($\sim$30 ) differs by $\sim$2  from that of core A ($\sim$32 ), which also differs by a similar amount from that ($\sim$34 ) obtained with lower angular resolution observations of other tracers (see e.g. GHLW). Discrepancies of this type are often found in high-mass star-forming regions (see, e.g., the case of AFGL5142 – Estalella et al. [@estalella1993]; Zhang et al. [@zhang2007]) and are likely related to the fragmentation process in molecular clumps. The interesting result is that both cores present velocity gradients, which may be due to expansion or rotation of the gas. Note that the sense of the gradient in core B is the opposite of that measured on larger scales (see Fig 5 of Little et al. [@little1985] and Figs. 8 and 9 of GHLW). However, the velocity range traced by the extended emission (30–41 ) differs significantly from that sampled by the  lines ($\sim$26–38 ), suggesting that we are observing cores whose velocity field is not tightly related to that of the largest, fragmenting structure where they are embedded. In the following, we will address this question for core B, which is better resolved and associated with the N–S jet/outflow.

The nature of core B: A Keplerian rotating disk? {#skep}
================================================

An elongated structure with a velocity gradient along its major axis, such as that observed in core B, can be interpreted as either a collimated bipolar outflow, or a rotating ring/disk seen close to edge on. We thus first investigate the possibility that we are observing the root of a large-scale bipolar outflow.

At first glance, both the orientation (see Fig. \[fcont\]) and sign of the velocity gradient seen in Fig. \[fvelo\] appear consistent with those of the N–S jet imaged by BFG (see the red-shifted emission at 40  in their Fig. 3). However, the CO emission in the northern jet lobe is detected by BFG also at 28  (see their Fig. 3), namely at blue-shifted velocities. Moreover, the direction of the  velocity gradient has a negative position angle, whereas the jet is slightly inclined to the east. Such a difference could be explained if the jet is precessing (as suggested by BFG) and lies close to the plane of the sky, because the jet direction would change from the small to the large scale. We believe that this explanation is not satisfactory, because the orientation of the  velocity gradient differs from that of the radio jet [*on the same scale*]{}. This can be seen in Fig. \[fvelo\], where the two free-free emission sources are aligned N–S, whereas the velocity gradient (and the core major axis) are inclined to the west (PA$\simeq$–20).

Based on these findings, we consider the alternative possibility that the velocity gradient is tracing a disk. To analyse the kinematics of the gas in better detail, we fitted the  emission in each 0.4  velocity channel with a 2-D Gaussian. This allows us to derive the peaks of the line emission at different velocities. Then all peaks can be plotted together to obtain a picture of the gas velocity field. This is done in Fig. \[fkepfit\]a, where for each peak we also plot the corresponding 50% contour level to give an idea of the size of the gas emitting at that velocity.

Two considerations are in order. The first is that the peaks outline a sort of elliptical pattern, roughly centered on the position of the continuum peak and oriented SSE–NNW. The second is that the most blue- and red-shifted peaks tend to converge towards the position of the sub-mm continuum peak. On this basis, one is tempted to hypothesize that the  emission is tracing Keplerian-like rotation, as this would explain why both the emission at systemic velocities and that at high (blue- and red-shifted) velocities are observed towards the rotation center.

![image](peaks_conts-kep_fit_port.ps){width="14cm"}

With this in mind, we have fitted the observed velocity pattern assuming Keplerian rotation about a point-like source. We stress that ours is a purely kinematical fit, where only the rotation velocity in a 2-D disk is calculated, and no estimate of the observed line intensity is computed. Also, no attempt is made to demonstrate the uniqueness of the model proposed and we envisage the possibility that other interpretations might be possible. Our simple model is sufficient to constrain a number of important parameters. The inputs of the model are: the LSR velocity of the system (), the central point mass () and its position ($x_0$,$y_0$), the inclination angle ($\theta$) of the rotation axis with respect to the plane of the sky, and the position angle ($\psi$) of the projected major axis of the disk. The best fit was obtained by varying all input parameters inside reasonable ranges and minimizing the expression $\sum_i{\left({V_i^{\rm observed}}-V_i^{\rm model}\right)^2}$ where $i$ indicates the peak in a generic velocity channel. For our calculations, we have used the $K$=2, 3, and 4 components of the (19–18) transition as well as other prominent lines of other species, such as ($7_{1,6}$–$6_{1,5}$) $v_{\rm t}$=1, ($14_{1,13}$–$14_{0,14}$), and HC$_3$N(37–36). All of these lines reveal a kinematical pattern very similar to that in Fig. \[fkepfit\]a, demonstrating that such a pattern does not depend on the tracer and mirrors a real physical structure. In Fig. \[fkepfit\]b we show an overlay of the velocity peaks of all the lines and a velocity map of the best-fit model. Clearly, the agreement between the computed and observed LSR velocities is remarkable. The best-fit parameters are the following: =30.0$\pm$0.3 , $x_0$=18$^{\rm h}$58$^{\rm m}$13027$\pm$0002, $y_0$=01403594$\pm$007, $\psi$=157$\pm$4, $\theta$=19$\pm$1, =18$\pm$3 .

A lower limit for the disk radius is $R_{\rm disk}\simeq2500$ AU, obtained from the maximum deprojected distance of the peaks from the center. One may wonder if such a big structure can undergo Keplerian rotation. The mass of the disk can be estimated from the continuum emission. The integrated flux density from core B at 350 GHz is 0.32 Jy. Assuming a dust absorption coefficient 0.5 cm$^2$g$^{-1}~(\nu/230.6~{\rm GHz})$ (Kramer et al. [@kramer2003]) and a gas-to-dust mass ratio of 100, we obtain $\sim$3  for a dust temperature of 100 K. Despite the large uncertainties on the dust opacity and temperature, we believe that the mass of core B is significantly less than the central mass ($\sim$18 ), which satisfies the condition for Keplerian rotation.

It is worth stressing that our findings are in good agreement with the recent study by Ilee et al. ([@ilee]). Through measurements of scattered light from 35, these authors find that the CO first overtone bandhead emission at 2.3  can be fitted with a Keplerian disk rotating about a 17.7  star.

The distribution of the molecular peaks in Fig. \[fkepfit\]b clearly shows that our observations detect only the NE side of the disk. We speculate that this could be an opacity effect. If the disk is optically thick in the relevant lines, flared, and inclined by 19, only part of the surface heated by the star is visible to the observer. This creates an asymmetry along the direction of the projected disk axis, with line emission being more prominent on the side (in our case the NE side) where the disk surface is visible. Clearly, radiative transfer calculations are needed to confirm this scenario, but we note that in our source the NE part of the disk axis should be pointing towards the observer, consistent with the orientation of the CO outflow (blue shifted to the NE and red shifted to the SW - see GHLW) and the obscuration seen to the SW in the IR images (see Fig. 2).

The stellar content of core B: A binary system? {#sdis}
===============================================

An issue that is worth discussing is whether an 18  YSO is compatible with the bolometric luminosity of the region. ($3\times10^4~L_\odot$; see Sect. \[sint\]). Depending on the adopted mass–luminosity relation, the luminosity expected for an 18  star ranges from $2.5\times10^4~L_\odot$ (Diaz-Miller et al. [@diazmiller1998]), to $6.6\times10^4~L_\odot$ (Martins et al. [@martins2005]). This means that the 18  star should be the main contributor to the luminosity of the whole star forming region. Such a possibility seems quite unlikely due to the presence of multiple cores (see Fig. \[fcont\]), one of which is an HMC possibly hosting at least another high-mass star (core A).

A possibility is that one is underestimating the true luminosity due to the “flashlight effect”, where part of the stellar photons are lost through the outflow cavities. According to the recent model by Zhang et al. ([@zhang2013]), when this effect is taken into account, the luminosity obtained assuming isotropic emission ($3.3\times10^4~L_\odot$) becomes as large as $7\times10^4$–$2.2\times10^5~L_\odot$ consistent with a single star of $\sim$$20$–$34~M_\odot$.

While the previous explanation is possible, the isotropic estimate appears more robust than a model-dependent value, and we thus consider another hypothesis, namely that one is dealing with a binary system. In this case, the luminosity is significantly reduced with respect to that of a single 18  object and may be as low as $\sim7\times10^3~L_\odot$ for equal members. The latter is much less than the bolometric luminosity, thus allowing for a significant contribution by other YSOs.

The existence of a binary system could also explain why the N–S jet associated with core B is not aligned with the disk rotation axis. The presence of a companion would in fact be sufficient to induce precession of the jet/outflow, as hypothesized by Shepherd et al. ([@shepherd2000]) to explain the observed precession of the bipolar outflow from the high-mass protostar . In this scenario, the outflow from core B would precess about an axis oriented NE–SW, i.e. along the bisector of the butterfly shaped IR nebula seen in Fig. \[flarge\]. The IR emission would arise from the cavity opened by the outflow itself during its precession, while the thermal radio jet would trace the current direction of the precessing axis.

The last question we address is the origin of the free-free emission from core B (see e.g. Fig. \[fvelo\]). Could this be tracing a (hypercompact)  region? According to GHLW, this source (n. 7 in their Table 1) has a spectral index $>$1.3 between 6 and 3.6 cm, compatible with free-free emission from an optically thick  region. Extrapolation of the 3.6 cm flux density (0.5 mJy) to 1.3 cm gives $>$1.9 mJy, in agreement with the 3 mJy flux measured, with lower angular resolution, by Codella et al. ([@codella2010]). It is hence possible that the emission is partially thick at 1.3 cm and the Lyman continuum estimate of $10^{45}~{\rm s}^{-1}$ obtained under the optically thin assumption is a lower limit. A binary system with a total mass of 18  has a Lyman continuum flux (see Diaz-Miller et al. [@diazmiller1998]) ranging from $5\times10^{44}$ s$^{-1}$ (for equal masses) to $1.8\times10^{47}$ s$^{-1}$. We conclude that we could be observing an  region ionized by a binary system at the center of a Keplerian disk. Should this be confirmed, 35 would represent a unique example of radio jet coexisting with an  region powered by the same YSO(s). Investigating such a short-lived transition phase in the evolution of an OB-type star could provide us with important clues on the formation process of these objects.

It is a pleasure to thank Göran Sandell for stimulating discussions on the 35 region and the anonymous referee for his/her constructive criticisms. We also acknowledge the support of the European ALMA Regional Center and the Italian ARC node. This paper makes use of the following ALMA data: ADS/JAO.ALMA\#2011.0.00275.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada) and NSC and ASIAA (Taiwan), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. This work also used observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. MSNK is supported by a Ciência 2007 contract, funded by FCT (Portugal) and POPH/FSE (EC).

Benjamin, R.A., Churchwell, E., Babler, B.L., et al. 2003, PASP, 115, 953 Birks, J. R., Fuller, G. A., & Gibb, A. G. 2006, , 458, 181 (BFG) Bonnell, I. A., & Bate, M. R. 2006, , 370, 488 Brebner, G. C., Heaton, B., Cohen, R. J., & Davies, S. R. 1987, , 229, 679 Cesaroni, R., Galli, D., Lodato, G., Walmsley, C.M., & Zhang, Q. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. of Arizona Press), 197 Codella, C., Cesaroni, R., L[ó]{}pez-Sepulcre, A., et al. 2010, , 510, A86 De Buizer, J. M. 2006, , 642, L57 Dent, W. R. F., Little, L. T., Kaifu, N., Ohishi, M., & Suzuki, S. 1985, , 146, 375 Dent, W. R. F., Little, L. T., Sato, S., Ohishi, M., & Yamashita, T. 1985, , 217, 217 Diaz-Miller, R. I., Franco, J., & Shore, S. N. 1998, , 501, 192 Dutrey, A., Guilloteau, S., & Ho, P. 2007, in Protostars and Planets V, ed. B. Reipurth, D. Jewitt, & K. Keil (Tucson: Univ. of Arizona Press), 495 Estalella, R., Mauersberger, R., Torrelles, J. M., et al.  1993, , 419, 698 Fuller, G. A., Zijlstra, A. A., & Williams, S. J. 2001, , 555, L125 Gibb, A. G., Hoare, M. G., Little, L. T., & Wright, M. C. H. 2003, , 339, 1011 (GHLW) Heaton, B. D., & Little, L. T. 1988, , 195, 193 Keto, E. 2007, , 666, 976 Kramer, C., Richer, J., Mookerjea, B., Alves, J., & Lada, C. 2003, , 399, 1073 Krumholz, M. R., Klein, R. I., McKee, C. F., Offner, S. S. R., & Cunningham, A. J. 2009, Science, 323, 754 Ilee, J.D., Wheelwright, H.E., Oudmaijer, R.D. et al. 2013, MNRAS, in press Lee, H.-T., Takami, M., Duan, H.-Y., et al. 2012, , 200, 2 Little, L. T., Dent, W. R. F., Heaton, B., Davies, S. R., & White, G. J. 1985, , 217, 227 Little, L. T., Kelly, M. L., & Murphy, B. T. 1998, , 294, 105 L[ó]{}pez-Sepulcre, A., Codella, C., Cesaroni, R., Marcelino, N., & Walmsley, C. M. 2009, , 499, 811 Martins, F., Schaerer, D., & Hillier, D. J. 2005, , 436, 1049 Shepherd, D. S., Yu, K. C., Bally, J., & Testi, L. 2000, , 535, 833 Velusamy, T., Marsh, K. A., Beichman, C. A., Backus, C. R., & Thompson, T. J. 2008, , 136, 197 Walther, D. M., Aspin, C., & McLean, I. S. 1990, , 356, 544 Wang, K.-S., van der Tak, F. F. S., & Hogerheijde, M. R. 2012, , 543, A22 Zhang, B., Zheng, X. W., Reid, M. J., et al. 2009, , 693, 419 Zhang, Y., Tan, J. C., De Buizer, J. M., et al. 2013, , in press Zhang, Q., Hunter, T. R., Beuther, H., et al. 2007, , 658, 1152 Zinnecker, H. & Yorke, H. W. 2007, ARA&A, 45, 481

[^1]: Estimate taken from the RMS database, available at http://www.ast.leeds.ac.uk/cgi-bin/RMS/RMS\_DATABASE.cgi
