next previous
Up: Spectroscopic and photometric

3. Results

3.1. Optical photometry

UBV photometry (Table 3 (click here)), carried out in January 1996, shows that Z CMa is presently in a low state at tex2html_wrap_inline1995. This value is in good agreement with the decreasing luminosity trend shown in the AAVSO database regarding Z CMa (Mattei 1996). We derived the following colors: B-V = + 1.00 and U-B = + 0.54.

 

 
Filter      J.D.      Mag Exp.
(s)
U 2450100.4780 11.78tex2html_wrap_inline2005.04 600
B 2450100.4910 11.24tex2html_wrap_inline2005.04 300
V 2450100.4976 10.24tex2html_wrap_inline2005.01 120

Table 3: UBV photometry obtained with BFOSC+CCD on 1996 January 17

3.2. Low-dispersion spectra

In the period March-April 1985 the star was in an active state, confirmed by the measured value V = 9.0, in comparison with V = 9.4 measured in April 1984 (Hessman et al. 1991). The variability trends of Htex2html_wrap_inline2019, tex2html_wrap_inline1859, tex2html_wrap_inline1857 and Fe II lines are shown in Fig. 1 (click here).

  figure353
Figure 1: Variation of the Z CMa spectrum from the 1984 quiescent state to the 1985 active state (IDS spectra). Transition of tex2html_wrap_inline1859 from pure absorption to P Cygni morphology, enhancement of tex2html_wrap_inline1857, and the appearance of FeII emission are seen

The main variations occurred from April 1984 to March-April 1985 are the following:

  1. An increase of about 100tex2html_wrap_inline2029 of the tex2html_wrap_inline1857 emission line and an increase of 100 kms-1 of the velocity of its blue-shifted absorption component (Fig. 1 (click here), Tables 4 (click here) and 5 (click here)),
  2. An increase of the FWHM width of the tex2html_wrap_inline1857 emission from 6 to 8.6 Å, corresponding to an expansion velocity increase of about 120 kms-1 (Fig. 1 (click here), Table 5 (click here)).
  3. The transition of the tex2html_wrap_inline1859 line from a pure absorption profile to a strong P Cygni profile and an increase of 160 kms-1 of its blue-shifted absorption component (Fig. 1 (click here), Table 5 (click here)).
  4. A strong increase of FeII emission number and the transient occurrence of strong P Cygni effects in some of the FeII lines (Fig. 1 (click here), Table 6 (click here)).

No significant variations are present in the NaID doublet.

The spectrum obtained on January 17 1996 shows a very prominent variation of the state of the star, characterized by a drastic change of the tex2html_wrap_inline1857 profile and by a 20tex2html_wrap_inline2029 increase of tex2html_wrap_inline1857 emission EW in comparison with the 1985 spectrum (Fig. 1 (click here), Table 4 (click here)). The presence of some slight P Cygni effects both in tex2html_wrap_inline1859 and in Htex2html_wrap_inline2019 seems to indicate a residual activity. A much more detailed analysis of the unusual tex2html_wrap_inline1857 profile and of other peculiar characteristics was performed on the medium dispersion BFOSC spectrum, the results of which are discussed in Sect. 3.3.

No clear evidence of variability of EW and profile on a monthly or daily time scale is present in any spectrum.

A study of very short-term variability (seconds to minutes), carried out with a time-sequence of 15 spectra on April 16, 1984 and with a time-sequence of 7 spectra on April 17, 1984, gave a negative result: the most important lines (tex2html_wrap_inline1857 and tex2html_wrap_inline1859) did not show any EW or morphological variation. This result can be expected because of the quiescent state of the star in this period, in which friction-driven fast variations in the gas of the accretion disk are unlikely to be observed using low-dispersion spectroscopy. A study of the short-term variability could give more significant results in the active states of the star. Unfortunately we could not carry out this monitoring in 1985 due to the limited available observing time.

 

 

Date

Absorption Emission

tex2html_wrap_inline2063tex2html_wrap_inline1859NaItex2html_wrap_inline1857 Htex2html_wrap_inline2019 tex2html_wrap_inline1859 tex2html_wrap_inline1857
Apr. 84 3.7 5.8 6.9 2.8 0.7 absent 14.6
Mar. 85 4.3 3.9 5.2 3.2 0.3 2.4 24.9
Apr. 85 4.4 3.5 5.7 2.9 0.2 2.6 28.3
Jan. 96 1.1 2.9 5.7 absent 2.5 0.7 32.5

Table 4: Equivalent widths in Ångstrom of the Balmer lines and the doublet NaI D from IDS spectra. Typical errors are tex2html_wrap_inline2061 20%

 

 
Date *RV tex2html_wrap_inline1857 (A) FWHM tex2html_wrap_inline1857 (E) *RV tex2html_wrap_inline1859 (A)
(km s-1) (Å, km s-1) (km s-1)
1984 Apr. 13 - 486.37 6.04, 276.21 - 334.97
1985 Mar. 15 - 582.03 8.60, 393.17 - 493.23

  A: Absorption, E: Emission

  * RV heliocentric.
Table 5: Outstanding dynamical parameters of tex2html_wrap_inline1857 and tex2html_wrap_inline1859 by IDS spectra

 

 

Present at Days of Months:
Identification Apr. 84 Mar. 85 Apr. 85 Jan. 96
FeII 4923.92 - 15, 16, 17 8, 9 17
FeII 5018.43 - 15, 16, 17 8, 9 17
FeII 5158.00 14 - - -
FeII 5169-71 - 15 9, 10 17
FeII 5199.20 - - 9, 10 -
FeII 5234.62 - 15, 17 9, 10 -
FeII 5269.54 - 15, 17 9, 10 -
FeII 5276.00 - 15, 17 9, 10 -
FeII 5284.09 13, 14, 15 15, 17 9, 10 -
FeII 5316.61 - 15, 17 9, 10 17
TiII 5336.80 - - 9, 10 -
FeII 5362.86 13, 14, 15 15, 17 9, 10 -
FeII 5534.90 - 15, 17 9, 10 -
FeII 6149.71 - 15, 17 9, 10 -
FeII 6238.38 - 15, 17 9, 10 -
FeII 6247.56 - 15, 17 9, 10 -
[OI] 6300.20 - 15, 17 (W) 9 (W) 17 (S)
FeII 6369.50 - - - 17
FeII 6432.65 - 15, 17 9, 10 17
FeII 6456.38 - 15, 17 9, 10 -
FeII 6516.05 - 15, 17 9, 10
FeII 6627.30 13 - - -
[SII] 6717.00 - - - 17
FeII 6729.90 - - - 17
FeII 7155.10 - - - 17
FeII 7376.50 - - - 17

W: Weak, S: Strong.

Table 6: Emission lines other than Balmer present in low-dispersion spectra

3.3. Medium and high-dispersion spectra

Our CASPEC spectrum on April 11, 1985 supplied the most detailed line analysis. As it is possible to see in Figs. 2 (click here), 3 (click here) and Table 7 (click here), the 15 reduced echelle orders present a great number of emission lines. The strongest P Cygni effect is present in tex2html_wrap_inline1857 (Fig. 4 (click here)), where we measured a heliocentric velocity tex2html_wrap_inline2119 from the absorption component. The velocity difference between the tex2html_wrap_inline1857 emission and absorption components is on the order of tex2html_wrap_inline2123, while the blue wing of the tex2html_wrap_inline1857 absorption component is characterized by a terminal velocity of about tex2html_wrap_inline2127. The doublet NaID presents many absorption components (Fig. 5 (click here)): the deepest one has tex2html_wrap_inline2129, the bluest one has tex2html_wrap_inline2131, extending up to a terminal velocity of tex2html_wrap_inline2123 (see also Table 7 (click here)). Similar values can be deduced from the CES spectrum. Other P Cygni features can be seen in all FeII lines but their absorption components are too weak to allow accurate equivalent width (EW) and radial velocity (RV) measurements. Much stronger P Cygni effects, as easily shown in the 1985 low-dispersion spectra, should have appeared in the range 4900-5200ÅÅ in FeII 4923Å, 5018Å and 5169-71ÅÅ lines, but unfortunately this range was not covered by CASPEC spectra. The average velocity of the FeII absorption components is - 90 kms-1 while the average velocity of the FeII emissions is about 38 kms-1, comparable with the average velocity of the tex2html_wrap_inline1857 emission component (tex2html_wrap_inline2147).

  figure495
Figure 2: CASPEC overall spectrum (5700-6250 Å) obtained during the 1985 active state. Identified lines are indicated

  figure500
Figure 3: CASPEC overall spectrum (6250-6750 Å) obtained during the 1985 active state. Identified lines are indicated

  figure505
Figure 4: tex2html_wrap_inline1857 line in 1985 CASPEC spectra (short dash), 1989 REOSC spectra (long dash) and 1996 BFOSC spectra (solid). Radial velocity (tex2html_wrap_inline2155) is heliocentric

  figure510
Figure 5: Absorption components of doublet NaID (CASPEC spectrum acquired in 1985 April 11). Heliocentric velocities (tex2html_wrap_inline2155) of absorption components of NaI 5889.9 Å are indicated

In addition to emission and P Cygni features, CASPEC spectra show clear evidence of double-peaked absorption lines at tex2html_wrap_inline2159 5915Å, 6142Å, 6192Å, 6440Å, 6451Å, 6496Å, 6644Å, 6678Å, 6664Å, 6709Å and 6719Å: we measured an average velocity of their peak separation tex2html_wrap_inline2161. The absorption line HeI 6678.15 Å, blended with the FeI disk doubled absorption line at 6678 Å (Fig. 6 (click here)), apparently is the only photospheric feature of interest in our CASPEC spectra. A blue-shift corresponding to tex2html_wrap_inline2163 was derived for the HeI 6678.15 Å absorption line.

Our CES spectra (Fig. 7 (click here), Table 8 (click here)) were mostly concentrated on tex2html_wrap_inline1857. No relevant monthly EW variation of the tex2html_wrap_inline1857 emission is detected between the spectra taken on April 8-9, 1985 and the spectrum taken on March 15, 1985. The EW's may be affected by an erroneous determination of the continuum level due to the small spectral range covered by CES spectra. At the same time, some differences are evident in absorption components between the March and April spectra.

  figure527
Figure 6: CASPEC spectrum (1985 April 11) showing HeI 6678.15 deep absorption line. Double-peaked absorption lines at 6644, 6679, 6709 and 6719 Å\ are also indicated

  figure532
Figure 7: CES spectra, obtained during the 1985 active state, showing variation of the tex2html_wrap_inline1857 absorption component within a period of approximately a month

From CES spectra we determine an average velocity of - 470 kms-1 for the tex2html_wrap_inline1857 absorption component and an average velocity of +44 kms-1 for the emission component. There is a velocity decrease of about 20 kms-1 from 15 March to 8-9 April, which we interpretate as a consequence of wind deceleration. Furthermore, CES spectra show also very strong CaII 8498Å and 8542Å emission, whose radial velocities seem to be in good agreement with the tex2html_wrap_inline1857 emission measurement. CaII lines are present with a (more or less) pronounced P Cygni profile. The derived velocities from absorption components are - 116 and - 196 kms-1 respectively. This difference can be interpreted as the effect of two equally dense and hot shells ejected at different velocities.

The small difference between the velocities measured with CES and the ones measured with CASPEC is probably due to the fact we used two different instruments. Because of this, we focus attention on measurements made with a single instrument and consequently to ratios obtainable between values of a given line parameter present in a specific spectrum.

From the medium-dispersion spectra taken with the RES spectrograph, we obtained tex2html_wrap_inline1857 (Fig. 4 (click here)) and tex2html_wrap_inline1859 profiles. The emissiony (EW = 15.48 Å) and absorption (EW = 6.42 Å) components of tex2html_wrap_inline1857 show a velocity difference of 656 kms-1, while the blue wing of the absorption component extends out to - 1150 kms-1. tex2html_wrap_inline1859 is present strongly in absorption (EW = 7.00 Å) with a weak emission component (EW = 0.10 Å), absorption and emission components show a velocity difference of 270 kms-1.

BFOSC spectra (Table 9 (click here) and Figs. 8 (click here) and 9 (click here)) obtained with a medium-dispersion of 22 Å/mm show the most interesting characteristics both in terms of variability and in terms of peculiar morphology. The following fundamental features are evident:

  figure550
Figure 8: BFOSC overall spectrum (5000-7500 Å) obtained in 1996 January 17. Identified lines are indicated

  figure555
Figure 9: BFOSC overall spectrum (7500-9800 Å) obtained in 1996 January 17. Identified lines are indicated

1. The tex2html_wrap_inline1857 profile is completely changed in comparison to the profile shown in the 1984-1985 spectra (Fig. 4 (click here)). In addition to the main emission a second blue shifted emission is present. The central absorption, shifted by - 516 kms-1 with respect to the main emission, corresponds perfectly to the absorption detected in the ESO spectra. The central wavelengths of the two tex2html_wrap_inline1857 emission components are separated by a velocity of over 800 kms-1.

2. The spectrum shows the usual FeII emissions which are normally present in the spectrum of Z CMa during the active states, but only the FeII 5169Å line shows a sharp P Cygni effect. All the other FeII lines do not show P Cygni features. The average velocity of the FeII emissions is + 77 kms-1.

3. The spectrum shows the forbidden [OI] 6300Å emission line (EW = 1.5 Å), much stronger than in the 1985 CASPEC spectrum, and the forbidden [SII] 6717Å (EW = 0.1 Å), [SII] 6731Å (EW = 0.5 Å) and [FeII] 7155Å (EW = 0.4 Å) emission lines. In particular, [OI] 6300Å and [SII] 6731Å\ present both an asymmetric profile with the red component much stronger than the blue, and with barycenters strongly blue-shifted.

4. The blue wing of the NaI 5889.9 Å absorption line corresponds to a velocity of over - 470 kms-1, about 80 kms-1 less than the values derived from April 1985 CASPEC and CES spectra.

5. In five cases, at 6192 Å, 6345 Å, 6496 Å, 6644 Å and 6664 Å, we identified double-peaked "disk absorption lines''. An average velocity of tex2html_wrap_inline2229 kms-1 was obtained from their peak separation. The shallow HeI 6678Å absorption line, present in CASPEC 1985 spectrum, is absent. At the same wavelength the BFOSC spectrum shows a strong double peaked "disk absorption line''.

3.4. IR photometry

As shown in Table 10 (click here), JHKLM photometry shows decreases of up to 0.2-0.3 mag from April 1984 to March-April 1985. In the same period, color indexes J-K, H-K and K-L are subject to a blueing effect consistent with a decrease of up to 0.2, 0.08 and 0.08 mag respectively. We notice that such IR luminosity increases and color blueing can be correlated with analogous variations detected in the optical range (Hessman et al. 1991). Similar variations in the IR range can be found in the literature (Kenyon & Hartmann 1991; Berrilli et al. 1992; Hamann & Persson 1992; Molinari et al. 1993; Noguchi et al. 1993). No luminosity or blue increases larger than 0.02-0.06 mag are recorded on a daily or monthly timescale.

3.5. Energy distribution

The energy distributions in the range 0.4-5 tex2html_wrap_inline1951m obtained in 1984 and 1985, are presented in Fig. 10 (click here). The optical components of these distributions are represented by low-dispersion spectra taken on March 17, 1985 and on April 17 1984. The IR components are given by an average of spectrophotometric data taken in March and April 1985, overlapped with IR photometric data taken in April 1984 and April 1985.

  figure589
Figure 10: Energy distribution of Z CMa in the range tex2html_wrap_inline1867 constructed from April 1984 (thick line) and March 1985 (dotted line) spectroscopic data merged with IR March-April 1985 spectrophotometric data (dashed line), IR April 1984 photometric data (filled circles) and IR April 1985 photometric data (open circles)

The following main features can be outlined from the energy distributions:

  1. The distributions peak around 1.6 tex2html_wrap_inline1951m.
  2. A deep absorption band, at about 1.9 tex2html_wrap_inline1951m.
  3. The increase of the 1985 IR flux with respect to that of 1984, by a factor of 1.4 in the J passband, 1.3 in H, 1.2 in K and 1.1 in L and M.

 

 
Identification Absorption Emission
EW RV** EW RV**
(Å) (km s-1) (Å) (km s-1)

NaI D

7.0 - 52
'' - 127
'' - 172
'' - 234
'' + 26*
FeII 5991.4 0.2 - 69 0.2 + 40
FeII 6238.4 weak - 66 0.5 + 37
FeII 6247.6 weak - 119 0.6 + 40
[OI] 6300.2 0.2 - 2.4
FeII 6369.5 weak - 74 0.1 + 25
FeII 6416.9 weak - 144 0.2 + 38
FeII 6432.6 weak - 105 0.4 + 45
FeII 6456.4 weak - 72 0.5 + 37
FeII 6516.0 weak - 84 0.7 + 38
Htex2html_wrap_inline2355 6.5 - 510 29.1 + 46
HeI 6678.1 1.0 - 207
FeII 6729.9 0.1 + 20

  * Interstellar line

  ** RV heliocentric.
Table 7: Equivalent widths (EW) and radial velocities (RV) of prominent lines present in CASPEC spectra acquired on April 11 1985. Error in EW is tex2html_wrap_inline2283, error in RV is tex2html_wrap_inline2287. Four values of RV are measured for NaID corresponding to different components as shown in Fig. 5 (click here), while the given EW is the total

 

 

Identification

Absorption Emission   Date
EW  RV** EW RV**
(Å) (km s-1) (Å)  (km s-1)
NaI D 3.5 - 51* Mar. 16
tex2html_wrap_inline1857 4.5 - 483 15.8 + 44 Mar. 16
'' 5.5 - 465 15.7 + 45 Apr. 8
'' 5.2 - 467 15.8 + 46 Apr. 9
CaII 8498.0 0.2 - 116 4.6 + 38 Apr. 10
CaII 8542.1 0.5 - 196 5.2 + 47 Apr. 10

  * Deepest component

  ** RV heliocentric.
Table 8: Equivalent widths (EW) and radial velocities (RV) of prominent lines present in CES spectra. Error in EW is tex2html_wrap_inline2283, error in RV is tex2html_wrap_inline2287

 

 

Identification

Absorption Emission
EW RV* EW RV*
(Å) (km s-1) (Å) (km s-1)
FeII 5169.0 0.3 0.3
FeII 5284.1 0.4
FeII 5455.6? 0.4 0.2
HgI 5460.7? 0.3 + 5.6
NaI D 5.4 - 40
[OI] 6300.2 1.5 - 214
FeII 6432.6 0.2 + 87
FeII 6456.4 0.2 + 97
FeII 6516.0 0.3 + 116
tex2html_wrap_inline1857 1.0 - 378 2.3 - 677
'' 32.5 + 138
[SII] 6717.0 0.1 - 6.0
[SII] 6731.3 0.5 - 102
[FeII] 7155.1 0.4 + 35
FeII 7376.5 0.5 + 48
OII 7938.1?
OI 8446.4 1.5 + 52
CaII 8498.0 5.6 + 68
CaII 8542.1 0.3 - 154 6.4 + 75
CaII 8662.1 0.1 - 173 5.4 + 60

  * RV heliocentric.

Table 9: Prominent lines present in BFOSC (tex2html_wrap_inline2413) spectra acquired on January 17 1996. Error in EW is tex2html_wrap_inline2283, error in RV is tex2html_wrap_inline2419. For lines with tex2html_wrap_inline2421 it was not possible to measure RV because of bad calibration due to echelle order overlapping

 

Date

    J.D.   J   H   K   L   M

1984 Apr. 17

2445807.558 6.11tex2html_wrap_inline2505 4.84tex2html_wrap_inline2507 3.68tex2html_wrap_inline2507 1.85tex2html_wrap_inline2511 0.93tex2html_wrap_inline2513
1984 Apr. 18 2445808.551 6.15tex2html_wrap_inline2515 4.88tex2html_wrap_inline2517 3.71tex2html_wrap_inline2517 1.90tex2html_wrap_inline2521 0.94tex2html_wrap_inline2523
1985 Mar. 17 2446141.507 5.82tex2html_wrap_inline2525 4.61tex2html_wrap_inline2515 3.52tex2html_wrap_inline2529 1.76tex2html_wrap_inline2517 0.80tex2html_wrap_inline2533
'' 2446141.666 5.80tex2html_wrap_inline2525 4.65tex2html_wrap_inline2515 3.52tex2html_wrap_inline2529 1.73tex2html_wrap_inline2517 0.77tex2html_wrap_inline2533
1985 Mar. 18 2446142.579 5.81tex2html_wrap_inline2525 4.63tex2html_wrap_inline2525 3.53tex2html_wrap_inline2529 1.75tex2html_wrap_inline2515 0.83tex2html_wrap_inline2555
1985 Apr. 9 2446164.546 5.84tex2html_wrap_inline2525 4.70tex2html_wrap_inline2515 3.60tex2html_wrap_inline2525 1.83tex2html_wrap_inline2515 0.85tex2html_wrap_inline2523
1985 Apr. 10 2446165.499 5.78tex2html_wrap_inline2517 4.66tex2html_wrap_inline2517 3.56tex2html_wrap_inline2525 1.81tex2html_wrap_inline2525 0.87tex2html_wrap_inline2575
1985 Apr. 11 2446166.490 5.83tex2html_wrap_inline2517 4.68tex2html_wrap_inline2515 3.59tex2html_wrap_inline2515 1.81tex2html_wrap_inline2515 0.83tex2html_wrap_inline2585
'' 2446166.497 5.82tex2html_wrap_inline2517 4.66tex2html_wrap_inline2515 3.56tex2html_wrap_inline2515 1.79tex2html_wrap_inline2515 0.80tex2html_wrap_inline2585
'' 2446166.516 5.82tex2html_wrap_inline2517 4.67tex2html_wrap_inline2515 3.56tex2html_wrap_inline2515 1.79tex2html_wrap_inline2515 0.80tex2html_wrap_inline2585
'' 2446166.523 5.81tex2html_wrap_inline2517 4.66tex2html_wrap_inline2515 3.57tex2html_wrap_inline2515 1.81tex2html_wrap_inline2515 0.79tex2html_wrap_inline2585

Table 10: JHKLM photometric data

 

The absorption band at 1.9 tex2html_wrap_inline1951m, due to the composite transitions of vibration-rotation bands of water vapor (Sato et al. 1992), is the only clear and seemingly constant characteristic of the IR spectrophotometric energy distribution. This absorption feature is the result of the heating effect on water-ice grains which are contained in the external region of the proto-stellar accretion disk. Such heating is caused by the ionization front developed during the 1985 phase.

Given the variability of Z CMa, different energy distributions are expected from data taken in different epochs. In particular, in our data this difference can be noticed in the general flux increase from 1984 to 1985, in a steep increase in the J and H fluxes and in the basic equality of H and K fluxes of 1985 (in comparison with 1984 IR data).

Previous authors (Berrilli et al. 1992; Natta et al. 1993) demonstrate that the envelope is the source of a strong IR excess starting at about 2 tex2html_wrap_inline1951m.


next previous
Up: Spectroscopic and photometric

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr