next previous
Up: Raman scattering in

5. Reference-model results

 

In order to examine the distribution of the Raman line formation, the position, frequency, weight, and polarization of each Raman-scattered photon packet were stored for the reference model, viewed at quadrature. The three spatial dimensions were then collapsed to produce two-dimensional images of the model outputs, which are shown in Figs. 2 (click here), 3 (click here), 4 (click here). The corresponding spectropolarimetric results are shown in Fig. 5 (click here).

The strongest Raman emission in this model originates close to the parent photon source, with the intensity dropping with increasing radial distance simply because of dilution. (This result must be treated with some caution, since the reference model neglects the HII region which will usually be associated with the hot source; within such a region the number of neutral-hydrogen scatterers will be negligible.) The intensity also falls off rapidly away from the red giant, because of the decreasing density of the wind. The bulk of the Raman scattering therefore occurs between the photon source and the red-giant photosphere, where the density of both Otex2html_wrap_inline3400 photons and Htex2html_wrap_inline3402 scatterers is high.

Multiple Rayleigh scatterings degrade the directionality of the OVI photon field, and hence the normalized polarized intensity tex2html_wrap_inline3406. Since there is no UV opacity in the reference model, the Raman photons emerge after three Rayleigh scatterings on average; we therefore expect the normalized polarized intensity to be greatest at moderate to large radial distances from the photon source, because the scatterers will see a more nearly radial radiation field (most scatterings occuring close to the source). Figure 2 (click here)d confirms this expectation, with normalized intensities approaching unity at tex2html_wrap_inline3408 from the parent source. Figures 2 (click here)a and b further confirm that the code generates the expected behaviour in Q and U.

A second effect of Rayleigh scatterings is broaden the velocity dispersion of the OVI photons. Local scatterings have little effect in this regard, as there is only a small velocity dispersion between the source atoms and scatterers, but more distant scatterings can introduce considerable shifts in the profile; scatterings across opposite hemispheres will introduce shifts corresponding to twice the wind velocity. These shifts are always redshifts, as every point in the wind sees every other point receding from it (the `expanding universe' paradigm).

  figure634
Figure 3: Results for the reference model (tex2html_wrap_inline34166825 Å line), viewed at quadrature (tex2html_wrap_inline34185 (click here)). Shown are a) the polarized intensity, tex2html_wrap_inline3420 (scaled logarithmically over 3 decades), and b) the line-of-sight component of the velocity (scaled linearly over the range -50 to +150 km stex2html_wrap_inline3426)

  figure644
Figure 4: The polarization of tex2html_wrap_inline34286825 Å Raman photons for the reference model viewed at quadrature, overlaid on a logarithmic intensity scale (scaled over 3 decades). 100% polarization corresponds to a vector of 1 tex2html_wrap_inline3430 in length

  figure650
Figure 5: The polarization spectrum of the reference model (Sect. 5 (click here)), viewed at quadrature. The error bars are one sigma. Note that the two Raman lines are plotted on different scales

  figure655
Figure 6: Observations of the tex2html_wrap_inline34326825 Å line in a) V455 Sco, b) Hen 1092, c) Hen 1242, and d) RR Tel, adapted from Paper I (Harries & Howarth 1996b). The data shown here have been corrected for the interstellar polarizations derived in Paper I, to facilitate comparison of the percentage polarizations with the models. The polarized-flux and position-angle data in the line are insensitive to such corrections, aside from zero-point shifts

5.1. The tex2html_wrap_inline34346825Å line

The spatially-integrated polarization spectrum (Fig. 5 (click here)) is basically composed of three separate peaks (at tex2html_wrap_inline34366821, 6831, and 6836 Å in the tex2html_wrap_inline34386825 Å line), most obvious in the polarized flux, but also identifiable in the intensity spectrum. The polarization in the central peak is orthogonal to the other two. The origin of this characteristic is most easily understood by noting that the polarization-vector directions are centrosymmetric about the parent source (Fig. 4 (click here)). Because of the cancellation of the U parameter in the spatially integrated flux, the position angle in the models must have a value of 90tex2html_wrap_inline3442 or 0/180tex2html_wrap_inline3444 (within the errors of the numerical `noise'). The 90tex2html_wrap_inline3446 features originate `above' and `below' the parent-photon source (i.e., along directions perpendicular to both the sight-line and the line-of-centres); thus the spatial origin of a given peak can be deduced from its position-angle dependence.

The bluemost peak is clearly produced in the region between the two stars (the only region of the red-giant wind seen as blue-shifted by the OVI source in the geometry of the reference model). The polarization spectrum is characterized by very high polarization in the blue line wing (as high as 80%), although there is very little line flux there. The highly polarized blue line-wing is observed in many symbiotics after the removal of continuum flux (cf. Schmid & Schild 1994; Paper I; Fig. 6 (click here)).

The central (redshifted) peak has polarization position-angle orthogonally polarized to the bluemost peak, and thus must originate in the two quadrants of the red-giant wind `above' and `below' the parent-photon source. The highest-velocity peak is displaced from rest by tex2html_wrap_inline3450100 km stex2html_wrap_inline3452 in `parent-velocity' space. The velocity image (Fig. 3 (click here)) confirms that this emission arises from scattering across hemispheres, from close to the parent-photon source to the `back' of the red giant. This third, broad peak is more obvious in the intensity spectrum than the polarized flux (Fig. 5 (click here)), showing that the `back' of the wind first sees most photons after they have already undergone multiple scatterings, which produce an effectively spatially extended source, resulting in reduced polarization (Sect. 11 (click here)).

Although in our models all photon packets are forced to Raman-scatter, the weighting scheme (Sect. 2.4 (click here)) ensures that the number ratio of input OVI photons to output Raman photons is correct. In the reference model this ratio is of order 10:1 for tex2html_wrap_inline34561032:6085 - that is, tex2html_wrap_inline345890% of OVI photons escape, in agreement with the observations reported by Espey et al. (1995).

5.2. The tex2html_wrap_inline34627082 Å line

The tex2html_wrap_inline34647082 Å simulation (Fig. 5 (click here)) shows broadly similar characteristics to the shorter-wavelength Raman line, in agreement with observations (Paper I). The most obvious differences between the lines are the factor tex2html_wrap_inline34666.7 difference in Stokes' I (for an assumed 2:1 photon ratio for tex2html_wrap_inline3470:tex2html_wrap_inline3472), and the change in relative strength of the central (tex2html_wrap_inline34747090 Å) and redshifted (tex2html_wrap_inline34767095 Å) polarized-intensity peaks.

The observed tex2html_wrap_inline3478 intensity ratios are in the range tex2html_wrap_inline34802-10, averaging tex2html_wrap_inline34824-5 (Allen 1980; Paper I). The reference-model result is slightly larger than this average value, but this mild discrepancy could easily be resolved if tex2html_wrap_inline3484, or if tex2html_wrap_inline3486:tex2html_wrap_inline3488. In our models, the only fixed parameters distinguishing the tex2html_wrap_inline34906825 and 7082 Å lines are the Rayleigh- and Raman-scattering cross-sections for the parent OVI photons, which are smaller for the longer-wavelength line. The difference in line intensity in the reference model is therefore directly attributable to the escape of a larger fraction of tex2html_wrap_inline34941038 Å than tex2html_wrap_inline34961032 Å photons (by a factor of tex2html_wrap_inline34983).

This also accounts for the differences in line structure: the red-giant wind is quite optically thin perpendicular to the line of centres, so there is little scattered flux from those quadrants, or from the quadrants furthest from the red giant. The blue-shifted polarized flux again arises from the region between the two stars, while the redmost peak arises from multiple scattering, and is slightly more intense (relative to the other two peaks) than the tex2html_wrap_inline3500 model because the ratio of Raman:Rayleigh scatterings is slightly higher, meaning that the average Raman photon has undergone fewer scatterings (which reduce the `coherence' of the polarization signal) than in the corresponding tex2html_wrap_inline3502 model.

Clearly, observed differences between the tex2html_wrap_inline35046825 and 7082 Å lines encode differences in the governing system parameters (e.g., the separate values of tex2html_wrap_inline3506 and tex2html_wrap_inline3508 appropriate to the formation of the two lines). Such differences have important diagnostic potential in principle, but the broad similarity of the two lines in the reference model (as well as in observations; Paper I) suggests that the model sensitivity of response of the Raman lines to those parameters can, to first order, be adequately demonstrated by investigating either line separately. Moreover, the tex2html_wrap_inline3510 ratio depends directly on the input tex2html_wrap_inline3512 ratio, which is not well known; a 2:1 ratio is expected under optically-thin conditions, but the resonance doublets of CIV, NV, and SiIV commonly show smaller ratios in practice (e.g., Shore & Aufdenberg 1993 and references therein). Henceforth we therefore concentrate on the stronger, and more reliably observed, tex2html_wrap_inline35206825 Å line.

5.3. Preliminary comparison with observations

Figure 6 (click here) shows polarization data for the S-type systems V455 Sco, Hen 1092, and Hen 1242, and for the D-type system RR Tel. These illustrate part of the diversity of polarimetric behaviour reported in Paper I.

The blueward rise in polarization seen in the model spectra is not apparent in the observations, as dilution by unpolarized continuum flux (which is not included in the model calculation) means that the polarization tends to zero in the line wings. It is possible to subtract off the diluting flux in the data (cf. Figure 3 in Paper I), but small uncertainties in continuum placement lead to large uncertainties in the polarization, particularly in the line wings. As commented in Paper I, we therefore believe that the polarized flux the more straightforward observable against which to compare the models.

The data for Hen 1092 are broadly consistent with the results from the reference model in that they show a central peak in PA, `flipped' at tex2html_wrap_inline3522 to the adjacent regions, although the relative strengths of the polarized-intensity peaks differ from those in the reference model. The data for Hen 1242, which are of better quality, also clearly show the tex2html_wrap_inline3524 flip between the bluemost and central peaks, although the polarized intensity of the bluemost peak is again much less than the reference model predicts. However, in the observations shown in Fig. 6 (click here), the central and redmost peaks show the same PA. In both V455 Sco and RR Tel the relative strengths of the intensity and polarized-intensity peaks are in satisfactory agreement with the model, but the central and bluemost peaks show the same PA, with a tex2html_wrap_inline3526 degree flip in the redmost peak. Thus these stars show the bluemost peak flipped (Hen 1242), the central peak flipped (Hen 1092), and the redmost peak flipped (V455 Sco, RR Tel) - the last configuration being the commonest, present in at least half the spectra reported in Paper I as showing a `flip'.

Given the range in observed spectropolarimetric characteristics noted above and reported on further in Paper I, as well as the number of free parameters in the model, it is not unexpected that there should be some disagreements between the data and the reference calculation. In order better to understand these disagreements, we have performed calculations which relax some important basic assumptions of the reference model.


next previous
Up: Raman scattering in

Copyright EDP Sciences
web@edpsciences.com