text:
# Comprehensive modelling of the planetary nebula LMC-SMP 61 and its [WC]-type central star

## 1 Introduction

Only a few studies have been devoted so far to a consistent modelling of a planetary nebula and of its central star. Such studies are useful to get a better insight into the relation between the nebula and its progenitor. Another, very important aspect is that this is the only way to test model atmosphere predictions in the Lyman continuum and thus to validate the model atmospheres. The works of Rauch, Köppen & Werner (1994, 1996), Peña et al. (1998), De Marco & Crowther (1998, 1999), De Marco et al. (2001) are examples of such studies, while Crowther et al. (1999) have performed a similar study on a Population I Wolf-Rayet ring nebula. Such investigations are particularly important in the case of planetary nebulae with Wolf-Rayet type central stars (which represent about 10% of all planetary nebulae), since recent work (e.g. Górny & Tylenda 2000, De Marco & Soker 2002) have completely changed previous views on the evolutionary status of these objects.

In general, Wolf-Rayet central stars of PNe belong to the \[WC\] sequence. In our galaxy most of these objects have been classified as \[WC-early\] or \[WC-late\] types, with only few objects of intermediate types (Tylenda, Acker & Stenholm 1993). In the Magellanic Clouds, the WR central stars are also of \[WC\] type (except for the extraordinary central star of LMC-N66, e.g. Peña et al. 1997b), but in this case they have been classified in the intermediate \[WC\] types. Peña et al. (1997a) suggested that such a difference might be a consequence of the differences in metallicity between the Milky Way and the Magellanic Clouds. In any case, \[WC\] central stars show spectral features identical to those of massive WC stars but at much lower luminosity, and they can be analyzed with the tools developed for massive WR stars (e.g. Hamann 1997).

In the present paper, we concentrate on the planetary nebula SMP 61 (also known as N203, WS 24 and LM1-37) which is among the brightest planetary nebulae in the Large Magellanic Cloud. This is a good case for a detailed study: The object is at a known distance modulus of 18.50 mag (Benedict et al. 2002), equivalent to 50.1 kpc. The central star is one of the brightest \[WC\] central stars in the LMC and therefore relatively easy to observe. It is of \[WC 5-6\] type (Monk, Barlow & Clegg 1988; Peña, Ruiz & Torres-Peimbert 1997a), therefore the nebular spectrum contains many lines of various excitation levels that allow refined diagnostics. In addition, the nebula appears to be spherical and integrated spectra are already available (Peña et al. 1997a).

Using HST STIS, we have secured high signal-to-noise spectra of SMP 61 in a wide spectral range. This provided strong constraints for our modelling of the central star atmosphere. The best fit model atmosphere was then used as an input to build a photoionization model for the nebula. The observations are described in Sect. 2, the stellar modelling technique is presented in Sect. 3, the spectral analysis of the central star of SMP 61 is discussed in Sect. 4. The nebular modelling strategy is exposed in Sect. 5 and the nebular model fitting of SMP 61 is presented in Sect. 6. The implications of our modelling are discussed in Sect. 7, and the main points of this study are summarized in Sect. 8.

## 2 New Observations and Previous Data

### 2.1 Observations

HST STIS spectroscopic data were obtained on 1998/12/30 (HST Cycle 7, program ID 7303). MAMA and CCD detectors were employed to cover the broadest spectral range possible. Prior to the spectroscopic observations, STIS images were obtained to acquire the target. The log of the observations is presented in Table 1.

Direct images with HST WFPC were searched in the archives. Apart from our own STIS images for target acquisition (see Table 1), we retrieved the FOC image W1ID0601T, obtained on 1993/11/08, with exposure time of 3 min through the filter F502N which isolated the strong \[O iii\] $`\lambda 5007`$ line (this image is shown in Vassiliadis et al. 1998b). We decided to use our target acquisition images to analyze the nebular morphology and surface brightness because they are better exposed and because the Vassiliadis et al. images were obtained before the COSTAR mission. Our images are obtained with a filter including the brightest emission lines: H$`\alpha `$, H$`\beta `$ and \[O iii\] lines.

We also considered ground-based spectrophotometric data for SMP 61 obtained by Monk, Barlow, & Clegg (1988), Meatheringham & Dopita (1991) and Peña et al. (1997a). Other available UV spectroscopic data are IUE data given in Peña et al. (1997) and HST FOS data from 1175 to 4800 Å published by Vassiliadis et al. (1998a).

All the available spectroscopic data are presented in Table 2, which lists the reddening corrected fluxes for the observed lines (in units of H$`\beta `$=100) obtained by the different authors. We also list the total H$`\beta `$ flux observed in the slit, $`F(\mathrm{H}\beta )`$, the logarithmic extinction at H$`\beta `$, $`C(\mathrm{H}\beta )`$, and the slit dimensions, when available. All the spectra mentioned above have good signal-to-noise, and the differences among observations by different authors are, in general, small. Note that observations by Peña et al. (1997a) were made both with a slit of 2″and a slit of 10″and, as noted by these authors, the line ratios are very similar in both slits. The H$`\beta `$ flux reported in column 2 of Table 2 is the one obtained with the 10″slit.

SMP 61 is one of the twelve planetary nebulae in the LMC detected by IRAS. The fluxes retrieved from the IRAS data base are: F(12$`\mu `$m) = 0.08 Jy (flux density quality = 2) F(25$`\mu `$m) = 0.13 Jy (flux density quality = 3), F(60$`\mu `$m)$``$0.16 Jy (flux density quality = 1).

### 2.2 Preliminary analysis of the observational data

In HST images the nebula appears perfectly spherical and shows no small-scale structure. The surface flux distribution, as derived from the image o57n02iuq, is shown in Fig. 1. The ionizing star dominates the emission in the central pixels. The nebula seems to have an inner radius of about 0.065″, showing the maximum surface brightness there. The surface brightness declines outwards and drops below 1/10 of its maximum at 0.26″.

The extinction of the nebula is modest. Throughout the paper, we use the value $`C(\mathrm{H}\beta )`$ = 0.18 (corresponding to $`E_{\mathrm{B}\mathrm{V}}=0.122\mathrm{mag}`$ derived by Peña et al. (1997a) from the observed Balmer decrement and we adopt the Seaton (1979) reddening law.

The STIS observations were essentially obtained to provide constraints on the central star. However, the 0.2″$`\times `$52″STIS slit includes a significant fraction of the nebula so that both stellar and nebular emission is seen. The nebular emission must of course be subtracted from the stellar emission in order to perform the stellar analysis (see section 4.1). Most of the emission lines detected in the UV are of stellar origin, except for the intense C iii\] 1909 Å, C ii 2321,26 Å, and \[O ii\] 2470 Å. Many bright nebular lines were detected in the optical region. The dereddened fluxes of all the nebular lines measured in the STIS calibrated spectra are presented in the first column of Table 2.

The STIS nebular lines together with the data presented by Peña et al. (1997a, column 2 in Table 2) have been used to derive first estimates of the physical conditions and of the chemical composition of the ionized gas. The atomic data are the same as used in the nebular modelling (see Sect. 5). The density derived from \[S ii\] $`\lambda 6731/6717`$ is 5400 cm<sup>-3</sup>, the temperature derived from \[O iii\] $`\lambda 4363/5007`$ is 10300 K, and the one derived from \[N ii\] $`\lambda 5755/6584`$ is 11400 K. These values were used to derive the ionic abundances, and the ionization correction factors from Kingsburgh & Barlow (1994) led to the following total abundance ratios : H : He : C : N : O : Ne : S : Ar = 1: .104 : $`\mathrm{7.35\hspace{0.17em}10}^4`$ : $`\mathrm{3.95\hspace{0.17em}10}^5`$ : $`\mathrm{2.65\hspace{0.17em}10}^4`$ : $`\mathrm{4.26\hspace{0.17em}10}^5`$ : $`\mathrm{6.0\hspace{0.17em}10}^6`$ : $`\mathrm{8.0\hspace{0.17em}10}^7`$. These abundances are employed in Sect. 6 as starting values for the photoionization modelling.

Regarding the IRAS data, Zijlstra et al. (1994) commented that the IRAS colors F(12)/F(25)$``$0.62 and F(25)/F(60)$``$0.81 of SMP 61 are very blue and typical of a young nebula.

## 3 Stellar modelling technique

### 3.1 Model atmospheres

The analysis of hot stars with strong winds like the \[WC\]-type central star of SMP 61 requires the application of model atmospheres that account for complex model atoms in non-LTE. Recent investigations revealed the necessity to include iron line-blanketing as well as the effect of density inhomogeneities (clumping) for a reliable determination of their stellar parameters and emergent flux distributions (Hamann & Koesterke 1998; Hillier & Miller 1999; Crowther et al. 2002; Gräfener et al. 2002).

In the present work, we utilize the Potsdam code for expanding atmospheres in non-LTE (Koesterke et al. 1992; Hamann et al. 1992; Leuenhagen & Hamann 1994; Koesterke & Hamann 1995; Leuenhagen et al. 1996; Hamann & Koesterke 1998; Koesterke et al. 2002; Gräfener et al. 2002). This code calculates the radiative transfer in the co-moving frame of reference for a spherically symmetric, stationary outflow. The atomic populations, the electron density, and the electron temperature throughout the extended atmosphere are determined from the equations of statistical- and radiative equilibrium. The system of statistical equations is solved in line with the radiation transport by application of the ALI formalism (accelerated lambda iteration, see Hamann 1985, 1986), whereas radiative equilibrium is obtained by a temperature correction procedure which is based on the method of Unsöld (1955) and Lucy (1964) (see Hamann & Gräfener 2003; Mihalas 1978, p.174). The statistical equations are solved for complex model atoms of He, C, O, Si, and the iron group, taking advantage of the concept of super levels for the inclusion of millions of iron-group line transitions (see Gräfener et al. 2002). Density inhomogeneities are accounted for by the assumption of small-scale clumps with a constant volume filling factor $`f_V`$ (Hamann & Koesterke 1998).

### 3.2 Model parameters

The model atmospheres are specified by the luminosity and radius of the stellar core, by the chemical composition of the envelope, and its density- and velocity structure. The basic parameters are: The stellar core radius $`R_{}`$ at Rosseland optical depth $`\tau _\mathrm{R}=20`$, the stellar temperature $`T_{}`$ (related to the luminosity $`L_{}`$ via the Stefan Boltzmann law), the chemical composition (given by mass fractions $`X_{\mathrm{He}}`$, $`X_\mathrm{C}`$, $`X_\mathrm{O}`$, $`X_{\mathrm{Si}}`$, and $`X_{\mathrm{Fe}}`$), the mass-loss rate $`\dot{M}`$, the terminal velocity $`v_{\mathrm{}}`$, and the “clumping factor” $`D=1/f_V`$ giving the density enhancement of the clumped matter with respect to the mean wind density. For the dependence of the velocity on radius, a relation of the form $`v(r)=v_{\mathrm{}}\left(1R_0/r\right)`$ is assumed. This outer velocity law is augmented by an exponentially decreasing density distribution in the hydrostatic domain. $`R_0`$ is suitably determined to connect both domains smoothly.

The parameters $`R_{}`$, $`\dot{M}`$, $`v_{\mathrm{}}`$, and $`D`$ are connected by the “transformed radius” $`R_\mathrm{t}R_{}(v_{\mathrm{}}/\sqrt{D}\dot{M})^{2/3}`$. Models with the same $`R_\mathrm{t}`$ show almost identical line equivalent widths (Schmutz et al. 1989). This invariance holds because the line emission in WR stars is dominated by recombination processes (see Hamann & Koesterke 1998).

## 4 Spectral analysis of the central star

In the present section, we perform a detailed spectral analysis of the central star of SMP 61 based on HST UV + optical observations. Owing to the reliable spectrophotometry of HST in combination with the known distance and the moderate interstellar reddening towards the LMC, the observations provide the absolute stellar flux distribution in the wavelength range from 1150 to 6000Å. The model comparisons in the present section are therefore performed with respect to the absolute flux, corrected for interstellar extinction applying the standard law of Seaton (1979) with $`E_{\mathrm{B}\mathrm{V}}=0.122\mathrm{mag}`$.

### 4.1 Subtraction of the nebular continuum

Due to the small angular diameter of the nebula ($`0.7\mathrm{}`$), the HST spectra ($`0.2\mathrm{}`$ slit width) include a significant part of the total nebular flux. Consequently, the observations clearly show nebular features superposed to the stellar spectrum. The most striking of these is the Balmer jump at 3646 Å in emission, indicating the dominance of the nebular continuum in the UV to blue part of the observation. A reliable subtraction of this continuum emission is essential for the spectral analysis of the central star.

A distinction between both contributions is in principle possible, because the nebula is spatially extended in contrast to the point-like star. To estimate the contribution of the nebula, two spectra with apertures of $`1.2\mathrm{}`$ and $`0.25\mathrm{}`$ are extracted from the UV part of the HST data. The difference between these two spectra is dominated by nebular emission. After scaling the difference spectrum to the flux of the nebular C ii line at 2327 Å, a coarse estimate for the pure nebular flux is obtained.

In Fig. 2, this spectrum is compared to the theoretical nebular continuum flux computed with the physical conditions relevant to the object. After applying the appropriate reddening the model flux is scaled by the ratio of observed to theoretical $`\mathrm{H}_\beta `$ line flux. Finally, the theoretical continuum is utilized for the nebular subtraction. Due to the good agreement between theory and observation, a possible stellar contribution to the observed Balmer jump can be excluded.

### 4.2 Spectral fit

The central star of SMP 61 shows a typical \[WC 5\]-type spectrum dominated by emission lines of He, C, and O. Apart from a narrower line width, the spectral appearance of SMP 61 is similar to the well-known massive WC 5 star WR 111. Consequently, the fit criteria and the fitting technique are similar to the analysis of WR 111 described by Gräfener et al. (2002). An important difference between both objects concerns the low iron abundance of SMP 61. The iron forest around 1300 Å – a combination of a large number of iron emission and absorption lines – is barely visible in the HST observations of this object.

The stellar analysis is based on fitting the absolute flux distribution including the observed wind emission lines of He ii, C iii, C iv, O iii, O iv, and O v. A hydrogen-free WC surface composition with a metallicity of $`0.3`$ solar is adopted (referring to the compilation of solar element abundances in Gray (1992, p. 318)). This value matches the derived nebular oxygen abundance. The iron mass fraction is correspondingly set to $`X_{\mathrm{Fe}}=\mathrm{4\hspace{0.17em}10}^4`$, and the silicon abundance to $`X_{\mathrm{Si}}=\mathrm{2\hspace{0.17em}10}^4`$. The other model parameters are kept free throughout the fitting procedure. The hydrogen abundance is set to zero. This assumption cannot be checked empirically, because the strong nebular Balmer emission masks any possible stellar contribution. However, in previous works no hydrogen is found on the surface of early-type \[WC\] stars, and even for late types only small amounts of $`X_\mathrm{H}<0.1`$ are detected (Koesterke 2001; De Marco & Barlow 2001).

The final model fit is presented in Fig. 3 together with the derived stellar parameters in Table 3. As expected, the stellar temperature, transformed radius, and surface composition are similar to the values obtained for WR 111 (Gräfener et al. 2002), while the luminosity is lower by 1.5 dex. A good fit quality is obtained with some exceptions like C iii 2297 Å (a transition that is mainly fed by dielectronic recombination), the classification line C iii 5696 Å, or the He ii/C iv blend at 5412/5470 Å. The latter reacts very sensitively to changes of the ratio of carbon to helium abundance. Additionally, our model with standard iron abundance produces too strong iron features in the UV.

Test calculations for different iron abundances are presented in Fig. 4. For a mass fraction around $`\mathrm{5\hspace{0.17em}10}^5`$ ($`0.04`$ solar) the overall strength of the iron forest is well reproduced, but the wavelengths of the single transitions still do not match the observation. For this reason we can only set an upper limit of $`X_{\mathrm{Fe}}<\mathrm{1\hspace{0.17em}10}^4`$ ($`0.07`$ solar) to the iron abundance. Compared to the standard LMC metallicity found in the nebula the central star shows a significant iron deficiency. In addition, the absence of N iv line emission at 1721 Å, which is relatively common among galactic \[WC\] stars, allows for the determination of an upper limit for the nitrogen surface abundance. Calculations with a full nitrogen model atom imply a value of $`X_\mathrm{N}<\mathrm{5\hspace{0.17em}10}^5`$.

For the subsequent nebular analysis, it is important to discuss the reliability of the derived ionizing fluxes. The main difficulty in the present analysis arises from the nebular continuum dominating the observed flux distribution from $`2500\text{Å}`$ to the visible wavelength range. At first sight, the situation seems to be well established, because the subtraction of the nebular continuum works without difficulties. In addition, the resulting stellar flux distribution is very well reproduced by our models, and the extinction parameter determined from the nebular spectrum is perfectly compensating the interstellar 2200 Å feature (see Fig. 5). However, slight uncertainties of the flux level in the optical wavelength range remain. For this reason, the determination of the stellar luminosity primarily relies on the UV, i.e. it becomes extremely sensitive to the applied reddening parameter $`E_{\mathrm{B}\mathrm{V}}`$.

Moreover, already small changes of the model parameters may significantly affect the derived values for the stellar luminosity and ionizing fluxes, because the observed flux distribution is still far away from the flux maximum at $`350`$ Å. Test calculations indicate that a variation of the model parameters, leading towards an inferior but still acceptable fit quality, can change the derived luminosity up to 0.1–0.2 dex. Including the observational uncertainties mentioned above, the possible error may reach values up to 0.3 dex.

In Fig. 5, the overall flux distribution as obtained from our model calculations is shown together with the observed flux distribution, and the ionization edges of O ii, He i and H i. The importance of testing the ionizing fluxes for a reliable determination of the overall flux distribution of the central star is obvious. A detailed analysis of the nebular spectrum is therefore capable to put substantial constraints on our atmosphere calculations.

## 5 Nebular modelling technique

### 5.1 The photoionization code

The photoionization models for the nebula were constructed using the code PHOTO which computes the radiative transfer of the ionizing photons in spherical symmetry with the “outward only” assumption for the diffuse radiation field and solves at each radius the equations of ionization equilibrium and thermal balance. The atomic data are the same as used by Stasińska & Leitherer (1996). The only modification is that low-temperature dielectronic recombination coefficients for S and Ar ions were introduced (assuming that they are equal to the computed values for the corresponding second-row elements, O and Ne respectively). The effects of dust are treated as described in Stasińska & Szczerba (1999, 2001).

### 5.2 The strategy

The aim is to reproduce all the observational constraints within a certain degree of tolerance (to be clarified later). The observational constraints are the emission line fluxes in the optical and the ultraviolet, the fluxes in the 12$`\mu `$m, 25$`\mu `$m and 60$`\mu `$m IRAS bands and the surface brightness distribution as measured by HST. Out of the zone affected by the central star, the surface brightness distribution is dominated by the \[O iii\] $`\lambda `$4959, 5007, H$`\beta `$ and H$`\alpha `$ lines (with the \[O iii\] lines giving the major contribution except, perhaps, in the outer parts of the nebula).

The models are computed using as an input the radiation field from the stellar atmospheres described in the previous section. The other input parameters are the inner nebular radius, the gas density distribution, the abundances of the elements and the characteristics of the dust grain population.

We assume that dust is composed of graphite grains with a classical power-law size distribution of exponent -3.5 between 0.003$`\mu `$m and 0.3$`\mu `$m. The total dust-to-gas mass ratio is chosen so as to reproduce the observed flux in the 25$`\mu `$m IRAS band (which includes nebular fine structure lines). As a matter of fact, in SMP 61, the observed IRAS fluxes correspond to a small amount of dust, which then does not affect significantly the nebular H$`\beta `$ luminosity or its ionization structure. We do not introduce an additional population of small organic grains, as in Stasińska & Szczerba (2001), since the effects of those grains are essentially to heat the gas and, as will be clear later, our initial models rather show an excess of heating.

We first determine the density distribution that reproduces the observed surface brightness distribution. By trial and error, we find that the observed surface brightness profile is well reproduced by a spherical model with a hydrogen density distribution described by

$$n=n_o\mathrm{exp}\left[\left(\frac{r\mathrm{5\hspace{0.17em}10}^{16}}{\mathrm{1.84\hspace{0.17em}10}^{17}}\right)^2\right],$$
(1)
where $`r`$ is the distance to the star in cm and $`n_o`$ is the density at the inner boundary.

The chemical composition with which we start is the one derived by classical empirical methods in Sect. 2.2. (the abundances of Mg, Si, Fe with respect to H are fixed at an arbitrary value of $`10^6`$ and do not affect our results). In a first step, we aim at reproducing the main characteristics of the nebula, which can be summarized by the following quantities: F(H$`\beta `$) the total nebular flux in H$`\beta `$, the nebular outer angular radius, $`\theta `$, defined as the radius of the nebular image at zero intensity, and several line ratios: \[S ii\] $`\lambda 6731/6717`$ to test the gas density, \[O iii\] $`\lambda 4363/5007`$ to test the gas temperature in the O<sup>++</sup> zone, \[O ii\] $`\lambda 3726,3729`$$`/`$\[O ii\] $`\lambda 7320,7330`$, to test the density and temperature in the O<sup>+</sup> zone; \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ to test the ionization structure; \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ to test the heating power of the central star, and \[O i\] $`\lambda 6300`$$`/`$H$`\beta `$ to test the importance of a warm neutral zone. It is only when we are reasonably satisfied with these constraints (which we will refer to as the restricted set) that we modify the elemental abundances, if needed, to obtain an acceptable model.

In order to easily visualize whether a model is acceptable, we compute for each constraint $`P`$ the quantity

$$\sigma _\mathrm{P}=(\mathrm{log}P_{\mathrm{mod}}\mathrm{log}P_{\mathrm{obs}})/t,$$
(2)
where $`P_{\mathrm{mod}}`$ is the value returned by the model, $`P_{\mathrm{obs}}`$ is the observed value, and $`t`$ the accepted tolerance in dex for this constraint. We determine the values of $`t`$ for each constraint from the estimated uncertainties in line fluxes, including the uncertainty in the dereddening process. Note that we use the outer nebular radius $`\theta `$ as a convenient way to visualize the constraint on the nebular size (but in fact we require the entire brightness profile to be reproduced by the model). Rather than using line fluxes with respect to H$`\beta `$, as is done usually, we prefer to use line ratios that are easier to interpret. Apart from the traditional temperature and density diagnostics mentioned above, we also use \[S ii\] $`\lambda 4069,4076`$$`/`$\[S ii\] $`\lambda 6717,6731`$ and O iii\] $`\lambda 1661`$$`/`$\[O iii\] $`\lambda 5007`$, which increase with nebular temperature and density, and C iii\] $`\lambda 1909`$$`/`$C ii$`\lambda 4267`$, which increases with nebular temperature. The line ratios \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$, C iii\] $`\lambda 1909`$$`/`$C ii\] $`\lambda 2326`$, N iii\] $`\lambda `$1750$`/`$\[N ii\] $`\lambda 6584`$, \[S iii\] $`\lambda `$6312$`/`$\[S ii\] $`\lambda 6717,6731`$, and \[Ar iv\] $`\lambda `$4740$`/`$\[Ar iii\] $`\lambda `$7136 reflect the ionization structure of the nebula. Finally, \[N ii\] $`\lambda 6584`$$`/`$\[O ii\] $`\lambda 3727`$, \[S ii\]$`\lambda `$6725$`/`$\[O ii\] $`\lambda 3727`$, C iii\] $`\lambda 1909`$$`/`$\[O iii\] $`\lambda 5007`$, \[Ne iii\] $`\lambda 3869`$$`/`$\[O iii\] $`\lambda 5007`$ and \[Ar iii\] $`\lambda `$7136$`/`$\[O iii\] $`\lambda 5007`$ are closely linked to the abundances of the elements with respect to oxygen while the IRAS fluxes reflect the abundance and temperature of dust. The entire list of constraints used in this study, together with the associated tolerances is given in Table 4. Note that we do not consider redundant lines, such as \[N ii\] $`\lambda `$6548 whose ratio with \[N ii\] $`\lambda 6584`$ is independent of the physical conditions in the nebula. The results of our models are presented in graphical form in the upcoming figures, where the values of $`\sigma _\mathrm{P}`$ are plotted for all the constraints $`P`$. The quality of a model is then easily judged by the location of the points with respect to the line of ordinate 0. Note that we consider a model as fully satisfying only if *each* of the values of $`\sigma _\mathrm{P}`$ is roughly found between $`1`$ and $`+1`$. *Any* strong deviation requires an explanation.

## 6 Results from the nebular model fitting

### 6.1 First models

We use the stellar atmosphere model described above as an input to the photoionization code. From the point of view of the nebular diagnostics available, there is no difference between using the model atmosphere with iron mass fraction $`X_{\mathrm{Fe}}=\mathrm{4\hspace{0.17em}10}^4`$ or $`X_{\mathrm{Fe}}=\mathrm{5\hspace{0.17em}10}^5`$. The hydrogen density distribution of the nebular model is given by Eq. (1). As will be made clear later on, one of the main shortcomings of the models is that the excitation, as measured by \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ is too high. We therefore discuss only models tailored to provide the lowest possible \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ ratio. For example, there is no point to present models in which the nebula would have an inner radius smaller than $`\mathrm{5\hspace{0.17em}10}^{16}`$ cm (which corresponds to 0.065″at the distance of the LMC). Also, while planetary nebulae are not necessarily ionization bounded, optically thin models are not discussed, since they lead to higher \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ ratios than ionization-bounded models. We tried various values of the inner density $`n_0`$. Note that the \[S ii\] $`\lambda 6731/6717`$ ratio is not very sensitive to $`n_0`$, due to the fact that the hydrogen density decreases outwards, and that, in addition, the electron density in the \[S ii\]$`\lambda `$6725 emitting zone is smaller than the hydrogen density. We show in Fig. 6 a series of models with $`n_0`$ equal to $`\mathrm{1.0\hspace{0.17em}10}^4`$ cm<sup>-3</sup> (circles), $`\mathrm{1.05\hspace{0.17em}10}^4`$ cm<sup>-3</sup> (squares), $`\mathrm{1.1\hspace{0.17em}10}^4`$ cm<sup>-3</sup> (triangles), $`\mathrm{1.15\hspace{0.17em}10}^4`$ cm<sup>-3</sup> (diamonds). All these models have a dust-to-gas mass ratio of $`\mathrm{3\hspace{0.17em}10}^4`$. We see that a value of $`n_0`$ between about $`1.02`$ and $`\mathrm{1.12\hspace{0.17em}10}^4`$ cm<sup>-3</sup> leads to an angular radius (and a surface brightness profile) compatible with the observations. But they give too large $`F(\mathrm{H}\beta )`$ (by a factor about 3 compared to the observed value after correction for extinction). We cannot help by adding more grains to absorb part of the Lyman continuum radiation, since the modelled flux in the 25$`\mu `$m IRAS band is already larger than observed. Thus, we conclude that, either the real number of stellar photons in the Lyman continuum is smaller than in the model atmosphere or that there is some leakage of ionizing photons. As seen in Fig. 6, the photoionization models also predict too large \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ and \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ (with a $`\sigma _\mathrm{P}`$ of 3–4), and slightly too large \[O iii\] $`\lambda 4363/5007`$. This cannot be improved by changing the oxygen abundance or by fine tuning the C/O and S/O ratios to fit the observed C iii\] $`\lambda 1909`$$`/`$\[O iii\] $`\lambda 5007`$ and \[S ii\]$`\lambda `$6725$`/`$\[O ii\] $`\lambda 3727`$ ratios. The models also have too low \[O ii\] $`\lambda 7320,7330`$$`/`$\[O ii\] $`\lambda 3726,3729`$ and \[S ii\] $`\lambda 4069,4076`$$`/`$\[S ii\] $`\lambda 6717,6731`$ ratios, the agreement with the observations being better for higher values of $`n_0`$, because the nebula is then more efficient in absorbing photons and the ionization front occurs at a higher value of the density.

In order to get away with the $`F(\mathrm{H}\beta )`$ problem, one could argue that the covering factor of the nebula is smaller than unity, so that some ionizing photons escape without interacting with the gas. But this would not help in solving the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$, \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ and \[O iii\] $`\lambda 4363/5007`$ problems.

We do not attempt two-sector models of different densities (such as those of Clegg et al. 1987 or Dudziak et al. 2000) for the following reasons. First, as mentioned in Sect. 2, the nebular image is extremely round. This, a priori does not mean that the nebula is spherical. Indeed, it could be elongated but seen pole on. But, in that case, the nebular extension along the line of sight would be larger than the one observed on the images. Replacing a portion of the model nebula with two sectors of lower density (to obtain a larger dimension) would considerably worsen the \[O ii\] $`\lambda 7320,7330`$$`/`$\[O ii\] $`\lambda 3726,3729`$ and \[S ii\] $`\lambda 4069,4076`$$`/`$\[S ii\] $`\lambda 6717,6731`$ ratios while not helping much for the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ ratio and worsening the \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ ratio, as we can judge from spherical models built at lower densities. Besides, if the low density sector is density bounded, this may indeed help with the H$`\beta `$ flux, but will worsen the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ ratio.

One way to lower the excitation of the nebula is to assume that it is clumpy. Models that differ in density and filling factor but have the same value of the product $`n_0^2ϵ`$, where $`ϵ`$ is the volume filling factor, will have about the same external radius but will have lower excitation for higher densities. Such a series of models is shown in Fig. 7, where the values of the couple ($`n_0`$, $`ϵ`$) are ($`\mathrm{1.05\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 1.0) (circles), ($`\mathrm{1.5\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.49) (squares), ($`\mathrm{2.0\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.27) (triangles), ($`\mathrm{3.0\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.122) (diamonds). As seen in the figure, models with higher densities considerably improve many of the line ratios, and especially the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ and \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ ratios as well as the \[O ii\] $`\lambda 7320,7330`$$`/`$\[O ii\] $`\lambda 3726,3729`$ and \[S ii\] $`\lambda 4069,4076`$$`/`$\[S ii\] $`\lambda 6717,6731`$ ratios and also the \[O i\] $`\lambda 6300`$$`/`$H$`\beta `$ ratio. It is then easy to fine tune the N, S and C abundances. A model with $`n_0`$ between $`\mathrm{2.0\hspace{0.17em}10}^4`$ and $`\mathrm{3.0\hspace{0.17em}10}^4`$ cm<sup>-3</sup> would then be in agreement with all the observational constraints (except C iii\] $`\lambda 1909`$$`/`$C ii$`\lambda 4267`$, but see later) if one assumed a covering factor of roughly 1/3.

The predicted emission in the IRAS bands deserves some comments. First, we note in Fig. 7 that models with lower excitation show lower fluxes in the IRAS bands. The reason is that the larger proportion of neutral hydrogen particles in the nebula reduces the amount of Lyman continuum photons absorbed by dust. Another point to make is that, in the hypothesis of incomplete coverage of the ionizing source by the nebula (and if dust is intimately mixed with the gas), the predicted IRAS fluxes should be larger than observed by roughly the same amount as $`F(\mathrm{H}\beta )`$ is larger than observed in our models computed for a covering factor of unity. Obviously, this is not the case of our predicted flux in the 12$`\mu `$m band (the observed flux in the 60$`\mu `$m band is only an upper limit so there is no problem there). Increasing the dust abundance only worsens the problem, since then the dust particles are more numerous to share the photons that heat them, and the temperature of the dust grains becomes smaller, decreasing the flux in the 12$`\mu `$m band with respect to that in the 25$`\mu `$m band. Actually, there may be a simple explanation to our under-prediction of the flux in the 12$`\mu `$m band. Our code does not handle PAHs, while PAHs are known to exist in planetary nebulae with Wolf-Rayet central stars (Szczerba et al. 2001, Cohen 2001) and they contribute to the emission in the 12$`\mu `$m band. Note also that our predictions concerning dust emission are given here only for the ionized part of the nebula. The presence of dust in the neutral zone would increase the fluxes in the $`25\mu `$m and $`60\mu `$m bands. We do not consider this in the modelling presented here, since it brings no useful information in the problem we are interested in. In the following, we therefore simply require from our models that they do not over-predict the 12$`\mu `$m and 25$`\mu `$m fluxes with respect to the observations.

We have found that with a model having an inner density of about $`2\mathrm{3\hspace{0.17em}10}^4`$ cm<sup>-3</sup> and a filling factor about $`0.30.1`$, we are able to reproduce many of the observational constraints, provided that we assume a covering factor of about 1/3. However, this model is still not entirely satisfactory. Mainly, the \[O iii\] temperature is still too high and so is the predicted \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ ratio. And a conspicuous shortcoming of the model is that the C ii$`\lambda 4267`$$`/`$C iii\] $`\lambda 1909`$ ratio is predicted too large by a factor of about 5. Note that, while the density structure of planetary nebulae is known generally not to be uniform on a small scale from high resolution images, the clumps or filaments of gas are actually embedded in a more diffuse medium. This diffuse medium is expected to be more highly ionized, and thus its presence will enhance the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ predicted by our simple model. We do not attempt to construct a two-phase model but we keep this in mind.

In the following, we look for more satisfactory solutions by releasing some of the assumptions of our models.

### 6.2 Models with modified stellar fluxes

Due to the compact and spherical appearance of the nebula, the covering factor of 1/3 which is needed to explain $`F(\mathrm{H}\beta )`$ seems to be very unrealistic. The obvious next step is therefore to assume that the discrepancy is at least partly due to an overestimation of the stellar ionizing radiation field in Sect. 4. As explained in Sect. 4.2, the error for the derived stellar luminosity may be as high as a factor of 2 under very pessimistic assumptions.

Let us therefore first consider the case where we assume that the number of photons emitted by the star at any energy larger than 13.6 eV is only 50% of what is predicted by the model atmosphere (from the point of view of the photoionization model, it is equivalent to assume that the star luminosity is lower than given in Sect. 4 by a factor 2). The results are shown in Fig. 8. Here, the values of the couple $`n_0`$, $`ϵ`$ are $`\mathrm{7.5\hspace{0.17em}10}^3`$ cm<sup>-3</sup>, 1.0 (circles), $`\mathrm{1.0\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.6 (squares), $`\mathrm{1.2\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.4 (triangles), $`\mathrm{1.4\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, 0.3 (diamonds) while the dust-to-gas mass ratio is still assumed to be $`\mathrm{3\hspace{0.17em}10}^4`$. Since the number of ionizing photons is smaller than in the previous case, we had to take smaller densities to match the observed size of the ionized nebula. Globally, the ionization parameter is nevertheless smaller, so that \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ can now be fitted even with a model with a filling factor of unity. Of course, the \[O ii\] $`\lambda 7320,7330`$$`/`$\[O ii\] $`\lambda 3726,3729`$ and \[S ii\] $`\lambda 4069,4076`$$`/`$\[S ii\] $`\lambda 6717,6731`$ ratios are better reproduced for the highest density models of this series, i.e. $`n_0`$ around $`\mathrm{1.5\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, which again require a small filling factor. The best models of this series are at least as good as the best models shown in Fig. 7 without requiring a small covering factor. Note that the \[O iii\] $`\lambda 4363/5007`$ ratio is still somewhat high in this series of models, and adjusting the S, C and Ne abundances (by lowering them) can only make the problem worse.

Another option to consider is a radiation field softened in the Lyman continuum. This would not only lower the \[O iii\] $`\lambda 5007`$$`/`$\[O ii\] $`\lambda 3727`$ ratio, but also reduce the heating. Fig. 9 shows a series of models identical to those of Fig. 8, but this time keeping 60% of the photons between 13.6 and 24.6 eV and 40% of the photons with energies larger than 24.6 eV. Clearly, the models from this series represented by the triangles are quite satisfactory also as regards the electron temperature and the \[O iii\] $`\lambda 5007`$$`/`$H$`\beta `$ ratio. \[S ii\]$`\lambda `$6725$`/`$\[O ii\] $`\lambda 3727`$ can easily be matched by reducing the original sulfur abundance by a factor $`1.32`$). The observed fluxes in the IRAS 25 $`\mu `$m and 60 $`\mu `$m bands are compatible with dust-to-gas mass ratios of up to $`\mathrm{3\hspace{0.17em}10}^3`$. The quality of the fit of the optical and UV data is not affected by the amount of dust in the allowed range. The only important mismatch of the model with the observations is the C iii\] $`\lambda 1909`$$`/`$C ii$`\lambda 4267`$ whose computed value is much higher than observed, by a factor of almost 4. This problem is discussed in the following subsection.

### 6.3 Models with inhomogeneous nebular composition

It is a known fact that recombination lines of C, and also of other elements like N, O, Ne generally indicate much higher abundances of the parent ions than collisionally excited lines (see e.g. Liu 2002 or Esteban 2002 and references therein). Many explanations have been proposed (see e.g. Stasińska 2002 for a review). The generally preferred one is to assume that the chemical composition of the nebulae is not uniform (e.g. Torres-Peimbert, Peimbert & Peña 1990, Liu et al. 2000). And indeed, there is  *direct* evidence at least in two nebulae (A 30 and A 78) of the presence of clumps of matter enriched in C (Jacoby & Ford 1983, Harrington & Feibelman 1984) and also N, O and Ne for A30 (Wesson, Liu & Barlow 2003). The central star of SMP 61 being a Wolf-Rayet, this nebula is a good candidate for carbon enrichment. If material seen in the stellar wind ends up in the nebula as cool clumps, these clumps will strongly emit in the C ii$`\lambda `$4267 line. We have mimicked such a situation by assuming that the inner zone of the nebula consists of helium- and carbon-rich material, maintaining the same density distribution as before. The number of combinations to explore is obviously very high, and we have run many simulations. In Fig. 10, we present the result of a series of simulations with $`n_0`$=$`\mathrm{1.4\hspace{0.17em}10}^4`$ cm<sup>-3</sup>, $`ϵ`$=0.3 (the best combination from Fig. 9) with the following abundances of He and C (relative to H by number) in the inner 3 10<sup>-3</sup>M of the nebula : 1.04 and 0.4 (circles), 0.78 and 0.3 (squares), 0.52 and 0.2 (triangles), 0.26 and 0.1 (diamonds). These ratios correspond to the mass ratio of He and C in the atmosphere of the central star, as derived in Sect. 4. Otherwise, the number abundance ratios with respect to hydrogen are: H : He : C : N : O : Ne : S : Ar = 1 : .09 : $`\mathrm{4.0\hspace{0.17em}10}^4`$ : $`\mathrm{3.95\hspace{0.17em}10}^5`$ : $`\mathrm{2.65\hspace{0.17em}10}^4`$ : $`\mathrm{4.26\hspace{0.17em}10}^5`$ : $`\mathrm{4.0\hspace{0.17em}10}^6`$ : $`\mathrm{8\hspace{0.17em}10}^7`$. We see that the last combination provides a very good fit to the data, including the C ii\] $`\lambda `$1909$`/`$C ii$`\lambda `$4267 line ratio. Probably, other combinations can be found that provide a similar fit.

## 7 Discussion

### 7.1 Stellar ionizing fluxes

The nebular analysis of SMP 61 indicates two possible problems concerning the theoretical flux distribution of the central star: The absolute number of Lyman photons is probably too large by a factor of 2, and the Lyman continuum radiation field is possibly too hard. The first problem may be discussed away by assuming a very pessimistic value for the uncertainty of the derived luminosity. However, the second point concerns the slope of the energy distribution in the flux maximum and cannot be easily resolved.

A simultaneous solution of both problems would be achieved by decreasing the effective temperature $`T_{\mathrm{eff}}`$ of the model atmosphere (i.e. increasing the radius where $`\tau _{\mathrm{Ross}}=2/3`$). In this way, photons would be distributed from the flux maximum to the observed wavelength range. In effect, the derived luminosity would be decreased, and a softer flux distribution would be obtained in the flux maximum. For a corresponding decrease of $`T_{\mathrm{eff}}`$, the stellar core temperature must be lowered significantly, or the mass-loss rate must be increased. Of course, both operations are strongly limited by the necessity to fit the observed spectrum of the central star. Test calculations show that the simultaneous fitting of the observed spectrum and the nebular ionizing fluxes is not possible with the model atmospheres applied in present work.

A possible solution concerns the treatment of dielectronic recombination in our stellar atmosphere code. In the present version dielectronic transitions are treated as optically thin, and their contribution to the rate equations is accounted for by the approach of Mihalas (1971). However, when C iv recombines to C iii in the outer part of the WC atmosphere, the ionization edge of C iii becomes optically thick. In this case also the corresponding dielectronic transitions become optically thick, and their recombination efficiency is reduced because recombination photons cannot escape. Indeed, first tests with an improved treatment show that the recombination from C iv to C iii is suppressed by this effect, leading to higher derived mass-loss rates and slightly lower stellar temperatures. A detailed investigation of this topic will be presented in a forthcoming paper.

### 7.2 The effect of He-rich clumps on the nebular spectrum

He-rich clumps, which we have invoked in Sect. 6.3, can in fact soften the radiation field available to the rest of the nebula under certain conditions, by selectively absorbing the photons above 24.6 eV. The total mass of helium $`M_{\mathrm{He}}^\mathrm{a}`$ required to absorb all the He<sup>0</sup>-ionizing photons emitted by the central star is of the order $`10^3`$/$`n_\mathrm{e}`$ M (where $`n_\mathrm{e}`$ is the electron density in cm<sup>-3</sup>). If the integrated mass of helium in the clumps is at least equal to this value and if the clumps are located close to the inner boundary of the nebula with an integrated covering factor of unity, then the He<sup>0</sup>-ionizing photons would be completely blocked by the He-rich material. Such an extreme situation does not correspond to the case of SMP 61, since \[O iii\] $`\lambda 5007`$ is emitted not only in the most central part but in the entire nebula. One can however imagine a less extreme situation where the high density He- and C-rich clumps are distributed over the entire volume of the nebula and soften the average radiation field available to the bulk of the nebular material. In such a case, a softening of the predicted stellar energy distribution in the Lyman continuum would perhaps not be necessary. Note however that such clumps are not efficient in absorbing radiation between 13.6 and 24.6 eV, unless they are very dusty, which for SMP 61 seems excluded by the observed IRAS fluxes. Obviously, the impact of these He- and C- rich clumps on the global spectrum of the nebula will strongly depend on their density, on the amount of mixing with nebular material and on their spatial distribution in the nebula.

In the case of SMP 61, we can estimate the maximum mass of helium contained in the clumps by considering that the He/C proportion is the same as in the stellar wind and by using the observed flux in the C ii$`\lambda 4267`$ line. The observed flux in this line corresponds to a carbon mass $`M_\mathrm{C}`$ of about $`10`$/$`n_\mathrm{e}`$ M, if most of the carbon is in C<sup>++</sup> form and if the line is emitted at a temperature of 8000 K, which is the case for the C-rich zone in our models. Since helium and carbon are roughly equal by mass in the stellar wind, this means that the total helium mass is much smaller than to $`M_{\mathrm{He}}^\mathrm{a}`$, implying that there is no significant blocking of He<sup>0</sup>-ionizing photons by the clumps. If the clumps are of much higher density than in our models, C<sup>++</sup> may be partly recombined in the clumps, increasing $`M_\mathrm{C}`$ with respect to $`M_{\mathrm{He}}^\mathrm{a}`$, but in that case the clumps will have a smaller cross section and will be less efficient to block the stellar radiation. Therefore, in the specific case of SMP 61, we do not think that the He-rich material will significantly soften the stellar radiation available to the nebula.

### 7.3 Stellar parameters and evolutionary status

The derived stellar parameters are in line with recent analyses of galactic \[WC\]-type central stars based on line-blanketed models (De Marco & Crowther 1998, 1999; De Marco et al. 2001; Crowther et al. 2003). In contrast to previous un-blanketed calculations (Koesterke & Hamann 1997a, b), the new models show a trend towards similar surface mass fractions for early and late \[WC\] subtypes, with $`X_\mathrm{O}10\%`$ and $`X_\mathrm{C}/X_{\mathrm{He}}1`$.

Due to the known distance to the LMC, SMP 61 offers the rare chance to determine the luminosity of a \[WC\]-type central star. From the analysis of the stellar spectrum alone we infer a value of $`10^{3.9}`$L. The nebular analysis reveals that this value is probably too large by up to 0.3 dex – dependent on the assumed nebular covering factor. The luminosity of SMP 61 is therefore in the range of $`10^{3.6}`$$`10^{3.9}`$L. For a second object, BD+30 3696, Li et al. (2002) derive a distance of 1.2 kpc from the angular expansion of the nebula. For this distance, Crowther et al. (2003) obtain a value of $`L_{}=10^{3.6}L_{}`$. Both values are in the range that is expected for the majority of all post-AGB stars, with a typical mass around 0.6 M (see Blöcker 1995; Herwig 2001b). Altogether, no evidence is found that \[WC\] central stars have different masses than H-rich objects although, admittedly, the number of cases with known $`L_{}`$ is small.

An interesting hint on the evolutionary status of SMP 61 is given by its low iron and nitrogen abundances. For a standard solar composition as given by Gray (1992, p. 318) or Grevesse & Sauval (1998) a ratio of $`\mathrm{log}(\mathrm{Fe}/\mathrm{O})=1.33`$ is expected. More recent investigations (Allende Prieto et al. 2001a, b) imply a lower solar oxygen abundance and give $`\mathrm{log}(\mathrm{Fe}/\mathrm{O})=1.19`$. The upper limit of $`X_{\mathrm{Fe}}<\mathrm{1\hspace{0.17em}10}^4`$ that is derived in Sect. 4 translates to $`\mathrm{log}(\mathrm{Fe}/\mathrm{H})<5.54`$ for a solar composition. In relation to the oxygen abundance of the nebula ($`\mathrm{log}(\mathrm{O}/\mathrm{H})=3.58`$) we obtain $`\mathrm{log}(\mathrm{Fe}/\mathrm{O})<1.96`$, i.e. iron is under-abundant by at least 0.63 dex in the atmosphere of the central star. In addition, no nitrogen is detected in our analysis. The upper limit of $`X_\mathrm{N}<\mathrm{5\hspace{0.17em}10}^5`$ lies 0.9 dex below the abundance derived for the nebula.

The detection of an iron deficiency is in line with recent results from analyses of PG 1159 stars, the probable descendants of \[WC\]-type central stars. Miksa et al. (2002) set upper limits from 0 to -1.5 dex solar for the iron abundance of 15 PG 1159 stars. Additionally, Werner et al. (2003) find an iron depletion of at least -1.5 dex for the central star of A78, a \[WC\]-PG 1159 transition object. Also Crowther et al. (2003) find evidence for an iron abundance of 0.3–0.7 dex below solar for BD+30 3696 and NGC 40. However, for these two objects the nebulae show oxygen abundances significantly below solar (see Peña et al. 2001), so that the iron deficiency may be attributed to a low initial metallicity.

Such an iron deficiency is expected for material which has been exposed to s-process nucleosynthesis in the He-intershell of thermally pulsating AGB or post-AGB stars (see Lugaro et al. 2003; Herwig et al. 2003). This region is expected to consist of partially He-burned material, i.e. mainly He, C, O and no N. The observed surface composition of SMP 61 therefore exactly resembles the abundance pattern that is expected for a He-intershell where the s-process has been active.

For the formation of central stars with \[WC\] surface composition, the H-rich layers above the He-intershell must be removed or mixed with intershell material. Herwig (2001a) and Blöcker (2001) demonstrated that the latter is possible by dredge-up of intershell material after an AGB-final thermal pulse (AFTP) or a late- or very late thermal pulse on the post-AGB (LTP, VLTP). In case of a LTP or VLTP very long evolutionary timescales on the post-AGB are expected, because the central star first evolves to the blue as H-burner, then back to the red during the thermal pulse, and again to the blue as He-burning \[WC\] star.

The observed absence of nitrogen for SMP 61 strongly favors the AFTP scenario because for a LTP or VLTP nitrogen is produced by the CNO-cycle in the H-burning phase. In the AFTP scenario the last thermal pulse occurs when the star is leaving the AGB, and the central star directly enters the post-AGB as a He-burner. However, to obtain a sufficiently high probability for the AFTP to occur at small enough H-envelope masses, a coupling of mass-loss to the thermal pulse cycle is required (Blöcker 2001). From the upper limit of $`X_\mathrm{N}<\mathrm{5\hspace{0.17em}10}^5`$ we can derive an upper limit for the mass of the remaining H-envelope after the final thermal pulse.

Under the assumption that nitrogen has nebular abundance in the H-envelope, and is completely destroyed in the He-intershell, it follows that the maximum mass of the H-envelope can be only 0.15 of the dredged-up mass from the intershell. Herwig (2001a) finds that the dredged-up mass is typically of the order of $`\mathrm{6\hspace{0.17em}10}^3`$M. Consequently, the remaining H-envelope mass at the final thermal pulse must have been smaller than $`\mathrm{9\hspace{0.17em}10}^4`$M. This value is significantly smaller than the corresponding envelope masses in the AFTP models of Herwig (2001a). For these models the mass-loss on the AGB has been tuned in such a way that envelope masses of $`\mathrm{3\hspace{0.17em}10}^2`$ and $`\mathrm{4\hspace{0.17em}10}^3`$M were obtained after the thermal pulse. The present work implies even lower values, and therefore also lower H-abundances for the resulting \[WC\]-type central stars.

Additional evidence against the LTP and VLTP scenarios is due to the very young appearance of the planetary nebula. Using the observed nebular expansion velocity of $`29.3\text{km}\text{s}^1`$ (Vassiliadis et al. 1998) and angular radius from Table 4 we find an expansion time of about 3000 yr, whereas evolutionary tracks of Blöcker (1995, 2001) for a 0.625 M central star undergoing a VLTP imply a timescale of at least 8000 yr to reach a temperature of $`85\mathrm{kK}`$ (dependent on the time where the zero point for the PN age is assumed). Moreover, Górny & Tylenda (2000) find general evidence against the LTP and VLTP scenarios due to the similarity of the \[WR\]PN and non-\[WR\]PN populations.

### 7.4 Stellar mass-loss history

We know from the analysis of the central star atmosphere that the presently observed mass-loss rate is $`\dot{M}`$ = $`10^{6.12}M_{}\text{yr}^1`$ and that the carbon mass fraction is 0.52 (see Table 3). As mentioned in Sect. 6.3, the total mass of carbon in the nebula required to produce the observed flux in the C ii$`\lambda 4267`$ line, $`M_\mathrm{C}`$, is about $`10`$/$`n_\mathrm{e}`$ M. In our model, the value of $`n_\mathrm{e}`$ in the zone where this line is mostly emitted is about $`\mathrm{2\hspace{0.17em}10}^4`$ cm<sup>-3</sup> which implies $`M_\mathrm{C}`$ $`\mathrm{5\hspace{0.17em}10}^4M_{}`$. Note that $`M_\mathrm{C}`$ is just slightly larger than the total mass of carbon in the C-rich zone, and that, if the density of the carbon-rich clumps is larger, the total carbon mass will be smaller. We thus infer that the star has spent about $`1300`$ yr in a similar state of mass-loss.

This timescale is significantly lower than the time expected for a PN central star to reach $`T_{\mathrm{eff}}=85\mathrm{kK}`$. Evolutionary tracks for He-burning stars of the corresponding luminosity from Vassiliadis & Wood (1994) give values around 3000 yr. Vassiliadis et al. (1998) show that this timescale fits the nebular age of SMP 61. H-burning tracks from Blöcker (1995) imply an age of 3000 yr for a post-AGB mass of 0.625 M (with $`L_{}=10^{3.9}L_{}`$) and 5000 yr for 0.605 M (with $`L_{}=10^{3.7}L_{}`$).

Our modelling then suggests that the mass-loss of SMP 61 may be intermittent, which would indicate that it is initiated by other processes than radiative acceleration. Interestingly, our atmosphere models would confirm this hypothesis, because the force due to radiation pressure is much too low to explain the observed mass loss: For the massive WC star WR 111, which has a very similar spectral appearance to SMP 61, the radiation pressure, as calculated in our models, supplies about one half of the energy necessary to drive the stellar wind (Gräfener et al. 2002). On the other hand, for SMP 61 the same models provide only 17% of the wind energy. However, one must be aware that the timescales for mass-loss and evolution derived above both suffer from large uncertainties. The spectroscopically derived mass-loss rate depends on the clumping factor ($`\dot{M}1/\sqrt{D}`$ for constant $`R_\mathrm{t}`$, see Sect. 3.2) which is only roughly known. The error margin for $`\dot{M}`$ may therefore be as large as $`\pm 50\%`$. The stellar evolutionary timescale is even more uncertain, because it depends on the stellar luminosity and the mass of the stellar envelope above the burning shell. The latter is strongly dependent on model assumptions. The nebular age is also uncertain, mostly because the nebular velocity changes during the course of evolution as shown by dynamical simulations (Mellema 1994, Villaver et al. 2002).

## 8 Summary

We have obtained HST STIS imaging and spectroscopic observations of the Magellanic Cloud planetary nebula SMP 61 and its central star. We have also searched the literature for other observational data concerning this object, in order to obtain a full set of observational constrains and to assess their uncertainties.

We have performed a detailed spectral analysis of the central star. For that, we used the Potsdam code for expanding atmospheres in non-LTE, which includes the statistical equations for many ions, including ions from the iron group. The observed stellar spectrum is well represented by a model with the following parameters: $`L_{}`$ = $`10^{3.96}L_{}`$, $`R_{}`$ = $`0.42R_{}`$, $`T_{}`$ = $`87.5\mathrm{kK}`$, $`\dot{M}`$ = $`10^{6.12}M_{}\text{yr}^1`$, $`v_{\mathrm{}}`$ = $`1400\text{km}\text{s}^1`$, clumping factor $`D`$ = $`4`$ and $`R_\mathrm{t}`$ = $`4.63R_{}`$. The abundances of the elements by mass are $`X_{\mathrm{He}}`$ = $`0.45`$, $`X_\mathrm{C}`$ = $`0.52`$, $`X_\mathrm{N}<\mathrm{5\hspace{0.17em}10}^5`$, $`X_\mathrm{O}`$ = $`0.03`$, $`X_{\mathrm{Fe}}`$ $`<\mathrm{1\hspace{0.17em}10}^4`$. The low iron and nitrogen abundances indicate that the star has undergone an AGB final thermal pulse with a very low H-envelope mass ($`<\mathrm{9\hspace{0.17em}10}^4M_{}`$). The iron depletion is probably due to s-process nucleosynthesis during previous thermal pulses.

The fluxes from the model stellar atmosphere were used as an input to our photoionization code to construct photoionization models of the nebula. We considered all the available observational constraints, within their error bars. These constraints are numerous and allow a detailed discussion of the models. We find that the fluxes from the stellar atmospheres allow to reproduce many features of the nebular emission. However, the observed nebular properties are better explained if the central star produces less photons in the Lyman continuum than the model atmosphere by about a factor 2. An even better fit is obtained if the Lyman continuum radiation is also softer than predicted by the stellar model. A possible solution involves the treatment of dielectronic recombination in the stellar atmosphere code and will be presented in a forthcoming paper.

We also find that the nebula must contain extremely carbon rich clumps in order to reproduce the observed C iii\] $`\lambda 1909`$$`/`$C ii$`\lambda 4267`$ ratio. The observational constraints do not allow us to describe these clumps in detail (density, chemical composition, location) but these clumps are likely produced by ejecta from the clumpy stellar wind, possibly partly mixed with the nebular material. Obviously, high signal-to noise and high spatial resolution observations of nearby planetary nebulae excited by a similar star would be extremely useful to better understand the process of enrichment of planetary nebulae by freshly made carbon (and helium).

From a comparison of the estimated total amount of carbon in the nebula with the carbon mass-loss rate measured in the central star, we infer that the strong WC mass-loss may have been active for only a limited period in the post-AGB evolution of SMP 61. However, this last result relies on uncertain timescales and is subject to revision.

###### Acknowledgements.

This work has been partially supported by DGAPA/UNAM (grant 114601), CONACYT/México (grant 32594-E), by the CNRS-CONACYT joint project (number 10385) and by the Deutsche Agentur für Raumfahrtangelegenheiten under grant DARA 50 OR 0008.


format: text

text:
# 1 Introduction

## 1 Introduction

Recent advances in high-energy astrophysics involve observations of extremely complex phenomena such as jets from active galactic nuclei (AGN), gamma-ray bursts (GRB), and ultrahigh-energy cosmic rays (UHECR). Observations of AGN jets, consisting of a highly collimated stream of material, show that the outflow expands at relativistic velocity and spans a distance scale of thousands of light-years. The collimation and production mechanism, most likely involving the dynamics of accretion disks around a black hole in the center of the AGN, are subjects of current research. Gamma-ray bursts, on the other hand, are some of the brightest observed light sources in the universe. The amount of electromagnetic energy output in a burst is equivalent to several times the solar mass released in a matter of seconds. The out-flowing materials of a GRB expand at relativistic velocity as well, and are possibly collimated, similar to an AGN jet. The nature of the progenitor and explanations for their observed characteristics are currently under debate. Much of our current understanding of these objects are inferred from the properties of the observed radiation.

A strong magnetic field is believed to exist in both the GRB and the AGN jets. The interaction of the relativisticly expanding material with the environment can lead to nonlinear plasma phenomena that result in the acceleration of particles to high energies. Ultrahigh-energy cosmic rays, with energies observed up to around $`10^{20}`$ eV, are believed to come from extra-galactic sources. The nature and origin of these cosmic rays as well as their acceleration mechanism are still a mystery.

The study of these extreme phenomena requires tremendous effort. So far, progress in our understanding has required a combination of observation, numerical simulations, and theoretical modeling. However, astrophysical observations must be carefully checked for instrumentation effects. And the complex numerical and theoretical calculations used to interpret these observations must be validated. Thus, it is important to calibrate the techniques used in the observations and to benchmark computer model calculations. Furthermore, since observational astrophysics deals with uncontrolled environments, laboratory experiments able to model the relevant extreme conditions would provide unique insight into the underlying physical mechanisms.

Laboratory studies, ranging from work on atomic spectroscopy, and the studies of hydrodynamics, radiation flow, and the equation-of-state using intense lasers, have been instrumental in astrophysics research. Recently, it has been suggested that accelerators can be used in the laboratory investigation of extreme astrophysical phenomena. In this paper we discuss possible experiments using intense particle and photon beams to verify astrophysical observations and to study relativistic plasma dynamics and ultrahigh-energy cosmic acceleration mechanisms. An overview of the accelerator facility at the Stanford Linear Accelerator Center (SLAC) is given in Section 2. In Section 3, we discuss calibration experiments, focusing on the discrepancy in the UHECR spectrum measured by two large-aperture cosmic ray experiments, and describe an experiment that may help resolve it. In Section 4, we discuss laboratory experiments that may improve our understanding of the underlying dynamics of high-energy astrophysics phenomena. We conclude with an outlook in Section 5.

## 2 An Overview of the SLAC Facility

The 3-km long linear accelerator is the backbone for SLAC’s high-energy physics research program. It is capable of delivering electrons and positrons with 50 GeV energy and 120 Hz repetition rate at $`10^{10}`$ particles per pulse. Currently it serves as the injector for the PEP-II storage ring to produce copious amount of B-meson particles for CP-violation measurements. It can also deliver beams to the fix-target experimental hall End-Station A (ESA) and the Final Focus Test Beam (FFTB). A schematic layout of the SLAC facility is shown in Figure 1. High intensity photon beams, tunable from X-ray to gamma-ray, can be derived using a variety of methods, such as undulators, laser-Compton back-scattering, and bremsstrahlung. Depending on the required wavelength, typical fluences of $`10^9`$ photons per pulse can be provided.

## 3 Calibration Experiments

### 3.1 Detector Calibration

High energy beams from the linac can be used to generate a variety of secondary and tertiary beams for calibration purposes. A secondary pion and positron beam with well-defined momentum can be generated using a combination of target and selection magnet system, with a beam intensity that can be set to below 1 particle per pulse. With the addition of a tagger magnet, this secondary beam can be converted into a photon beam with known energy up to 20 GeV.

This test beam setup in the ESA has been used for the GLAST satellite mission, whose objective is to study energetic astrophysical gamma rays with energies in the 20 MeV to TeV range. The GLAST detector package consists of sophisticated silicon tracker and CsI calorimeter. It is important to calibrate its response and understand the various analysis algorithms in a controlled test beam environment before its space launch. Details on this experiment can be found elsewhere.

### 3.2 X-ray Spectroscopy

Recent X-ray observations of AGN galaxies have revealed features in the iron emission lines that are characteristic of Doppler shifts and gravitational redshifts expected from accretion disk models. The emission lines can be thought of as “clocks” moving in various circular orbits around the black hole. To further probe the spacetime structure in the accretion disk, high resolution imaging and broad-band spectroscopy, such as those planned for the Constellation-X and MAXIM missions, are needed. A detailed laboratory measurement of heavy element atomic transitions and associated polarization effects will also be required for a proper interpretation of the observational data.

For this purpose, an intense X-ray source, such as those available at a synchrotron light-source facility would be valuable. The next generation of linac-based light-source, with peak brilliance in the range of $`10^{32}`$-$`10^{33}`$ photons/sec/mrad2/mm2/0.1% bandwidth at 1-10 keV, could also play a role in this study.

### 3.3 Air Fluorescence Efficiency Measurement

The study of ultrahigh-energy cosmic rays has been based on observations of the secondary shower particles resulting from interactions in the atmosphere. For cosmic ray energies above $``$1014 eV, the shower particles can reach ground level and extend over a large area. One observation technique uses an array of sparsely spaced ground detectors to measure the density of these shower particles, which is related to the energy of the primary cosmic ray. The Akeno Giant Air Shower Array (AGASA) near Tokyo, Japan, for example, covers an area of approximately 100 km2, with 100 detector units separated by about 1 km from each other.

The cosmic ray shower also generates a trail of fluorescent light. The fluorescence is emitted nearly isotropically, mostly by the nitrogen molecules in air excited by shower secondaries. Instead of studying the transverse profile of the shower, as in the ground array approach, fluorescence-based detectors use a system of mirrors and photomultipliers to image the shower’s longitudinal development. The fluorescence luminosity is related to the primary’s energy; and the shape of the longitudinal profile provides information on the primary’s composition. This technique is used by the Fly’s Eye detector and its upgraded version, the High Resolution Fly’s Eye (HiRes).

Studies of UHECR events showed that they are not related to any known galactic sources. If they originated in extra-galactic sources, interactions with the cosmic microwave background radiation would result in the attenuation of their energy. The flux above 1019 eV is expected to drop significantly due to the production of pions – the so-called GZK cutoff. However, the Fly’s Eye/HiRes and AGASA experiments have observed events greater than 1020 eV, well above the GZK cutoff. The two experiments have now accumulated similar exposure at the highest energies. With increased statistics, differences between the two measurements have become apparent. In particular, the flux measured by HiRes is systematically lower than that reported by AGASA above $`4\times 10^{18}`$ eV; there is also a difference in the energy at which the observed power-law spectrum changes slope, the so-called “ankle” structure. This can be due to tails in the energy resolution function or other systematic errors, and is currently being investigated by both experiments.

One possible contribution to this discrepancy is the air fluorescence yield. Current understanding of air fluorescence, based on previous measurements, is incomplete. Many issues still remain: the detailed shape of the fluorescence spectrum, the pressure and atmospheric impurities dependences, and the dependence of fluorescence yield on shower particle energy. The associated systematic uncertainty is estimated to be 15%. A more precise measurement is desired as improvements are being made to other systematics in the observation.

Recently, it has been suggested that the high intensity electron beams at SLAC can be used for such study. At the relevant energies, air showers produced by a cosmic ray hadron is a superposition of electromagnetic sub-showers. At the shower’s maximum, it consists of mostly electrons with energies dissipated to the 100 MeV level, near the critical energy of air. Further shower development is dominated by energy loss through ionization and excitation rather than shower particle creation. SLAC’s electron beams interacting in an air-equivalent alumina target produce similar secondary electron energy distributions – see Figure 2. The SLAC E-165 experiment – Fluorescence in Air from Showers (FLASH) – has been proposed to study in detail the fluorescence yield in an air shower. It aims to make precision measurements of the total air fluorescence yield, as well as the spectral, pressure, composition, and energy dependencies. Details on this experiment have been presented elsewhere at this Workshop.

Other examples of accelerator-based experiments that support astrophysical investigations are measurement of the Landau-Pomeranchuk-Migdal (LPM) effect, which has implication for photon/hadron identification at high energies, and observation of the Askarian effect, which can be used to detect UHE neutrinos. These experiments have been carried out using SLAC beams.

## 4 Relativistic Plasma Experiments

While important issues remain to be resolved in the observational results of super-GZK events, the existence of extra-galactic UHECR above $`10^{18}`$ eV is well established. The nature of these cosmic rays and their acceleration mechanism are still a mystery, and various models have been proposed as solutions. In the so-called “top-down” approach, the decay products of massive particles produced in the early universe could account for the observed UHECRs, especially those above the GZK cutoff. Certain “grand-unification” theories predict the existence of particles with mass around $`10^{16}`$ GeV. Particles more massive than this, if they were to explain super-GZK events, would have to be produced continuously since their lifetimes would be extremely short. In some theories these particles can be emitted from topological defects created between causally disconnected regions during early epochs of cosmological phase transitions.

In the “bottom-up” approach, conventional particles accelerated in powerful astrophysical systems are thought to be responsible for the observed UHECR spectrum. The acceleration mechanisms are complex, involving strong magnetic fields and nonlinear plasma effects. Diffusive shock acceleration has been the generally accepted model. More recent ideas include unipolar induction acceleration and high gradient plasma acceleration in wakefields created by Alfvén shocks. Possible acceleration sites are AGNs and gamma-ray bursts where typically relativistic plasma outflows are present. The key observational feature of UHECR is the power-law spectrum. The appropriate spectral index is predicted by existing models. Our goal is to experimentally test some of these models in the laboratory.

Typical beams delivered for experimentation in the FFTB are in short pulses pico-seconds long, $`10\mu `$m in radius, and consists of $`10^{10}`$ particles. Thus, the pulse power is approximately 40 Petawatts, and the intensity is $``$1020 W/cm2. The energy density in the bunch is on the order of $`10^{13}`$ J/m3. For comparison, the threshold for high-energy-density conditions, the energy density in a hydrogen molecule or the bulk moduli of solid-state materials, is $`10^{11}`$ J/m3 The strong nonlinear and collective responses of a bunched relativistic particle beam to external stimuli are some of the important characteristics of a high-energy-density system relevant for astrophysical studies. Here we discuss possible relativistic plasma experiments. In particular, we explore the possibility of merging electron and positron beams to form a kinetically relativistic plasma, allowing the laboratory investigation of cosmic high-energy acceleration and radiation production phenomena.

### 4.1 e+e- Beams as Relativistic Plasma

Neutral co-moving e+e- beams have been investigated in an effort to improve the luminosity limit in high energy e+e- storage ring colliders. The disruptive effect of one beam’s electromagnetic fields on the other can be compensated, in principle, by colliding neutral beams. This idea had been studied using two pairs of 0.8 GeV beams. The experiment demonstrated beam-charge compensation with improved luminosity.

For our purpose, the 1-GeV electron and positron beams emerging from the damping rings at the beginning of the SLAC linac (see Figure 1) could be combined, forming an e+e- plasma streaming at relativistic velocity. The transverse positions of the two beams would be aligned to micron precision using high resolution beam position monitors. The temporal locations would be synchronized using the damping rings’ RF phase control, which is stable at the sub-picosecond level. This level of precision beam control has been demonstrated in measurements of wake-fields in accelerator structures. The concept is illustrated in Figure 3.

For a relativistic bunched beam, temperature can be defined in terms of its emittance. Analogous to entropy, emittance is a measure of disorder. The discussions here follow those in Lawson. The beam’s transverse temperature is given by

$$kT_{}=\frac{\beta ^2\gamma mc^2ϵ^2}{4\sigma _r^2}$$
(1)
where $`ϵ`$ is the beam’s emittance, $`\gamma `$ the Lorentz factor, and $`\sigma _r`$ the transverse beam size. The longitudinal temperature due to energy spread is negligible for relativistic beams.

The other plasma parameters can now be calculated. Results are shown in Table 1, using typical SLAC parameters at the exit of the damping rings. Plasma parameters are given in the frame co-moving with the beams. As can be seen, the number of particles inside a “Debye sphere” ($`N_D`$) is much greater than one, so that the effects of individual particles on each other are negligible compared to the collective effects, and the plasma description is indeed appropriate.

For typical AGN jet parameters, the plasma length scales are much smaller than the jet dimensions. Thus, the AGN jet plasma is usually treated as having infinite extend. For typical relativistic bunched beams, however, the Debye radius ($`\lambda _D`$) is smaller than the bunch length but larger than the transverse beam size. As a consequence, the perpendicular plasma waves (involving particle motion in the transverse direction) have different properties compared to those excited in an infinite plasma. However, properties of the parallel propagating waves remain the same as those in an infinite plasma. The laboratory e+e- plasma discussed here can thus model the parallel propagating waves in an infinite plasma. As discussed below, this mode is most relevant for AGN jet dynamics.

So far our discussion have concentrated on neutral plasmas. The composition of astrophysical jets is, however, far from being understood. Magnetic confinement is generally accepted as the collimation mechanism, but it is also highly unstable. Models of current-carrying jets provide a possible alternative mechanism where the self-magnetic fields create a pinching force. This is very similar to the plasma-lens effect familiar to the beam-plasma physics community. Non-neutral plasma instabilities relevant for AGN jets could be studied using charged beams readily available at a facility such as SLAC. Possible experiments are under study.

### 4.2 Scaling Laws and Relevance to Astrophysics

The challenge for laboratory astrophysics is to create a terrestrial setting which can be scaled to the astrophysical environment. Magnetohydrodynamic (MHD) models have been used to describe many astrophysical processes such as bow-shock excitation in AGN jets or supernova explosions. The MHD equations have the property that they are invariant under the appropriate scale transformations. This has been the basis, for example, for designing laser experiments to simulate supernova remnants.

The MHD models are applicable when certain assumptions are satisfied. These, however, may not be applicable to the astrophysical conditions of interest here. In the following, we discuss a more general approach based on kinetic plasma theory. In particular, we concentrate on astrophysical plasma processes that might be investigated using high-energy-density particle beams.

The observed non-thermal radiation spectrum from AGNs is the subject of many recent studies. In some models, broad-band Blazar emission has been attributed to synchrotron radiation and/or various forms of Compton processes. While in other models, it is described by the production of photon-pairs from the decay of mesons produced via the interaction of energetic protons with ambient photon and/or matter. These models successfully describe various features in the observed spectrum, and thus are useful for understanding the radiation processes. But such phenomenological approach does not describe details of the underlying micro-physical dynamics of AGN jets. In particular, it does not address the issue of how the relativistic jet gives rise to energetic electrons and/or protons which subsequently produce the radiation. For example, these models typically assume that diffusive shock acceleration produces the required power-law spectrum.

In the plasma physics approach, details of the underlying dynamics for transferring kinetic energy in the relativistic jet into radiation are described. In the model proposed by Schlickeiser et al., the jet is described by a one-dimensional outflow consisting of electron and positron pairs with bulk relativistic velocity, directed parallel to a uniform background magnetic field. The pairs have non-relativistic temperature in the co-moving frame. The e+e- jet propagates into an interstellar medium consisting of cold protons and electrons.

This two-stream multi-fluid system is studied in the jet rest frame. The analysis starts with a general phase space distribution, and the calculations then give the dispersion relations of the parallel propagating electrostatic (longitudinally polarized) and low-frequency transverse (Alfven-type) plasma waves. These waves are excited via a two-stream instability in the pair plasma. For typical AGN parameters, the calculations show that the jet kinetic energy is transferred via plasma turbulence to the initially cold interstellar protons and electrons, which then reach a plateau distribution in momenta. The resulting radiation spectrum is consistent with observation.

These kinetic plasma calculations also show that the instability build-up times and growth rates scale with the densities and the bulk relativistic factor, while the damping rates scales also with temperature. For example, the time it takes to build up the transverse instability in the protons is given by $`t_{t,p}(1/\omega _{p,e})(n_j/n_i)(m_p\mathrm{\Gamma }/m_e)^{4/3}`$, where $`\omega _{p,e}`$ is the electron plasma frequency, $`n_j`$ and $`n_i`$ are the jet and interstellar plasma densities, and $`\mathrm{\Gamma }`$ is the bulk Lorentz factor. The Landau damping rate is found to scale with $`\mathrm{\Theta }^{3/2}\omega _{p,e}\mathrm{\Gamma }^2\mathrm{exp}[(\mathrm{\Gamma }1)/\mathrm{\Theta }]`$, where $`\mathrm{\Theta }=kT/m_ec^2`$ is the dimensionless temperature parameter.

### 4.3 Parameters for Laboratory Experiments

To determine whether the parameters of the relativistic plasma created by merging electron and positron bunches are relevant for an experimental investigation of AGN dynamics, the various dynamical time scales are calculated. The results are shown in Figure 4 for the parameters given in Table 1. In the setup being considered here, the pair plasma in the co-moving frame appears to be $``$1-m long to the ambient plasma traveling through it. As can be seen from Figure 4, all dynamical time scales are shorter than the plasma traversal time: the time during which the relativistic plasma and the ambient plasma interact with each other. Typical plasma time scales are shown as the inverse plasma frequency. The build-up of the electrostatic waves is rather quick, for both the electrons and the protons, even with a fairly thin ambient plasma. The build-up of the transverse waves takes much more time, particularly for the protons, in which case an ambient plasma density of 1015 cm-3 in the laboratory is required.

Also, the maximum growth rate of the electrostatic turbulence is much greater than the Landau damping rate; similarly, the transverse turbulence growth rate is much larger than the cyclotron damping rate. Thus, this set of experimental parameters is in a regime where strong nonlinear plasma turbulence similar to those excited in AGN jets can be created and studied in detail experimentally in the laboratory. Further theoretical calculations are needed to guide the design of the experiment. A detailed numerical simulation using particle-in-cell techniques is needed as the next step.

The transverse magneto-hydrodynamic (Alfven-type) wave is especially interesting for testing various cosmic acceleration mechanisms. This type of turbulence is crucial in the formation of collisionless shocks and for efficient particle deflection in the diffusive shock acceleration process. The Alfven wave is also expected to excite plasma wakefields, which can provide high gradient particle acceleration. The spectrum and polarization properties of the radiation produced in the interaction of this e+e- “jet” with an ambient plasma can be measured and compared with astrophysical observations. Detailed simulation studies for these experiments are underway.

The laboratory experiments described here could have applications beyond the understanding of AGN jet dynamics. The dynamics in the polar caps of a spinning neutron star have been studied in the context of relativistic streaming electron-positron plasma. Also, if GRB radiation is beamed, its dynamics would be similar to those found in AGN jets. Thus, our laboratory experiments would also shed light on these systems.

## 5 Summary and outlook

The field of laboratory astrophysics holds promise to the understanding of some of the most exciting astrophysical observations today. We have shown that particle accelerators are excellent tools for laboratory astrophysics, providing calibration data for observations and bench-marking computer models, as well as creating extreme conditions that make possible investigation of astrophysical dynamics in a terrestrial laboratory. SLAC, with the existing expertise and infrastructure, is well-positioned to contribute to this rapidly growing field. The proposed ORION facility for advanced accelerator research and beam physics will also be able to support dedicated laboratory astrophysics experiments with its unique combination of high quality electron beams and diagnostic lasers.

## Acknowledgments

I would like to thank the Workshop organizers for their kind hospitality. I would also like to thank P. Chen, R. Noble, K. Reil, and P. Sokolsky for many fruitful discussions.