Issue |
A&A
Volume 509, January 2010
|
|
---|---|---|
Article Number | L7 | |
Number of page(s) | 5 | |
Section | Letters | |
DOI | https://doi.org/10.1051/0004-6361/200913704 | |
Published online | 19 January 2010 |
LETTER TO THE EDITOR
The nature of the recent extreme outburst of the Herbig Be/FU Orionis binary
Z Canis Majoris
,
T. Szeifert1 - S. Hubrig2 - M. Schöller3 - O. Schütz1 - B. Stelzer4 - Z. Mikulásek5
1 - European Southern Observatory, Alonso de Cordova 3107, Casilla
19001, Santiago 19, Chile
2 -
Astronomisches Institut Potsdam, An der Sternwarte 16, 14482 Potsdam,
Germany
3 -
European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching,
Germany
4 -
INAF-Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134
Palermo, Italy
5 -
Department of Theoretical Physics and Astrophysics, Masaryk University,
Brno, Czech Republic
Received 19 November 2009 / Accepted 17 December 2009
Abstract
Context. Z CMa is a binary system which consists of two
young stars: a Herbig AeBe component Z CMa NW embedded in a
dust cocoon and a less massive component Z CMa SE, which is
classified as a FU Orionis type star. Associated to the binary system
is a giant parsec-size jet. Past spectropolarimetric observations
showed that the position angle of the linear optical polarization is
perpendicular to the jet axis, indicating that the visual light escapes
the cocoon via cavities aligned with the jet axis and is then scattered
back into the line of sight of the observer. Recently the system showed
the largest outburst reported during the almost 90 years of
available observations.
Aims. We present new spectrophotometric and spectropolarimetric data obtained in 2008 during the recent outburst phase.
Methods. The data obtained in the visual spectral range at
medium spectral resolution were used to study the geometry of the
system from the linear polarization spectra as well as its magnetic
field from the circular polarization spectra.
Results. During the recent outburst we detected that the
Z CMa system is polarized by 2.6% in the continuum and
emission line spectrum, with a position angle still perpendicular to
the jet. From the high level of polarization we concluded that the
outburst is associated with the dust-embedded Herbig AeBe
NW component. The deep absorption components of the Balmer lines
in the velocity frame which extend from zero velocity and reach a wind
velocity of 700 km s-1,
together with the absence of a red-shifted broad emission at similar
velocities, indicate a bi-polar wind. We did not detect a significant
mean longitudinal magnetic field during the outburst, but in the data
obtained in 2004 we detected the possible presence of a rather
strong magnetic field of the order of
1 kG.
However, we critically review the applied method of magnetic field
measurements in the presence of a strong stellar wind. The main result
of our studies is that the bolometric luminosity of Z CMa remained
surprisingly constant during the recent outburst. We conclude that
either the geometry of the cavity through which the light escapes from
the cocoon has opened a new path or that the screen of dust, which
reflects the light toward the observer became more efficient, causing
the observed increase of the visual brightness by about 2
5.
Key words: stars: individual: Z Canis Majoris - stars: pre-main sequence - stars: winds, outflows - stars: magnetic field - stars: variables: general - binaries: close
1 Introduction
Z CMa (HD 51370) is a young visual binary system consisting of a
12
Herbig AeBe star and a 3
FU Ori companion of an age
of 0.3 Myr (Alonso-Albi et al. 2009). Both components
are surrounded by active accreting disks. Associated with the binary system is
a giant parsec-size jet (Poetzel et al. 1989). First evidence
for the binarity was an elongation of the K-band image of this star toward the
northwest (Leinert & Haas 1987;
Koresko et al. 1989). The binary was resolved by
Koresko et al. (1991) with near-IR speckle interferometry.
The angular separation of the binary was
only
0
1 with a position angle of 120
.
At that time the
south-eastern (Z CMa SE) component was dominant in the visual spectral range
and near-infrared J and H bands, while the north-western component was found brighter at
2.2
m and longer wavelengths. Thiébaut et al. (1995)
resolved the binary system with speckle interferometry at 730 nm and 656 nm.
From the data obtained in the quiescent state in 1989.84 the authors deduced
L = 1300
,
R = 13
and
K for the Z CMa SE FU Ori component
and L = 2400
,
R = 1690
,
and
K
for the near-sinfrared Z CMa NW component.
The authors attributed the excess flux associated with Z CMa NW at
optical wavelengths to scattered light, arguing that photons originating from
the hidden central source are scattered back into the line of sight.
Scattering of the photons is assumed to take place along the walls of a
cavity (usually ascribed to a strong bipolar outflow) evacuated within the
envelope. Radio observations in 2005 and 2006 in the cm and mm range suggest
that we see Z CMa NW at an inclination of 30
(Alonso-Albi et al.
2009). This means that we do not observe the inner walls of the
bipolar cavity within the dusty envelope, where the photons are scattered.
The scattered light scenario is well-consistent with optical
spectropolarimetric observations of Whitney et al. (1993),
who detected in 1991-1992 a significant linear polarization of 1% to 2% in
the continuum and up to 6% in the strong emission lines. The different
polarization of continuum and lines have led Whitney et al.
(1993) to suggest that the emission lines have their origin in
the Z CMa NW component. At that time the NW component contributed about 20%
of the flux of the whole system at visible wavelength. The Z CMa SE
component is assumed to be unpolarized, the NW component is assumed to
show 6% linear polarization and the composite spectrum of the binary would
then depend on the flux ratio at the respective wavelength. The position
angle of the polarization of 150
in November 1991 and 154
in March 1992 was roughly perpendicular to the 3.6 parsec scale jet of
Z CMa at a position angle of 60
,
discussed by
Poetzel et al. (1989).
The same authors reported several velocity components in the jet, of which
the fastest knots reaches a velocity of about -620 km s-1.
This position angle of the polarization indicates that the
light from the Z CMa NW component primarily escapes from the envelope
by scattering off the walls of a jet-blown cavity. A smaller scale jet with a
similar position angle was recently discovered in X-rays by Stelzer et al.
(2009).
![]() |
Figure 1: Light curve of Z CMa over the last decades. Open symbols denote the observations from the American Associations of Variable Stars Observers (AAVSO). The filled symbols correspond largely to the ASAS data and during the peak of the light curve to our measurements (blue asterisks). Epochs of NACO observations, the high resolution images as summarized by Thiébaut et al. (1995) and the epoch of spectropolarimetric observations (arrows) are indicated on the top. |
Open with DEXTER |
According to photometric and visual observations available from various
archives, Z CMa recently exhibited the largest outburst reported
during the almost 90 years of available observations.
The light curves presented in Fig. 1 were extracted from
the archives of the All Sky Automated Survey (ASAS) rejecting the data
points with quality flags C and D (Pojmanski 2002).
More data over a larger time span were retrieved from the AAVSO data base.
For the later visual estimates the data points taken for a given month
were averaged after clipping data points beyond a 2
threshold.
Irregular variability was observed in the 1930s. Later
Z CMa steadily increased its brightness, and reached a phase of constant
brightness until about 1973. Since then the star was fading from almost
9
in 1973 to 10
5 in 2007. Superimposed to the fading light curve we
can identify in Fig. 1 several more recent outbursts, the
most prominent ones were observed in 1987 and 2000 with a duration similar to
that of eruptions of EX Lupi-type (EXors) young stars. Between
2002 and 2006 the light curve appeared very irregularly.
A gradual increase of the visual brightness started at the end of 2007,
and reached the highest brightness between November 2008 and May 2009. This
bright phase has now ended, with the latest observations reporting a visual
brightness fainter than 10
after a very rapid decline.
In this letter we present new optical spectropolarimetric observations of linear and circular polarization in the Z CMa system obtained during the recent outburst and discuss the cause for this outburst.
![]() |
Figure 2: Spectrophotometric flux of Z CMa in 2004 (lower spectrum) and 2008 (upper spectrum). |
Open with DEXTER |
![]() |
Figure 3:
Spectropolarimetry of Z CMa in H |
Open with DEXTER |
2 Results and discussion
2.1 Spectroscopy and linear polarization
The spectropolarimetric data sets provide very important insights into the
nature of the at visual wavelengths unresolved binary star. We found
significant changes of the line shapes in the linearly polarized light with
respect to the data published by Whitney et al. (1993). In our
data obtained during the extreme outburst, we see a flat polarization spectrum
with an average polarization P=2.6% % at 4350 Å at a position
angle of 160
.
This polarization angle is roughly perpendicular to the parsec scale jet
associated with Z CMa and similar to observations presented
by Whitney et al. (1993), who measured in November 1991
P = 2.0%, PA = 150
and four months later
P = 1.4%, PA = 154
.
However, in the deep absorption troughs
of the Balmer lines and other strong wind absorption lines, we detect
depolarization signatures (see Fig. 3). These line
components caused by the accelerated wind show a composite spectrum:
we observe a broad blue trough typically associated with an accelerated
wind, visible to about -600 km s-1, with a less deep absorption up to
about -800 km s-1. As is shown in more detail in Fig. 4,
this is not mirrored by emission on the red side of the lines as would be
expected for spherical symmetric winds. The trough is rather deep, indicating
that a large fraction of the stellar disk is obscured by the wind acceleration
zone.
![]() |
Figure 4:
Logarithmic spectral shape of the Balmer lines H |
Open with DEXTER |

We note that the depth and the width of the absorption at the blue-shifted
side of the Balmer lines indicates a bipolar geometry of the wind in absence
of a broad 600 km s-1 red-shifted emission. If no such broad red
emission is observed, it can be argued that either there is no gas
streaming out except out of the polar region, or that the recombination
radiation in the line is not visible through the cavities, but blocked by the
disk or the cocoon. For the H line though we find a broad
emission component beyond the velocity of the wind absorption feature. The
H
absorption trough is not as deep as for other Balmer
lines. This broad component is caused by electron
scattering as is often seen in luminous hot stars and therefore not
related to high velocity outflows.
The spectral slope as indicated in Fig. 2 appears redder during the recent outburst, which contradicts an increase of temperature as source of the luminosity change as proposed by van den Ancker et al. (2004). Similarly, the very insignificant change of the depth of the Balmer jump would argue against a dramatic change of stellar parameters.
2.2 Circular polarization
The analysis of circular polarization measurements obtained during the
outburst reveals no significant detection. In the data obtained in May 2008
we detect a weak signal of circular polarization in the trough of all Balmer
lines, which does however not match the unpolarized line profile in wavelength.
While reviewing the archive data already published by Wade et al.
(2007) we do detect a clear circular polarization signal in
particular at the blue shifted side of the high-n Balmer lines. For the
low excitation lines from H
to H
we measure
G,
while from the higher excitation Balmer lines a very significant field of
G (see Figs. 5 and 6) was
obtained. Still, these results have to be taken with some caution.
To properly interpret the circular polarization signal, it is necessary to
consider the composite nature of the binary system spectrum,
the superimposed effects of the stellar wind, and the variability of the
target. The B-field diagnostics as applied above correlate the s-shaped
profile in V/I with the s-shaped profile of the spectral gradient
dI/d
of the Voigt profiles.
The B-field diagnostics may for instance be disturbed by the
light of the companions and in particular by the stellar wind features due
to the spectral gradient dI/d
along the x-axis of the diagnostic
diagram (see Figs. 5 and 6).
It should also be considered that the target is variable on very short time
scales, such that spectra may have changed between the respective exposures at
.
It is furthermore unclear if these fields were detected in the dust-embedded
Herbig Ae star or in the SE companion. In the first case
the strong field is only seen in the fainter Balmer lines because
the high opacity of the expanding wind at the blue-shifted part
of the strong lines would not allow light from the photosphere
to escape from the NW component. Alternatively, without any indication
of the contribution of the respective two stars to the UV and
visible spectra, one can equally explain the data if the
magnetic field was detected in the SE component. This component may have
dominated the spectra of 2004 in the wavelength range of the Balmer
lines where we measured the Stokes V features, whereas later in May 2008
the contribution of Z CMa SE became insignificant on the blue-shifted
side of the P Cygni line of Z CMa NW in a way that we detect a hint
of circular polarization in the red-shifted part of the hydrogen lines only,
but no S-shaped feature, which would lead to a magnetic field detection.
Nor would the wind-dominated spectral line gradients dI/d
linearly
correlate with V/I. With the few available spectra it is not yet possible
to reliably conclude on the origin of the features in the circular
polarization spectrum.
![]() |
Figure 5:
Circular polarization in the Balmer line spectra observed in
November 2004 ( on the bottom) together with the spectral gradients
dI /d |
Open with DEXTER |
2.3 Origin of the outbursts
Despite the different origin of the NIR Ks band flux from the
surrounding dust cocoon and the optical flux with an origin from
the pole of the star, we presume that a strongly changed optical/UV
source would heat up the dust, and the increase of the luminosity of
the stellar core would be re-emitted by thermal radiation in the
-band and beyond. There is however no significant change of the
-band flux as far as we can tell from the K-band flux ratios given
in Table 2 or from comparing the most recent report of infrared
observations by Antoniucci et al. (2009) of
at the peak of the light curve in March 2009,
with
reported by Koresko et al. (1991).
Similarly there is no change in the X-ray properties reported by
Stelzer et al. (2009) with respect to earlier measurement
at a ``quiescent'' phase. If the short-term optical brightening was an
EXor-like accretion event, X-ray variability would have been expected.
Our polarimetric data lead to a new interpretation in which
the X-ray emission associated with the binary may be attributed to the
lower-mass SE component, which is not involved in the outburst.
The apparently insignificant change of the K-band thermal emission
rather argues against the strongly increased luminosity of
during short outbursts, as estimated by
van den Ancker et al. (2004) from the visual SEDs with respect
to only L = 2400
estimated by Thiébaut et al.
(1995) from the SED including the thermal emission.
It rather underlines our interpretation that the bolometric luminosity
during the outbursts remains constant, but the amount of visible light
transmitted through the envelope has changed.
The derived high luminosity was primarily a consequence of the
early spectral type estimated for the NW source from the He I and
O I lines, which lead to very high bolometric and extinction
corrections. We detect these lines with blue shifted broad absorption profiles
and with P Cygni profiles respectively, both indicating an origin of the
lines in the accelerated wind rather than in the photosphere of the NW
component. The wind however, with an indication of bipolarity and eventually
with indications for a magnetic field, may well be excited by other
sources instead of by the thermal radiation from the star.
The short-term outbursts are nevertheless superimposed to long-term
FU Ori like changes in the visual brightness, with a slow brightening by
one magnitude until about 1970 and a subsequence fading by a similar amount.
It is well possible that any of the binary components contributes to these
changes.
3 Summary
We have shown that the NW component of Z CMa was undergoing a historically
bright outburst in 2008/2009. Our polarimetric data support a variable
scattering geometry of the NW envelope as the cause for the optical outburst,
rather than FU Ori or EXor events from the SE component.
From the changes in the linear polarization spectra
taken at the outburst phase we conclude that the scattering path in which the
light of Z CMa NW escapes from the dust cocoon must have changed. The
linear polarization from the now dominant NW component was found to be smaller
than the 6% estimated at the fainter phase in 1991.
We do not find changes in the K-band magnitude in contrast to the recent
change of 2
5 in the visual. A significant fraction of the variability
may be explained by the changed, more direct line of sight in which we now
see the Z CMa NW:
more light from the central source would reach the observer through a
widened cavity, while the central source directly observed in the K-band
did not significantly change in luminosity.
We acknowledge with thanks the variable star observations from the AAVSO International Database contributed by observers worldwide and used in this research. This work profited greatly from data extracted from the ASAS data base. Based on observations made with the European Southern Observatory telescopes obtained from the ESO/ST-ECF Science Archive Facility.
References
- Antoniucci, S., Arkharov, A. A., Giannini, T., et al. 2009, ATel, #2024 [Google Scholar]
- Appenzeller, I., Fricke, W., Fürtig, W., et al. 1998, Messenger, 94, 1 [Google Scholar]
- Alonso-Albi, T., Fuente, A., Bachiller, R., et al. 2009, A&A, 497, 117 [Google Scholar]
- Bagnulo, S., Szeifert, T., Wade, G. A., et al. 2002, A&A, 389, 191 [Google Scholar]
- Garcia, P. J. V., Thiébaut, E., & Bacon, R. 1999, A&A, 346, 892 [Google Scholar]
- Koresko, C. D., Beckwith, S. V. W., & Sargent, A. I. 1989, AJ, 98, 1394 [NASA ADS] [CrossRef] [Google Scholar]
- Koresko, C. D., Beckwith, S. V. W., Ghez, A. M., et al. 1991, AJ, 102, 2073 [NASA ADS] [CrossRef] [Google Scholar]
- Leinert, C., & Haas, M. 1987, A&A, 182, L47 [Google Scholar]
- Lenzen, R., Hartung, M., Brandner, W., et al. 2003, SPIE, 4841, 944 [Google Scholar]
- Millan-Gabet, R., & Monnier, J. D. 2002, ApJ, 580, L167 [NASA ADS] [CrossRef] [Google Scholar]
- Poetzel, R., Mundt, R., & Ray, T. P. 1989, A&A, 224, L13 [Google Scholar]
- Pojmanski, G. 2002, Acta Astron., 52, 397 [NASA ADS] [Google Scholar]
- Stelzer, B., Hubrig, S., Orlando, S., et al. 2009, A&A, 499, 529 [Google Scholar]
- Thiébaut, E., Bouvier, J., Blazit, A., et al. 1995, A&A, 303, 795 [Google Scholar]
- van den Ancker, M. E., Blondel, P. F. C., Tjin A Djie, H. R. E., et al. 2004, MNRAS, 349, 1516 [NASA ADS] [CrossRef] [Google Scholar]
- Wade, G. A., Bagnulo, S., Drouin, D., et al. 2007, MNRAS, 376, 1145 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Whitney, B. A., Clayton, G. C., & Schulte-Ladbeck, R. E. 1993, ApJ, 417, 687 [NASA ADS] [CrossRef] [Google Scholar]
Online Material
Table 1: Journal of spectroscopic observations.
Table 2: Journal of NACO observations.
![]() |
Figure 6:
Magnetic field measurements using the data obtained in November 2004.
The longitudinal magnetic field is derived following the equation
|
Open with DEXTER |
Appendix A: Observations and data reduction
The observational material discussed in this section includes FORS 1
spectropolarimetric data (Appenzeller et al. 1998) obtained
in the framework of our programs carried out in May 2008 and in December 2008.
Additionally, we used ESO archive data from 2004 obtained in circular
spectropolarimetric mode.These data were used by Wade et al.
(2007) to search for the presence of longitudinal magnetic fields
in a larger number of Herbig AeBe stars. These older observations (see
Table 1) were carried out with a slit width of 0
8
and grism 600B (
). The flux spectra in the ordinary and
extraordinary ray of the Wollaston prism were extracted for all position
angle settings of the retarder waveplate. Since there were no available flux
calibrators observed for this setup in the respective night, we used flux
calibrators taken in normal spectroscopic mode without the Wollaston prism.
The new spectropolarimetric data from 2008 were obtained with the same
instrument, but now equipped with a higher UV response detector system with a
smaller pixel size.
During the December run we have used the higher dispersion
grism 1200B in addition to the data with the 600B grism. In all our
observations the slit width was set to 0
4 to optimize the resolution
for the polarimetric measurements (
and
,
respectively).
The spectrophotometric calibrators retrieved from the archive were observed
without polarization optics on different nights. One source of uncertainty
in the flux calibration is therefore the unclear spectral response of the
wave retarder plates, for which we presume a high throughput over the full
spectral range. However, significant slit losses are expected with this
setup due to the narrow slit widths and the difficult acquisition of the
target on the slit considering the extreme brightness of the target for
8 m-class telescopes. To remedy this difficulty related to the absolute
calibration, we re-scaled the flux at 5500 Å to the monthly average of the
AAVSO visual observations of 9
4 in November 2004 and 8.4 in
December 2008.
In Fig. 2 we present spectrophotometric
fluxes for the Z CMa system in 2004 and 2008, respectively.
Further we corrected the result with the chromatic zero angle of the
super-achromatic half wave retarder plate
.
for the proper interpretation of the measured polarization angle.
To perform linear polarization measurements, a Wollaston prism and a half-wave
retarder plate rotated in 22.5
steps between 0 and 67.5
were
used. For circular polarization measurements the quarter-wave retarder plate
was used at the positions +45
and -45
.
Our FORS 1 data were supplemented by three high-resolution
FEROS spectra (R=48 000) obtained in February and March 2009 to study the
behavior of hydrogen lines. The journal of spectroscopic observations is
presented in Table 1.
To better understand the parameters of the resolved binary system we retrieved all available NACO adaptive optics near-infrared data (Lenzen et al. 2003) from the ESO archive. As for the FORS 1 data sets, the target is rather bright for imaging-mode observations with 8 m-class telescopes, and accordingly the various data sets were taken with a mixed choice of neutral density and narrow band filters. As a consequence we cannot provide absolute photometry for NACO data with the required level of confidence, but only the flux ratio of the two components (see Table 2) and the astrometric parameters of the binary. In case of normal jittered observation sequences we have used the Eclipse data reduction package provided by ESO to NACO users. In other cases we had only two acquisition exposures or a few jitter positions.
Footnotes
- ... Majoris
- Based on observations made with ESO Telescopes at the Paranal Observatories under programme IDs: 074.C-0442, 081.C-0410, and 282.C-5041.
- ...
- Appendix A, Tables 1 and 2, and Fig. 6 are only available at http://www.aanda.org
- ... plate
- The data is available on the instrument web pages: http://www.eso.org/sci/facilities/paranal/instruments/fors/
All Tables
Table 1: Journal of spectroscopic observations.
Table 2: Journal of NACO observations.
All Figures
![]() |
Figure 1: Light curve of Z CMa over the last decades. Open symbols denote the observations from the American Associations of Variable Stars Observers (AAVSO). The filled symbols correspond largely to the ASAS data and during the peak of the light curve to our measurements (blue asterisks). Epochs of NACO observations, the high resolution images as summarized by Thiébaut et al. (1995) and the epoch of spectropolarimetric observations (arrows) are indicated on the top. |
Open with DEXTER | |
In the text |
![]() |
Figure 2: Spectrophotometric flux of Z CMa in 2004 (lower spectrum) and 2008 (upper spectrum). |
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Spectropolarimetry of Z CMa in H |
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Logarithmic spectral shape of the Balmer lines H |
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Circular polarization in the Balmer line spectra observed in
November 2004 ( on the bottom) together with the spectral gradients
dI /d |
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Magnetic field measurements using the data obtained in November 2004.
The longitudinal magnetic field is derived following the equation
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.