A&A 450, 667-680 (2006)
DOI: 10.1051/0004-6361:20054255
F. Herpin1 - A. Baudry1 - C. Thum2 - D. Morris2, - H. Wiesemeyer2
1 -
Observatoire Aquitain des Sciences de l'Univers, Laboratoire d'Astrodynamique, d'Astrophysique et d'Aéronomie de Bordeaux, CNRS/INSU UMR N
5804, BP 89, 33270 Floirac, France
2 -
IRAM, 300 rue de la Piscine, Domaine Universitaire, 38406 Saint Martin d'Hères, France
Received 25 September 2005 / Accepted 29 December 2005
Abstract
Aims. We study the polarization of the SiO maser emission in a representative sample of evolved stars in order to derive an estimate of the strength of the magnetic field, and thus determine the influence of this magnetic field on evolved stars.
Methods. We made simultaneous spectroscopic measurements of the 4 Stokes parameters, from which we derived the circular and linear polarization levels. The observations were made with the IF polarimeter installed at the IRAM 30 m telescope.
Results. A discussion of the existing SiO maser models is developed in the light of our observations. Under the Zeeman splitting hypothesis, we derive an estimate of the strength of the magnetic field. The averaged magnetic field varies between 0 and 20 Gauss, with a mean value of 3.5 Gauss, and follows a 1/r law throughout the circumstellar envelope. As a consequence, the magnetic field may play the role of a shaping, or perhaps collimating, agent of the circumstellar envelopes in evolved objects.
Key words: masers - polarization- surveys - stars: late-type - stars: evolution - stars: magnetic field
The prodigious mass loss observed in numerous and widespread evolved stars make these objects the main recycling agents of the interstellar medium, and thus one of the most important objects in the Universe. Even though our knowledge of evolved stars has considerably improved over recent years, some of their main characteristics remain insufficiently understood (see the review by Herwig 2003): which mechanisms are responsible for their drastic change of geometry when evolving to the Planetary Nebula (hereafter PN) stage? What is powering the mass-loss so efficiently and could the magnetic field be playing a major role?
Important information about the physics and chemistry prevailing in the circumstellar envelope (hereafter CSE) of evolved stars can be retrieved from the radiowave line emission of molecules, especially from maser emission (see the review by Bujarrabal 2003). These envelopes can be probed at different depths through the study of three masing molecules, OH, H2O, and SiO. Our current knowledge indicates that:
Research on astronomical masers polarization is very active but is made difficult both by the lack of specific instrumental facilities and by the excitation and propagation of the masers themselves. Until now, numerous polarimetric observations of OH masers have been done, several of H2O masers, but few of SiO maser emission. Few SiO polarimetric observations have been done with VLBI giving the very first images of the magnetic field in some objects (e.g. Kemball & Diamond 1997, in TX Cam). Most of the early SiO studies were done in linear polarization. The first complete SiO polarimetric observations were performed by Johnson & Clark (1975), then by Troland et al. (1979); emission was found to be typically 15-30% linearly polarized and to exhibit no circular polarization. Barvainis et al. (1987) and McIntosh et al. (1989) measured circular polarization of 1-9% in several stars. Circular (0-4%) and linear (3.7-9.7%) polarizations were measured in VY CMa by McIntosh et al. (1994). Later, Kemball & Diamond (1997) made the first image of the magnetic field in the atmosphere of TX Cam, measuring a circular polarization level of 5% with some features showing polarization up to 30-40%.
It must be stressed that SiO, like H2O, is a non-paramagnetic species. Zeeman splitting exists but the sublevels overlap; the effect is thus undetectable, so only net polarization can be used to trace the magnetic field. The current status of our knowledge of the magnetic field strength can be summarized as follows:
An electromagnetic plane wave is defined by two components (horizontal and vertical):
![]() |
(1) |
![]() |
(2) |
Its energy flux is described by the 4 Stokes parameters:
![]() |
(3) |
![]() |
(4) |
![]() |
(5) |
![]() |
(6) |
![]() |
Figure 1:
V, U, Q, and I Stokes parameters for R Leo. I, V, U,
and Q are given in Kelvins (
![]() |
Open with DEXTER |
We present here spectroscopic measurements of the 4 Stokes parameters (see Fig. 2). The observations were made with the IF polarimeter installed at the IRAM 30 m telescope on Pico Veleta, Spain (Thum et al. 2003). Simultaneous measurements of I, U, Q, V allow us to calculate I, ,
and
for each velocity channel. The polarization angle calibration (i.e. the sign of Stokes U) was verified by observations of the Crab Nebula. Moreover, planets (polarization of planets is negligible at our frequency) have been used to check the instrumental polarization on the optical axis.
The instrumental beam polarization is known to be stronger in Stokes Q and U than in Stokes V (known to be 2-3%, see Thum et al. 2003, comparable to our sensitivity as stated elsewhere). If some detections from sources with weak
are from a bad or uncertain pointing, they naturally induce a value of
that is weaker than
.
A strong instrumental polarization in Stokes V would instead be due to a bad phase tracking (the IF polarimeter works in a manner quite similar to that of an adding interferometer, and good phase tracking is essential). From several tests (Thum et al. 2003; Wiesemeyer et al. 2004), we know that the polarization seen for weak SiO components with (Q, U, V) = (+ - -), (- - +) or (- + -) is instrumental polarization. We see that signature for only 3 objects (R Crt, R UMi and
RT Vir). Some instrumental polarization may thus contaminate the observations of these objects. All instrumental parameters were carefully calibrated through specific procedures described in Thum et al. (2003). The error on
is
%.
The SiO (v=1, J=2-1) line observations at 86.243442 GHz were carried out towards 57 stars in August and November 1999 with the IRAM 30 m radiotelescope. The pointing was regularly checked directly on the star itself (for the vast majority of objects). In order to obtain flat baselines, we used the wobbler switching mode. The system temperature of the SIS receiver ranged from 110 to 170 K. The front-ends were the facility receivers A100 and B100, and the back-end was the autocorrelator. The lines were observed with a spectral resolution of 0.3 kms-1. The integration times were 4-10 min using the wobbler switching. The forward and main beam efficiencies were 0.92 and 0.77 at 3 mm, respectively. (Additional SiO (v=1, J=5-4) line observations at 215.596 GHz were also performed in most stars studied here; results will be reported elsewhere.)
Our source sample (see Table 1) consists of 43 Miras, 7 Semi-Regular stars (hereafter SR), 2 IR late-type stars, 1 irregular variable, 3 supergiants, and 1 Planetary Nebula (QX Pup), all selected from our SiO maser master catalogue (Herpin & Baudry, private communication). Coordinates and the main characteristics of the objects are given in Table 1. Nearly 60% of stars in this table have been observed with the HIPPARCOS satellite and thus have excellent optical positions; such positions were adopted in our work.
Table 1: Stars observed in this work. The stellar type is derived from the literature as are the mass loss rates (e.g. Loup et al. 1993) and the period (e.g. AAVSO data).
The values of the polarization level presented here (see Table 2) are those measured for the different components within the SiO maser emission profile for each star. Examples are given in Fig. for a few stars. We only considered the well-identified components (distinct peaks or strong wing emission separated from the bulk emission well, according to the noise). Some interesting cases are briefly presented below.
Some profiles show isolated emission red/blue-shifted from the main emission, and are more strongly circularly polarized (e.g. IRAS 18204-1344). These peculiar characteristics imply a different spatial origin for the main and higher/lower velocity components. Sometimes the circular polarization varies regularly across the profile (e.g. R Leo), but sometimes not.
In T Lep, the SiO emission shows two peaks linked by a plateau; the circular polarization varies linearly across the profile from -11 to 10% (see Fig. ). Several objects (e.g.
S UMi, IK Tau) show the same
pattern.
The IR late-type source IRAS18158-1527 exhibits a complex profile with several well defined components, each of them differently polarized indicating a complex maser structure probably with different maser spots contributing to the whole emission. The red wing emission is highly polarized (43%). Such a complex multi-component maser line profile and "semi-circle'', convex,
pattern appear to be characteristic of SR objects (see other similar objects in our sample and R Crt in Fig.
). Nevertheless, the circular polarization pattern observed in the Mira star U Lyn is a convex profile like those encountered in SR objects.
One of the most studied Mira stars is R Leo. Its profile is made of a strong emission with a blue broad line wing. Main and linewing emissions are strongly polarized, negative and positive
%, respectively. R Leo is a very well studied object exhibiting a bipolar jet throughout its envelope. The clear symmetry observed between the positive and negative circular polarization patterns in the main and wing line emissions suggests that the maser emission comes from the jet lobes. We note that the Mira star RS Lib exhibits an emission and polarization pattern similar to that observed in R Leo, so they may have the same spatial structure.
The circular polarization level in several of our objects has already been measured by Kemball & Diamond (1997) or Barvainis et al. (1987). For TX Cam and W Hya, our results are consistent with previous observations:
We emphasize that with a single dish radiotelescope we cannot spatially distinguish between the different maser spots producing the SiO profile (various masers spots contribute in the various features observed at a given velocity). The whole SiO maser emission region, hence all the maser cells, lie within the 29 arcsec of the 30 m (but not necessarily with a uniform distribution) while the SiO emission covers less than 40 milliarcsec in TX Cam (Kemball & Diamond 1997) so about everything is beam averaged. This means that any conclusion on the geometry of the objects observed here would be very uncertain. Only global trends or global geometry can be discussed. One of the consequences of this spatial resolution problem is that, if the polarization vectors are distributed isotropically around the object, the average polarization level that we measure is zero, even if the maser emission produced in each SiO cell is polarized well.
A global analysis of our data in Table 2 shows the following. We find that
varies between 0 and 70%, and
between 0 and
43%; hence, polarization vectors are not distributed isotropically. Emission from Mira-type objects clearly tends to have a relatively high linear (
%,
%) and circular polarization
(
%,
%). Note that the emission from the PN QX Pup is highly polarized, and, in contrast, maser emission from supergiants shows very weak polarization (
%,
%), with the exception of one maser component in S Per. Moreover, all observations show that the polarization level varies across the maser line profile (see Fig.
), i.e. the different spectral components of the maser emission producing the profile come from different localizations in the SiO shell and have different polarization levels. The highest polarization level for one object can be encountered either in the main peak or in the other components.
Semi-Regular objects (RU Aur, R Crt, W Hya, S Per, AH Sco, R UMi, RT Vir) have a common circular polarization pattern with the central main emission unpolarized and other peak emission or wings strongly polarized: a characteristic "semi-circle'' (i.e. convex shape) pattern for
is observed (see R Crt in Fig.
). The infrared late-type star IRAS18158-1527 exhibits a similar pattern, thus suggesting that this star is a semi-regular object.
A group of objects (W And, NV Aur, T Cas, R Com, T Lep, IK Tau, S Ser, S UMi) shows approximately the same
pattern (see T Lep in Fig.
); the circular polarization varies linearly across the line profile from a positive value to a negative one (or the contrary). The only common spectral characteristic of the SiO emission from these stars is the presence of an plateau-like emission on top of which the narrow emission peaks are located.
![]() |
Figure 2:
Position angle of polarization (![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
In this section, we first discuss our source sample in the frame of the 2-color diagram. Then we briefly summarize the existing SiO maser polarization theories. Finally, we discuss our data set in this context and estimate the stellar magnetic field strength.
Stars of our sample can be plotted in a [12] - [25], [25] - [60] color-color diagram (van der Veen & Habing 1988; [12], [25], and [60] stand respectively for 12, 25, and 60 microns IRAS-fluxes). This diagram is partitioned into several regions (see Fig. 3) defined by van der Veen & Habing as follows: Region I, oxygen-rich non-variable stars without circumstellar shells; Region II, variable stars with young O-rich circumstellar shells; Region IIIa, variable stars with more evolved O-rich circumstellar shells; Region IIIb, variable stars with thick O-rich circumstellar shells; Region IV, variable stars with very thick O-rich circumstellar shells; Region V, Planetary Nebulae and non-variable stars with very cool envelopes; Region VIa, non-variable stars with relatively cold dust at a large distance; Region VIb, variable stars with relatively hot dust close to the star and relatively cold dust at a large distance; Region VII, variable stars with more evolved C-rich circumstellar shells.
Figure 3 represents the linear and circular polarization level for the main SiO emission component from each star. Most of the objects in our sample fall in regions II and IIIa and do not show particular characteristics, except for S Cas (an S-type star) where the circular polarization is high. Mira-type stars are in regions I, II, IIIa, and VII. The IR late-type objects are in VIb and VII. The SRa semi-regular variable W Hya is in I, while the SRb stars lie in IIIa (RU Aur, RT Vir, R Crt), II (R UMi), and Src in IIIa (S Per) and VII (AH Sco). The Red Supergiants, VY CMa and VX Sgr, and the IR supergiant IRAS 18204-1344 lie in VII. Objects in Region VII do not exhibit strong polarization compared to other objects, perhaps because of their more C-rich circumstellar shells (e.g. AU Aur) or because of the presence of hot dust close to the star, implying less SiO abundance and thus weaker emission, making the polarization measurement less significant. The presence of hot dust may also influence the pumping of the SiO molecules and thus the polarization level; the optically thick, hence isotropic, radiation field of hot dust can assist the collisional pumping. This could apply to UX Cyg (an irregular variable), IRAS 18055-1433, and IRAS 18158-1527 in region VIb. In region V QX Pup is a PN that exhibits strong polarized emission. Note that IRAS 18055-1433 and IRAS 18158-1527 show very strong circular polarization in their line wings. We may conjecture here that wing emission comes from more outer layers than those where the main line is excited (Herpin et al. 1998); as a consequence, the SiO cells giving rise to wing emission are less influenced by the presence of hot dust (preferably lies in the inner layers).
For more than one decade, two opposing schools have grown up to explain SiO maser emission. SiO polarization theory is described in: (i) Watson (e.g. Watson & Wyld 2001; Wiebe & Watson 1998; Nedoluha & Watson 1994); (ii) Elitzur (2002, 1998, 1996, 1994).
![]() |
Figure 3: Two-color diagram for our sample. The size of the empty and filled circles for each star are proportional to the linear and absolute circular polarization level, respectively, measured for the main peak of SiO emission. |
Open with DEXTER |
The main difference between the two approaches rests in the pumping mechanisms. While anisotropic pumping associated with a weak field produces high
and quite significant
in Watson's model, a strong magnetic field is necessary with the more classical pumping mechanisms used in Elitzur's model. Details about both models can be found in the Elitzur review (2002).
We may summarize the main characteristics of Watson's model as follows:
![]() |
(7) |
![]() |
(8) |
![]() |
(9) |
![]() |
(10) |
![]() |
(11) |
![]() |
(12) |
Table 2:
Derived parameters of the different components of the SiO
maser emission profile for each star. Only the well-identified components are given (distinct peak or strong wing emission separated from the bulk emission). Note that the polarization is fractional. The
is the rms derived from the
plot.
We may look for any correlation between
and
.
As shown in Figs.
and 3, it first appears that in all cases
is larger than
as predicted by all models. More precisely,
is noticeable if and only if
is high. If we plot the values of
and
derived for all maser components in this work (Fig. 4) and make a regression fit to our data, we obtain
.
The circular polarization level tends to vary approximately linearly with
,
in agreement with the fact that no
is detected towards sources with a marginal
detection. Note that in Fig. 4, four objects, S Per, IRAS18055-1433, S Cas, and
Hya, do not follow the same general trend observed for the rest of the sample. (The case of the Supergiant S Per, however, is an exception as it exhibits substantial
while
.) This observation may in fact favor Watson's model. It must be stressed again that the beam-averaged polarization that we measure makes any conclusion uncertain. In fact, due to this averaging, we should observe no correlation at all, even if one would exist!
As mentioned earlier (see Sect. 3.2) we observe relatively high circular polarization rates in several stars (
%,
%). These values are larger than those predicted from Watson's model (e.g. Nedoluha & Watson 1994). We also have not been able to find any correlation of
with total intensity. Finally, although we adopt the Zeeman case to derive the magnetic field strength (see Sect. 4.5), we cannot conclude firmly from present observations which maser theory prevails for SiO emission.
![]() |
Figure 4:
Absolute fractional polarization ![]() ![]() ![]() |
Open with DEXTER |
Some recent studies tend to demonstrate the importance of the magnetic field in evolved objects. Bujarrabal et al. (2001) show that for 80% of the PPNe in their sample the fast molecular flows have momenta that are too high to be powered by radiation pressure alone (1000 times larger in some cases), which may be explained by magnetic field. Moreover, X-ray emission found in evolved stars (e.g. Hünch et al. 1998) may indicate the presence of a hot corona that possibly results from magnetic activity. Very recently, a magnetic field was discovered for the first time in central stars of PN (Jordan et al. 2005) and estimated to kiloGauss, much stronger than what we find here from our SiO data in QX Pup.
New models involving the magnetic field have been developed to try to explain the
morphological changes of an object during its transition from the AGB stage to the PN stage; B plays the role of a catalyst and of a collimating agent. The simplest models are based on a moderately weak magnetic field alone (
Gauss at the stellar surface, a few 1013 cm, i.e. at a radiis of a few AU, Soker 1998). The influence of B is stressed by the work of Smith et al. (2001) and Greaves (2002) in VY CMa. But the role of B can only be decisive when its energy density is greater than the radiative pressure, i.e. when B is greater than around 10 G close
to the stellar surface in the SiO region (see Soker & Zaobi 2002).
Arguing that such a strong field may be very unusual, Balick & Frank (2002) explain that B alone cannot produce the observed structures, and a combination of several factors thus has to be considered (rotation, magnetic field and presence of a companion).
Soker & Harpaz (1992) first proposed a model with a weak magnetic field ( G) and included slow rotation together with the presence of a companion to transform the envelope (and lead for example to the peculiar geometry observed in NGC 6826 or
NGC 6543). Even if the star is not binary, B can have a great influence locally (Palen & Fix 2000). A significant magnetic field can form cold spots on the star's surface and a slow rotation of the star can then increase the field strength to build up a dipolar magnetic field varying as 1/r3 (Matt et al. 2000); such a field is stronger at the equator and may thus lead to an axisymmetrical mass loss.
The main argument against the dominant influence of the magnetic field on the shaping of the circumstellar envelope is that a strong field seems to be necessary to dominate the dynamics of the gas. However, several authors (Pascoli 1985, 1992, 1997; Chevalier & Luo 1994; García-Segura 1997; Gurzadyan 1997;
Delamarter 2000) have demonstrated the strong influence of a reasonable toroïdal magnetic field embedded in the normal radiation-driven stellar wind (magnetic wind bubble theory, hereafter MWB). This field has a strength between a few Gauss and a few 10 Gauss at a few stellar radii (the SiO region is believed to be at
cm or
AU), varies as 1/r2, then as 1/r at larger radii; therefore
mG at
1016,17 cm or 700-7000 AU. These results are confirmed by the simulations of García-Segura et al. (2001) for the PN He 2-90. Even if the origin of the wind is not explained by these models, it seems clear that a magnetic field is essential for generating fast collimated outflows (Kastner et al. 2003).
There are many models of magnetic jet production and collimation and some, or all of them, are applicable to various star geometries. One most interesting study was performed by Blackman et al. (2001) in which the magnetic field emerges from the AGB stellar core and the resulting 1 G field helps to collimate the radiation-driven wind or a stronger, more anisotropic, magnetically driven wind.
The exact interpretation, in terms of magnetic field, of our observations depends on the adopted specific SiO maser model (see Sects. 4.2 and 4.3). From the current knowledge of the strength of the magnetic field in the OH and H2O layers, we expect B// of a few Gauss at least in the SiO maser region (see also Kemball & Diamond 1997), i.e. at 5-10 AU from the central object. This tends to invalidate Watson's model, and furthermore tends to agree with a field varying in r-1,2 as predicted by Elitzur's model. However, Vlemmings et al. (2005) measured the circular polarization of the H2O maser emission in a few evolved stars with the VLBA observations and showed that the magnetic field is either a solar-type field (with a r-2 field strength dependence) or a dipole magnetic field (with a r-3 dependence) in their sample.
Table 3:
Average magnetic field strengths derived from .
In the following, we decide to use Elitzur's theory (Zeeman case) to infer magnetic field strength from the circular polarization levels. From Eq. (12), we thus calculate the mean value of the magnetic field B// for each star and give results in Table 3.
For our sample, B// is between 0 and 20 Gauss with a mean value of 3.5 G. This value, combined with the strength of the field in more outer layers of the envelope (OH and H2O masers), agrees with a B-field variation law in 1/r, closer to Elitzur's model. As explained in the Introduction and in Sect. 4.4, B alone can be the main agent to shape the circumstellar envelope if its value is higher than around 10 Gauss. This means that only S Cas, RU Aur, IRAS 18055-1433, and IRAS 18158-1527 may have a magnetic-field ruled geometry (B> 10 G in these objects). The rest of our sample shows that B is strong enough to be dominant at this stage of the AGB star evolution, but it should be associated with rotation and the presence of a companion as suggested in models mentioned in Sect. 4.4. Even though our B measurements are beam-averaged, they suggest in many cases that they are not too much (not orders of magnitude) below the critical value; local B values may exceed the critical value in many cases, and therefore participate in the shaping of the AGB envelopes. Our estimated values of B are consistent with the MWB theory (toroïdal magnetic field) or the model of Blackman et al. (2001).
From Elitzur (1996) and our measurements
%,
%, we can estimate that
is around 0.017 and
,
respectively, for Mira-type objects and semi-regular variables (see Eq. (9) in Sect. 4.2). Moreover, according to Elitzur (1996, see Fig.
), there are no stationary physical solutions for propagation at
;
i.e. at
the radiation is not polarized. As xB is small (<0.02), from Fig.
of Elitzur (1996) we can estimate the volume of phase space in which propagation of linear polarization in a maser is possible or not (
or
), and then calculate that the probability for a random magnetic axis to be aligned with a given direction (our line of sight) is higher than
.
Our present estimate is around 18%. Therefore, Elitzur's model predicts that 18.4% of the SiO 86 GHz masers should not be linearly polarized, because such polarized masers cannot propagate if the magnetic field, although weak, is closer than
to the line of sight (propagation direction). Hence non-polarized maser emissions do not imply that there is no magnetic field or only a weak one. In our sample, roughly 13% of the SiO maser components have no detectable or very weak (
)
polarization.
We looked without success for a possible correlation between the polarization rates and physical parameters, such as the known envelope asymmetry, the presence of SiO-maser high-velocity linewings (see Herpin et al. 1998), or the mass loss rate. If the magnetic field plays an important role in the shaping of the object, one may expect to find a relationship between the strength of B// (thus )
and the geometry of the object. Unfortunately, no trend is clearly found in our data. Nevertheless it is known that radiative pressure is driving the wind in AGB objects, so it is not surprising that we find no correlation between the B strength and a known asymmetry in our sample. Of course, our stellar sample would require new observations with sufficient spatial resolution (VLBI) to confirm the present results; the same type of study should also be conducted toward several Proto-PN and PN objects.
We have made a study of the SiO maser polarization in a representative sample of evolved stars, simultaneously measuring, for the first time, the 4 Stokes parameters. From our measurements we derive the circular and linear polarization levels and show that, due to the beam averaging of our polarization measurements, we cannot firmly discriminate between the two dominant theories of SiO maser emission. In particular, VLBI observations of our source sample are absolutely necessary to distinguish between Zeeman or non-Zeeman theories. Nevertheless, the magnetic field strength was derived assuming Elitzur's model. The value
varies between 0 and 20 Gauss, with a mean value of 3.5 G. As a consequence, we suggest that the magnetic field plays a significant role in the evolution of these objects. Within the frame of the Zeeman theory, the magnetic field could shape or even collimate the gas layers surrounding the AGB objects. Emission from Mira-type objects clearly tends to have a higher linear (
%,
%) and circular polarization (
%,
%). Basically, if there is a real correlation
between
and the strength of the magnetic field, this trend may indicate that the magnetic field may be stronger in Mira objects than in Semi-Regular variables, at least in the inner layers of the circumstellar envelope.
To better understand the mechanisms at work with the magnetic field, complementary studies have to be conducted, and in particular, the presence of a companion has to be investigated in a large sample of objects. Of course VLBI maps of the magnetic field in these stars are essential. Another important objective is to investigate the evolution of the magnetic field and its influence during the transition from the AGB star phase to the PN stage.
Acknowledgements
The authors are grateful to M. Elitzur for reading and commenting on this paper. We also thank W. D. Watson for his useful comments and suggestions. The authors are indebted to the staff of the IRAM 30 m telescope who most efficiently helped during the observations and to R. Mauersberger who closely followed part of these observations. Finally, we also thank the referee for several useful comments.
![]() |
Figure 2:
Position angle of polarization (![]() ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |