Sub-0.1″ optical imaging of the Z CMa jets with SPHERE/ZIMPOL
Key Words.:stars: variables: Herbig Ae/Be, FU Ori – stars: winds, jets – stars: individual: Z CMa – circumstellar matter – stars: pre-main sequence – techniques: high angular resolution imaging
Context:Crucial information on the mass accretion-ejection connection in young stars can be obtained from high spatial resolution images of jets in sources with known recurrent accretion outbursts.
Aims:Using the VLT/SPHERE ZIMPOL instrument, we observed the young binary Z CMa that is composed of a Herbig Be star and a FUor object, both driving a jet. We aim to analyse the structure of the two jets, their relation with the properties of the driving sources, and their connection with previous accretion events observed in this target.
Methods:We obtained optical images in the H and [O i] 6300 lines at the unprecedented angular resolution of 0.03 arcsec, on which we have performed both continuum subtraction and deconvolution, thereby deriving results that are consistent with each other.
Results:Our images reveal extended emission from both sources: a fairly compact and poorly collimated emission SW of the Herbig component and an extended collimated and precessing jet from the FUor component. The compact emission from the Herbig star is compatible with a wide-angle wind and is possibly connected to the recent outburst events shown by this component. The FUor jet is traced down to 70 mas (80 AU) from the source and is highly collimated with a width of 26-48 AU at distances 100-200 AU, which is similar to the width of jets from T Tauri stars. This strongly suggests that the same magneto-centrifugal jet-launching mechanism also operates in FUors. The observed jet wiggle can be modelled as originating from an orbital motion with a period of 4.2 yr around an unseen companion with mass between 0.48 and 1 M. The jet mass loss rate was derived from the [O i] luminosity and comprises of between and M yr. This is the first direct measurement from a jet in a FUor. If we assume previous mass accretion rate estimates obtained through modelling of the accretion disk, the derived range of would imply a very low mass-ejection efficiency (/ ), which is lower than that typical of T Tauri stars.
Star formation occurs through a mechanism that couples accretion from a disk and ejection of matter in the form of winds and collimated jets, which concur in removing angular momentum from the system. The accretion process is considered to be quasi-stationary, but time-variable accretion bursts caused by a quick rise of the mass accretion rate from the disk are sometimes observed. Studies of jets from sources that show episodic accretion events allow one to trace the past history of the accretion bursts, as one may expect to detect spatial structures within the jet that can be connected with particularly intense mass-ejection events. Fundamental insights come from the inner regions close to the source, where the jets are less contaminated by interactions with the ambient medium, so that their morphology can be tightly connected with variations of the ejection efficiency.
The Z CMa binary is a peculiar system that is suited to investigate the occurrence of mass ejection in young sources undergoing episodic accretion. This binary is composed of a Herbig Be star (the NW component), which shows EXor-like outbursts (Szeifert et al. 2010), and a second component (SE) that exhibits broad double-peaked optical absorption lines that are typical of a FU Ori object (Hartmann et al. 1989; Hartmann & Kenyon 1996). The two sources, located at a distance of 1150 pc (Herbst et al. 1978), are separated by 0 1 and only recently high angular resolution observations with adaptive optics (AO) systems and interferometers have allowed the separatele determination of the properties of the two sources (Bonnefoy et al. 2016, submitted; Hinkley et al. 2013). The Herbig star, which is responsible for the variability of the binary (Bonnefoy et al. 2016, submitted), has shown recurrent variations in the visual magnitude of 2-3 mag in the last ten years with recent major events recorded in 2008, 2011, and 2015 (see light curves in Canovas et al. 2012; Maehara & Ukita 2015). These variations have been imputed to both variable accretion and to changes of the extinction along the line of sight (Benisty et al. 2010; Szeifert et al. 2010).
An extended optical jet from Z CMa was originally detected by Poetzel et al. (1989). Recently, Whelan et al. (2010, hereafter W10) discovered that both sources drive their own jet, at slightly different position angles (P.A.): 245 (east of north) for the Herbig and 235 for the FU Ori. The jet from the NW source has a very high radial velocity reaching up to 600km s while the jet from the SE component has a velocity of 200/300 km s. The large-scale optical jet observed by Poetzel et al. (1989) is usually associated with the Herbig component (e.g. Garcia et al. 1999, W10). Additional emission that is not collimated well is observed at intermediate velocities, which cannot however be unambiguously associated with one of the two components (W10). The analysis of these jets offers us the opportunity to obtain direct information on the mass loss and indirect information on the accretion rate. This is especially interesting in the case of the FUor component, since no direct mass loss measurements from jets are available for such objects. Wind mass loss rates as high as M yrwere estimated on the prototype source FU Ori from line profiles (Croswell et al. 1987; Calvet et al. 1993). Accretion rates estimates for FUor objects are also very scarce and range from to M yr(Audard et al. 2014).
We present here VLT/SPHERE (Beuzit et al. 2008) high angular resolution optical observations of Z CMa in the [O i]6300Å and H lines. The extreme AO performance of SPHERE allows us to have direct images of the Z CMa system at unprecedented high contrast and spatial resolution (0 03).
2 Observations and data reduction
Observations of Z CMa in the H and [O i] lines were performed using the ZIMPOL instrument
(Thalmann et al. 2008; Schmid et al. 2012) of SPHERE on the night of 31 March 2015.
Simultaneous images of Z CMa in H and in the adjacent continuum were acquired with the two ZIMPOL cameras using
the B_Ha (=655.6 nm, =5.5 nm) and Cnt_Ha (=644.9 nm, =4.1 nm) filters. Exposures were taken in field-stabilized mode using two different field rotation angles and swapping the filters in front of each camera to optimize artefact removal during the reduction. This resulted in 60 frames per filter (15 for each setup) for a total exposure time of 30 minutes.
Images in the [O i] line with the OI_630 filter
The raw data were processed using the SPHERE-ZIMPOL IDL pipeline (version 1.1) developed at ETH Zürich, which performs bias, dark, flat-field correction, re-centring of the dithered frames, and allows for derotation, selection, subtraction, and average or median combination of the frames. The Cnt_Ha exposures were used to remove the continuum emission from the B_Ha and OI_630 filter images, so as to produce continuum-subtracted H and [O i] images in which the line-emitting structures around the components are evidenced. Details about the reduction procedure for the continuum subtraction are given in Appendix A. From the FWHM of the PSF of the two stars measured in the Cnt_Ha frames we estimate an effective angular resolution of 30 mas for our images. The binary separation and position angle are analysed in Appendix B.
We also performed an alternative processing of the ZIMPOL images with the multi-component Richardson-Lucy (MC-RL) deconvolution method developed by La Camera et al. (2014), which is optimized for the reconstruction of high dynamical range images, such as those of jets from young stars (e.g. Antoniucci et al. 2014). The MC-RL method is able to separately reconstruct the star(s) and the underlying diffuse emission in the image. We extracted the PSF from the Cnt_Ha image via the blind deconvolution algorithm described in Prato et al. (2013, 2015) and employed to reconstruct both the [O i] and H images.
3 Results and discussion
The final continuum-subtracted [O i] and H images are shown in the upper part of Fig. 1a,b. We can recognise two major structures: an extended collimated jet from the FUor component and a more compact and poorly collimated emission SW of the Herbig star. These structures are especially evident in the [O i] image, which is much less affected by subtraction residuals than the H image. In our [O i] continuum-subtracted image we can trace the flows down to about 90 mas from the sources (corresponding to 105 AU). The morphology of both the FUor jet and the compact emission from the Herbig is thoroughly confirmed by the deconvolved images (lower part of Fig. 1a,b), which represent the diffuse emission in the image as reconstructed by the MC-RL algorithm. The [O i] deconvolved image shows details of the FUor jet even closer to driving source, down to 70 mas (80 AU). The deconvolved H image provides evidence of the non-collimated emission SW of the Herbig star, which was not fully visible in the continuum-subtracted image owing to the large residuals.
3.1 Outflow from the Herbig Be star
The extended [O i] and H emission detected SW of the Herbig star is fairly compact and poorly collimated, suggesting an origin in a wide-angle wind more than in a jet. The flow has an extent of 012, corresponding to 140 AU, while the observed opening angle is with a median P.A. around 230. In [Fe ii] KECK observations by W10, collimated emission travelling at velocities up to 600 km sand with a P.A.245 is detected up to distances of 04 from the source, along with a more diffuse emission at lower velocities (100-200 km s). The wide-angle compact emission that we observe is likely connected to this diffuse emission, while the high-velocity collimated jet is not visible in our images. Although the Herbig jet is visible at high velocities in the [Fe ii] data reported by W10 (which are not flux calibrated), the flux contrast between the FUor and Herbig jet is indeed greater than a factor 20 for velocities between 400 and 90 km s. This also agrees with the non-detection of the jet in the SINFONI data shown by W10. As we can estimate a mean signal-to-noise ratio of about 12 on the FUor jet, we conclude that we did not observe the collimated jet from the Herbig because of sensitivity limits.
Interferometric measurements of the Br line during the 2008 outburst (Benisty et al. 2010) suggested a bipolar (non fully spherical) wind from the star on mas scales, roughly aligned in the direction of the large-scale jet. As the Br emission disappeared outside the outburst phase, the authors concluded that the wind emission is connected to the event of increased mass accretion driving the outburst. To understand if the emission detected in our images is related to such bipolar wind, we evaluated whether the signal may be attributed to an ejection event from the recent Herbig star bursts, which occurred in January 2008 and February 2011 (Canovas et al. 2012), and more recently in January 2015 (see photometry from the Kamogata Wide-field Survey Maehara & Ukita 2015). Given the observed extent of the flow, we find that it should travel at a tangential velocity of 90, 170, and 2700 km sto have been ejected from the 2008, 2011, or 2015 events, respectively. Ruling out the latter event, which implies an exceedingly large velocity, a connection with the 2008 and 2011 events is possible.
3.2 Jet from the FUor
The collimated jet from the FUor object can be identified with the Jet B observed in [Fe ii] by W10, which travels at radial velocities between 100 and 400 km s. The morphology of the jet in our [O i] image is however much more structured than in the [Fe ii] images. In the continuum-subtracted image, we can identify three peaks on the brightest part of the jet, which we name A, B, and C, while we denote with D a more diffuse nebulosity located farther out (Fig.1a,b). The section of the jet between knots A and C shows a very apparent wiggling, which we analyse in detail in Sect. 3.2.1. Outside the ring of enhanced speckle noise (with radius 03), we also detect a weak emission knot that we identify as the knot K1 of W10. As in the [Fe ii] images, this knot is not aligned with the rest of the jet although it is probably connected with it. The knot shows an appreciable proper motion with respect to the W10 images. In the assumption that [Fe ii] and [O i] trace the same gas, we measure an increase of 186 mas for the knot distance from the source in the two images, which were taken 5.3 years apart. This corresponds to a tangential velocity of 192 km s. Hereafter we assume that the rest of the jet has the same tangential velocity. With this assumption, we can infer the times at which the A-D knots were ejected, obtaining dynamical timescales of 3.8, 5.4, 6.0, and 9.3 yr for knots A, B, C, and D, respectively. We therefore expect that only knot D should be visible in the W10 [Fe ii] image. Indeed, only one knot is clearly identifiable in that image and its dynamical timescale is consistent with it being knot D of our image. The derived tangential velocity combined with the observed radial velocities between 100 and 400 km s implies a jet inclination angle with respect to the plane of the sky in a range between 28° and 64°.
From the observed wiggling we can estimate a mean P.A. of the jet of 236.5, which approximately intersects the knot B, and a variation of about 3 around such direction. We analysed the jet wiggling in our [O i] image by considering 12 contiguous 3-pixel-wide slices orthogonal to the jet axis covering the section between knots A and C, where we have the best signal to noise. We averaged the counts within each slice to obtain the spatial profiles of the jet at various distances. Using a linear combination of a second degree polynomial (to fit the local continuum) and a Gaussian function, we computed the profile peak positions as a function of the distance from the star, which are reported in Fig. 1c for both the subtracted and deconvolved [O i] image.
The jet profiles are barely resolved with fitted full width and half maxima (FWHM) of the Gaussians in the continuum-subtracted image spanning between 35 and 50 mas with a general trend of increasing width with distance (Fig. 1c). When corrected for the instrumental profile, these indicate a jet intrinsic FWHM in the range 26-48 AU at distances 100-200 AU from the star; these valuesare in substantial agreement with widths observed in Class II jets and also in less evolved sources (Cabrit et al. 2007) and testifies a high degree of collimation of the jet. This finding is an important indication that the same magneto-centrifugal jet-launching mechanism is probably at work in objects with high accretion rates such as FUors with their massive peculiar disks.
With a jet tangential velocity of 192 km s, the observed wiggling occurs on timescales of a few years, which suggests an origin due to the interaction with a close undetected companion. We therefore explored the hypothesis that the wiggling arises from the orbital motion of the jet-emitting source around this companion. We adopt the formulation of Anglada et al. (2007), which considers a ballistic jet from a star in a circular orbit with a constant ejection velocity orthogonal to the orbital plane. Using Eq. 8 of Anglada et al. (2007), which gives the observed jet shape on sky as a function of the jet half-opening angle and of the projected angular period of the wiggles , we fit the peak positions measured in the continuum-subtracted image (Fig. 1c) and obtain the best-fit parameters and mas. These values indicate an absolute orbital radius of the jet source AU and, for a jet tangential velocity of 192 km s, an orbital velocity km sand a orbital period yr. The fitted wiggling is also roughly compatible with the position of the more external and diffuse knot D, whereas the location of the farther K1 knot deviates by more than 50 mas from the expected position. However, at larger distances from the sources we expect that the interaction of the jet with the ambient medium may produce significant deviations of the jet stream. In the adopted formulation the total mass of the binary system is given by MAUyr , where is the ratio of the mass of the companion to the total mass. Previous estimates of the FUor mass range between 1 and 3 M (Hartmann et al. 1989, 1986, Van den Ancker 2004). If we consider this as the possible range for the total mass of the binary, we obtain (mass of the jet source)=0.52 M =0.48 M and =2.0 M =1.0 M for the lower and upper limit, respectively. The expected separation would be between 3.8 and 2.6 AU, corresponding to 3.3-2.3 mas, which is well below our angular resolution. This scenario is compatible with the binary hypothesis considered by Millan-Gabet et al. (2006) to explain the low interferometric visibility measured in the band. No signature of such a binary was visible in the optical spectra of Hartmann et al. (1989), which had a resolution of 12 km sand were however dominated by disk emission.
Jet mass ejection rates
Further clues to the launch of the FUor jet and the accretion of its driving source can be gathered by measuring the jet mass loss rate (). In principle, both the [O i] and H flux in the jet could be used for this purpose; however, the H emissivity is subject to a much higher uncertainty owing to the contribution of both recombination and collisional excitation (see Bacciotti & Eislöffel 1999). In addition, H could be optically thick at the jet base, therefore not tracing all the jet mass. Therefore, we chose to consider only the [O i] line for determining the . To derive the [O i] flux we first flux-calibrated the images following the procedure described in Appendix C. The jet [O i] flux was then measured in a rectangular box of 160 mas in length, which covers the part of the jet between knots A-C, and 85 mas in width. Four equivalent boxes on adjacent areas were used to evaluate the background contribution. Based on the standard deviation of the background counts in these areas we estimated a 10% uncertainty on the jet count measurement. The jet flux eventually derived is erg s cm. Estimates of the visual extinction towards the target span between 2.4 and 4.6 mag (Stelzer et al. 2009, and references therein), which suggest a relatively low circumstellar on source. Hence, we considered a purely interstellar extinction towards the jet and assumed a value = mag based on the 1 kpc distance of Z CMa (Savage & Mathis 1979). The final extinction-corrected flux is erg s cm.
From the line luminosity, we can compute by considering the projected jet length,
the measured tangential velocity of 192 km s, and the oxygen emissivities (see e.g. Antoniucci et al. 2008),
which depend on the gas temperature and electron density .
For the computation, we considered an oxygen abundance of 4.6 (Asplund et al. 2005) and assumed that all oxygen is in neutral state.
The emissivities were taken from the grids of the the NEBULIO database
Appendix A: Continuum-subtraction procedure
For the H observations, we performed a frame-by-frame subtraction of the B_Ha and Cnt_Ha images to remove the continuum emission and evidence the line-emitting structures around the components. Prior the subtraction, each Cnt_Ha frame must be normalized to the corresponding B_Ha frame to account for the different filter band-passes. The normalization factor is empirically computed from the ratio of the total counts measured in the same region of the two images. The comparison between the photon count ratios between the two Z CMa components in the continuum and in the H frames clearly reveals a strong H emission from the NW component, in agreement with previous observations of the separate spectra of the two stars (Hinkley et al. 2013). This makes it difficult to find a normalization that efficiently removes the intense continuum from both stars at the same time. Indeed, if we normalize to the signal of the SE component, the point-like H emission from the Herbig star creates a very strong residual PSF signal that hampers the detection of the structures around both objects. Conversely, by normalizing to the Herbig component, we oversubtract the signal from the FUor and eventually create a deep residual negative region around the binary where no information can be retrieved. The best compromise is found with a normalization factor calculated over a region including both components. After de-rotating the observations taken at a different field position angle, we finally performed a median of all the subtracted frames.
For the [O i] data, we produced a median [O i] image from all the frames, after de-rotating those obtained at a different field position angle. Analogously, we created a median image of the Cnt_Ha frames, which we normalized to the signal of the Herbig star and then subtracted from the median [O i] image. The residuals from continuum subtraction are more compact for [O i], probably because the point-like [O i] emission from the NW source is much weaker than H.
Appendix B: Binary separation and orbital motion
We derived the separation and P.A. of the Z CMa components by computing the average centroids of the two stars from the 60 frames acquired with the Cnt_Ha filter. Considering the ZIMPOL pixel scale of 3.6006 0.0051 mas/pix and the vertical axis position angle of 358.39 0.11 deg (Ginski et al., in prep.), we obtain a separation of 114.23.1 mas and a position angle of 136.91.5 deg, which corresponds to a projected length of about 131 AU. This measurement and others collected in a time span of 26 years (Koresko et al. 1991; Haas et al. 1993; Thiebaut et al. 1995; Barth et al. 1994; Millan-Gabet et al. 2002; Bonnefoy et al. 2016, submitted) are reported in Fig.2. The plot provides evidence of the variation of the secondary component position relative to the Herbig star. Our measurement indicates that the separation is still increasing with the position angle, suggesting a (projected) elliptical shape. If we consider masses of 16 and 3M for the Herbig and FUor (van den Ancker et al. 2004) and assume an inclined circular orbit, the still increasing separation indicates an orbital period longer than 345 yr. A more detailed analysis of the orbit is presented in (Bonnefoy et al. 2016).
Appendix C: Flux calibration of the [O i] image
To perform the flux calibration of the [O i] image, we first calibrated the Cnt_Ha image using the conversion coefficient reported by Schmid et al. (2016, in prep.), assuming a 10% uncertainty on the provided factor. To check the goodness of this calibration, we derived an estimate of the binary integrated -band magnitude by measuring the counts in a circular aperture including both components. We thus obtain mag for the binary. The measured magnitude is consistent with the integrated magnitudes of and reported by the Kamogata Wide-field Survey on 30 Mar 2015 (Maehara & Ukita 2015), i.e. the night before our SPHERE observations, also considering that the typical photometric variability of Z CMa on a timescale of one day is below 0.1 mag.
To flux calibrate the [O i] image, we then made the assumption that the continuum flux density in the Cnt_Ha image (at nm) is the same for the [O i] image (at nm). On this basis, we were able to derive a count-to-flux conversion factor for the [O i] image by considering a circular area centred on the SE component on both images. We selected the FUor because we expect a weaker on-source [O i] line emission on this component than on the Herbig star. Using different radii for the circular aperture we derive conversion factors that differ for less than 2%. To be conservative, also considering the previous continuum flux density assumption, we associate a relative uncertainty of 10% on the computed factor. This, in conjunction with the initial uncertainty on the Cnt_Ha image conversion coefficient, allows us to evaluate an accuracy of 20% for the flux calibration of the [O i] image.
- Simultaneous observations on the continuum are not possible with this filter, which is mounted on the common wheel of the two cameras.
- Anglada, G., López, R., Estalella, R., et al. 2007, AJ, 133, 2799
- Antoniucci, S., La Camera, A., Nisini, B., et al. 2014, A&A, 566, A129
- Antoniucci, S., Nisini, B., Giannini, T., & Lorenzetti, D. 2008, A&A, 479, 503
- Asplund, M., Grevesse, N., & Sauval, A. J. 2005, in ASP Conf. Ser. 336: Cosmic Abundances as Records of Stellar Evolution and Nucleosynthesis, ed. T. G. Barnes, III & F. N. Bash, 25–+
- Audard, M., Ábrahám, P., Dunham, M. M., et al. 2014, Protostars and Planets VI, 387
- Bacciotti, F. & Eislöffel, J. 1999, A&A, 342, 717
- Benisty, M., Malbet, F., Dougados, C., et al. 2010, A&A, 517, L3
- Beuzit, J.-L., Feldt, M., Dohlen, K., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 701418
- Bonnefoy, M., Chauvin, G., Dougados, C., et al. 2016, A&A, in press
- Cabrit, S., Codella, C., Gueth, F., et al. 2007, A&A, 468, L29
- Calvet, N., Hartmann, L., & Kenyon, S. J. 1993, ApJ, 402, 623
- Canovas, H., Min, M., Jeffers, S. V., Rodenhuis, M., & Keller, C. U. 2012, A&A, 543, A70
- Croswell, K., Hartmann, L., & Avrett, E. H. 1987, ApJ, 312, 227
- Ellerbroek, L. E., Podio, L., Kaper, L., et al. 2013, A&A, 551, A5
- Garcia, P. J. V., Thiébaut, E., & Bacon, R. 1999, A&A, 346, 892
- Giannini, T., Antoniucci, S., Nisini, B., Bacciotti, F., & Podio, L. 2015, ApJ, 814, 52
- Hartmann, L. & Kenyon, S. J. 1996, ARA&A, 34, 207
- Hartmann, L., Kenyon, S. J., Hewett, R., et al. 1989, ApJ, 338, 1001
- Herbst, W., Racine, R., & Warner, J. W. 1978, ApJ, 223, 471
- Hinkley, S., Hillenbrand, L., Oppenheimer, B. R., et al. 2013, ApJ, 763, L9
- La Camera, A., Antoniucci, S., Bertero, M., et al. 2014, PASP, 126, 180
- Machida, M. N. 2014, ApJ, 796, L17
- Maehara, H. & Ukita, N. 2015, The Astronomer’s Telegram, 6874
- Maurri, L., Bacciotti, F., Podio, L., et al. 2014, A&A, 565, A110
- Millan-Gabet, R., Monnier, J. D., Akeson, R. L., et al. 2006, ApJ, 641, 547
- Podio, L., Eislöffel, J., Melnikov, S., Hodapp, K. W., & Bacciotti, F. 2011, A&A, 527, A13
- Poetzel, R., Mundt, R., & Ray, T. P. 1989, A&A, 224, L13
- Prato, M., La Camera, A., Bonettini, S., & Bertero, M. 2013, Inverse Problems, 29, 5017
- Prato, M., La Camera, A., Bonettini, S., et al. 2015, New Astronomy, 40, 1
- Savage, B. D. & Mathis, J. S. 1979, ARA&A, 17, 73
- Schmid, H.-M., Downing, M., Roelfsema, R., et al. 2012, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 8446, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, 84468Y
- Stelzer, B., Hubrig, S., Orlando, S., et al. 2009, A&A, 499, 529
- Szeifert, T., Hubrig, S., Schöller, M., et al. 2010, A&A, 509, L7
- Thalmann, C., Schmid, H. M., Boccaletti, A., et al. 2008, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, Vol. 7014, Ground-based and Airborne Instrumentation for Astronomy II, 70143F
- van den Ancker, M. E., Blondel, P. F. C., Tjin A Djie, H. R. E., et al. 2004, MNRAS, 349, 1516
- Whelan, E. T., Dougados, C., Perrin, M. D., et al. 2010, ApJ, 720, L119