Issue |
A&A
Volume 635, March 2020
|
|
---|---|---|
Article Number | A128 | |
Number of page(s) | 13 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201937118 | |
Published online | 24 March 2020 |
Planetary nebulae seen with TESS: Discovery of new binary central star candidates from Cycle 1
1
Departamento de Astrofísica, Centro de Astrobiología (INTA-CSIC), ESAC Campus, Camino Bajo del Castillo s/n, 28692 Villanueva de la Cañada, Madrid, Spain
e-mail: alba.aller@cab.inta-csic.es
2
Spanish Virtual Observatory, Spain
3
Instituto de Astrofísica de Canarias, 38205 La Laguna, Tenerife, Spain
4
Departamento de Astrofísica, Universidad de La Laguna, 38206 La Laguna, Tenerife, Spain
5
Instituto de Astrofísica de Andalucía – CSIC, C/Glorieta de la Astronomía s/n, 18008 Granada, Spain
Received:
15
November
2019
Accepted:
12
February
2020
Context. It has become clear in recent years that binarity plays a crucial role in many aspects of planetary nebulae (PNe), particularly with regard to the striking morphologies they exhibit. To date, there are nearly 60 known binary central stars of PNe (bCSPNe). However, both theory and observation indicate that this figure represents only the tip of the iceberg, with the Galactic PN population hosting orders of magnitude more stars.
Aims. We are involved in a search for new bCSPNe with the aim of enhancing the statistical validation of the key role of binarity in the formation and shaping of PNe. New discoveries of bCSPNe and their characterization carry important implications not only for understanding PN evolution, but also for studying binary evolution and the common-envelope phase, which is still poorly understood.
Methods. We used data from the TESS satellite to search for variability in the eight CSPNe that belong to the two-minute cadence of preselected targets in Cycle 1, with their available pipeline-extracted light curves. We identified strong periodicities and analysed them in the context of the binary scenario.
Results. All the CSPNe but one (Abell 15) show clear signs of periodic variability in TESS. The cause of this variability can be attributed to different effects, some of them requiring the presence of a companion star. We find simple sinusoidal modulations in several of the systems, compatible with irradiation effects. In addition, two of the central stars (PG 1034+001 and NGC 5189) also show photometric variations due to ellipsoidal variations and other signs of variability that are probably caused by star spots or relativistic Doppler-beaming. The case of the well-studied Helix Nebula is of particular interest; here we constructed a series of binary models to explain the modulations we see in the light curve. We find that the variability constrains the possible companion to be very low-mass main-sequence star or sub-stellar object. We also identify, in substantial detail, the individual pulsation frequencies of NGC 246.
Key words: planetary nebulae: general / techniques: photometric / binaries: general
© ESO 2020
1. Introduction
In recent years, a broad consensus has been reached on the role of central star binarity in explaining the departure from sphericity in planetary nebulae (PNe). Single-star scenarios can no longer explain the wide variety of striking morphologies shown by these stellar end products. There are more than 3000 PN identified in the Milky Way (Frew & Parker 2010) and roughly 80% of them have complex morphologies which depart from sphericity; this includes, for example, elliptical, bipolar, and multipolar PNe, as well as other structures like jets, disks, and knots present in many of the PNe (see e.g. Manchado et al. 1996a,b; Sahai & Trauger 1998; Parker et al. 2006). Since the discovery of the first binary central star of a planetary nebula (bCSPN) by Bond (1976), the total sample has been steadily growing. To date, there are fewer than 60 binary central stars known (see Boffin & Jones 2019). However, the amazing variety of morphologies observed makes it highly likely that binaries are mainly responsible for the formation of a significant fraction (if not the majority) of non-spherical PNe. This means that the true population of bCSPNe should be far larger than what is currently known.
The orbital periods of the known sample of bCSPNe cluster between two hours and four days, approximately (Boffin & Jones 2019). Just a few systems populate the longer period zone, although this is probably an observational bias since they are really difficult to detect through photometric variability. This is why most of the long period binary central stars have been discovered through radial velocity observations (see e.g. Van Winckel et al. 2014).
With regard to the close binary fraction (i.e. the fraction of PNe hosting binaries with periods between a few hours and a few days), recent surveys have shown that they represent about 12–21% (Miszalski et al. 2009) of the CSPNe, although this number may be a lower limit, especially given the observational challenges involved in constraining this fraction. Most of the known bCSPNe have been detected through modulations in the light curve (see e.g. Bond 2000; Miszalski et al. 2009), revealing the presence of close companions with very short orbital periods (Miszalski et al. 2009, 2011). Apart from eclipsing binaries, typical photometric modulations include irradiation effects (i.e. the reflection of the light of the hot primary by the less massive companion due to the large temperature difference between them), ellipsoidal modulations (i.e. tidal distortions due to the proximity of both stars) and Doppler beaming effects (which are the photometric imprints of the Doppler effect). These short-period binary systems must have undergone a common-envelope (CE) phase (Iben & Livio 1993), during which the primary fills its Roche lobe and transfers mass to its companion in an unstable way, forming a CE within which the two components spiral-in towards each other, resulting in a much closer binary system and the ejection of the CE. This process culminates in the formation of a (most likely aspherical) PN.
Space missions like NASA’s Transiting Exoplanet Survey Satellite (TESS; Ricker et al. 2015) provide a unique opportunity to detect new binary central stars. TESS was launched on 18 April 2018 and its science operations were begun after a three-month commissioning period on 25 July 2018. It is a two-year, all-sky photometric survey charged with the main goal of obtaining high-precision photometry of more than 200 000 selected bright stars with a cadence of about two minutes, observing in 26 different sky sectors. With four cameras observing each TESS data sector, the resulting field-of-view is 24 × 96°, which allows it to obtain photometry of any target within this wide field. These full frame images (FFIs) are read out approximately every 30 min, establishing the cadence of the photometry for the targets in the FFIs. As with Kepler, the high-precision, high-cadence obeservations and the continuous photometric data that TESS provides in a long time span (in comparison with the limitations that we have from the ground) make TESS an ideal machine for searching for bCSPNe.
In this paper, we use TESS observations to find new candidate short-period binaries in the nuclei of PNe. Even though the initial sample is not very large (only eight central stars were observed in the short-cadence mode in Cycle 1), we detected significant photometric variability for seven targets, that is, more than 80% of the sample. In Sect. 2, we describe the PNe observed by TESS, the TESS observations, and the analysis of the data. In Sect. 3, we discuss the variability and periodicities of each central star and we conclude in Sect. 4 with some final remarks.
2. TESS photometry
2.1. The PN sample
In its first year of observations, TESS observed the southern ecliptic hemisphere in a total of 13 sectors. This includes roughly half of the CSPNe in our Galaxy. However, only eight of them are among the targets observed at the two-minute cadence and, therefore, their light curves have already been published in the archive. These eight PNe are summarised in Table 1. The rest of the southern PNe observed by TESS are available in the FFIs, with a longer cadence of 30 min. They do not have light curves published in the archive, but their photometry can be obtained directly from the images. We are currently in the process of extracting the light curves from the FFIs in order to obtain a complete dataset for all the CSPNe, the results of which will be presented in a future publication.
Log of the TESS observations for our CSPNe sample.
2.2. Description of the TESS data
As mentioned above, TESS is equipped with four identical cameras covering each sector. Each of these four cameras has four 2k × 2k CCDs with a pixel scale of 21 arcsec pixel−1. The detector bandpass covers the spectrum from 600 to 1000 nm and is centred on 786.5 nm, just as for the traditional Cousins I-band. We used the light curves processed by the Science Processing Operations Center (SPOC) pipeline (Jenkins et al. 2016) and retrieved them from the Mikulski Archive for Space Telescopes (MAST1). SPOC light curves contain both the SAP (Simple Aperture Photometry) flux and the Pre-search Data Conditioning SAP (PDCSAP) flux. In the PDCSAP, long-term trends have been removed from the fluxes using the co-trending basis vectors. For the purpose of this work, we used the PDC light curves without any additional treatment of the data, only removing those data points with a non-zero “quality” flag since they could be affected by some anomalies, such as cosmic ray events or instrumental issues (see Sect. 9 in the TESS Science Data Products Description Document2).
We downloaded the target pixel file (TPF) of each object to inspect the aperture and possible contamination sources. This check of the TPFs is especially important because of the faintness of our sample in comparison with most of the two-minute cadences targets observed by TESS. Figure 1 shows the TPF of each of the eight central stars in the sample, marked with a white cross. The red circles represent the Gaia sources from the DR2 catalogue (Gaia Collaboration 2018), scaled by magnitude contrast against the target source. We over-plotted all the sources with a magnitude contrast up to Δm = 3 (i.e. three magnitudes fainter than our target, which corresponds to ∼6% of contamination if entirely inside the aperture). The aperture mask used by the pipeline to extract the photometry was also plotted over the TPF. In five cases (NGC 7293, PG 1034+001, Abell 7, Abell 15, and RWT 152), there is no contamination from other sources in the aperture mask. In the case of NGC 246 there is minimal contamination by an object of Gaia magnitude 14.2 (2.4 magnitudes fainter than our target). Finally, NGC 5189 and NGC 2867 have severe contamination by several sources of similar magnitudes to the central star, which makes the photometry difficult to interpret. We will discuss these cases in detail in Sect. 3.
![]() |
Fig. 1. Target pixel files (TPFs) of each central star in the sample. The red circles are the sources of the Gaia DR2 catalogue in the field with scaled magnitudes (see legend). White crosses indicate the location of the targets. The aperture mask used by the pipeline to extract the photometry is also marked. Pixel scale is 21 arcsec pixel−1. These plots have been created with the tpfplotter code, publicly available at https://github.com/jlillo/tpfplotter. |
2.3. Data analysis and modelling
A Lomb-Scargle analysis was performed for each object in order to identify possible periodic variability in the data. This first inspection resulted in the detection of strong periodicities in all but one of the targets (Abell 15). We analysed the data in the time-flux domain with a simple sinusoidal model that includes the amplitude (A), the zero level (Z0), period (P), and the time of conjunction (T0) as parameters. In this first analysis, we assumed uniform priors on all of the parameters, namely P ∈ 𝒰(0.01, 10) days, A ∈ 𝒰(0, 40) parts per thousand (ppt), T0 ∈ 𝒰(tmin, tmax) days, and Z0 ∈ 𝒰(0.9, 1.0).
In order to correctly sample the posterior probability distribution of each of these parameters, we used the implementation of Goodman & Weare’s affine invariant Markov chain Monte Carlo (MCMC) ensemble sampler emcee3, developed by Foreman-Mackey et al. (2013). We used a first phase with 30 walkers and 5000 steps to find a first guess for the solution. Then we re-sampled the position of the walkers in a small N-dimensional ball around the best walker and ran a second phase with the same number of walkers and 2500 steps. All steps from the latter phase were used to compute the marginalised posterior distributions of each parameter (no thinning was applied). The final chains are thus composed of 75 000 steps that were used to compute the marginalised posterior probabilities for each parameter.
In the case of Abell 7, two significant periodicities appear in the periodogram. The MCMC is not able to distinguish between these two similar periods. So we ran two separate MCMCs with Gaussian priors centered on each of the two periodicities (P1 = 3.0 d, and P2 = 2.6 d, see Sect. 3.7). The periodogram of the residuals indicate that after removing one of the sinusoidal signals, the other periodicity still remains above the 0.1% FAP (False Alarm Probability). Hence, we conclude that both frequencies are independent and probably come from different sources.
In two cases (PG 1034+001 and NGC 5189), the simple sinusoidal model did not produce a successful fit because of the presence of a two uneven maxima in the phase-folded light curves. We interpret this as a combination of ellipsoidal variations, Doppler beaming, and reflection. We then proceeded with a new model including all these three effects with the shape:
where θ = 2π(t − T0)/P). The amplitudes of the ellipsoidal, Doppler beaming, and reflection effects are represented by the parameters Aell, Abea, and Aref, respectively. It is important to note, however, that the model aims at representing the phase dependency for the three effects (ellipsoidal beaming and reflection) but does not put physical constraints on the physical parameters. Hence, it must be taken as the combination of sinusoidal functions with maxima at phases 0.25, 0.5, and 0.75. We proceeded in the same manner as explained above to obtain the posterior distributions of the parameters involved in this new model and obtained successful evidence against the simple sinusoidal model.
After applying the appropriate model to the each system, all the chains converged to the parameters shown in Table 3 (for the objects showing only a simple sinusoidal) and Table 4 (for the cases of PG 1034+001 and NGC 5189, for which the three different effects have been modelled). These results are discussed in detail in Sect. 3.
It is important to note that applying a purely sinusoidal model to the effects described can only be considered a rough estimation as the true light curve morphologies are a complex function of the stellar and systemic parameters (gravity and limb darkening, orbital inclination, etc.).
The precision of the TESS data also allows us to estimate the radial velocity amplitude of the stars in which a Doppler beaming factor is (a priori) present (as explained by van Kerkwijk et al. 2010, for the case of KOI-74). Thus, following the Eqs. (2) and (3) from van Kerkwijk et al. (2010), we can calculate the variation of the photon rate nγ observed by a telescope as follows:
where fDB is the Doppler boost pre-factor, λ = 7865 Å from TESS bandpass and T is the temperature of the primary. As such, we can obtain an estimation for the binary mass function. We note that the functional form of the equations of van Kerkwijk et al. (2010) assumes that one component of the binary completely dominates the observed flux from the system (as is the case for KOI-74). In systems where the stars have similar luminosities, the Doppler beaming effect of both stars leads to a sum that is destructive.
Figures 2–6 show the time-flux light curve (top panels), the Lomb-Scargle periodograms (right panels) as well as the phase-folded light curves of each central star with the period derived in the fitting (left panels). The FAPs at 10%, 1%, and 0.1% are indicated in each periodogram with grey horizontal dotted lines. Also, we have over-plotted in blue the periodogram of the residuals. For the purpose of visualisation, weighted means are employed to obtain the binned fluxes presented in the phase folded light curves. Bin sizes were chosen to ensure 20 and 100 datapoints per orbit for the red and grey datapoints shown, respectively.
![]() |
Fig. 2. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 7293 and RWT 152. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3. |
![]() |
Fig. 3. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 246 and PG 1034+001. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3 (for NGC 246) and Table 4 (for PG 1034+001). |
![]() |
Fig. 4. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 5189 and NGC 2867. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3 (for NGC 2867) and Table 4 (for NGC 5189). |
![]() |
Fig. 5. Time-flux light curve (top panel), phase-folded light curves (left panels) and periodograms (right panels) of Abell 7. Two different bin sizes are shown with red and grey circles. The periods derived from the corresponding fittings can be found in Table 3. |
![]() |
Fig. 6. Time-flux light curve (top panels), phase-folded light curve (left panel) and periodogram (right panel) of Abell 15. Two different bin sizes are shown with red and grey circles. |
3. Variability analysis and discussion
In the following section we describe the analysis of each central star individually and discuss the phase variations.
3.1. NGC 7293
Widely known as Helix nebula, NGC 7293 is the nearest PN to us and possibly one of the most well-studied in the literature, with more than 800 references (most of these related to its amazing and complex morphology). It has been well-characterised thanks to numerous observations from a wide range of wavelengths and telescopes. It is a double ring PN, with an inner thick disk and multitude of cometary knots present in the nebular centre (see e.g. Meaburn et al. 1998; Henry et al. 1999; O’Dell et al. 2005; Van de Steene et al. 2015, and references therein). In addition, the whole nebula is surrounded by a massive molecular envelope (Young et al. 1999). The central nucleus is a hot (Teff = 100 000 K; Napiwotzki 1999) DAO white dwarf with unexpected hard X-ray emission (see e.g. Leahy et al. 1994; Guerrero et al. 2001). The origin of this emission is suggested to be an unseen dMe (dwarf M star with emission lines) companion, which would be also consistent with the variable Hα emission line profile found by Gruendl et al. (2001). However, there has beeen no evidence (until now) for such a companion in spectroscopic or photometric observations, nor in the HST imaging. Also, Su et al. (2007) detected a dust disk around the central star by means of Spitzer observations. This disk was interpreted as a debris disk, although it differs substantially from the typical dust disc found around other white dwarfs.
The TESS light curve (see Fig. 2, first panel) shows for the first time a clear photometric variability with 2.77 d period. The signal in the periodogram is very strong and above the 0.1% FAP. As mentioned in Sect. 2.2, there is no contamination in the aperture photometry by other sources in the field (within a magnitude limit of Δm = 3) which means that the variability is real and definitively comes from NGC 7293. The light curve phase-folded with the 2.77 d period (Fig. 2, first left panel) shows an amplitude of 1.7 ppt (see also Table 3; 1 ppt equals to 1.086 mmag). If the variability were attributable to ellipsoidal modulations, then the true period would be double and the light curve would display two minima per cycle (at each opposition). However, this interpretation is unlikely given that the period would imply that the size of the tidally distorted star would be rather large, inconsistent with the measured surface gravity of the hot central star (Napiwotzki 1999) or a secondary with low enough luminosity to have escaped detection in previous studies. Perhaps a more reasonable interpretation of the variability is that it originates from the irradiation of a cool companion (e.g. Horvat et al. 2019). In order to explore this hypothesis, we constructed a series of models using the PHOEBE 2 code (Prša et al. 2016; Jones et al. 2019) assuming that the mass, temperature, and radius of the primary corresponds to what is determined in Napiwotzki (1999) and that the secondary is a main sequence star with mass, temperature, and radius relying on the empirical relations of Eker et al. (2018). The primary was assumed to be a blackbody with limb-darkening coefficients extrapolated from Gianninas et al. (2013), while the main sequence secondaries were modelled using PHOENIX atmospheres (Husser et al. 2013) and interpolated limb-darkening as described in Prša et al. (2016). The inclination of the binary was fixed based on the inclination of the nebula i = 28° (Henry et al. 1999) as is expected and observed in all systems for which both inclinations are known (Hillwig et al. 2016).
We did not encounter any satisfactory fits for secondaries with masses in the range of 0.16 M⊙ ≤ M2 ≤ 2.5 M⊙ (the range covered by the model atmospheres and empirical mass-radius-temperature relations used, which do not imply a secondary so bright that it would likely have been discovered by previous studies), with all models producing significantly greater variability amplitudes than those observed. Even adjusting the parameters of the system, for example using the more modern evolutionary tracks of Miller Bertolami (2016), which imply a lower mass and smaller radius than those used by Napiwotzki (1999) or the lower inclination of O’dell (1998, i ∼ 21°), a satisfactory fit could not be found. This leaves a number of possible solutions: (1) The variability may be due to irradiation of a smaller body (i.e. a sub-stellar or planetary mass object) or (2) the variability may not be due to irradiation. In this case, the variability could be due to Doppler beaming (van Kerkwijk et al. 2010, see Sect. 3.4 for more details), or perhaps due to spots on a rotating star. However, the scenario in which the Doppler beaming is the cause of the variability could be practically discarded, since the radial velocity amplitude derived following Eq. (2) (vr ≈ 460 km s−1) would involve a mass function of f(m) = (Porb/2πG) vr = (m2 sin i)3/(m1 + m2)2 = 28 M⊙, which is highly unlikely. The clearest test of these possibilities will be high-resolution time-series spectroscopy, capable of differentiating the radial velocity signals that would be associated with these hypotheses.
3.2. RWT 152
RWT 152 is a very faint, bipolar planetary nebula with an sdO in its nucleus (Aller et al. 2016). The morphology and the displacement of the central star with respect to the nebula’s centre may indicate the presence of a possible binary system. From the high-resolution, long-slit spectra, Aller et al. (2015) derived an inclination of 84 degrees for the equatorial plane of the nebula, which implies that the possible bCSPN (if any) would be almost edge-on.
The periodogram of the TESS light curve is quite different from the rest of the central stars in the sample (see Fig. 2, second panel). We find a forest of significant periodicities above the 0.1% FAP between 0.1 and 10 days approximately. The MCMC analysis found the most prominent peak (corresponding to ∼1.7 d period) as the most probable, with an amplitude of 1.6 ppt (similar to that found in NGC 7293). The nature of the most significant periodicity as well as the other significant peaks in the periodogram is not easily interpretable and radial velocity observations are essential for unveiling it.
One possibility is that the observed variability is caused by pulsations or stellar rotation. Thanks to a frequency analysis in the time domain (Barceló Forteza et al. 2015), we found the most significant peak at 6.2 μHz (or 1.9 days) and its closest peak of similar amplitude at 5.6 μHz (or 2.1 days). Both peaks may correspond to the rotation period of a star with differential rotation. We also see multiples of this signature such as twice (∼11 μHz) and thrice (∼17 μHz) its value (see Table 2) with a non-significant deviation in frequency, lower than the resolution (7%). These signatures may be produced by the presence of spots. In addition, we used the wavelet transform (Mathur et al. 2010) to find a proper value of rotation, obtaining similar results: 2.0 ± 0.6 days (5.8 ± 1.6 μHz).
Multiples of rotational peaks of the power spectrum of RWT 152.
3.3. NGC 246
NGC 246 exhibits an elliptical, shell-like appearance with a clumpy structure. The western edge is considerably brighter than the rest of the shell. This, and the displacement of the PN nucleus from the geometrical center of the nebula support the idea of the interaction with the ISM (for more details, see e.g. Muthu et al. 2000; Szentgyorgyi et al. 2003). The central star of this PN (HIP 3678) is a known hierarchical triple stellar system with an early to mid-K dwarf (K2–K5, named HIP 3678 B) and a mid M dwarf (M5–M6, named HIP 3678 C) as companions of the primary white dwarf (see Adam & Mugrauer 2014). This hot, hydrogen-deficient planetary nebula nucleus (PNN) was also identified as a low-amplitude pulsator with periods in the range of 24–31 min (i.e. 543–683 μHz) and amplitude of ∼2 mmag (Ciardullo & Bond 1996). Later, González Pérez et al. (2006) found significant daily changes in these pulsations, with the appearance of a new peak in the power spectra, which is peculiar in comparison with other observed PNNi, which have some variations in the amplitudes of certain peaks but never the appearance of a new peak.
In the TESS light curve we clearly identify the periodic photometric variability due to the aforementioned pulsations (with several significant peaks in the periodogram from 550 to 700 μHz approximately, see Fig. 7). For the first time, we are able to resolve all the individual frequencies due to the different pulsations modes. We can now conclude that what previous authors observed was, in fact, the broad envelope corresponding to all the “forest” of frequencies in the periodogram, which are concentrated in two different distributions in the periodogram. However, González Pérez et al. (2006) and Ciardullo & Bond (1996) did not have enough precision to resolve all the individual frequencies. We reserve a more detailed analysis of these pulsations for a future publication (Sowicka et al., in prep.).
![]() |
Fig. 7. Periodogram of NGC 246 in the low frequencies. The different peaks are due to the pulsations. |
In addition to the frequencies associated with pulsations, we found a new periodic photometric variability with a period of 6.8 days and an amplitude of nearly 0.6 ppt (see Fig. 3, first panel). This variability may be ascribed to a new source in the system. It is important to note that, as mentioned above, the photometry of this central star could be contaminated by another source that is inside the aperture mask. This source is HIP 3678 B and has a Gaia magnitude of G = 14.2 (2.4 magnitudes fainter than the central star HIP 3678 A). Both stars are in the same pixel of the TESS camera, which makes impossible to separate the contribution of both sources. The difference in magnitudes of the two stars corresponds to a ∼11% of contamination of the fainter star over the brighter. Hence, the origin of the variability is not clear and dedicated ground-based photometry on the two sources is necessary for unveiling its origin.
3.4. PG 1034+001
The PN surrounding this hot hydrogen-deficient DO white dwarf is probably the largest known PN on the sky. It was discovered by using the Sloan Digital Sky Survey spectra and confirmed with narrow-band Hα and [O III] images, which revealed a complex and non-defined nebulosity with an angular size of about 2° of diameter (Hewett et al. 2003). Later, Rauch et al. (2004) reported the detection of other much more extended emission structures around this central star.
In the Lomb-Scargle periodogram of the TESS dataset (see Fig. 3, second panel), we find two prominent peaks above the 0.1% FAP, at 1.85 d and 0.93 d. We can definitively state that the signal at 0.93 d must be an alias, since when we recalculate the residual Lomb-Scargle periodogram after subtracting the signal at 1.85 d, the other peak at 0.93 d is no longer present. Figure 3 shows the light curve phase-folded with the strongest period of 1.85 d. At first glance, it is possible to interpret the variability as the result of at least two different effects. On the one hand, we can identify a sinusoidal curve with two maxima per orbital period, which points to ellipsoidal modulation variability. This is the imprint in the light curve of the tidal distortions caused by the gravitational interaction between the two components of a binary system: the secondary star acts as a gravitational perturber to the external layers of the primary, inducing changes in the spherical shape of the star resulting in variability in the light curve along the orbital path. On the other hand, the asymmetry in the ellipsoidal modulation could indicate the presence of relativistic Doppler-beaming effects (Zucker et al. 2007). This effect is basically the photometric signature of the radial velocity variations. The shifted spectral energy distribution due to the velocity variations produces flux variability as seen by a wavelength-fixed bandpass. Its detection from ground based telescopes is challenging because of the low amplitudes that this effect produces (≤0.1%; Zucker et al. 2007) and the high precision and continuous observations along a long-timespan needed to detect it. To date, the only bCSPN exhibiting Doppler-beaming effects is the central star of the multi-shell PN AMU 1 (Aller et al. 2013), found in the analysis of the Kepler light curves by De Marco et al. (2015). First Kepler and now TESS have opened the window to this new detection technique of stellar variability, in particular in the field of bCSPNe. Finally, irradiation effects (i.e. the reflection of the light of the hot primary by a less massive companion) could also be present. In the phase-folded light curve, the final model is a combination of the equations for the three effects (reflection, ellipsoidal variability, and Doppler beaming). The amplitudes of each effect are listed in Table 3. We can conclude that the predominant effect would be Doppler beaming and that reflection is the least significant. However, some incongruities can be derived from these values, especially from the amplitude of the Doppler beaming, which we will discuss in the following.
Results from the light curve study for the five central stars analyzed with a simple sinusoidal model.
Results from the light curve study for the two central stars analyzed with models including reflection, ellipsoidal and doppler beaming effects.
It is interesting to compare the light curve of PG 1034+001 to that of the binary star KOI-74, which consists of an A-type star and a white dwarf (Bloemen et al. 2012; van Kerkwijk et al. 2010). The main difference between both light curves is that apart from reflection, ellipsoidal modulation, and Doppler beaming, KOI-74 also shows eclipses in its light curve. In the case of PG 1034+001, the amplitude of the possible Doppler beaming variations is a factor ∼30 larger than the same amplitude in KOI-74. Following Eq. (2) we can estimate the radial velocity amplitude of the star from the Doppler beaming (as we did in the case of NGC 7293). Assuming a temperature of Teff = 115 000 K for the primary (Werner et al. 2017), we derive fDB = 1.08, which corresponds to a radial-velocity semi-amplitude of vr ≈ 1095 km s−1. This yields a mass function of f(m) = 252 M⊙. With one of the components being a white dwarf (with 0.62 M⊙ determined by Werner et al. 2017) this would imply a minimum mass for the companion of ≈252 M⊙, essentially ruling out Doppler beaming.
Another more plausible scenario to account for such a huge amplitude is that the variability we are detecting is the consequence of the O’Connell effect (O’Connell 1951; Milone 1968; Johnston et al. 2019), whereby ellipsoidally-modulated binaries are frequently found to present uneven maxima in their light curves. It had previously been suggested that this could be the result of a so-called “periastron effect” in which the binaries are somewhat eccentric and experience an increased irradiation effect close to periastron (Roberts 1906); however it quickly became clear that in most cases the brighter maxima does not coincide with periastron (Dugan 1916) or that their orbits were essentially circular (O’Connell 1936), thus ruling out this hypothesis. It is now generally thought that the effect is likely due to a hot or cool spot (or spots) on one of the component stars, perhaps via accretion streams which would naturally lead to hot spots (Wilsey & Beaky 2009). However, introducing such spots into a model leads to significant degeneracy, making the drawing of reasonable conclusions nigh on impossible without employing Doppler tomography techniques (Bell et al. 1990). As such, we make no attempt to characterise the nature of PG 1034+001, which appears to display a pronounced O’Connell effect.
While a binary origin is the most plausible for the variability shown by PG 1034+001, it is important to note that spots on a single star could also reproduce the observed light curve, particularly given the recent observational evidence for the existence of long-lived spots on white dwarfs and similar objects (see e.g. Dupuis et al. 2000; Kilic et al. 2015).
3.5. NGC 5189
NGC 5189 exhibits a very complex morphology, with multiple filamentary and knotty structures distributed over the nebula (see e.g. Danehkar et al. 2018). Sabin et al. (2012) classified it as a quadrupolar PN and also reported the detection of a dense and cold infrared torus. The nucleus of NGC 5189 is a [WO1] Wolf-Rayet-type CS and was recently identified as a binary by Manick et al. (2015) with a period of 4.04 d. They also found other significant peaks at 3.54, 1.32, and 0.80 d, although they concluded that these are most likely to be aliases. In that paper, the authors proposed as the most plausible companion a massive white dwarf with a mass of m2 ≥ 0.5 M⊙.
In Fig. 4, we present the periodogram of the TESS data, which shows several peaks above the 0.1% FAP. Before computing the MCMC fitting, we visually inspected the most prominent peaks individually, finding several sources of periodicities. On the one hand, the strongest peak at 0.125 d shows a clear variability with ∼40 ppt amplitude (see Fig. A.1), which we were unable to properly fit with a simple sinusoidal function. We also tried the three-component model mentioned previously, but the results provided an unphysical solution for the amplitude of the beaming effect. The alias of this signal is also significant in the periodogram at 0.062 d and it does not persist in the residual Lomb-Scargle periodogram once the first periodicity is subtracted (with the three-component model). Neither of the periodicities identified result in a phase-folded light curve which can be easily associated with typical binary variability. However, in spite of the strong deviations from the sinusoidal, the light curve when folded on the strongest peak in the periodogram (see Fig. A.1) could feasibly be the result of irradiation in a high inclination system (which can result in strong asymmetries between minima and maxima).
The second most significant and interesting peak is 0.858 d, which shows a sinusoidal function when we phase-folded the light curve with such period. However, after inspecting the signal at 1.71 d (twice the period of 0.858 d), we realised that this was the true period since we have two minima per cycle. As the case of PG 1034+001, we identify three possible effects in the light curve of NGC 5189: reflection, ellipsoidal variability or Doppler beaming. In Fig. 4, we have plotted the light curve phase folded with this period and the final model accounting for the three effects. In the residual periodogram can be seen that the strong peak at 0.858 d is an alias. However, the 0.125 d period (and its alias) persist in the periodogram, which means that the nature of these two periodicities are different. As in the case of PG 1034+001, the amplitude of the Doppler beaming derived from the fitting of the light curve is too large and would imply a very massive companion in the system (see Sect. 3.4 for more details), which is somewhat difficult to explain. Again, the most plausible and likely explanation is that the asymmetry in the ellipsoidal modulation is due to a star spot and so, we are indeed detecting the O’Connell effect and not the Doppler beaming. We speculate that a plausible scenario would be a system with two white dwarfs, as proposed by Manick et al. (2015), but with different effective temperatures. However, the lack of information about the system does not allow us to do a more detailed model of the binary and radial velocity observations are required to constrain the parameters of the two components.
Surprisingly, we do not see any variability at the 4.04-day period previously published by Manick et al. (2015) to be the strongest in their radial velocity data. Among the other periodicities mentioned by the authors, we find photometric counterparts close to their 3.54 d and 0.80 d periods (the corresponding peaks in the TESS periodogram are located at 3. and 0.85, respectively, both are above the FAP 0.1%). Finally, we find no evidence for the 1.32 d periodicity, which has also been reported by Manick et al. (2015).
It is really important to consider that, as mentioned in Sect. 2.2, the aperture mask used to extract the photometry of NGC 5189 is subject to severe contamination by several sources of similar magnitudes as the central star (see Fig. 1). This makes it very difficult to discuss where both periodicities might come from. Another possibility for overcoming the apparent discrepancy between the radial velocity data from Manick et al. (2015) and the photometric data presented here is the presence of stellar winds. It is known that Wolf-Rayet-type CSPNe have strong stellar winds with mass-loss rate up to ∼10−7 Myr−1. Stellar wind variability may introduce some apparent velocity variations in the absorption line profiles (see e.g. De Marco et al. 2004). Therefore, a contribution of a variable stellar wind could be present in the radial velocity curve of NGC 5189 leading to the derivation of a (perhaps) spurious period where there is no sufficient data to separate the variability due to wind and orbital motion.
3.6. NGC 2867
This PN is catalogued as round or elliptical with multiple shells or haloes by Stanghellini et al. (1993). As the case of NGC 5189, NGC 2867 is a planetary nebula with a [WO1-2] spectral type as a central star. It was identified as a pulsating star by González Pérez et al. (2006), with a period of 769 s (∼12 min). For this object, we found two different light curves in the TESS archive, taken in sectors 9 and 10, respectively. We inspected the target pixel file (TPF) of both and checked the Tess Input Catalogue (TIC) name and the coordinates to be sure that both light curves correspond to the same object. For the analysis, we combined both data sets into a single light curve.
We detect a significant periodic variability at 4.5 days and amplitude of 0.4 ppt (see Fig. 4). However, as in the case of NGC 5189, the aperture mask includes several other stars, some of them with similar magnitude to our target, which makes it challenging to identify the origin of the variability. In this case, photometric ground-based follow-up, with better spatial resolution, of each target within the aperture is required in order to discern the source of the variability. We note that detections of photometric variations of 0.4 ppt in amplitude are feasible from ground-based telescopes.
3.7. Abell 7
Abell 7 is a very faint, slightly elongated PN (with a size of 15′×13′) and a clumpy structure inside the main shell (see deep image by D. Goldman4). The central star of this PN was classified as a possible physical binary by Ciardullo et al. (1999) after resolving a very faint companion (∼5.3 mag fainter than the central star) 0.91 arcsec away in the HST frames. De Marco et al. (2013) estimated an M5V spectral type for this wide companion based on the measurement of I- and J-band flux excess.
In the analysis of the TESS light curve, we identify two clear strong peaks at approximately 2.6 days and three days. The periodogram of the residuals (in blue in Fig. 5) indicates that both signals are significant since after removing one of the signals, the other still remains above the 0.1% FAP. Hence, both signals are independent and most probably the origin of each is different. In any case, the periods of the detected signals could not be compatible with an M star at 0.91 arcsec and, as consequence, both photometric periodicities would be associated with other companions (if they are indeed of binary origin).
3.8. Abell 15
This old PN was discovered by Abell (1966) on the Palomar Observatory Sky Survey (POSS) plates. The nebula has been classified as round or ring-shaped PN (Phillips 2003; Hromov & Kohoutek 1968). However, the lack of deep images or high-resolution, long-slit spectra of the nebula prevent us from doing a more detailed description of the morphology. Its central star is the only one in our sample that does not have significant peaks in the periodogram (see Fig. 6). The light curve does not show variability above ∼8 mmag at the 95% confidence level. This non-detection of variability is not surprising since Abell 15 is, apparently, the only round PN in our sample (only Abell 7 would be close to round, although it shows a slightly deviation from sphericity; see above). This result is compatible with the idea that round PN are formed from single stars.
4. Final remarks and conclusions
We analysed the TESS light curves available in the MAST archive of eight CSPNe in order to identify new binary candidates based on high-precision photometry. In all cases except one, we detect clear signs of variability, which means a positive detection rate of more than 80%. In most of the cases, these photometric variations can be attributed to the effect of a companion star or stars to the nebular progenitor. In seven out of the eight objects, the light curves show a variability that is consistent with an irradiation effect. Two of the targets (PG 1034+001 and NGC 5189) also exhibit ellipsoidal variability. We also discuss the possible detection of relativistic Doppler beaming effects and the O’Connell effect in several objects. In the case of NGC 7293 (Helix Nebula), the most well-studied PN in the sample, it was possible to construct a series of binary models to probe the possible companion properties which would lead to the observed variability in its light curve. These models lead us to conclude that if the variability is, indeed, due to binarity then the companion is likely a very low-mass main-sequence star or sub-stellar object.
The orbital periods found for the seven central stars range from 1.6 to 6.8 days, approximately. These periods would imply that the systems are post-common-envelope and while they are not exceptionally long, they do lie in the longer period tail of the observed post-common-envelope bCSPN period distribution.
In the case of the central star of NGC 246, an already known low- amplitude pulsator, the signals attributed to the pulsations are clearly identified in the periodogram. For the first time, all the frequency modes are resolved though we reserve a full analysis of these modes for a future paper (Sowicka et al., in prep.).
The only central star for which no significant variability is found is Abell 15. This is the only PN in the sample with an apparently round morphology, which is consistent with the theory that round PNe are formed via a single star evolution. Intriguingly, the other seven detections are PNe with more complex morphologies, which are believed to be typical of central star binarity.
Three central stars in the sample exhibit a possible contamination by other sources in the TESS aperture within a magnitude limit of Δm = 3. For these, ground-based photometry is necessary in order to discern the source of the variability. Also, radial velocity follow-up is essential in all cases to constrain the orbital parameters of the binary candidates found in this work.
See https://emcee.readthedocs.io/en/stable/ for further documentation.
Acknowledgments
We thank our anonymous referee for their useful comments that have improved the interpretation and discussion of the data. AA acknowledges support from Government of Comunidad Autónoma de Madrid (Spain) through postdoctoral Grant “Atracción de Talento Investigador” 2018-T2/TIC-11697. JLB and SBF acknowledge support from the Spanish State Research Agency (AEI) through projects No. ESP2017-87676-C5-1-R and No. MDM-2017-0737 Unidad de Excelencia “María de Maeztu” – Centro de Astrobiología (INTA-CSIC). DJ acknowledges support from the State Research Agency (AEI) of the Spanish Ministry of Science, Innovation and Universities (MCIU) and the European Regional Development Fund (FEDER) under grant AYA2017-83383-P. DJ also acknowledges support under grant P/308614 financed by funds transferred from the Spanish Ministry of Science, Innovation and Universities, charged to the General State Budgets and with funds transferred from the General Budgets of the Autonomous Community of the Canary Islands by the Ministry of Economy, Industry, Trade and Knowledge. LFM acknowledges partial support by grant AYA2017-84390-C2-1-R, co-funded with FEDER funds, and financial support from the State Agency for Research of the Spanish MCIU through the “Center of Excellence Severo Ochoa” award for the Instituto de Astrofísica de Andalucía (SEV-2017-0709). Authors are very grateful to Cristina Rodríguez López for guiding us with the analysis of RWT 152 and NGC 246. This research has made use of the SIMBAD database, operated at the CDS, Strasbourg (France), Aladin, NASA’s Astrophysics Data System Bibliographic Services, and the Spanish Virtual Observatory (http://svo.cab.inta-csic.es) supported from the Spanish MICINN/FEDER through grant AyA2017-84089. This publication makes use of data collected by TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. We acknowledge the use of pipelines at the TESS Science Office and at the TESS Science Processing Operations Center. This paper makes use of data products from the Mikulski Archive for Space Telescope.
References
- Abell, G. O. 1966, ApJ, 144, 259 [NASA ADS] [CrossRef] [Google Scholar]
- Adam, C., & Mugrauer, M. 2014, MNRAS, 444, 3459 [NASA ADS] [CrossRef] [Google Scholar]
- Aller, A., Miranda, L. F., Ulla, A., et al. 2013, A&A, 552, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Aller, A., Miranda, L. F., Olguín, L., et al. 2015, MNRAS, 446, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Aller, A., Miranda, L. F., Olguín, L., Solano, E., & Ulla, A. 2016, MNRAS, 462, 3945 [NASA ADS] [CrossRef] [Google Scholar]
- Barceló Forteza, S., Michel, E., Roca Cortés, T., & García, R. A. 2015, A&A, 579, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bell, S. A., Rainger, P. P., & Hilditch, R. W. 1990, MNRAS, 247, 632 [NASA ADS] [Google Scholar]
- Bloemen, S., Marsh, T. R., Degroote, P., et al. 2012, MNRAS, 422, 2600 [NASA ADS] [CrossRef] [Google Scholar]
- Boffin, H. M. J., & Jones, D. 2019, The Importance of Binaries in the Formation and Evolution of Planetary Nebulae [CrossRef] [Google Scholar]
- Bond, H. E. 1976, PASP, 88, 192 [NASA ADS] [CrossRef] [Google Scholar]
- Bond, H. E. 2000, ASP Conf. Ser., 199, 115 [NASA ADS] [Google Scholar]
- Ciardullo, R., & Bond, H. E. 1996, AJ, 111, 2332 [NASA ADS] [CrossRef] [Google Scholar]
- Ciardullo, R., Bond, H. E., Sipior, M. S., et al. 1999, AJ, 118, 488 [NASA ADS] [CrossRef] [Google Scholar]
- Danehkar, A., Karovska, M., Maksym, W. P., & Montez, Jr., R. 2018, ApJ, 852, 87 [NASA ADS] [CrossRef] [Google Scholar]
- De Marco, O., Bond, H. E., Harmer, D., & Fleming, A. J. 2004, ApJ, 602, L93 [NASA ADS] [CrossRef] [Google Scholar]
- De Marco, O., Passy, J.-C., Frew, D. J., Moe, M., & Jacoby, G. H. 2013, MNRAS, 428, 2118 [NASA ADS] [CrossRef] [Google Scholar]
- De Marco, O., Long, J., Jacoby, G. H., et al. 2015, MNRAS, 448, 3587 [NASA ADS] [CrossRef] [Google Scholar]
- Dugan, R. S. 1916, ApJ, 43, 130 [NASA ADS] [CrossRef] [Google Scholar]
- Dupuis, J., Chayer, P., Vennes, S., Christian, D. J., & Kruk, J. W. 2000, ApJ, 537, 977 [NASA ADS] [CrossRef] [Google Scholar]
- Eker, Z., Bakış, V., Bilir, S., et al. 2018, MNRAS, 479, 5491 [NASA ADS] [CrossRef] [Google Scholar]
- Foreman-Mackey, D., Hogg, D. W., Lang, D., & Goodman, J. 2013, PASP, 125, 306 [CrossRef] [Google Scholar]
- Frew, D. J., & Parker, Q. A. 2010, PASA, 27, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Gaia Collaboration (Brown, A. G. A., et al.) 2018, A&A, 616, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gianninas, A., Strickland, B. D., Kilic, M., & Bergeron, P. 2013, ApJ, 766, 3 [NASA ADS] [CrossRef] [Google Scholar]
- González Pérez, J. M., Solheim, J. E., & Kamben, R. 2006, A&A, 454, 527 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gruendl, R. A., Chu, Y.-H., O’Dwyer, I. J., & Guerrero, M. A. 2001, AJ, 122, 308 [NASA ADS] [CrossRef] [Google Scholar]
- Guerrero, M. A., Chu, Y.-H., Gruendl, R. A., Williams, R. M., & Kaler, J. B. 2001, ApJ, 553, L55 [NASA ADS] [CrossRef] [Google Scholar]
- Henry, R. B. C., Kwitter, K. B., & Dufour, R. J. 1999, ApJ, 517, 782 [NASA ADS] [CrossRef] [Google Scholar]
- Hewett, P. C., Irwin, M. J., Skillman, E. D., et al. 2003, ApJ, 599, L37 [NASA ADS] [CrossRef] [Google Scholar]
- Hillwig, T. C., Jones, D., De Marco, O., et al. 2016, ApJ, 832, 125 [NASA ADS] [CrossRef] [Google Scholar]
- Horvat, M., Conroy, K. E., Jones, D., & Prša, A. 2019, ApJS, 240, 36 [NASA ADS] [CrossRef] [Google Scholar]
- Hromov, G. S., & Kohoutek, L. 1968, Bull. Astron. Inst. Czech., 19, 1 [NASA ADS] [Google Scholar]
- Husser, T. O., Wende-von Berg, S., Dreizler, S., et al. 2013, A&A, 553, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Iben, Jr., I., & Livio, M. 1993, PASP, 105, 1373 [Google Scholar]
- Jenkins, J. M., Twicken, J. D., McCauliff, S., et al. 2016, SPIE Conf. Ser., 9913, 99133E [Google Scholar]
- Johnston, K. B., Haber, R., Caballero-Nieves, S. M., et al. 2019, Comput. Astrophys., 6, 4 [CrossRef] [Google Scholar]
- Jones, D., Conroy, K. E., Horvat, M., et al. 2019, ApJS, submitted [arXiv:1912.09474] [Google Scholar]
- Kilic, M., Gianninas, A., Bell, J., et al. 2015, ApJ, 814, 31 [NASA ADS] [CrossRef] [Google Scholar]
- Leahy, D. A., Zhang, C. Y., & Kwok, S. 1994, ApJ, 422, 205 [NASA ADS] [CrossRef] [Google Scholar]
- Manchado, A., Guerrero, M. A., Stanghellini, L., & Serra-Ricart, M. 1996a, The IAC Morphological Catalog of Northern Galactic Planetary Nebulae (La Laguna, Spain: Instituto de Astrofisica de Canarias (IAC)) [Google Scholar]
- Manchado, A., Stanghellini, L., & Guerrero, M. A. 1996b, ApJ, 466, L95 [NASA ADS] [CrossRef] [Google Scholar]
- Manick, R., Miszalski, B., & McBride, V. 2015, MNRAS, 448, 1789 [NASA ADS] [CrossRef] [Google Scholar]
- Mathur, S., García, R. A., Catala, C., et al. 2010, A&A, 518, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Meaburn, J., Clayton, C. A., Bryce, M., et al. 1998, MNRAS, 294, 201 [NASA ADS] [CrossRef] [Google Scholar]
- Miller Bertolami, M. M. 2016, A&A, 588, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Milone, E. E. 1968, AJ, 73, 708 [NASA ADS] [CrossRef] [Google Scholar]
- Miszalski, B., Acker, A., Moffat, A. F. J., Parker, Q. A., & Udalski, A. 2009, A&A, 496, 813 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Miszalski, B., Corradi, R. L. M., Boffin, H. M. J., et al. 2011, MNRAS, 413, 1264 [NASA ADS] [CrossRef] [Google Scholar]
- Muthu, C., Anandarao, B. G., & Pottasch, S. R. 2000, A&A, 355, 1098 [NASA ADS] [Google Scholar]
- Napiwotzki, R. 1999, A&A, 350, 101 [Google Scholar]
- O’Connell, D. J. K. 1936, Publ. Riverv. Coll. Obs., 3, 33 [NASA ADS] [Google Scholar]
- O’Connell, D. J. K. 1951, Publ. Riverv. Coll. Obs., 2, 85 [NASA ADS] [Google Scholar]
- O’dell, C. R. 1998, AJ, 116, 1346 [NASA ADS] [CrossRef] [Google Scholar]
- O’Dell, C. R., Henney, W. J., & Ferland, G. J. 2005, AJ, 130, 172 [NASA ADS] [CrossRef] [Google Scholar]
- Parker, Q. A., Acker, A., Frew, D. J., et al. 2006, MNRAS, 373, 79 [NASA ADS] [CrossRef] [Google Scholar]
- Phillips, J. P. 2003, MNRAS, 344, 501 [NASA ADS] [CrossRef] [Google Scholar]
- Prša, A., Conroy, K. E., Horvat, M., et al. 2016, ApJS, 227, 29 [NASA ADS] [CrossRef] [Google Scholar]
- Rauch, T., Kerber, F., & Pauli, E. M. 2004, A&A, 417, 647 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ricker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, J. Astron. Telesc. Instrum. Syst., 1, 014003 [NASA ADS] [CrossRef] [Google Scholar]
- Roberts, A. W. 1906, MNRAS, 66, 123 [NASA ADS] [CrossRef] [Google Scholar]
- Sabin, L., Vázquez, R., López, J. A., García-Díaz, M. T., & Ramos-Larios, G. 2012, Rev. Mex. Astron. Astrofis., 48, 165 [NASA ADS] [Google Scholar]
- Sahai, R., & Trauger, J. T. 1998, AJ, 116, 1357 [NASA ADS] [CrossRef] [Google Scholar]
- Stanghellini, L., Corradi, R. L. M., & Schwarz, H. E. 1993, A&A, 279, 521 [NASA ADS] [Google Scholar]
- Su, K. Y. L., Chu, Y. H., Rieke, G. H., et al. 2007, ApJ, 657, L41 [NASA ADS] [CrossRef] [Google Scholar]
- Szentgyorgyi, A., Raymond, J., Franco, J., Villaver, E., & López-Martín, L. 2003, ApJ, 594, 874 [NASA ADS] [CrossRef] [Google Scholar]
- Van de Steene, G. C., van Hoof, P. A. M., Exter, K. M., et al. 2015, A&A, 574, A134 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- van Kerkwijk, M. H., Rappaport, S. A., Breton, R. P., et al. 2010, ApJ, 715, 51 [NASA ADS] [CrossRef] [Google Scholar]
- Van Winckel, H., Jorissen, A., Exter, K., et al. 2014, A&A, 563, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Werner, K., Rauch, T., & Kruk, J. W. 2017, A&A, 601, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wilsey, N. J., & Beaky, M. M. 2009, Soc. Astron. Sci. Annu. Symp., 28, 107 [NASA ADS] [Google Scholar]
- Young, K., Cox, P., Huggins, P. J., Forveille, T., & Bachiller, R. 1999, ApJ, 522, 387 [NASA ADS] [CrossRef] [Google Scholar]
- Zucker, S., Mazeh, T., & Alexander, T. 2007, ApJ, 670, 1326 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Additional figure
![]() |
Fig. A.1. Phase-folded light curve (left panel) and periodogram (right panel) of NGC 5189. Two different bin sizes are shown with red and grey circles. |
All Tables
Results from the light curve study for the five central stars analyzed with a simple sinusoidal model.
Results from the light curve study for the two central stars analyzed with models including reflection, ellipsoidal and doppler beaming effects.
All Figures
![]() |
Fig. 1. Target pixel files (TPFs) of each central star in the sample. The red circles are the sources of the Gaia DR2 catalogue in the field with scaled magnitudes (see legend). White crosses indicate the location of the targets. The aperture mask used by the pipeline to extract the photometry is also marked. Pixel scale is 21 arcsec pixel−1. These plots have been created with the tpfplotter code, publicly available at https://github.com/jlillo/tpfplotter. |
In the text |
![]() |
Fig. 2. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 7293 and RWT 152. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3. |
In the text |
![]() |
Fig. 3. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 246 and PG 1034+001. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3 (for NGC 246) and Table 4 (for PG 1034+001). |
In the text |
![]() |
Fig. 4. Time-flux light curves (top panels), phase-folded light curves (left panels) and periodograms (right panels) of NGC 5189 and NGC 2867. Two different bin sizes are shown with red and grey circles (see legend for each individual object). The periods derived from the corresponding fittings can be found in Table 3 (for NGC 2867) and Table 4 (for NGC 5189). |
In the text |
![]() |
Fig. 5. Time-flux light curve (top panel), phase-folded light curves (left panels) and periodograms (right panels) of Abell 7. Two different bin sizes are shown with red and grey circles. The periods derived from the corresponding fittings can be found in Table 3. |
In the text |
![]() |
Fig. 6. Time-flux light curve (top panels), phase-folded light curve (left panel) and periodogram (right panel) of Abell 15. Two different bin sizes are shown with red and grey circles. |
In the text |
![]() |
Fig. 7. Periodogram of NGC 246 in the low frequencies. The different peaks are due to the pulsations. |
In the text |
![]() |
Fig. A.1. Phase-folded light curve (left panel) and periodogram (right panel) of NGC 5189. Two different bin sizes are shown with red and grey circles. |
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.