Issue |
A&A
Volume 642, October 2020
|
|
---|---|---|
Article Number | A71 | |
Number of page(s) | 10 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/202038049 | |
Published online | 06 October 2020 |
Rapid evolution of [WC] stars in the Magellanic Clouds⋆
Space Radio-Diagnostics Research Centre, University of Warmia and Mazury, ul.Oczapowskiego 2, 10-719 Olsztyn, Poland
e-mail: marcin.hajduk@uwm.edu.pl
Received:
30
March
2020
Accepted:
30
June
2020
We obtained new spectra of fourteen Magellanic Cloud planetary nebulae with the South African Large Telescope to determine heating rates of their central stars and to verify evolutionary models of post-asymptotic giant branch stars. We compared new spectra with observations made in previous years. Five planetary nebulae showed an increase in excitation over time. Four of their central stars exhibit [WC] features in their spectra, including three new detections. This raises the total number of [WC] central stars of PNe in the Magellanic Clouds to ten. We compared determined heating rates of the four [WC] central stars with the He-burning post-asymptotic giant branch evolutionary tracks and the remaining star with the H-burning tracks. Determined heating rates are consistent with the evolutionary models for both H and He-burning post-asymptotic giant branch stars. The central stars of the PNe that show the fastest increase of excitation are also the most luminous in the sample. This indicates that [WC] central stars in the Magellanic Clouds evolve faster than H-burning central stars, and they originate from more massive progenitors.
Key words: stars: AGB and post-AGB / planetary nebulae: general / stars: evolution / Magellanic Clouds
The reduced spectra are only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/cat/J/A+A/642/A71
© ESO 2020
1. Introduction
Stars with initial masses in the range of 1 − 8 M⊙ lose most of their envelope on the asymptotic giant branch (AGB). During an early AGB phase, helium burning proceeds in a thin shell on the top of the carbon-oxygen core. When the helium shell reaches the hydrogen shell, the nuclear burning becomes unstable. A series of helium shell flashes punctuate quiescent hydrogen burning.
The large energy generation due to a helium flash induces a convective instability between the shells (Iben & Renzini 1983; Herwig 2005). When the H shell becomes active again, part of the intershell material is dredged up to the surface.
The envelope of the star is rapidly expelled due to the intensive mass loss experienced during the AGB phase. This reduces the number of helium shell flashes. The mass ejected by AGB stars significantly contributes to the chemical evolution of galaxies and stellar systems (Karakas 2010).
After a star ejects most of its envelope, it terminates the AGB phase and its evolution accelerates. Heating rate increases from Ṫ∗ ≲ 0.1 kK yr−1 to 1 K yr−1 ≲ Ṫ∗ ≲ 10 kK yr−1 (Miller Bertolami et al. 2019). The evolutionary timescale of post-AGB stars is a very strong function of their envelope mass and the nuclear burning rate. Post-AGB stars occupy a relatively narrow range in mass between 0.5 and 0.9 M⊙, but evolutionary timescales differ by a few orders of magnitude within this range.
Post-AGB stars cross the H-R diagram from low to high temperatures almost at constant luminosity, as long as the shell sources remain active (Paczyński 1970). Once the star is hot enough, the ionization front expands throughout the shell ejected during the AGB phase. The shell forms a planetary nebula (PN). The excitation of the PN increases with the increasing temperature of the central star.
A central star of a PN may experience a final thermal pulse after it leaves AGB (Schoenberner 1979). A final thermal pulse can punctuate evolution of a post-AGB star immediately after it leaves the AGB (called AFTP), during its horizontal evolution, when H-burning shell is active (late thermal pulse – LTP), or during cooling (very late thermal pulse – VLTP; Koesterke 2001).
In the result of the VLTP, the pulse driven convecting zone eventually reaches out into the H-rich envelope (Herwig 2005). The LTP and AFTP lead to the diluted shell and significantly reduced H abundance. An AFTP can lead to both a considerable enrichment with carbon and oxygen and to the dilution of hydrogen. This depends on the envelope mass. The smaller envelope mass, the less hydrogen is left in the envelope (Blöcker 2001).
H-deficient central stars constitute at least 30% of all the central stars in the Galaxy (Weidmann & Gamen 2011). More than half of them are [WC] stars, which show emission line spectra typical for WR stars of carbon sequence (Acker & Neiner 2003). Hydrogen-deficient central stars are believed to expose their intershell region in their atmospheres due to a late helium-shell flash (Werner & Herwig 2006).
The evolution of post-AGB stars remains one of the most difficult topics in stellar evolution (Miller Bertolami et al. 2019). The models must rely on assumptions and approximations such as mixing length theory or stellar wind theory. Observational verification is necessary for the validation of the evolutionary models.
Stellar luminosities and temperatures have been traditionally used to study the evolution of post-AGB stars. However, post-AGB stars occupy a narrow range in luminosities in the H-R diagram. Their distances are little known and the stars emit most of their energy in the ultraviolet range. Observational errors propagate large uncertainties in the determination of the luminosities.
An alternative set of parameters useful for studying the evolution of post-AGB stars is the stellar temperature and heating rate. The heating rate is a sensitive measure of the stellar mass. The pace of the temperature evolution is determined by the ratio of the stellar luminosity and remaining envelope mass.
Gesicki et al. (2014) derived heating rates from stellar temperatures and nebular ages for a sample of Galactic Bulge PNe. They needed to accelerate the evolutionary models by a factor of three to fit the white dwarf mass distribution and asteroseismological masses of central stars of PNe.
New models by Miller Bertolami (2016) confirmed the results obtained by Gesicki et al. (2014). He obtained three to ten shorter timescales of post-AGB evolution than Bloecker (1995) and Vassiliadis & Wood (1994). His models were also brighter by ∼0.1 − 0.3 dex for the same mass. He took advantage of an updated treatment of the constitutive microphysics and included an updated description of the mixing processes and winds.
An alternative method to derive heating rate is to measure the temperature evolution of the central star. This could be derived indirectly by comparison of the nebular line fluxes in different epochs. The change of the [O III] 5007 Å-to-H β flux ratio in time is sensitive to the temperature evolution of cool central stars (Hajduk et al. 2015).
In this paper, we derived heating rates for a sample of central stars in the Magellanic Clouds and compared them with the existing evolutionary models. An advantage of Magellanic Cloud objects over Galactic PNe is the precisely determined distance. Moreover, Magellanic Cloud PNe are often compact and spatially unresolved by ground based telescopes. In such cases, the ion stratification in the nebula would not affect the measured flux ratios.
2. Observations
We observed the fourteen PNe with the South African Large Telescope (SALT) equipped with the Robert Stobie Spectrograph (Kobulnicky et al. 2003; Burgh et al. 2003). The SALT telescope has a fixed azimuthal angle. The Magellanic Clouds were accessible with the SALT near the meridian (hour angle from −2 to 2), so that the observing track length was sufficiently long. Under these conditions, atmospheric refraction should not play a significant role.
Three chips recorded each spectrum in different wavelength ranges: 4345–5347 Å, 5403–6393 Å, and 6448–7397 Å. We used longslit with the projected width of two arcsec. The resulting resolution was about 800 in the center of the spectral range. The slit was oriented north-south. The range of air masses for the targets and standards was 1.21 − 1.36 (Tables 1 and 2). The log of observations is given in Table 1.
Log of the PNe observations.
Log of the standard star observations.
The spectra were extracted and wavelength calibrated. We calibrated the flux scale using a standard observed during the same night. On three nights: November 10, December 21, and December 22, flux standards were not observed. For these three nights we used the sensitivity curve derived for the closest night.
We measured the emission line fluxes with the Gaussian fit. We computed the errors of individual lines according to the formulae provided by Lenz & Ayres (1992). We added the derived error (root mean-square) of the sensitivity curve. We computed an error using at least two closest sensitivity curves if the standard was not observed during the given night. The fluxes were dereddened using the extinction law given by Fitzpatrick (1999) and the total to selective extinction ratio RV = 3.1 (Table A.1). We derived the extinction from the observed Hα-to-Hβ line flux ratio.
3. New [WC] central stars
Hydrogen- and helium-burning post-AGB stars have different evolutionary tracks. Helium-burning tracks show slower evolution by approximately a factor of three in comparison to H-burning tracks for the same remnant mass.
In order to compare the evolution of the central stars of PNe with the evolutionary tracks, one needs to know which shell is active. This can be inferred from the photospheric abundances of a star. A hydrogen-poor atmosphere indicates that the star is a helium burner.
Most of the H-deficient central stars show [WC] spectra (Weidmann & Gamen 2011). They are relatively easy do identify, even if the stellar continuum is unobservable, since they show wide and prominent stellar wind lines imposed on the nebular spectrum (Acker & Neiner 2003). A hydrogen-rich atmosphere does not unambiguously indicate which shell is the main source of the energy in the star. However, we assumed the stars that do not show emission lines to be H burners.
One of the stars in the sample (MGPN SMC 8) is a known [WC] central star (Pena et al. 1997). We detected three more [WC] stars in our sample: SMP LMC 31, SMP LMC 55, and SMP SMC 1. They show emission lines of C III, C IV, and He II. The He II 4686 Å line in SMP LMC 31 has a full width at half maximum (FWHM) of 13 Å, while SMP LMC 55 and SMP SMC 1 have FWHMs of C IV 4650 Å lines of 12 Å and 8 Å, respectively. Neither of the new [WC] stars show the C IV 5806 Å line. Thus, we classify them as [WC 11] type (Acker & Neiner 2003). However, faint C II 7231 and 7236 Å lines would suggest an earlier type.
4. Heating rates
We combined new and archival spectra of PNe to study the temperature evolution of their central stars. We used the [O III] 5007 Å-to-Hβ line ratio. This ratio is sensitive to the changes of stellar temperature for cool central stars due to a much higher ionization potential of O+ than hydrogen. Both lines are relatively strong and located close to each other, so wavelength-dependent calibration errors are reduced.
We fit a linear function to the [O III] 5007 Å-to-Hβ line flux ratio observed in different epochs for each PN (Fig. 1). We discarded one point by Reid & Parker (2006) for SMP LMC 64 from the fitting, which exceeded all other points by a factor of two. The fit slope was larger than one, at least by 3σ in five PNe, indicating an increasing stellar temperature.
![]() |
Fig. 1. Evolution of [O III] 5007 Å-to-Hβ line flux ratio in the Magellanic Cloud PNe. The dashed line fits the observed trend. The horizontal error bars correspond to the flux uncertainties. The vertical error bars in the time domain correspond to the observations, of which the dates were not precisely specified. Meatheringham & Dopita (1991) did not specify dates of their observations precisely, thus the error bars in time domain are large. |
Among the PNe showing the flux evolution, one (MGPN SMC 8) was observed five times, three PNe were observed four times, and one PN has three observations. The observations show a linear increase of the [O III] 5007 Å-to-Hβ flux ratio in time. An exception is SMP LMC 55. The third observation by Reid & Parker (2006) is significantly higher than the last one. However, the [O III] 5007 Å-to-Hβ flux ratio determined by Reid & Parker (2006) appears to be systematically higher than fluxes measured by other authors.
We performed spherically symmetric photoionization models of these five PNe with the Cloudy v. 17.01 code. All the models were constrained by the observed absolute Hβ flux. We fit the nebular line flux ratios observed in 2013 (Table A.1).
The nebular sizes and stellar luminosities are available for three objects observed with the Hubble Space Telescope by Villaver et al. (2003) and Villaver et al. (2004). However, for the PN MGPN SMC 8 temperature and luminosity, determination is most likely affected by the stellar contribution to the He II 4686 Å line flux (Shaw et al. 2010). For this and two other PNe, stellar luminosity was a free parameter. The adopted luminosities, or those derived from our photoionization models L*, are given in Table 3.
Input parameters and results of the photoionization modeling.
We varied stellar temperature to obtain the best fit. We used atmosphere models by Castelli et al. (2003). The grid covers the temperature range of 3.5 − 50 kK.
We used nebular abundances by Leisy & Dennefeld (2006) as an input parameter, but varied it to better fit the observed emission line fluxes. For some cases (e.g., MGPN LMC 39, SMP LMC 31, and SMP SMC 1) we were not able to achieve the observed [O III] 5007 Å-to-Hβ line flux ratio with the O/H abundances derived by Leisy & Dennefeld (2006). At low metallicities, the [O III] 5007 Å-to-Hβ line flux ratio depends on the nebular O/H abundance and stellar temperature (Stasińska 2007). For the average SMC oxygen abundance, the [O III] 5007 Å-to-Hβ flux ratio reaches a maximum of five at a temperature of about 100 kK (Stanghellini et al. 2003). For the LMC, the models reach maximum of ten at similar stellar temperatures.
Table 3 shows the input parameters and results of the photoionization modeling of the five objects. The modeled line fluxes are compared to the fluxes measured in 2013 in Table A.4.
We varied stellar temperature to fit the [O III] 5007 Å-to-Hβ line flux ratio observed in other epochs. After we obtained the temperatures in different years, we used them to compute heating rates Ṫ∗. The uncertainty of the measured [O III] 5007 Å-to-Hβ line flux ratio propagates to the measured temperatures in different epochs and subsequently to the heating rates presented in Table 3.
We computed heating rates from the evolutionary tracks as a function of stellar temperature. We used models by Miller Bertolami (2016) for the star that does not show emission lines, and Vassiliadis & Wood (1994) for four [WC] central stars. Miller Bertolami (2016) computed only three He-burning models spanning a very low range of masses.
Heating rates derived from the evolutionary tracks were interpolated to fit observed heating rates. We determined model-dependent luminosities LMOD and corresponding model-dependent masses MMOD. The errors of LMOD and MMOD shown in Table 3 result from the uncertainties in Ṫ∗.
We did not include systematic errors in measured heating rates in Table 3. These depend mostly on the determination of the O/H abundance. For an adopted uncertainty of the log (O/H) of 0.3 in dex, the model-dependent stellar mass MMOD is accurate to about 0.04 M⊙, and log of luminosity to 0.15 dex. Moreover, the model-dependent masses for [WC] stars derived using tracks by Vassiliadis & Wood (1994) may be overestimated by 0.038 M⊙ (Gesicki et al. 2014).
Heating rates derived from the evolutionary tracks do not vary strongly with stellar temperature. The uncertainty of stellar temperature determination affects the derived model-dependent parameters to a much lower extent than the uncertainty in the determination of O/H abundance.
5. Discussion
We attempted to select PNe with relatively cool central stars to study their evolution. Our sample fulfils the criterion of [O III] 5007 Å-to-Hβ ratio ≲3 for all but one object. However, this criterion appeared to be insufficient in some cases. The [O III] 5007 Å line may be weak in hot PNe, due to photoionization of O++ to O+ + + (Marigo et al. 2001; Stasińska 1989). Also, the [O III] 5007 Å line can be weak in PNe with low O abundances, despite a high temperature. This appears to be the case for three PNe in our sample: Jacoby LMC 17, MGPN LMC 35, and SMP SMC 28, which show prominent He II 4686 Å lines. These PNe are not expected to show rapid evolution of the [O III] 5007 Å flux.
Young PNe can mimic symbiotic stars due to high nebular density and strong dust emission. Iłkiewicz & Mikołajewska (2017) listed two PNe from our sample, SMP LMC 31 and SMP LMC 104A, as symbiotic star candidates. However, we did not find any indication of the cool giant in their spectra. In addition, SMP LMC 31 show stellar wind lines consistent with a [WC] star.
Eleven out of fourteen PNe in our sample may be considered as young PNe. Five of them show an evolution of the [O III] 5007 Å line flux. The five central stars showing temperature evolution are most likely the most massive in our sample, and thus the most luminous. An independent luminosity determination confirms that these stars show the highest luminosities (Table 4).
Luminosities of the central stars of PNe obtained using different methods: model-dependent luminosities determined from heating rates, luminisities derived from Zanstra temperatures and stellar magnitudes, from the integrated spectral energy distribution, and constrained from abolute Hβ flux and nebular modeling.
The remaining six central stars may evolve too slowly to show significant evolution of the [O III] 5007 Å-to-Hβ flux ratio on a timescale of decades. The collected observations are sensitive down to heating rates of about 5 − 10 K yr−1. Lower heating rates, corresponding to the initial masses below 1.25 M⊙ (final masses below 0.532 M⊙) (Miller Bertolami et al. 2019), would not be detected.
Figure 2 shows the modeled evolution of nebular [O III] 5007 Å-to-Hβ flux ratio with increasing stellar temperature for two different O/H abundance ratios. The [O III] 5007 Å-to-Hβ flux ratio scales approximately with the O abundance. The flux ratio does not increase linearly throughout the whole temperature range of the central star, but it can be approximated as linear in sufficiently small temperature intervals. The heating rate is highest between 35 and 38 kK. Then it slows down. The five central stars of which the nebulae show flux evolution have stellar temperatures in this range.
![]() |
Fig. 2. Two models of the evolution of the [O III] 5007 Å-to-Hβ flux ratio with the temperature of the central star in a PN. The model parameters are identical to MGPN LMC 39, and the second model the (O/H) abundance is lowered by 0.3 in dex. |
The central star of MGPN LMC 39 is probably the most massive H-burning central star in our sample. Its final mass of 0.57 M⊙ is the same as the peak of the mass distribution for the central stars of Galactic PNe. The remaining H-burning central stars in our sample are less massive. This suggests that central stars in the Magellanic Clouds are less massive than the stars in the Galaxy. Hajduk et al. (2014) derived heating rates of 45 ± 7 K yr−1 for Galactic PN Hen 2-260, significantly higher than for Magellanic Cloud central stars.
We compared the model-dependent luminosities LMOD with luminosities determined independently using different methods: using Zanstra temperatures and a black-body fit (log T* (H I)), by fitting spectral energy distribution (SED) (LSED) and by photoionization modeling (L*) (Table 4). Unfortunately, the luminosities yielded from observations using different methods show a relatively large scatter.
The model-dependent luminosity agrees with the observed luminosity for SMP SMC 1 and is within the observed range of stellar luminosities for other objects. For SMP LMC 55, the model-dependent luminosity is close to the luminosity derived from photoionization modeling by Dopita & Meatheringham (1991), but lower by 0.56 dex from the luminosity of our photoionization model. However, the absolute Hβ flux used by Dopita & Meatheringham (1991) was 0.46 lower in log than that used by us from Reid & Parker (2006). If this discrepancy were real, this would indicate an unexpected change in luminosity of this object.
We are not able to discriminate between Miller Bertolami (2016) and Vassiliadis & Wood (1994) or Bloecker (1995). Miller Bertolami (2016) models were faster and more luminous for the same masses. However, these models all show similar evolutionary speeds for the same luminosity of the central star (but different masses). We would be able to discriminate between these models if we obtained an independent mass determination.
The [WC] central stars constitute about 15% of Galactic central stars of PNe. This ratio was suggested to be as low as 5% in the Magellanic Clouds (Pena et al. 1997). Pena et al. (1997) reported only two [WC] central stars in the Small Magellanic Cloud (SMC) and six in the Large Magellanic Cloud (LMC). Five stars are intermediate types [WC 4] or [WC 4-5], and only one is [WC 9] (Monk et al. 1988). Margon et al. (2020a) and Margon et al. (2020b) announced the discovery of two more [WC 11] stars in the LMC. They suggested that a large population of late-type [WC] stars may exist, undetected due to observational selection. Late-type [WC] stars show faint emission lines, which could be missed in observations performed with smaller telescopes. We discovered three new late-type [WC] stars.
Interestingly, the four [WC] PNe observed by us show rapid evolution. This implies relatively high stellar mass. Hajduk et al. (2015) observed a relatively rapid change of the [O III] 5007 Å-to-Hβ flux ratio for the Galactic PNe M 1-11 and M 1-12 containing [WC] central stars. This suggests that the evolution of Galactic [WC] stars also is faster than H-rich central stars and originate from more massive progenitors. This was independently confirmed by Weidmann & Gamen (2011), who found that H-poor central stars are more concentrated toward the Galactic center and Galactic plane than the H-rich group. The DB white dwarfs, which are the possible progeny of [WC] central stars, show an average mass of 0.651 M⊙, higher than the 0.598 M⊙ derived for DA white dwarfs, which originate from H-rich central stars (Kleinman et al. 2013).
Acknowledgments
MH thanks the Ministry of Science and Higher Education (MSHE) of the Republic of Poland for granting funds for the Polish contribution to the International LOFAR Telescope (MSHE decision no. DIR/WK/2016/2017/05-1) and for maintenance of the LOFAR PL-612 Baldy (MSHE decision no. 59/E-383/SPUB/SP/2019.1). Polish participation in SALT is funded by the MSHE Grant No. DIR/WK/2016/07. MH acknowledges financial support from National Science Centre, Poland, grant No. 2016/23/B/ST9/01653. Some of the observations reported in this paper were obtained with the Southern African Large Telescope (SALT) under programme 2013-2-POL_OTH-001 (PI M. Hajduk).
References
- Acker, A., & Neiner, C. 2003, A&A, 403, 659 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Blöcker, T. 2001, Ap&SS, 275, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Bloecker, T. 1995, A&A, 299, 755 [NASA ADS] [Google Scholar]
- Boroson, T. A., & Liebert, J. 1989, ApJ, 339, 844 [NASA ADS] [CrossRef] [Google Scholar]
- Burgh, E. B., Nordsieck, K. H., Kobulnicky, H. A., et al. 2003, in Prime Focus Imaging Spectrograph for the Southern African Large Telescope: Optical Design, eds. M. Iye, & A. F. M. Moorwood, SPIE Conf. Ser., 4841, 1463 [Google Scholar]
- Castelli, F., & Kurucz, R. L. 2003, in Modelling of Stellar Atmospheres, eds. N. Piskunov, W. W. Weiss, & D. F. Gray, IAU Symp., 210, A20 [Google Scholar]
- Dopita, M. A., & Meatheringham, S. J. 1991, ApJ, 367, 115 [NASA ADS] [CrossRef] [Google Scholar]
- Fitzpatrick, E. L. 1999, PASP, 111, 63 [NASA ADS] [CrossRef] [Google Scholar]
- Gesicki, K., Zijlstra, A. A., Hajduk, M., & Szyszka, C. 2014, A&A, 566, A48 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hajduk, M., van Hoof, P. A. M., Gesicki, K., et al. 2014, A&A, 567, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hajduk, M., van Hoof, P. A. M., & Zijlstra, A. A. 2015, A&A, 573, A65 [CrossRef] [EDP Sciences] [Google Scholar]
- Herwig, F. 2005, ARA&A, 43, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Iben, I. Jr., & Renzini, A. 1983, ARA&A, 21, 271 [NASA ADS] [CrossRef] [Google Scholar]
- Iłkiewicz, K., & Mikołajewska, J. 2017, A&A, 606, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kamath, D., Wood, P. R., & Van Winckel, H. 2014, MNRAS, 439, 2211 [NASA ADS] [CrossRef] [Google Scholar]
- Karakas, A. I. 2010, MNRAS, 403, 1413 [NASA ADS] [CrossRef] [Google Scholar]
- Kleinman, S. J., Kepler, S. O., Koester, D., et al. 2013, ApJS, 204, 5 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Kobulnicky, H. A., Nordsieck, K. H., Burgh, E. B., et al. 2003, in Prime Focus Imaging Spectrograph for the Southern African Large Telescope: Operational Modes, eds. M. Iye, & A. F. M. Moorwood, SPIE Conf. Ser., 4841, 1634 [Google Scholar]
- Koesterke, L. 2001, Ap&SS, 275, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Leisy, P., & Dennefeld, M. 2006, A&A, 456, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lenz, D. D., & Ayres, T. R. 1992, PASP, 104, 1104 [NASA ADS] [CrossRef] [Google Scholar]
- Margon, B., Manea, C., Williams, R., et al. 2020a, ApJ, 888, 54 [CrossRef] [Google Scholar]
- Margon, B., Massey, P., Neugent, K. F., & Morrell, N. 2020b, ApJ, 898, 85 [CrossRef] [Google Scholar]
- Marigo, P., Girardi, L., Groenewegen, M. A. T., & Weiss, A. 2001, A&A, 378, 958 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meatheringham, S. J., & Dopita, M. A. 1991, ApJS, 75, 407 [NASA ADS] [CrossRef] [Google Scholar]
- Miller Bertolami, M. M. 2016, A&A, 588, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miller Bertolami, M. M. 2019, in IAU Symposium, eds. F. Kerschbaum, M. Groenewegen, & H. Olofsson, IAU Symp., 343, 36 [Google Scholar]
- Monk, D. J., Barlow, M. J., & Clegg, R. E. S. 1988, MNRAS, 234, 583 [NASA ADS] [CrossRef] [Google Scholar]
- Paczyński, B. 1970, Acta Astron., 20, 47 [NASA ADS] [Google Scholar]
- Pena, M., Ruiz, M. T., & Torres-Peimbert, S. 1997, A&A, 324, 674 [Google Scholar]
- Reid, W. A., & Parker, Q. A. 2006, MNRAS, 373, 521 [NASA ADS] [CrossRef] [Google Scholar]
- Schoenberner, D. 1979, A&A, 79, 108 [NASA ADS] [Google Scholar]
- Shaw, R. A., Stanghellini, L., Villaver, E., & Mutchler, M. 2006, ApJS, 167, 201 [NASA ADS] [CrossRef] [Google Scholar]
- Shaw, R. A., Lee, T.-H., Stanghellini, L., et al. 2010, ApJ, 717, 562 [NASA ADS] [CrossRef] [Google Scholar]
- Stanghellini, L., Shaw, R. A., Mutchler, M., et al. 2002, ApJ, 575, 178 [NASA ADS] [CrossRef] [Google Scholar]
- Stanghellini, L., Shaw, R. A., Balick, B., et al. 2003, ApJ, 596, 997 [NASA ADS] [CrossRef] [Google Scholar]
- Stasińska, G. 1989, A&A, 213, 274 [Google Scholar]
- Stasińska, G. 2007, ArXiv e-prints [arXiv:0704.0348] [Google Scholar]
- van Aarle, E., van Winckel, H., Lloyd Evans, T., et al. 2011, A&A, 530, A90 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vassiliadis, E., & Wood, P. R. 1994, ApJS, 92, 125 [NASA ADS] [CrossRef] [Google Scholar]
- Vassiliadis, E., Dopita, M. A., Morgan, D. H., & Bell, J. F. 1992, ApJS, 83, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Villaver, E., Stanghellini, L., & Shaw, R. A. 2003, ApJ, 597, 298 [NASA ADS] [CrossRef] [Google Scholar]
- Villaver, E., Stanghellini, L., & Shaw, R. A. 2004, ApJ, 614, 716 [NASA ADS] [CrossRef] [Google Scholar]
- Weidmann, W. A., & Gamen, R. 2011, A&A, 526, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Werner, K., & Herwig, F. 2006, PASP, 118, 183 [Google Scholar]
Appendix A: Line fluxes
Measured and dereddened line fluxes in the LMC PNe with the relevant uncertainties.
Measured and dereddened line fluxes in the MC PNe with the relevant uncertainties.
Measured and dereddened line fluxes in the SMC PNe with the relevant uncertainties.
Modeled and observed line-flux ratios in the five MC PNe showing the evolution of the [O III] 5007 Å line flux.
All Tables
Luminosities of the central stars of PNe obtained using different methods: model-dependent luminosities determined from heating rates, luminisities derived from Zanstra temperatures and stellar magnitudes, from the integrated spectral energy distribution, and constrained from abolute Hβ flux and nebular modeling.
Measured and dereddened line fluxes in the LMC PNe with the relevant uncertainties.
Measured and dereddened line fluxes in the MC PNe with the relevant uncertainties.
Measured and dereddened line fluxes in the SMC PNe with the relevant uncertainties.
Modeled and observed line-flux ratios in the five MC PNe showing the evolution of the [O III] 5007 Å line flux.
All Figures
![]() |
Fig. 1. Evolution of [O III] 5007 Å-to-Hβ line flux ratio in the Magellanic Cloud PNe. The dashed line fits the observed trend. The horizontal error bars correspond to the flux uncertainties. The vertical error bars in the time domain correspond to the observations, of which the dates were not precisely specified. Meatheringham & Dopita (1991) did not specify dates of their observations precisely, thus the error bars in time domain are large. |
In the text |
![]() |
Fig. 2. Two models of the evolution of the [O III] 5007 Å-to-Hβ flux ratio with the temperature of the central star in a PN. The model parameters are identical to MGPN LMC 39, and the second model the (O/H) abundance is lowered by 0.3 in dex. |
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.