EDP Sciences logo
Subscriber Authentication Point
Search
Menu
Home All issues Volume 620 (December 2018) A&A, 620 (2018) A84 Full HTML
Free Access
Issue
A&A
Volume 620, December 2018
Article Number A84
Number of page(s) 12
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/201833981
Published online 30 November 2018

© ESO 2018

1 Introduction

Although the interacting wind model for the formation of planetary nebulae (PNe) has been accepted for a long time (Kwok et al. 1978; Kwok 1982, 2002) and morphological sequences were drawn in detail (Balick 1987), the exact mechanism is not yet clearly understood up to now. However, significant progress have been made in the past years. Several numerical simulations starting with Okorokov et al. (1985) followed by Frank & Mellema (1994); Mellema & Frank (1995), and by the Schönberner group (e.g. Perinotto et al. 1998, 2004; Steffen et al. 1998; Corradi et al. 2000; Schönberner et al. 2014), have been carried out. The latter performed a one-dimensional hydrodynamic simulation based on the evolutionary tracks and the simulations of the circumstellar shell around asymptotic giant branch (AGB) stars. The results of their simulation showed that a shock-bounded ionized main PN shell, moving supersonically into the AGB-material, compresses the inner parts of the matter into a dense but thin shell. The unaffected AGB-material becomes ionized as well and forms a rapidly expanding halo. Because of the drop in luminosity of the central star of the PN (CSPN) toward the white dwarf regime, the outer part of the PN shell recombines quickly forming a second inner halo. The disadvantage of this simulation is that it just considers one dimension. Thus, more complex interactions of the winds such as Rayleigh–Taylor instabilities, which quickly build a clumpy environment as we see in HST images, are not taken into account. For this reason, the filling factor (ϵ) has usually been set to unity. Moreover, comparing these models with observations of older evolved nebulae suffers from interaction to the environment and the interstellar matter (ISM). Thus, the study of nebulae with higher galactic latitudes is of special interest. However, there might be an influence on other properties changing the wind physics, mainly due to the lower metallicity. Photoionizing radiative transfer codes such as the steadily evolving Cloudy1 (Ferland et al. 1998, 2013, 2017; Ferland 2000) are often therefore used to characterize such nebulae and to identify their membership to the stellar populations (Emprechtinger et al. 2004, 2005; Öttl et al. 2014).

The PN IC 5148 (PN G002.7-52.4) is a nebula which is to date not well investigated in detail. First listed in the Second Index Catalogue of Nebulae and Clusters of Stars (Dreyer 1910) with two independent entries as numbers IC 5148 and IC 5150 discovered by Swift (1899) and Gale (1897) independently, it was finally discovered to be the same object by Hoffmeister (1961). Morphologically it is declared as a round nebula in all catalogs and Chu et al. (1987) classified it as a multiple shell planetary nebula (MSPN) due to a small step in the H α image with a radius ratio of only 1:1.2 between the two structures. In addition, the intensity decrease was found smaller than that for typical MSPNe. A few years later the authors searched with a larger field systematically around many nebulae for extended emission features without detecting the very low surface brightness halo we investigate here (Hua et al. 1998). While earlier spectroscopic studies took only a small fraction of one or two spectral lines to obtain radial velocity and expansion of the nebula (e.g. Meatheringham et al. 1988), to our knowledge up to now the only, more extended spectroscopic analysis were performed by Kaler et al. (1990) and by Kingsburgh & Barlow (1994) within two surveys of 75 and 80 southern PNe, respectively. Both studies used pre-CCD electronic spectral scanning devices, taking a very small aperture region of the nebula and detected only a hand full of lines. They end up coinciding in the result that the nebula has about galactic disk abundance or only slight underabundance despite its large galactic latitude (b approximately −52°). Further, thesurvey spectra used in the PN catalog of Acker et al. (1992) gives intensities of only four lines. The inspection of the original data file of this survey, provided now at the Hong Kong/AAO/Strasbourg PN data base (thereafter HASH2; Parker et al. 2016; Bojičić et al. 2017), does not recover more usable lines above the noise level.

The CSPN of IC 5148 was first classified by Acker et al. (1982), and studied in detail by using the International Ultraviolet Explorer (IUE) by Kaler & Feibelman (1985). They revealed that it is a very evolved hot central star. The IUE color index was not sufficient to derive an exact value of the temperature, but they give values based on their three methods from approximately 100–220 kK. As later stated by Patriarchi & Perinotto (1991b,a), using the same data, no strong wind features were found. Moreover, a later investigation shows that this is common for CSPN of nebulae which have passed the knee in the Hertzsprung–Russell diagram (HRD) toward the beginning white dwarf cooling track. The only published optical study of the CSPN is the classification of Mendez (1991), giving it the type of a hgO(H) star (a hydrogen rich high gravity star with very broad Balmer absorptions).

Previous distance estimations obtained from so-called statistical distance scales span from 0.53 kpc using the CSPN evolution (Phillips 2005), 1.06 kpc using the ionized mass versus surface brightness relation (Cahn et al. 1992; Stanghellini et al. 2008), 1.1 kpc using the IRAS dust temperature (Tajitsu & Tamura 1998), 1.37 kpc by using the Hα surface brightness-radius relation (Frew et al. 2016), 1.8 kpc using a non-linear calibration of the radio brightness vs. diameter (using Buckley & Schneider 1995, calculated by us), 2.1 kpc using the radio temperature vs. radius relation (using van de Steene & Zijlstra 1994, calculated by us) and even 3.59 kpc using the mass-radius distance (Zhang 1995).

In this paper we present the recent discovery of extended emission features around the hitherto known nebula IC 5148 and a detailed medium-high resolution spectroscopic study of its central star with the ESO VLT X-shooter instrument. We also performed a detailed spectroscopic study with long-slit spectroscopy from ESO FORS2 in various directions of the nebula and the newly discovered halo. The latter is, to our knowledge, one of the first and deepest yet found of its type in a halo. Finally, using the now available distance estimation from the second data release DR2 of Gaia (Gaia Collaboration 2018) and the derived parameters from our CSPN spectroscopy, we modeled the main nebula by using the photoionizing radiative transfer code Cloudy (Ferland et al. 2017).

thumbnail Fig. 1

Co-added direct images in [O III] (top panel) and H α (bottom panel). The slit positions of the spectroscopy are indicated in the top panel. While in the inner part, a linear intensity scale covering the whole range from zero to peak flux of the nebula was used, for regions with flux below 5% a steep asinh between 0.01 and 1% of the peak of the linear part was used to enhance the by orders of magnitude fainter structures.

2 Data

The narrow band imaging was derived at the CDK20 Planewave telescope of ITelescope.net in Siding springs (31°16′24′′ south, 149°03′52′′ east) during several nights in July 2013. An Apogee U16MD9 camera equipped with a Kodak 16803 CCD without focal reducer and without binning was used. A total exposure time of 11.0 h with the [O III] filter by co-adding 22 frames and 4.5 h with the H α filter by co-adding nine exposures were obtained. The resulting images, shown in Fig. 1, were used to identify the features of interest and the position angles for the follow up investigations. To obtain a correct line ratio [O III]/H α image, the resulting frames were differentially calibrated by means of the long-slit spectroscopy derived at ESO. The astrometry was determined by 40 surrounding stars of the Gaia DR2 catalog (Gaia Collaboration 2018) resulting in a plate scale of 0.′′ 5447 per pixel. Stars near the nebula on the co-added images have a full width at half maximum (FWHM) of 1.′′ 82 and 1.′′ 96 in H α and [O III], respectively.

Long-slit spectra were obtained with the FORS2 spectrograph (Appenzeller et al. 1998) mounted on the Cassegrain focus of ESO VLT UT1 (Antu) in the nights October 5, 2016 from 3:15 to 4:03 UT; October 6, 2016 from 0:15 to 0:30 UT and, October 10, 2016 from 1:17 to 1:45 UT in service mode as filler program with minor sky quality. In total, seven spectra were collected, of which, three have a position angle PA = 30° and four have a position angle PA = 150° (from N over E) on the sky through the center of the PN. The slit positions were selected to cover the features in the wide halo of the PN. We used the long-slit mode of FORS2 since this enabled us to cover the entire halo using the full slit length of 6.′8 with the standard collimator (slit width = 0.′′7). All spectra were taken with 14 min. exposure time. We used GRISM 300V and the GG435 order separation filter, covering a wavelength range from 455 to 889 nm. The MIT/LL CCD mosaic and the standard focal reducer collimator result in a 0.′′ 2518 pixel−1 spatial scale. This setup leads to a final wavelength dispersion of 0.33 nm pixel−1. The night sky lines were measured with a resolution R = λ∕Δλ of 200 and R = 360 at the blue and the red ends of the spectrum, respectively. We derived a FWHM of the stellar sources along the slit of about 1.′′ 17 at the blue end and 0.′′95 at the red part of the spectrum. This corresponds well to the reported DIMM seeing of the ESO meteo monitor of 1.′′ 1 at 500 nm. The data were reduced incorporating the standard calibration mode using the ESO FORS2 pipeline v5.3.113 (Freudling et al. 2013). The resulting flux calibration was compared with the expected continuum flux of the CSPN, finding a scatter of less than 5% between the individual observations obtained at clear, but non-photometric sky conditions. Two spectra at position angle 30° were taken at slightly minor sky conditions with some cirrus clouds on October 5, 2016. These were scaled with a constant wavelength independent factor using the central star observations. We used the software package molecfit (Smette et al. 2015; Kausch et al. 2015) to create a model of the telluric absorption lines by applying it on the high signal-to-noise ratio (S/N) central star spectrum. As the PN and its extended halo covered the whole slit length, we used the software skycorr (Noll et al. 2014) to correct for the sky emission lines. This approach also enabled us to derive a secondary correction of the wavelength calibration. We averaged the spectra for each sky direction to achieve the final spectra of the faint sources used for this study. Finally, a correction by photometric values of 20% was applied (see Sect. 3).

The medium-high resolution echelle spectra were taken with the X-shooter instrument (Vernet et al. 2011) mounted at the ESO VLT UT2 (Kueyen) in the night of October 5, 2016 from 0:23 to 0:57 UT at slightly degraded weather conditions as filler program in service mode. This instrument covers the entire wavelength range from 0.3 to 2.5 μm in three spectral arms (UVB arm from the ultraviolet to the B band with an E2V CCD44-82; VIS arm in the visual regime with a MIT/LL CCD; NIR arm in the near-infrared regime with an Hawaii2RG chip), at medium resolution simultaneously. While the NIR and UVB frames had an exposure time of 125 s each, the VIS arm exposures were 213 s each (all 0.8′′ slit). A total of 29 frames were obtained and co-added during data reduction with the standard ESO X-shooter pipeline v2.84 (Modigliani et al. 2010). The spectra were normalized to continuum to derive line profiles for the spectral analysis. Thus, the slightly degraded weather conditions did not harm. We did not apply a telluric absorption correction since all lines of interest were in undisturbed regions. The S/N was about 45 in the UVB and 25 in the VIS arm. That was lower than expected from exposure time estimators and the DSS magnitude of the central star. It could be caused by a combination of cirrus clouds and the 1.6′′ seeing at that time. The S/N of the NIR part of the spectrum was insufficient to derive any parameters and so was not used further.

3 Analysis

3.1 The central star

The X-shooter spectrum was normalized to unity flux and searched for absorption lines. The only one appearing in the VIS arm spectrum was the H α line. Nevertheless, it is heavily contaminated by the emission of the Hydrogen line by the nebula and by the two strong [N II] lines. As shown below in the following section, these lines are asymmetric and doubled due to nebular expansion. Moreover, they vary significantly along the slit. Thus, a reliable removal was not possible and the H α line was not used for the further analysis. The NIR arm has a too-low S/N value (<10 in the J band and <5 in the H and K bands) to derive a useful spectrum. Therefore, the analysis of the CSPN parameters is focused only on the UVB arm of the X-shooter instrument.

We overlaid in our analysis the high resolution non-LTE spectra for central stars of planetary nebulae of the Tübingen stellar atmosphere (thereafter TMAP, Rauch 2003) as offered at Rauch & Reindl (2014) to the identified lines in the UVB arm spectra (Fig. 2). The clear appearance of the He ii 468.6 nm line in absorption, not seen by Mendez (1991), and the He ii 420.0+454.1 nm lines without major absorption, places the CSPN in their classification as O(H) and not as hgO(H). The number of lines and the S/N does not allow for a sophisticated determination of the log (g) value. Thus, the evolution along a standard model by Miller Bertolami (2016) was assumed. As Napiwotzki (2001) showed, the result is not sensitive on the selection of log(g) for the optical flux analysis. A likelihood analysis of all lines, weighted by equivalent width of the line to consider the weakness of the lines like He ii 454.1 nm, 420.0 nm or H ϵ, was used to derive a (log T = 5.146) for the central star.

The calibrated FORS2 spectra of the CSPN were first compared individually and then combined to derive the final flux. The derived flux curve was folded with the filter curves of Bessell (1990) and multiplied with the zero point given in Bessell et al. (1998)5. This results in V = 16. m34. The DENIS IC = 16. m476 ± 0. m10 lead with a typical color of such very hot white dwarfs (V − IC) ≈−0. m35 ± 0. m05 from Landolt & Uomoto (2007) for a VDENIS = 16. m13 ± 0. m15, but the target is very much near the survey limit. Kaler & Feibelman (1985) derived magnitudes of V (UV) = 16. m16 ± 0. m17 using the IUE spectra between 130 and 185 nm and assuming a 130 000 K blackbody spectrum and zero extinction due to the high galactic latitude. The flux calibration allows us to directly estimate the interstellar foreground extinction to be very low. The Balmer decrement of the nebular lines in our FORS2 spectra (only H α and H β are covered by these spectra) give us an E(B − V) = 0. m031 ± 0.010. The central star continuum, using stellar models from the Tübingen Model Atmosphere Package TMAP (Rauch 2003) accessed via the German Astrophysical Virtual Observatory (GAVO) service TheoSSA (Rauch & Ringat 2011), confirm the almost extinction-free line of sight. The Galaxy Evolution Explorer (GALEX) data base give the similar value of E(B − V) = 0. m021 based on the NUV flux of mNUV = 13. m480 ± 0. m004. Using the GALEX FUV data in the small 7′′ aperture and applying that to the TMAP NLTE model together with the effective temperature derived by our spectroscopy with X-shooter and the interstellar extinction obtained from the Balmer lines, with the extrapolation as used in Kaler & Feibelman (1985) for IUE data, we achieve a visual magnitude of 16. m 35. The spectral energy distribution (SED) from 130 to 800 nm fits very well using a low extinction value of E(B − V) = 0. m03 (see Fig. 3). The overlaid TMAP spectra also show how marginally the SED differ for these stars even at IUE and GALEX wavelengths, and thus why Kaler & Feibelman (1985) were unable to derive a more accurate CSPN temperature. We further adopted this low value for the interstellar extinction. This value is in contradiction to the spectroscopy by Kaler et al. (1990) and Kingsburgh & Barlow (1994), who found, using several lines of the Balmer decrement, values as high as c(H β) of 0.4 and 0.38, respectively (corresponding to E(B − V) = 0. m26). However, using their published line intensities from H α, H β, H γ and H δ individually give the large scatter of values of c from even slightly negative values of −0.03 up to 0.85. An extinction value as high as c = 0.4 does neither match with our spectra, nor with the IUE and GALEX photometry. Moreover, it is also unlikely with respect to the high galactic latitude of the target6. On the other hand, the H α/H β = 1.7 in Acker et al. (1992) is nonphysical. We thus take these values no further into account.

There are some issues with Gaia DR2 for the extraction of nebular contamination (see Kimeswenger & Barría 2018). However, since the central star of IC 5148 appears isolated from its main nebula, we expect for a mostly reliable photometry from Gaia. We thus compared isolated blue compact white dwarfs in the lists of Landolt photometric standards (Landolt 1992) and spectrophotometric standards (Landolt & Uomoto 2007) to derive color equations dedicated to this type of stars for the photometry into standard system (see Fig. 4). With this we derive a VGaia = 16. m14 ± 0. m03, using GGaia = 16. m097 and Gaia color (BP−RP) = −0.635. We finally adopted the latter value as reliable, and corrected the FORS2 fluxes by 20% to take into account the filler program minor sky conditions.

To test the result and to approve for the photometric flux calibration, the FORS2 spectrum was weighted with λ4 (giving for the Raleigh Jeans tail of an exact blackbody a flat spectrum) and compared with the TMAP model downgraded to the FORS2 resolution (see Fig. 5). The deviation of the continuum is in most parts on the order of a few percent. The four stellar lines in the region, namely, Hα, Hβ, He I 540.6 nm, and He II 468.5 nm (taking into account the slight error in the continuum calibration at the blue end) also fit perfectly in their depths. The latter together to the fact that the He I/He II ratio is very sensitive to the stellar temperature, independently confirms our X-shooter result on Teff.

thumbnail Fig. 2

Spectral lines from the medium-high resolution spectroscopy of the UVB Arm of X-shooter. Overlaid are the model grid spectra along the evolutionary track log (g) values of a CSPN according to Miller Bertolami (2016) from the high resolution stellar spectra grid from the German Astrophysics Virtual Observatory (GAVO) grid resources (Rauch & Reindl 2014).
thumbnail Fig. 3

UV photometry from IUE and GALEX with our calibrated FORS2 spectrum. The model spectral energy distributions from TMAP are overlaid, normalized to the FORS2 data.

3.2 Nebular expansion, systemic velocity and astrometry

The X-shooter spectra have only a 9′′ long slit and thus do not cover the nebula as a whole. But its high resolution allows the determination of the line splitting along the line of sight in the center of the nebula, thus giving the direct value for the nebular expansion velocity (see Fig. 6). Moreover, the midpoint of the two lines gives a good value for the systemic velocity. Low resolution spectroscopy often suffer from the fact that the two unresolved lobes are not identical in intensity. Thus, for example, in the case of unresolved lines the result wouldbe biased toward to red values. The measured expansions were independent of the element type within the error. This is in fact astonishing, as the high ionized species often give higher expansion velocities than the low ionized species (Weinberger 1989; Gesicki et al. 2014). We measured a peak to peak expansion velocity of 47.9 ± 1.5 km s−1. This is in very good agreement with the value of 53.4 km s−1 given by Meatheringham et al. (1988) if taking into account the much lower spectral resolution of 11.5 km s−1 used by the authors. The nebula thus has a fairly large expansion compared to other roundish nebula in the catalogs (e.g. Weinberger 1989, and citations therein). The center between the two well separated peaks was used as system velocity. The measured system blueshift was corrected for the local standard of rest (LSRK) using the IRAF7 procedure rvcorrect. This result in RVLSRK = −36.7 ± 1.7 km s−1. The systematic difference from the hitherto published value of −23 km s−1 (Meatheringham et al. 1988) most likely originates from the asymmetry of intensities of the redshifted and of the blueshifted part of the shell. This moves the center toward the brighter peak in the case of low resolution spectroscopy.

The astrometry was derived by 40 surrounding stars from the Gaia DR2 catalog (Gaia Collaboration 2018) which resulted in a rms of about 105 mas. We found only one old plate, which was scanned with high resolution from the ESO red survey plate scanned at MAMA and provided via CDS Aladin. It dates from September 5 1985. This was calibrated in the same homogeneous way with the same Gaia stars. The rms of this calibration was 80 mas. The resulting positions using our astrometric solution and that of Gaia (Lindegren et al. 2018) are given by

While the ESO R plate corresponds perfectly using the Gaia proper motion, the stacked CDK20 images have a offset slightly bigger than the error bars, most likely having its origin in a remaining distortion correction on the stacked images. But as it is still far below the FWHM caused by the seeing, we are still able to trust the division of images to derive line ratios. Using the Gaia DR2 distance of 1.3 kpc (see Sect. 3.3) we derive a proper motion velocity in the plane of sight of 62 ± 5 km s−1 toward southeast.

thumbnail Fig. 4

Colorequations of hot white dwarfs from the Landolt photometric standards (Landolt 1992) and from the spectrophotometric standards (Landolt & Uomoto 2007) versus Gaia DR2. The target position used for calibrations is given in the lower panel. While all other stars are hot nearby compact white dwarfs without PNe, the encircled target is the central star of the Helix nebula.
thumbnail Fig. 5

λ4 weighted FORS2 spectrum with the TMAP model. We intentionally did not remove the telluric O2 and H2 O absorptions at the red end of the spectrum. The blue horizontal lines show the required correction of the line depths at the blue end due to the deviation to the continuum.
thumbnail Fig. 6

Expansion of the PN and system velocity in the Galaxy from X-shooter spectra after correction for LSRK.
Gaia DR2: Epoch 2015.5, Equinox J2000.0
α = 21h59m35. s095392053 ± 0. s000000001
μα = −0.00652 arcsec yr−1
δ = −39°23′ 08.′′35833212 ± 0.′′00000002
μδ = −0.′′00779 arcsec yr−1
CDK20: Epoch 2013.52, Equinox J2000.0
α = 21h59m35. s1066 ± 0. s0080
δ = −39°23′ 08.′′439 ± 0.′′12
ESO R: Epoch 1985.68.5, Equinox J2000.0
α = 21h59m35. s1084 ± 0. s0053
δ = −39°23′ 08.′′714 ± 0.′′08

3.3 Distance, luminosity, and position in the Galaxy

The most commonly used statistical distance scales are that one of Cahn et al. (1992) and Stanghellini et al. (2008) giving a distance to IC 5148 of 1.06 kpc. Using the same radio data, but the radio brightness temperature method of van de Steene & Zijlstra (1994) we obtain a distance of 2.1 kpc. More recently, Frew et al. (2016) established a method based on the optical Hα surface brightness giving 1.37 kpc. Using a recent calibration with the Gaia DR1 data for the statistical methods mentioned above by Stanghellini et al. (2017), we then get a distance of 1.88 kpc. In a new investigation, containing a large sample of sources, Kimeswenger & Barría (2018) showed for some issues concerning Gaia DR2 (Gaia Collaboration 2018; Luri et al. 2018) parallaxes for PNe with high surface brightness around the central stars. However, in this case the blue Gaia color (BP–RP) suggests for no contamination, and so we believe we are able to use the value provided by Gaia. Thus, we get a Gaia distance of . For the further analysis a final distance of 1.3 kpc is applied.

The radial velocity, accurate proper motion and the Gaia distance allows us to derive the spatial movement vector. The target moves in galactic coordinates toward lower galactic longitude at 22 km s−1 and marginally increases its distance below the galactic plane with −7 km s−1 while it moves to larger galactocentric distance with +68 km s−1. Thus, whilethe peculiar fraction from normal galactic rotation and change in z is small, it is significantly moving outwards in the Galaxy.

Using our spectra and photometry of the central star we are also able to derive a spectroscopic distance (Napiwotzki 1999, 2001) by means of the dereddened apparent magnitude, the mass of the central star and, the surface gravity log (g). With a temperature derived from the X-shooter data we obtained a range of distances. We compared these to the recent evolutionary tracks (Fig. 7) by Miller Bertolami (2016). We derived with the evolutionary tracks, the spectroscopic model distances and with the expansion velocity, a dynamical age τdyn for each, according to (1)

where θ = 67.′′2 is the angular radius averaged over all directions, vexp = 47.9 km s−1 the expansion velocity, and Dphot is a distance as function of the effective temperature T (formula 1 in Napiwotzki 2001) using the final mass Mfin of the track and the photometry VGaia. This was then compared to the model age τmod. Even if we take into account an uncertainty factor of two in the dynamical age, the intersection within the given temperature range derived for the CSPN fits only for the following tracks: initial mass Mini = 1.25 M⊙, metallicity Z = 0.01 giving a distance of D = 1.8 kpc and, Mini = 1.5 M⊙, metallicity Z = 0.02 ending up in D = 1.3 kpc. Extrapolating the track for Mini = 1.25 M⊙, metallicity Z = 0.02 (which stopped due to a late helium flash in the calculation of the stellar model grid) we derive an intermediate position at D = 1.5 kpc. We obtain similar results within 20% using the tracks of Vassiliadis & Wood (1994) and assuming an initial to final mass relation (Vassiliadis & Wood 1993) to derive mass and log(g). The tracks of Blöcker (1995) do not cover such low mass stars. This intersection lowers slightly, but within the error by 0.8 σ the effective temperature of the star as derived above from the X-shooter spectra. Any other intersection of dynamical age with model age do not fit by more than a factor of three. Hence, these evolutionary tracks are in good agreement with the distances and temperatures previously derived from the nebular parameters and lead us to the input parameter set for the best fitting model with a initial mass of Mini = 1.5 M⊙ and Z = 0.02. The resulting parameters are listed in Table 2.

thumbnail Fig. 7

Upper panel: evolutionary tracks of Miller Bertolami (2016) as log(g) vs. log (T) for initial masses of 1.00, 1.25, 1.5 and 2.0 M⊙ (color coded) for metallicities of Z = 0.01 and Z = 0.02 (line thickness coded). The corresponding spectroscopic distances are indicated in light gray. The vertical dark gray lines indicate the range of temperatures as derived from our spectroscopy. Dotted lines mark extrapolation of the track, which ended earlier. Lower panel: computed model ages from the tracks above (solid lines) vs. the dynamical ages (dashed lines). The circles indicate the solutions by intersections of the model age with the dynamical age. The corresponding solutions found in the lower panel are marked in the upper panel as well. Crosses mark the estimated solution of the extrapolated track. The intersections of the 1.0 M⊙ (red) fall faroutside the graph.
Table 1

Derived input parameters for the CSPN.

3.4 Line nebular analysis

While the lines of H α, H β, [O III] and He i follow well the same spatial behavior of the main nebula rim (Fig. 8, upper panel), there is a clear concentration of the low ionized species [N II] and [S II] toward the outer edge (Fig. 8, middle). Moreover, the central concentration of the high ionized species He ii and [Ar IV] relative to the lower ionization lines of He i (Fig. 8, lower panel) is more pronounced. The four cuts given by the two slits from the central star outwards in radial distance show fairly the same shapes, hence the lines were averaged as radial profiles for the Cloudy modeling below. The line ratio [O III] 500.7 nm/H β appears constant at high values of seven to eight. On the other hand, the line ratio RHe = (He ii 468.6 nm/He i 587.5 nm) which is normally used to define for an excitation class of the nebula as a total (Ratag et al. 1997; Martins & Viegas 2002), appears flat at ≈3.1 for an inner radius of 30″ and declines near perfectly linear to 0.5 out at 65″. We derived the H α/H β ratio for both slit positions to test for a significant fraction or variations of the internal extinction within the nebula. Despite the differences in line intensities and some fluctuations between the two slit positions, we are unable to detect any of such variations (Fig. 9).

The RN II = ([N II] 654.8 nm + [N II] 658.4 nm)/[N II] 575.5 nm line ratio was used to derive an electron temperature as we shown in Fig. 10. The blue nitrogen line has a fairly low S/N, therefore the individual features observed at the line ratio should not be over-interpreted as intrinsic electron temperature variations, otherwise a spatial smoothing would be required. The temperatures lead to an average value of around Te = 11 500 K using the calibrations given in Osterbrock & Ferland (2006). There is a slight tendency to asymmetry in the northern part of the nebula, but in general terms it appears homogeneous with an electron temperature well above 104 K. Unfortunately, there were no other temperature sensitive line triples like e.g. [O III], [Ar III] and [S III] covered in our spectral range. Temperatures derived from higher ionized states outside the colder clumps would have been of certain interest (see shock analysis Sect. 3.5).

Further we used the RS II = ([S II] 671.6 nm/[S II] 673.1 nm) line ratio to derive the electron density ne in the nebula. We used the average electron temperature Te derived above and the two calibration ϱ(RS II, Te) provided by Proxauf et al. (2014). The values are very constant at the lowest limit of the diagnostic diagram as shown in Fig. 11. Thus, the values must be interpreted mainly as upper limits. As the flux varies significantly and thus the integrated volume emissivity which is a function of the average density ⟨ϱ⟩, we have to assume taking into account this nearly constant values for ϱ, that the filling factor varies over the nebula. This was observed in the NIR molecular analysis for the nearby Helix nebula (Matsuura et al. 2009) and was also a result of the modeling of NGC 2438 in Öttl et al. (2014).

thumbnail Fig. 8

Spatial distribution of the major spectral lines. Upper panel: brightest lines of oxygen, hydrogen and helium. Middle panel: the low ionized states of nitrogen and sulfur. The [N II] 654.8 nm line was skipped for clarityas it caries the same information as the [N II] 658.4 nm line. Lower panel: fainter lines of the two different ionization stages of helium. As the line around 471 nm shows the shape of the higher ionized He ii lines, we have to assume that it is dominated by the [Ar IV] 471.1 nm emission and fewer contribution of the He i 471.3 nm emission. In all panels, solid lines denote to the slit at position angle 30°, while the dashed lines give the slit at 150°. The region contaminated by the central star was eliminated. Around 30″ south of the nebula center, the FORS2 instrument has a chip gap.
thumbnail Fig. 9

H α/H β normalized to a Case B line ratio (Osterbrock & Ferland 2006) to test for internal extinction.
thumbnail Fig. 10

Electron temperature derived from [N II] lines. Despite the large variation in intensity itself, the temperatures are fairly homogeneous. Line coding for the two slit directions is as in Fig. 9.

3.5 Gas ionization tendency and shock heating

The line ratios log (Hα/[N II]) versus log (Hα/[S II]) have been used in the past to identify photoionized PNe from that excited by additional mechanisms, such as found in supernova remnants orin H ii regions. First introduced by Garcia Lario et al. (1991), this scheme is refined in Magrini et al. (2003) and extensively reviewed by Frew & Parker (2010). Öttl et al. (2014) extended this analysis to two dimensions to use it for the study of NGC 2438. We recently carried out a more detailed analysis using the multiple shell PNe NGC 3242, NGC 6826 and NGC 7662 (Barría & Kimeswenger 2018). After investigating on their so called extended 2D diagnostic diagrams (E2DD) we concluded that major deviations from the inside to outside ionization tendency appeared only at regions were FLIERs or LISs are located. Previous studies on excitation properties of FLIERs and LISs (see e.g. Akras & Gonçalves 2016; Akras et al. 2017), support for this scenario. Nevertheless, no FLIERs or LISs are visible in IC 5148, we followed the same strategy and studied the nebular ionization tendency through the main fragmentation structures in the nebula by means of its E2DD. As is shown in Fig. 12, where colors denote different nebular regions, we do not claim for regions at the main nebula showing significant deviations from the general ionization tendency. The regression line gives for a 1.218 ± 0.016 inclination. The fit does not include lines ratios from the halo structure. At the small PNe sample investigated in Barría & Kimeswenger (2018) we noted a possible correlation between the slope of the fit in the E2DD and the CSPNe temperature. However, the previously studied nebula exhibit similar properties such as subsolar abundances, low expansion velocities, high excitation ranges, similar kinetic ages, and the presence of FLIERs/LISs. IC 5148, on the other hand, shows for a higher velocity expansion, no FLIERs/LISs, a solar abundance, and is highly evolved. Despite the fact that IC 5148 follow the tendency we found in the previously studied nebulae, this result has to be treated with caution, until a further study on the E2DD slope values considering a bigger and diversified sample of PNe could be conducted.

Following the ideas presented in Guerrero et al. (2013) to identify for shock heated regions, we derived the [O III]/H α in the FORS2 spectra as well as in the CDK20 images. As the seeing is a function of wavelength, we used the CSPN to derive a variation of the PSF. Considering the S/N on the central star was very high, we were able to identify not only the central component, but also the wide wings, known to happen as a result of rare scintillation events of the air. The combined PSF, giving for an error of far below 1%, consists of a combined Gaussian with a component covering 94% of the peak intensity (84% of the integrated flux) and a wing with a FWHM that is 3.5 times wider with 6% of the peak intensity (16% of the total flux). The width varied with wavelength as expected being about 7% larger at 500 nm (1.′′ 05) than that we found at 660 nm (0.′′98). Nevertheless, the wide wing reaches out to 3.′′0, containing asignificant contribution to the total flux. This is of special importance as the line flux drops in the region of interest by more than a order of magnitude. While the main nebula shows a very homogeneous line ratio of around 2.5 (see Fig. 13), both, the latter spectral analysis as well as the image analysis displayed at Fig. 14, claim for shock features revealed by a thin shell of enhanced oxygen emission with line ratios significantly higher than ten. According to Guerrero et al. (2013), this skin of shock-excited material develops when a fast shell expands into a lower density medium. Moreover, here also the Mach number changes. To test the reliability of this result, the data were also folded with two and three times the derived PSF. As the feature in the line ratio gets only slightly weaker, respectively, by 20 and 60% in the ratio, but still giving a line ratio of five, we can be sure that it is not an artifact of the data analysis. Moreover, we find a slight steady increase already in the last 8–12′′ inside the front, where the intensity itself resembles still a nearly constantly high flux. We assume that this is caused by the afterglow of the heating in the post-shock region. The S/N of H α images as well as for the spectra in the halo is not sufficient to derive reasonable numbers for the regions outside the main nebula.

thumbnail Fig. 11

Density sensitive line ratio of [S II] lines. Despite the large variation in intensity itself, the temperatures are fairly homogeneous. Line coding for the two slit directions is as in Fig. 9.
thumbnail Fig. 12

Extended 2D diagnostic diagram of the nebula along the two slits positions and out to the halo in the regions with S/N >3 for all lines included in the graph.
thumbnail Fig. 13

Shock front analysis along the slit in the two FORS2 spectra using ≤[O III]/H α. The spectra were folded with a PSF derived from the central star to obtain same resolution. Solid and dashed lines indicate for the 30° and 150° slit spectra, respectively.
thumbnail Fig. 14

Shock front analysis as in Fig. 13 but using the CDK20 narrow band images. The color coding is linear for 0 ≤ [O III]/H α ≤ 5.

3.6 Cloudy model

Physical parameters of the different regions along the main nebula, were investigated by means of Cloudy modeling. The model was performed with version 17.00 of Cloudy (Ferland et al. 2017) and covering only the main nebula. The Cloudy code makes use of an extensive atomic and molecular database together to the possibility to include grains and incorporates a filling factor to simulate clumpiness.

The previously derived CSPN parameters such as the effective temperature, luminosity and surface gravity listed in Table 2, were used in the model as input parameters for the ionizing source. Thus, we modeled for a central initial mass star of Mini = 1.5 M⊙, a metallicity of Z = 0.02 dex and a modeled age of 8500 yr. The latter measured from the point where the post-AGB object has a log (Teff) = 3.85 (Miller Bertolami 2016). The spectral energy distribution of the CSPN was obtained via state-of-the-art H-Ni NLTE model stellar atmospheres provided by Rauch (2003). An expanding spherical geometry was assumed for the nebula. The net line emission provided by Cloudy per volume-unit (erg cm-3s-1), was used to create a model 2D projected emissivity profile along the line of sight by an in-house IDL code. In this way, we were able to make a direct comparison of the observed emissivity profiles to the modeled ones. Abundances were set in the Cloudy model to that of a typical metal rich PN as incorporated in the model code from Aller & Czyzak (1983) and Khromov (1989). An assumed Gaia DR2 distance of 1.3 kpc was used. Although, we also attempt to model the nebula by using the two additional intersection distance solutions found in Sect. 3.3 (1.8 and 1.5 kpc). However, the use of these higher distances did not reach convergence for a good model.

The modeling was carried out setting an initial tabulated density profile derived from the H α emissivity (see details of this procedure in Barría & Kimeswenger 2018). As was pointed out in Sect. 3.4, despite the large intensity variation of the [S II] lines, the density sensitive line ratio RSII is fairly constant. Hence, we have to assume that the filling factor varies through the nebula. Variations in clumpiness and filling factors were also observed in NGC 2438 in the same way (Öttl et al. 2014) and resolved even in molecular lines in other nearby nebulae (e.g. Matsuura et al. 2009, and references therein). The estimated line ratio RSII ~ 1.4 (see Sect. 3.4) give us an upper limit for the electron density of ne ~ 200 cm-3 (Osterbrock & Ferland 2006; Proxauf et al. 2014). It has to be noted here, that the hydrogen density profile in Cloudy is representative for the density within the filled fraction of the nebulosity. In this sense, the integrated nebular mass for such a clumpy model can only be achieved by multiplying the density by the corresponding filling factor in each region. Based on this initial density profile as a starting point, we manually iterated to convergence for the final model using the observed line profiles of [N II] 658.3 nm, [S II] 673.0 nm, He i 587.5 nm, He ii 468.5 nm, [O III] 500.7 nm and Hα 656.3 nm. The observed profiles were obtained by averaging the four probes taken radially by the two slit directions. Convergence was derived only by dividing the main nebula into three independent spherical shell regions with independently filling factors as free parameters. The final model was achieved after added the resulting best models for each region (see Fig. 15). Final derived Cloudy density profile (called hden) is also shown in Fig. 15. The resulting parameters for the best nebula model are listed in Table 3.

While we notice very good agreement between the observed and modeled [N II], [S II], He i, He ii and Hα emissivity profiles, we observed moderate deviations in the [O III] profile. A modeled underestimation on the [O III] emission becomes more evident at the outer regions of the main nebula. Considering this, we suggest here for two possible scenarios: one involving emission from a thin hot interclump medium and, a second one related to the residual heating of the post-shock afterglow. The Cloudy code definition of the filling factor assumes that the mass is distributed solely in the clumps. Thus, a possible additional contribution from ionized [O III] gas coming from the interclump medium, cannot be considered for our Cloudy clumpy model. Under the second proposed scenario, spectroscopic and image analysis gives evidence for a thin layer of enhanced emission due to a shock front. We might consider then for a possible enhancement of the [O III] emission, as a result of the heated post-shock region. This corresponds well also to the previous line analysis shown in Fig. 13.

thumbnail Fig. 15

Best fit Cloudy model. Solid lines are giving the model, dashed lines the data averaged in the four directions taken in the slit. Lowest panel: derived Cloudy hden density profile.
Table 2

Main derived nebular properties and nebular parameters obtained from the best Cloudy model.

4 Nebular halo

A very special set of features in this nebula is occupied by the outer halo. While many nebulae in this evolutionary stage do show one or even more halos, this one has a few very special features to be mentioned in detail.

4.1 Innerhalo

The inner halo is well represented by an exponential decline of the intensity into nearly all four slit directions equally (Fig. 16). It can be detected in Hα and [N II] (658.4 nm) up to a radius of about 100′′ (corresponding to 0.85 pc at a distance of 1.3 kpc) with a constant line ratio of Hα/[N II](658.4 nm) ≈ 4. A similar (or only marginally flatter) decline is seen for the [O III] line, which is nearly an order of magnitude brighter. Hence, this one can be followed up to the very outskirts (Fig. 16). All four slit positions give about the same [O III] decline until a radius of 110′′. Beyond that, only the two northern ones, not influenced by the outer halo bow (see Sect. 4.3), seem to continue straight as before and after the halo rings (between 110 and 135′′; see Sect. 4.2). At 175′′ (corresponding to 1.5 pc at D = 1.3 kpc) both signals vanishes in the noise. The high [O III]/Hα ≈ 8 (even with a slight steady increase toward larger radius) implies that the halo is ionized to a high degree (see modeling and discussion of the halo in NGC 2438 by Öttl et al. 2014). This ratio and its steady increase is caused bythe hardening of the radiation due to the rapid decline of the absorption efficiency of hydrogen κH ∝ χ−3 for photon energies χ >1 Ryd (Preite-Martinez & Pottasch 1983). The effect for the halos of MSPNe is discussed in detail in Guerrero & Manchado (1999). As we cannot measure the density of the halo directly, a low density ρ ∝ r−2 was used to test the continuum radiation transmitted by the CLOUDY model. The incident radiation to the halo is depleted between 1.1 and 2.5 Ryd to more than 85%, while the very hot CSPN having its radiation maximum above 3 Ryd is still very strong at ionization energies of He ii and O iii. The result fit well to this assumption of a photoionized material. As we do not have decent density information tracers, the results can used only to follow the tendency and not can be used in absolute quantities.

thumbnail Fig. 16

Brightest emission lines in the halo. Upper panel: two northern and one southeast section of the halo. Lower panel: southwest section. The location of the two halo rings are marked by gray lines.
thumbnail Fig. 17

[O III] image of the halo region. Left panel: after masking the main nebula and removinga Gaussian shaped decline representing the normal halo intensity distribution. Right panel: same image with two circles centered exactly at the CSPN at 115 and 128′′.

4.2 Halo rings

Multiple outer faint rings and arcs have been found at several planetary nebulae among different morphological types (see e.g. Terzian & Hajian 2000; Corradi et al. 2004; Ramos-Larios et al. 2016). While several different formation mechanisms have been proposed to explain these structures, their origin is still not that clear, but strongly connected to the mass-loss processes at the latest phase of the AGB. While these features appears scarcely among roundish PNe (~ 7%, Ramos-Larios et al. 2016), we found in the [O III] images of IC 5148, at a distance of 115 and 128′′, corresponding to 1.00 and 1.08 pc at D = 1.3 kpc, respectively, two exact concentric rings centered perfectly respect to the central star. Moreover, the images are deep enough to exclude the typically found multiple structures down by nearly an order of magnitude in brightness. To enhance visibility a power law shaped decline of the regular halo was subtracted from the images (see Fig. 17). The striking is not only the perfect circular shape, but also the fact, that they increase the [O III]/Hα ratio to at least 30 (lower limit), as these rings are invisible in the other lines. According to Corradi et al. (2004), the spacing between rings correlates with the post-AGB age of the nebula. Although IC 5148 appears older than the nebulae sample studied by these authors, the large spacing between the rings found here (Δ θ= 13′′, Δ r = 0.08 pc at 1.3kpc) reinforces the likelihood of this being a quite-evolved nebula.

To date, the most plausible scenario for the formation of such structures is correlated to periodic or quasi-periodic enhancements in the mass-loss rate at the end of the AGB phase. As Ramos-Larios et al. (2016) reviewed, several mechanisms have been suggested as being responsible for these periodic or quasi-periodic fluctuations such as solar-like magnetic activity in AGB stars (Soker 2000), long-period oscillations in AGB stars (Icke et al. 1992), viscous momentum transfer between grains and gas (Simis et al. 2001) or even binary interactions (Kim & Taam 2012a,b). Some other mechanisms as thermal pulses or AGB pulsations have been dismissed basically due to an inconsistency of the time lapses between successive observed rings. But these scenarios produce normally either only a single (or widely spaced) rings or a large number of them. Especially the dust-based mechanisms (e.g., Simis et al. 2001) are dominated by reflected light of the main nebula, which would not change the line ratio so much. Certainly a deep study on the origin of the rings found in IC 5148 is far beyond the scope of this paper, however as an alternative scenario and following the discussion about shock signatures in Guerrero et al. (2013), we proposed that these structures might be caused by shock fronts. As the other lines, especially also [S II] are not detected anymore out there, these features require deeper observations for further verification.

4.3 Outer halo bow

The “bow” is in the outermost part of the halo, peaking into a diffuse clump in the SW region and having a decent visible tail until the SE slit position. It extends from 145 to 200′′ from the central star, which corresponds to 0.92–1.26 pc at a distance of 1.3 kpc. It is at a completely different excitation state than the ring features described above. As shown in Fig. 16 it has an Hα/[N II] (658.4 nm) line ratio around unity and [O III](500.7 nm)/Hα of about 3. In this respect, it almost exactly resembles the state of the central nebula and gives an excitation below that found in the innermost halo region (from 70 to 100′′) inside the rings. It nicely continues the excitation direction of the main nebula from inside to outside as shown in Fig. 12. Thus, it is most likely photoionized material. Weather it was an unusual asymmetric mass loss event on the AGB, or interaction with ISM remains unclear. Although the latter is fairly unlikely at this large distance from the galactic plane. Also, the low foreground extinction and the homogeneity of the Balmer line ratio over the whole nebula indicate that a patchy cloudy environment is very unlikely.

5 Conclusions

We performed a detailed photometric and spectroscopic analysis of the high galactic latitude roundish nebula IC 5148, an older evolved nebulae, undisturbed from interaction with the interstellar matter. We find that the distance derived from Gaia DR2 corresponds well to those obtained by statistical methods. Moreover, the very blue color in the Gaia photometry indicates, as shown in Kimeswenger & Barría (2018), a good separation of the star from the nebula and give us confidence in the low formal error given by the trigonometric solution of Gaia which results in a distance of kpc. Furthermore, we were able to derive a more reliable interstellar foreground extinction leading to E(B − V) = 0. m03 with different independent methods. This is significantly lower than the hitherto often-used values of Kaler et al. (1990) and Kingsburgh & Barlow (1994). Medium-high resolution spectroscopy was used to derive a central star temperature of 140 kK. However, after comparison to state-of-the-art evolutionary tracks (Miller Bertolami 2016) we obtain a real value to be slightly lower and thus we adopted 130 kK and a gravity of log(g) ≈ 7 on the white dwarf cooling track. The progenitor therefore must have been a low mass star of initially 1.5 M⊙. This fits well to the expectation that very massive stars should not be found so far from the galactic plane. We find from our observationsof the central star, as well as from the nebular investigation, that IC 5148, despite its high galactic z, is not, similar to other nebulae in these regions (e.g. A 15, A 20, MeWe 1-3; Emprechtinger et al. 2004), underabundant in metals. This was somewhat unexpected and in contradiction to the postulation of Soker (2002), that spherical PNe might form only in metal poor environment.

IC 5148shows a perfect shock structure at its outermost edge of the main nebula, as predicted in the hydrodynamic models (e.g. Perinotto et al. 2004). It is only visible by using the high ionized species as already suggested and discussed in Guerrero et al. (2013). The shock sensitive diagrams, using low ionized species following Magrini et al. (2003), are not sensitive to this effect in the outermost regions.

The presented CLOUDY photoionization model shows generally a self-consistent model of the nebula. However, the observed excess of the forbidden [O III] line could not be modeled in that way. We speculate that this might be an afterglow of the after-shock region, responsible for driving up the less frequent oxygen atoms even by some small rest of overheated helium. A dedicated time-dependent radiative transfer code with source terms of shock heating and afterglows would be required for such a model.

The reported discovery of structured halo emission is enigmatic. The inner halo region resembles the perfect circular shape as that of the main nebula, while further out a bow like structure is seen at one side of the nebula only. While the latter showsline ratios as seen from normal photoionized material, the spherical rings have an extremely large [O III]/H α. Spectroscopy even deeper than ours would be required to resolve this puzzle.

Acknowledgements

We would like to thank the referee Martin Guerrero for his helpful comments. The study is based on observationsmade with ESO Telescopes at the La Silla Paranal Observatory under program ID 098.D-0332. This work has made use of data from the European Space Agency (ESA) mission Gaia (https://www.cosmos.esa.int/gaia), processed by the Gaia Data Processing and Analysis Consortium (DPAC, https://www.cosmos.esa.int/web/gaia/dpac/consortium). Furtheron this research has made use of the SIMBAD database, operated at CDS, Strasbourg, France (Wenger et al. 2000), the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration and has made use of “Aladin Sky Atlas” developed at CDS, Strasbourg Observatory, France (Bonnarel et al. 2000). D.B. and this research were financed by the FONDO ALMA-Conicyt Programa de Astronomía/PCI 31150001 and W.K. is funded by the Hochschulraumstrukturmittel provided by the Austrian Federal Ministry of Education, Science and Research (BMBWF).

References

  1. Acker, A., Gleizes, F., Chopinet, M., et al. 1982, Publication Speciale du Centre de Donnees Stellaires, 3 [Google Scholar]
  2. Acker, A., Marcout, J., Ochsenbein, F., et al. 1992, The Strasbourg-ESO Catalogue of Galactic Planetary Nebulae. Parts I, II. (European Southern Observatory) [Google Scholar]
  3. Akras, S., & Gonçalves, D. R. 2016, MNRAS, 455, 930 [NASA ADS] [CrossRef] [Google Scholar]
  4. Akras, S., Gonçalves, D. R., & Ramos-Larios, G. 2017, MNRAS, 465, 1289 [NASA ADS] [CrossRef] [Google Scholar]
  5. Aller, L. H., & Czyzak, S. J. 1983, ApJS, 51, 211 [NASA ADS] [CrossRef] [Google Scholar]
  6. Appenzeller, I., Fricke, K., Fürtig, W., et al. 1998, The Messenger, 94, 1 [NASA ADS] [Google Scholar]
  7. Balick, B. 1987, AJ, 94, 671 [NASA ADS] [CrossRef] [Google Scholar]
  8. Barría, D., & Kimeswenger, S. 2018, MNRAS, 480, 1626 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bessell, M. S. 1990, PASP, 102, 1181 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
  11. Blöcker, T. 1995, A&A, 299, 755 [NASA ADS] [Google Scholar]
  12. Bojičić, I. S., Parker, Q. A., & Frew, D. J. 2017, in Planetary Nebulae: Multi-Wavelength Probes of Stellar and Galactic Evolution, eds. X. Liu, L. Stanghellini, & A. Karakas, IAU Symp., 323, 327 [NASA ADS] [Google Scholar]
  13. Bonnarel, F., Fernique, P., Bienaymé, O., et al. 2000, A&AS, 143, 33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  14. Buckley, D., & Schneider, S. E. 1995, ApJ, 446, 279 [NASA ADS] [CrossRef] [Google Scholar]
  15. Cahn, J. H., Kaler, J. B., & Stanghellini, L. 1992, A&AS, 94, 399 [NASA ADS] [Google Scholar]
  16. Chu, Y.-H., Jacoby, G. H., & Arendt, R. 1987, ApJS, 64, 529 [NASA ADS] [CrossRef] [Google Scholar]
  17. Corradi, R. L. M., Schönberner, D., Steffen, M., & Perinotto, M. 2000, A&A, 354, 1071 [Google Scholar]
  18. Corradi, R. L. M., Sánchez-Blázquez, P., Mellema, G., Gianmanco, C., & Schwarz, H. E. 2004, A&A, 417, 637 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Dreyer, J. L. E. 1910, MNRAS, 59, 105 [Google Scholar]
  20. Emprechtinger, M., Forveille, T., & Kimeswenger, S. 2004, A&A, 423, 1017 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  21. Emprechtinger, M., Rauch, T., & Kimeswenger, S. 2005, A&A, 431, 215 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  22. Ferland, G. J. 2000, Rev. Mex. Astron. Astrofis., 27, 153 [Google Scholar]
  23. Ferland, G. J., Korista, K. T., Verner, D. A., et al. 1998, PASP, 110, 761 [Google Scholar]
  24. Ferland, G. J., Porter, R. L., van Hoof, P. A. M., et al. 2013, Rev. Mex. Astron. Astrofis, 49, 137 [Google Scholar]
  25. Ferland, G. J., Chatzikos, M., Guzmán, F., et al. 2017, Rev. Mex. Astron. Astrofis., 53, 385 [NASA ADS] [Google Scholar]
  26. Frank, A., & Mellema, G. 1994, A&A, 289, 937 [NASA ADS] [Google Scholar]
  27. Freudling, W., Romaniello, M., Bramich, D. M., et al. 2013, A&A, 559, A96 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  28. Frew, D. J., & Parker, Q. A. 2010, PASA, 27, 129 [NASA ADS] [CrossRef] [Google Scholar]
  29. Frew, D. J., Parker, Q. A., & Bojičić I. S. 2016, MNRAS, 455, 1459 [NASA ADS] [CrossRef] [Google Scholar]
  30. Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  31. Gale, W. F. 1897, Astron. Nachr., 143, 293 [NASA ADS] [CrossRef] [Google Scholar]
  32. Garcia Lario, P., Manchado, A., Riera, A., Mampaso, A., & Pottasch, S. R. 1991, A&A, 249, 223 [NASA ADS] [Google Scholar]
  33. Gesicki, K., Zijlstra, A. A., Hajduk, M., & Szyszka, C. 2014, A&A, 566, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  34. Guerrero, M. A., & Manchado, A. 1999, ApJ, 522, 378 [NASA ADS] [CrossRef] [Google Scholar]
  35. Guerrero, M. A., Toalá, J. A., Medina, J. J., et al. 2013, A&A, 557, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Hoffmeister, C. 1961, Die Sterne, 37, 204 [NASA ADS] [Google Scholar]
  37. Hua, C. T., Dopita, M. A., & Martinis, J. 1998, A&AS, 133, 361 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  38. Icke, V., Frank, A., & Heske, A. 1992, A&A, 258, 341 [NASA ADS] [Google Scholar]
  39. Kaler, J. B., & Feibelman, W. A. 1985, ApJ, 297, 724 [NASA ADS] [CrossRef] [Google Scholar]
  40. Kaler, J. B., Shaw, R. A., & Kwitter, K. B. 1990, ApJ, 359, 392 [NASA ADS] [CrossRef] [Google Scholar]
  41. Kausch, W., Noll, S., Smette, A., et al. 2015, A&A, 576, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  42. Khromov, G. S. 1989, Space Sci. Rev., 51, 339 [NASA ADS] [CrossRef] [Google Scholar]
  43. Kim, H., & Taam, R. E. 2012a, ApJ, 744, 136 [NASA ADS] [CrossRef] [Google Scholar]
  44. Kim, H., & Taam, R. E. 2012b, ApJ, 759, 59 [NASA ADS] [CrossRef] [Google Scholar]
  45. Kimeswenger, S., & Barría, D. 2018, A&A, 616, L2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Kingsburgh, R. L., & Barlow, M. J. 1994, MNRAS, 271, 257 [NASA ADS] [CrossRef] [Google Scholar]
  47. Kwok, S. 1982, ApJ, 258, 280 [NASA ADS] [CrossRef] [Google Scholar]
  48. Kwok, S. 2002, in Interacting Winds from Massive Stars, eds. A. F. J. Moffat & N. St-Louis, ASP Conf. Ser., 260, 245 [NASA ADS] [Google Scholar]
  49. Kwok, S., Purton, C. R., & Fitzgerald, P. M. 1978, ApJ, 219, L125 [NASA ADS] [CrossRef] [Google Scholar]
  50. Landolt, A. U. 1992, AJ, 104, 340 [NASA ADS] [CrossRef] [Google Scholar]
  51. Landolt, A. U., & Uomoto, A. K. 2007, AJ, 133, 768 [NASA ADS] [CrossRef] [Google Scholar]
  52. Lindegren, L., Hernández, J., Bombrun, A., et al. 2018, A&A, 616, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  53. Luri, X., Brown, A. G. A., Sarro, L. M., et al. 2018, A&A, 616, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  54. Magrini, L., Perinotto, M., Corradi, R. L. M., & Mampaso, A. 2003, A&A, 400, 511 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  55. Martins, L. P., & Viegas, S. M. 2002, A&A, 387, 1074 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  56. Matsuura, M., Speck, A. K., McHunu, B. M., et al. 2009, ApJ, 700, 1067 [NASA ADS] [CrossRef] [Google Scholar]
  57. Meatheringham, S. J., Wood, P. R., & Faulkner, D. J. 1988, ApJ, 334, 862 [NASA ADS] [CrossRef] [Google Scholar]
  58. Mellema, G., & Frank, A. 1995, MNRAS, 273, 401 [NASA ADS] [Google Scholar]
  59. Mendez, R. H. 1991, Evolution of Stars: the Photospheric Abundance Connection (the Netherlands: Springer), 145, 375 [Google Scholar]
  60. Miller Bertolami M. M. 2016, A&A, 588, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Modigliani, A., Goldoni, P., Royer, F., et al. 2010, in Observatory Operations: Strategies, Processes, and Systems III, Proc. SPIE, 7737, 773728 [CrossRef] [Google Scholar]
  62. Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
  63. Napiwotzki, R. 2001, A&A, 367, 973 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  64. Noll, S., Kausch, W., Kimeswenger, S., et al. 2014, A&A, 567, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  65. Okorokov, V. A., Shustov, B. M., Tutukov, A. V., & Yorke, H. W. 1985, A&A, 142, 441 [NASA ADS] [Google Scholar]
  66. Osterbrock, D. E., & Ferland, G. J. 2006, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Sausalito, CA: University Science Books) [Google Scholar]
  67. Öttl, S., Kimeswenger, S., & Zijlstra, A. A. 2014, A&A, 565, A87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. Parker, Q. A., Bojičić, I. S., & Frew, D. J. 2016, J. Phys. Conf. Ser., 728, 032008 [CrossRef] [Google Scholar]
  69. Patriarchi, P., & Perinotto, M. 1991a, A&AS, 91, 325 [NASA ADS] [Google Scholar]
  70. Patriarchi, P., & Perinotto, M. 1991b, Proc. Astron. Soc. Aust., 9, 309 [NASA ADS] [CrossRef] [Google Scholar]
  71. Perinotto, M., Kifonidis, K., Schoenberner, D., & Marten, H. 1998, A&A, 332, 1044 [NASA ADS] [Google Scholar]
  72. Perinotto, M., Schönberner, D., Steffen, M., & Calonaci, C. 2004, A&A, 414, 993 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  73. Phillips, J. P. 2005, MNRAS, 357, 619 [NASA ADS] [CrossRef] [Google Scholar]
  74. Preite-Martinez, A., & Pottasch, S. R. 1983, A&A, 126, 31 [NASA ADS] [Google Scholar]
  75. Proxauf, B., Öttl, S., & Kimeswenger, S. 2014, A&A, 561, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  76. Ramos-Larios, G., Santamaría, E., Guerrero, M. A., et al. 2016, MNRAS, 462, 610 [NASA ADS] [CrossRef] [Google Scholar]
  77. Ratag, M. A., Pottasch, S. R., Dennefeld, M., & Menzies, J. 1997, A&AS, 126, 297 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  78. Rauch, T. 2003, A&A, 403, 709 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  79. Rauch, T., & Ringat, E. 2011, Astronomical Data Analysis Software and Systems XX, eds. I. N. Evans, A. Accomazzi, D. J. Mink, & A. H. Rots, ASP Conf. Ser., 442, 563 [NASA ADS] [Google Scholar]
  80. Rauch, T., & Reindl, N. 2014, Asymmetrical Planetary Nebulae VI Conf., 77 [Google Scholar]
  81. Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103 [NASA ADS] [CrossRef] [Google Scholar]
  82. Schönberner, D., Jacob, R., Lehmann, H., et al. 2014, Astron. Nachr., 335, 378 [NASA ADS] [CrossRef] [Google Scholar]
  83. Simis, Y. J. W., Icke, V., & Dominik, C. 2001, A&A, 371, 205 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  84. Smette, A., Sana, H., Noll, S., et al. 2015, A&A, 576, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  85. Soker, N. 2000, ApJ, 540, 436 [NASA ADS] [CrossRef] [Google Scholar]
  86. Soker, N. 2002, A&A, 386, 885 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  87. Stanghellini, L., Shaw, R. A., & Villaver, E. 2008, ApJ, 689, 194 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  88. Stanghellini, L., Bucciarelli, B., Lattanzi, M. G., & Morbidelli, R. 2017, New A, 57, 6 [NASA ADS] [CrossRef] [Google Scholar]
  89. Steffen, M., Szczerba, R., & Schoenberner, D. 1998, A&A, 337, 149 [NASA ADS] [Google Scholar]
  90. Swift, L. 1899, MNRAS, 59, 568 [NASA ADS] [CrossRef] [Google Scholar]
  91. Tajitsu, A., & Tamura, S. 1998, AJ, 115, 1989 [NASA ADS] [CrossRef] [Google Scholar]
  92. Terzian, Y., & Hajian, A. R. 2000, in Asymmetrical Planetary Nebulae II: From Origins to Microstructures, eds. J. H. Kastner, N. Soker, & S. Rappaport, ASP Conf. Ser., 199, 33 [NASA ADS] [Google Scholar]
  93. van de Steene, G. C., & Zijlstra, A. A. 1994, A&AS, 108, 485 [NASA ADS] [Google Scholar]
  94. Vassiliadis, E.,& Wood, P. R. 1993, ApJ, 413, 641 [NASA ADS] [CrossRef] [Google Scholar]
  95. Vassiliadis, E.,& Wood, P. R. 1994, ApJS, 92, 125 [NASA ADS] [CrossRef] [Google Scholar]
  96. Vernet, J., Dekker, H., D’Odorico, S., et al. 2011, A&A, 536, A105 [Google Scholar]
  97. Weinberger, R. 1989, A&AS, 78, 301 [NASA ADS] [Google Scholar]
  98. Wenger, M., Ochsenbein, F., Egret, D., et al. 2000, A&AS, 143, 9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  99. Zhang, C. Y. 1995, ApJS, 98, 659 [NASA ADS] [CrossRef] [Google Scholar]
5

In Bessell et al. (1998) Table A.2, the zp(fλ) and the zp(fν) are reversed in the formulas for the magnitudes – confirmed by author in private communication.

7

Image Reduction and Analysis Facility available at http://iraf.noao.edu/

All Tables

Table 1

Derived input parameters for the CSPN.

In the text
Table 2

Main derived nebular properties and nebular parameters obtained from the best Cloudy model.

In the text

All Figures

thumbnail Fig. 1

Co-added direct images in [O III] (top panel) and H α (bottom panel). The slit positions of the spectroscopy are indicated in the top panel. While in the inner part, a linear intensity scale covering the whole range from zero to peak flux of the nebula was used, for regions with flux below 5% a steep asinh between 0.01 and 1% of the peak of the linear part was used to enhance the by orders of magnitude fainter structures.
In the text
thumbnail Fig. 2

Spectral lines from the medium-high resolution spectroscopy of the UVB Arm of X-shooter. Overlaid are the model grid spectra along the evolutionary track log (g) values of a CSPN according to Miller Bertolami (2016) from the high resolution stellar spectra grid from the German Astrophysics Virtual Observatory (GAVO) grid resources (Rauch & Reindl 2014).
In the text
thumbnail Fig. 3

UV photometry from IUE and GALEX with our calibrated FORS2 spectrum. The model spectral energy distributions from TMAP are overlaid, normalized to the FORS2 data.
In the text
thumbnail Fig. 4

Colorequations of hot white dwarfs from the Landolt photometric standards (Landolt 1992) and from the spectrophotometric standards (Landolt & Uomoto 2007) versus Gaia DR2. The target position used for calibrations is given in the lower panel. While all other stars are hot nearby compact white dwarfs without PNe, the encircled target is the central star of the Helix nebula.
In the text
thumbnail Fig. 5

λ4 weighted FORS2 spectrum with the TMAP model. We intentionally did not remove the telluric O2 and H2 O absorptions at the red end of the spectrum. The blue horizontal lines show the required correction of the line depths at the blue end due to the deviation to the continuum.
In the text
thumbnail Fig. 6

Expansion of the PN and system velocity in the Galaxy from X-shooter spectra after correction for LSRK.
In the text
thumbnail Fig. 7

Upper panel: evolutionary tracks of Miller Bertolami (2016) as log(g) vs. log (T) for initial masses of 1.00, 1.25, 1.5 and 2.0 M⊙ (color coded) for metallicities of Z = 0.01 and Z = 0.02 (line thickness coded). The corresponding spectroscopic distances are indicated in light gray. The vertical dark gray lines indicate the range of temperatures as derived from our spectroscopy. Dotted lines mark extrapolation of the track, which ended earlier. Lower panel: computed model ages from the tracks above (solid lines) vs. the dynamical ages (dashed lines). The circles indicate the solutions by intersections of the model age with the dynamical age. The corresponding solutions found in the lower panel are marked in the upper panel as well. Crosses mark the estimated solution of the extrapolated track. The intersections of the 1.0 M⊙ (red) fall faroutside the graph.
In the text
thumbnail Fig. 8

Spatial distribution of the major spectral lines. Upper panel: brightest lines of oxygen, hydrogen and helium. Middle panel: the low ionized states of nitrogen and sulfur. The [N II] 654.8 nm line was skipped for clarityas it caries the same information as the [N II] 658.4 nm line. Lower panel: fainter lines of the two different ionization stages of helium. As the line around 471 nm shows the shape of the higher ionized He ii lines, we have to assume that it is dominated by the [Ar IV] 471.1 nm emission and fewer contribution of the He i 471.3 nm emission. In all panels, solid lines denote to the slit at position angle 30°, while the dashed lines give the slit at 150°. The region contaminated by the central star was eliminated. Around 30″ south of the nebula center, the FORS2 instrument has a chip gap.
In the text
thumbnail Fig. 9

H α/H β normalized to a Case B line ratio (Osterbrock & Ferland 2006) to test for internal extinction.
In the text
thumbnail Fig. 10

Electron temperature derived from [N II] lines. Despite the large variation in intensity itself, the temperatures are fairly homogeneous. Line coding for the two slit directions is as in Fig. 9.
In the text
thumbnail Fig. 11

Density sensitive line ratio of [S II] lines. Despite the large variation in intensity itself, the temperatures are fairly homogeneous. Line coding for the two slit directions is as in Fig. 9.
In the text
thumbnail Fig. 12

Extended 2D diagnostic diagram of the nebula along the two slits positions and out to the halo in the regions with S/N >3 for all lines included in the graph.
In the text
thumbnail Fig. 13

Shock front analysis along the slit in the two FORS2 spectra using ≤[O III]/H α. The spectra were folded with a PSF derived from the central star to obtain same resolution. Solid and dashed lines indicate for the 30° and 150° slit spectra, respectively.
In the text
thumbnail Fig. 14

Shock front analysis as in Fig. 13 but using the CDK20 narrow band images. The color coding is linear for 0 ≤ [O III]/H α ≤ 5.
In the text
thumbnail Fig. 15

Best fit Cloudy model. Solid lines are giving the model, dashed lines the data averaged in the four directions taken in the slit. Lowest panel: derived Cloudy hden density profile.
In the text
thumbnail Fig. 16

Brightest emission lines in the halo. Upper panel: two northern and one southeast section of the halo. Lower panel: southwest section. The location of the two halo rings are marked by gray lines.
In the text
thumbnail Fig. 17

[O III] image of the halo region. Left panel: after masking the main nebula and removinga Gaussian shaped decline representing the normal halo intensity distribution. Right panel: same image with two circles centered exactly at the CSPN at 115 and 128′′.
In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.