Issue |
A&A
Volume 567, July 2014
|
|
---|---|---|
Article Number | A30 | |
Number of page(s) | 12 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201423800 | |
Published online | 08 July 2014 |
IP Eridani: A surprising long-period binary system hosting a He white dwarf⋆,⋆⋆
1
Institut d’Astronomie et d’Astrophysique, Université Libre de
Bruxelles,
CP 226, boulevard du Triomphe,
1050
Brussels,
Belgium
e-mail:
tmerle@ulb.ac.be
2
Instituut voor Sterrenkunde, Katholieke Universiteit Leuven,
Celestijnenlaan
200D, 3200
Heverlee,
Belgium
Received: 13 March 2014
Accepted: 15 May 2014
We determine the orbital elements for the K0 IV + white dwarf (WD) system IP Eri, which appears to have a surprisingly long period of 1071 d and a significant eccentricity of 0.25. Previous spectroscopic analyses of the WD, based on a distance of 101 pc inferred from its Hipparcos parallax, yielded a mass of only 0.43 M⊙, implying it to be a helium-core WD. The orbital properties of IP Eri are similar to those of the newly discovered long-period subdwarf B star (sdB) binaries, which involve stars with He-burning cores surrounded by extremely thin H envelopes, and are therefore close relatives to He WDs. We performed a spectroscopic analysis of high-resolution spectra from the HERMES/Mercator spectrograph and concluded that the atmospheric parameters of the K0 component are Teff = 4960 K, log g = 3.3, [Fe/H] = 0.09 and ξ = 1.5 km s-1. The detailed abundance analysis focuses on C, N, O abundances, carbon isotopic ratio, light (Na, Mg, Al, Si, Ca, Ti) and s-process (Sr, Y, Zr, Ba, La, Ce, Nd) elements. We conclude that IP Eri abundances agree with those of normal field stars of the same metallicity. The long period and non-null eccentricity indicate that this system cannot be the end product of a common-envelope phase; it calls instead for another less catastrophic binary-evolution channel presented in detail in a companion paper.
Key words: stars: abundances / binaries: spectroscopic / stars: evolution / stars: individual: IP Eridani / subdwarfs / white dwarfs
Appendix A is available in electronic form at http://www.aanda.org
HERMES 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/qcat?J/A+A/567/A30
© ESO, 2014
1. Introduction
IP Eri is a very interesting system consisting of a K0 (sub)giant and a He white dwarf (WD) and it is tempting to relate this system to the family of subdwarf B (sdB) binaries where the hot component is a bare He-burning core surrounded by an extremely thin H envelope (Heber 2009). These systems owe their properties to envelope ejection, likely due to binary interaction, as they evolve along the red giant branch. Their evolution shares some similarities with that of IP Eri and such systems are important benchmarks for binary evolution.
Several long-period eccentric systems (with P ~ 103 d) were recently discovered among sdB stars (Østensen & Van Winckel 2011, 2012; Vos et al. 2012, 2013; Deca et al. 2012; Barlow et al. 2012, 2013). Using new radial-velocity data collected with the HERMES/Mercator spectrograph (Raskin et al. 2011), we show in this paper (Sect. 3) that IP Eri adds to this new class of long-period eccentric systems.
Since the He-WD progenitor did not evolve along the asymptotic giant branch (AGB), it had no chance to produce s-process elements and to pollute its companion (the present K0 subgiant), so that the latter should not appear as a barium star (a family of K giants with enhanced abundances of s-process elements; Bidelman & Keenan 1951). It is therefore of interest to perform a chemical analysis of the K0 subgiant in the IP Eri system to confirm the absence of overabundances of s-process elements. This is the second objective of the present paper, which is organized as follows: Sect. 2 gives an overview of the properties of the IP Eri system. Sect. 3 presents the radial-velocity data and the ensuing orbit. After deriving the atmospheric parameters of IP Eri in Sect. 4, the abundance analysis is presented in Sect. 5, with emphasis on s-process elements. These abundances are then compared with expectations for barium stars and for non-s-process-polluted stars, as derived from the abundance trends observed in large samples of field stars reflecting the chemical evolution of the Galaxy (Sect. 6). Section 7 briefly confronts our results with predictions from binary-evolution scenarios presented in detail in a companion paper (Siess et al. 2014). Section 8 summarizes our results.
2. The IP Eri system
IP Eri (HD 18131, HIP 13558, WD 0252-055, EUVE J0254-053) has attracted attention since its discovery as an extreme UV (EUV) source both by ROSAT (Pounds et al. 1993; Pye et al. 1995) and EUVE (Bowyer et al. 1994, 1996; Malina et al. 1994). In a subsequent analysis, Vennes et al. (1995) find that an older International Ultraviolet Explorer (IUE) spectrum already revealed that a WD dominates the IUE spectrum below 200 nm while a K0 spectrum prevails at longer wavelengths. Their model-atmosphere analysis of the EUV photometry reveals a hot, hydrogen-rich (DA) WD (with an effective temperature of about 30 000 K) that is the most likely source of the EUV emission. In the grid of models fitted to the WD spectrum by Burleigh et al. (1997), the model with a temperature of 29 290 K, a gravity of log g = 7.5 and a mass of 0.43 M⊙ locates the WD at a distance of 93 pc, consistent with the Hipparcos parallax of the K star, as discussed below. The corresponding age for the WD is then 7 Myr (Burleigh et al. 1997). The WD parameters inferred by Burleigh et al. (1997) locate it among He WDs in the log g − log Teff diagram of Driebe et al. (1998). Its mass of 0.43 M⊙ is incompatible with a CO WD since the minimum CO core mass at the base of the AGB (just at the end of core He-burning) is 0.51 M⊙ (for a star of initial mass 0.9 M⊙; e.g., Eq. (66) of Hurley et al. 2000). An independent study by Vennes et al. (1998) concludes at a somewhat higher mass (0.48−0.52 M⊙) for the IP Eri WD, which is still, for the most part of this range, compatible with the He nature of the WD.
![]() |
Fig. 1 The Ca ii H and K lines, from the HERMES spectrum on HJD 2 455 066.748 (2009, August 23). The emission present in the cores is somewhat more intense than the one presented in Fig. 5 of Vennes et al. (1997). |
The revised Hipparcos parallax of 9.82 ± 0.94 mas corresponds to a distance of 101 ± 11 pc and a distance modulus of 5.03 ± 0.2 (van Leeuwen 2007). With V = 7.32 (Cutispoto et al. 1995), one gets an absolute magnitude of 2.29 for the K star, which indicates that it is a subgiant. The photometric data for the K0 subgiant are V = 7.32, U − B = 0.74, B − V = 0.98, V − Rc = 0.52, V − Ic = 0.97 (Cutispoto et al. 1995) and J = 5.709,H = 5.263, K = 5.090 (2MASS; Cutri et al. 2003). The corresponding V − K index of 2.23 implies an effective temperature of 4900 K (Bessell et al. 1998). Using a Bayesian method, Bailer-Jones (2011) obtains log Teff = 3.70 ± 0.01 (5012 K) for IP Eri.
Although neither the ASAS lightcurve (Pojmanski 1997) nor the monitoring performed by Cutispoto et al. (1999) reveals variability at the 0.1 mag level, a more accurate monitoring by Strassmeier et al. (2000) has uncovered a 0.045 mag variability over 30 d in the Strömgren y band, hence its classification as a BY Dra variable in the General Catalogue of Variable Stars (with name IP Eri). The star is slowly rotating (V sin i< 5 km s-1; Cutispoto et al. 1999). Should the 30 d photometric variability be due to rotation, a value of ~1 km s-1 for the rotational velocity would result from the radius of 3.8 R⊙ derived in (Siess et al. 2014) for IP Eri. Nevertheless, this star exhibits moderate Mg II h and k emission in the IUE spectrum and Ca ii H & K emission (Fig. 5 of Vennes et al. 1997; Strassmeier et al. 2000). Figure 1 shows the emission cores in the Ca ii lines, as seen on our HERMES/Mercator spectrum obtained on HJD 2 455 066.748 (2009, August 23). This emission core is somewhat more intense than the one presented by Vennes et al. (1997).
3. Orbital elements
The 18 high-resolution spectra used to compute the spectroscopic orbital elements were obtained with the HERMES/Mercator spectrograph (Raskin et al. 2011), operating at an average resolution of 85 000 in high-resolution mode and with a spectral range of [λ370−900 nm]. The spectra were reduced with the HERMES pipeline, and the radial velocities, computed by cross-correlating the observed spectra with an Arcturus template, are on the IAU wavelength system defined by the standards from Udry et al. (1999a,b). The individual radial velocities Vr are listed in Table 1. The errors are dominated by the drift of the air refractive index caused by the atmopsheric pressure variations in the spectrograph room (see Fig. 9 in Raskin et al. 2011). The long-term stability during the 4 years of operations of the HERMES spectrograph turns out to be ~40 m s-1 as derived from the standard deviation of the radial velocities of the monitored IAU standards. This may be considered as the precision on the radial velocities produced by HERMES.
The corresponding orbital solution is listed in Table 2 and displayed in Fig. 2. The uncertainties of about 30 m s-1 on the orbital parameters Vγ and K1 are consistent with the precision of the spectrograph. Our orbital solution does not include the radial-velocity measurement Vr = 14.94 ± 0.10 km s-1 obtained by Chubak et al. (2012) on JD 2 455 261, but we checked a posteriori that it falls on the orbital solution within 0.1 km s-1. It is worth mentioning that IP Eri was imaged using the HST by Barstow et al. (2001) and the binary was not resolved, setting an upper limit on the orbital period of about 19 yr, consistent with the 2.9-yr period found here. No orbit could be adjusted to the O−C residuals of the IP Eri system (see lower panel of Fig. 2) when testing for the presence of a hypothetical third component.
The current mass function of the system M⊙ constrains the primary mass M1 to be lower than 4.27 M⊙, if M2 = 0.43M⊙ for the He-WD (see Sect. 2).
Radial velocities used for computing the orbit of the IP Eri system.
Spectroscopic orbital elements of IP Eri.
4. Atmospheric parameters
To derive the atmospheric parameters and the detailed abundances, we used two spectra with high signal-to-noise ratios from the list of Table 1 (namely those obtained on August 22 and September 5, 2009, respectively HJD 2 455 066.748 and 2 455 080.698). The first guesses for the atmospheric parameters were determined from the photometry. The color indices J − K = 0.62 and V − K = 2.23 yield a first estimate of Teff= 4900 K for the effective temperature, using the calibrations of Bessell et al. (1998), in perfect agreement with previous estimates. From the calibration of MK spectral types provided by Cox (2000), we adopted log g = 2.1 as initial guess for the surface gravity of a giant star of spectral type K0 (despite the fact that the luminosity class suggests a higher gravity). We initially assumed the metallicity to be solar.
The atmospheric parameters were then determined iteratively using the BACCHUS pipeline developed by one of the author (TMa; see also Jofre et al. 2014) in the context of the Gaia-ESO survey (Gilmore et al. 2012). This pipeline is based on the 1D LTE spectrum-synthesis code Turbospectrum (Alvarez & Plez 1998; Plez 2012) and allows an automated determination of effective temperature Teff, surface gravity log g, metallicity [Fe/H] and microturbulent velocity ξ. We used MARCS model atmospheres (Gustafsson et al. 2008) along with a selection of neutral and singly ionized Fe lines that have been selected for the analysis of stellar spectra in the framework of the Gaia-ESO survey. Oscillator strengths are from the VALD database (Kupka et al. 2000). The classical method to obtain the atmospheric parameters consists in avoiding trends in the [Fe/H] vs. χ and [Fe/H] vs. W/λ relations (where χ is the lower excitation energy of the considered line, W its measured equivalent width and λ its wavelength) but also in forcing lines of Fe i and Fe ii to yield the same abundance (see Fig. 3).
The equivalent widths are automatically measured using spectrum synthesis with the atmospheric parameters determined at the previous iteration. The synthetic spectra were convolved with a Gaussian function with full width at half maximum of 6.5 km s-1. Only iron lines having reduced equivalent widths (W/λ) lower than 0.025 mÅ/Å were kept in the analysis. We thus obtain the following atmospheric parameters: Teff= 4960 ± 100 K, log g = 3.3 ± 0.3, [Fe/H] = + 0.09 ± 0.08 and ξ = 1.5 ± 0.1 km s-1. We used the solar reference values from Grevesse et al. (2007) where A⊙(Fe) = 7.45.
![]() |
Fig. 2 The orbital solution for IP Eri. The lower panel shows the observed (O) minus the calculated (C) orbits. |
![]() |
Fig. 3 Line by line abundance analysis for Fe i and Fe ii lines as a function of the excitation potential χ and the reduced equivalent width W/λ. The grey areas represent the standard deviations around the mean abundances of Fe i and Fe ii. |
5. Abundances
The detailed abundance analysis was performed (in the framework of Local Thermodynamical Equilibrium – LTE) using the abundance module of the BACCHUS pipeline. The selection of atomic and molecular lines was performed over the whole wavelength range covered by the HERMES/Mercator spectrograph. The atomic line list used for the detailed analysis is given in Appendix A. It includes the isotopic shifts for Ba ii (with an update for isotopes 130 and 132) and the hyperfine structure for La ii from Masseron (2006). The CH molecular line list is from Masseron et al. (2014). The references for the other molecular line lists (TiO, SiO, VO, C2, CN, NH, OH, MgH, SiH, CaH and FeH) can be found in Gustafsson et al. (2008). Line fitting is essentially based on a least-square minimization method and all lines are visually inspected to check for possible bad fits (due to, e.g., line blends, cosmic hits, ...). The results of the detailed abundance analysis are presented in Table 3.
5.1. C, N, O and 12C / 13C
The numerous neutral carbon atomic lines lead to an abundance of [ C / Fe ] = 0.08 ± 0.16 dex. The nitrogen abundance is determined from numerous molecular CN lines selected over a large spectral range [λ640−890 nm]. An example of the fit of CN lines in the [λ812.2−813.2 nm] region is shown in Fig. 4: synthetic spectra with [ N / Fe ] = −0.12 ± 0.3 dex (see Table 3) are compared with the observed spectrum of IP Eri. The oxygen abundance is derived from the sole [O i] λ630.03 nm line which is supposed to be free from NLTE effects (Asplund 2005). Another forbidden line at λ636.38 nm is in the red wing of a Ca i autoionization line which is difficult to fit in our spectrum. The oxygen triplet at λ777.19, 777.42 and 777.54 nm gives an abundance A(O) ≈ 9.46 ± 0.07 dex, in disagreement by about 0.7 dex with respect to the λ630.03 nm line. This discrepancy is mainly due to the NLTE effect of diffusion in the triplet as clearly explained by Asplund (2005). No OH lines are available in the HERMES spectrum to better constrain the oxygen abundance. The resulting C/O ratio is 0.55, in agreement with the solar value.
Results of the chemical abundance analysis.
The carbon isotopic ratio 12C / 13C is deduced from the 12CN and 13CN molecular lines in the [λ799.5–801.5 nm] range (see, e.g., Barbuy et al. 1992; Drake & Pereira 2008). Specifically, we can compare the 12CN triplet between λ800.3 and λ800.4 nm with the 13CN feature at λ800.45 nm, and the 12CN weak line at λ801.0 nm with the 13CN very weak line at λ801.05 nm. The latter is in the far red wing of the former. The high resolution of the HERMES spectrograph is able to separate the two components. The best fit of the entire [λ799.5−801.5 nm] spectral range gives a carbon abundance of A(C) = 8.56 ± 0.10 dex, which confirms the abundance derived from atomic carbon lines. With this value, we can only deduce a lower limit for the carbon isotopic ratio of 12C / 13C ≥ 20, in accordance with the weakness of the 13CN features.
![]() |
Fig. 4 Example of CN-line fitting. |
![]() |
Fig. 5 Example of Ba ii line fitting (Ba ii λ585.37 nm). |
![]() |
Fig. 6 Example of Ce ii line fitting (Ce ii λ533.01 nm). |
5.2. Light elements
Neutral lines from Na, Mg, Al, Si and Ca provide reliable abundances with a standard deviation lower than, or of the order of, 0.1 dex. Ca is the most enriched among the investigated α-elements, and has the lowest dispersion. The Ti abundance derived from Ti i lines is consistent with that of Mg and Ca, whereas the Ti abundance derived from Ti ii is 0.1 dex lower but still within the statistical uncertainties. IP Eri is slightly enriched in α-elements ([α/Fe] = 0.17 ± 0.06 when considering Mg, Ca and Ti). The iron abundance determination is illustrated in Fig. 3 and is consistent with the derived metallicity.
5.3. s-process elements
Abundances for elements from the first two s-process peaks are measurable in IP Eri and their values are listed in Table 3. Examples of synthetic spectra are shown around two second-peak s-process elements (see Figs. 5 and 6) with solar ([Ba,Ce/Fe] = 0.00), enhanced ([Ba,Ce/Fe] = 1.00) and actual abundances ([Ba/Fe] = 0.18 and [Ce/Fe] = −0.22). For elements from the first s-process peak, lines of neutral and singly ionized Y and Zr are available. Abundances from neutral and ionized species agree within the statistical uncertainties, excepted for Y, with abundances from neutral lines being somewhat larger than those derived from ionized lines. This trend, if real, could be explained by the NLTE mechanism of overionization of the dominant neutral species due to UV radiation of non-local origin as shown, e.g., for Mg by Merle et al. (2011) or for Fe by Lind et al. (2012). These NLTE effects altering the ionization equilibrium of s-process elements should be investigated, but such an analysis is beyond the scope of this paper.
An average light-s-process abundance of [ls/Fe] = 0.06 ± 0.04 is obtained, based on Sr i, Y ii, and Zr ii abundances, as compared to [hs/Fe] = −0.03 ± 0.05, based on Ba ii, La ii, and Ce ii abundances. The subsolar Ce ii abundance is supported by a similarly subsolar Nd abundance, as derived from Nd ii lines. Thus, there does not seem to be a significant s-process enrichment in IP Eri. This issue is discussed further in Sect. 6, where the IP Eri abundances are confronted with abundances in barium stars and with expectations from the chemical evolution of the Galaxy.
![]() |
Fig. 7 Comparison of light-element abundances in IP Eri (black dot), in field stars (Edvardsson et al. 1993, grey plusses), and in mild barium stars (Rojas et al. 2013, black crosses). |
6. Confrontation with normal field stars and barium stars
We have compared the abundances of IP Eri with those of field stars and with those from a sample of mild barium stars, to confirm that the abundances of s-process elements in IP Eri derived in Sect. 5.3 are by no means peculiar.
We first consider the light-element abundances [X/Fe] in IP Eri which are compared in Fig. 7 with abundances in a large sample of F and G stars from the galactic disc (Edvardsson et al. 1993). We corrected for the zero-point solar abundance offset, since Edvardsson et al. (1993) used a different reference value of A⊙(Fe) = 7.51. Figure 7 shows that the light-element abundances in IP Eri are in relatively good agreement with those of moderately metal-rich stars. The Na i, Mg i, Al i and Ti i abundances of IP Eri are located within the abundance distribution of the galactic sample. O i (not shown in Fig. 7) and Ca i are slightly more abundant as compared to the disc stars. The oxygen overabundance may be trusted though, since Edvardsson et al. (1993) used a scaling relation to transform abundances from the high excitation lines that they used (λ615.8 nm as well as the λ777.3 nm triplet) to the abundance from the [O i] λ630.0 nm line that we used. For Ca, there are no lines in common with the Edvardsson et al. analysis. We used saturated and strong lines which are affected by NLTE effects, as shown by Mashonkina et al. (2007). These authors predict positive NLTE abundance corrections (between 0 and 0.07 dex for a model with Teff= 5000 K, log g = 3 and [Fe/H] = 0), which, if accounted for, would further strengthen the Ca enrichment as compared to the Edvardsson et al. (1993) values for field stars. We have also compared IP Eri light-element abundances with those of a sample of mild barium stars from Rojas et al. (2013). Their light-element abundances are similar to those of field stars.
On the contrary, the Si i abundance of IP Eri (for which most of the lines used are common to the two studies) is lower than the value derived for the disc-star distribution. The IP Eri abundance of Ti i matches the average abundance of the disc stars of similar metallicities, but our value suffers from large error bars due to a larger set of lines in our analysis (9 against 4 for Edvardsson et al. 1993). This large dispersion may be partially explained by the fact that all the lines have equivalent widths larger than 150 mÅ and are consequently less sensitive to abundance changes.
The s-process element abundances in IP Eri are compared with those of disc stars in Fig. 8. The comparison sample is from Edvardsson et al. (1993) for Y ii, Reddy et al. (2003) for Sr i and Ce ii, Pompéia et al. (2011) for Zr ii, Ba ii and La ii. The slight enrichment in Zr ii and Ba ii observed in IP Eri and is typical of the enrichment trend observed for disc stars. Y ii and Ce ii appear peculiar in that they are underabundant in IP Eri as compared to disc stars of similar metallicities.
![]() |
Fig. 8 Comparison of s-process element abundances of IP Eri with field stars (Edvardsson et al. 1993, grey plusses; Reddy et al. 2003, grey squares; Pompéia et al. 2011, grey asterisks) and mild barium stars (Rojas et al. 2013, black crosses). |
For Sr i, the only abundances available for comparison are from Reddy et al. (2003). They are based on a single line, and unfortunately, the Reddy et al. sample includes only stars with solar and sub-solar metallicities, giving the false impression that the IP Eri abundances are discrepant. For La ii, data from Pompéia et al. (2011) are the only ones available for comparison, and the slight La ii underabundance (with respect to the Sun) observed in IP Eri matches the trend observed among disc stars.
The comparison with the sample of mild barium stars of Rojas et al. (2013) is not discriminating for first-peak s-process elements (top panels of Fig. 8). But abundances from second-peak s-process elements in IP Eri are clearly lower than those of mild barium stars (bottom panels of Fig. 8). From these comparisons, we conclude that the slight enrichment in α and s-process elements in IP Eri is largely consistent with the chemical evolution of the Galaxy, so that there is no obvious signature from a chemical pollution resulting from mass transfer.
7. The evolutionary context of IP Eri
![]() |
Fig. 9 Evolutionary channels for the formation of a He WD. The dashed lines refer to channels where the eccentricity can be preserved (see text for details). |
He WDs form when a star loses its hydrogen-rich envelope before it ignites helium. For a single star, this is not possible within a Hubble time since only stars with masses M ≲ 0.45−0.5M⊙ can avoid helium ignition. A binary scenario is therefore required.
As illustrated in Fig. 9, different evolutionary channels can account for these objects. The first one involves mass transfer by Roche lobe overflow (RLOF). The long period of IP Eri (1071 d) imposes that mass transfer starts while the star is already on the red-giant branch (late case B). Because of the presence of a deep convective envelope in the donor, two outcomes are possible depending on the mass ratio q = Mdonor/Mgainer. If q> 1.3−1.5 (Soberman et al. 1997; Hurley et al. 2002), the mass transfer is dynamically unstable; after a rapid stage of common-envelope evolution, a short-period system forms (channel Ia in Fig. 9). In the alternative configuration (q< 1.3−1.5), soon after RLOF starts, the mass ratio reverses and subsequent mass transfer leads to the expansion of the orbit. The outcome is then a long-period system similar to IP Eri (channel Ib). However, in this RLOF scenario, tidal interactions are very strong because of the extended convective envelope of the Roche-filling donor star and the orbit always circularizes. This channel thus cannot explain the high eccentricity of IP Eri.
A solution to the eccentricity problem has been described in a companion paper (Siess et al. 2014). Based on binary-evolution calculations with the code BINSTAR (Siess et al. 2013; Davis et al. 2013; Deschamps et al. 2013), we showed that if the envelope of the He-WD progenitor is lost via tidally-enhanced winds (channel II), the circularization can be avoided. Tout & Eggleton (1988) suggested that the presence of a companion star can substantially increase the mass-loss rate of the evolved component. In this situation, the donor star loses its envelope while remaining inside its Roche potential and tidal forces are significantly reduced. Moreover, if the system has an initial eccentricity, the orbital wind mass-transfer modulation (Soker 2000) provides an eccentricity-pumping mechanism that counteracts the tidal circularization. We showed that such a scenario is able to account for all the orbital properties of a system like IP Eri.
![]() |
Fig. 10 The positions of IP Eri and its twin system HR 1608 (black circled crosses; Beavers & Eitter 1988, quoted by Landsman et al. 1993) in the eccentricity – period diagram are compared with those of the long-period sdB binaries, with (mostly) pre-mass-transfer binaries (normal K and M giants), and with post-mass-transfer binaries (Barium and S stars). Symbols are as indicated in the figure label. |
In Fig. 10 we compare the orbital properties of IP Eri with those of pre- (K giants in open clusters from Mermilliod et al. 2007 and binary M giants from Famaey et al. 2009; Jorissen et al. 2009) and post-mass-transfer binaries (barium and S stars from Jorissen et al. 1998, updated with some recently published orbits from Gorlova et al. 2014). We also show in the e-log P diagram the location of the G5 IV + WD system HR 1608 (=63 Eri). This system has an eccentricity (e = 0.30 ± 0.06) and a period (P = 903 ± 5 d) very similar to those of IP Eri (from Beavers & Eitter 1988,quoted by Landsman et al. 1993). The revised Hipparcos parallax (18.53 ± 0.84 mas; van Leeuwen 2007) yields a distance range 51−56 pc for HR 1608, which implies a WD mass around 0.4 M⊙ (Landsman et al. 1993), and thus a He-WD. A more recent analysis of the WD parameters by Vennes et al. (1998) suggests instead a higher mass range (0.51−0.67 M⊙). We also include in Fig. 10 the sdB binaries with long periods (black crosses) from Østensen & Van Winckel (2011, 2012), Vos et al. (2012, 2013), Deca et al. (2012), and Barlow et al. (2012, 2013).
IP Eri and HR 1608 have the largest eccentricities when compared to the long-period sdB systems, a fact whose significance is difficult to assess with so few He-WD systems. What seems significant, however, is the fact that (long-period) sdB and He-WD systems occupy a rather restricted period range around 103 d. In that respect, they differ from the barium and S binaries, two other families of post-mass-transfer systems, which spread over a more extended period range. This difference might be related to the fact that sdB and He-WD binaries are the end products of mass transfer occurring on the first red giant branch, whereas barium and S systems involved mass transfer on the AGB.
8. Summary
IP Eri is an unusual long-period binary system with a high eccentricity. Our abundance analysis reveals that it is not enriched in s-process elements, giving additional support to He-WD nature of the hot component that avoided evolution along the AGB. The abundances of the light elements do not reveal any peculiarities and within the error bars, the abundances are very close to solar. What is more surprising is that IP Eri alike its twin HR 1608 have the largest eccentricities among their closely related sdB systems and are very concentrated around a thousand-day period, which is likely related to the fact that the donor stars lose their envelope on the red-giant branch. However, the statistics remains small and additional data are highly desirable in order to confirm the formation channel of these systems.
Online material
Appendix A: Linelist
Line list used for determining atmospheric parameters and chemical composition of atomic species.
Acknowledgments
This work has been partly funded by an Action de recherche concertée (ARC) from the Direction générale de l’Enseignement non obligatoire et de la Recherche scientifique – Direction de la recherche scientifique – Communauté française de Belgique. T.M. is supported by the FNRS-F.R.S. as temporary post-doctoral researcher under grant No. 2.4513.11. The Mercator telescope is operated thanks to grant number G.0C31.13 of the FWO under the Big Science initiative of the Flemish governement. Based on observations obtained with the HERMES spectrograph, supported by the Fund for Scientific Research of Flanders (FWO), the Research Council of K.U. Leuven, the Fonds National de la Recherche Scientifique (F.R.S.-FNRS), Belgium, the Royal Observatory of Belgium, the Observatoire de Genève, Switzerland and the Thüringer Landessternwarte Tautenburg, Germany.
References
- Alvarez, R., & Plez, B. 1998, A&A, 330, 1109 [NASA ADS] [Google Scholar]
- Asplund, M. 2005, ARA&A, 43, 481 [NASA ADS] [CrossRef] [Google Scholar]
- Bailer-Jones, C. A. L. 2011, MNRAS, 411, 435 [NASA ADS] [CrossRef] [Google Scholar]
- Barbuy, B., Jorissen, A., Rossi, S. C. F., & Arnould, M. 1992, A&A, 262, 216 [NASA ADS] [Google Scholar]
- Barlow, B. N., Wade, R. A., Liss, S. E., Østensen, R. H., & Van Winckel, H. 2012, ApJ, 758, 58 [NASA ADS] [CrossRef] [Google Scholar]
- Barlow, B. N., Liss, S. E., Wade, R. A., & Green, E. M. 2013, ApJ, 771, 23 [NASA ADS] [CrossRef] [Google Scholar]
- Barstow, M. A., Bond, H. E., Burleigh, M. R., & Holberg, J. B. 2001, MNRAS, 322, 891 [NASA ADS] [CrossRef] [Google Scholar]
- Beavers, W. I., & Eitter, J. J. 1988, in BAAS, 20, 737 [Google Scholar]
- Bessell, M. S., Castelli, F., & Plez, B. 1998, A&A, 333, 231 [NASA ADS] [Google Scholar]
- Bidelman, W. P., & Keenan, P. C. 1951, ApJ, 114, 473 [NASA ADS] [CrossRef] [Google Scholar]
- Bowyer, S., Lieu, R., Lampton, M., et al. 1994, ApJS, 93, 569 [NASA ADS] [CrossRef] [Google Scholar]
- Bowyer, S., Lampton, M., Lewis, J., et al. 1996, ApJS, 102, 129 [NASA ADS] [CrossRef] [Google Scholar]
- Burleigh, M. R., Barstow, M. A., & Fleming, T. A. 1997, MNRAS, 287, 381 [NASA ADS] [CrossRef] [Google Scholar]
- Chubak, C., Marcy, G., Fischer, D. A., et al. 2012 [arXiv:1207.6212] [Google Scholar]
- Cox, A. N. 2000, Allen’s astrophysical quantities (Springer) [Google Scholar]
- Cutispoto, G., Pallavicini, R., Kuerster, M., & Rodono, M. 1995, A&A, 297, 764 [NASA ADS] [Google Scholar]
- Cutispoto, G., Pastori, L., Tagliaferri, G., Messina, S., & Pallavicini, R. 1999, A&AS, 138, 87 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, VizieR Online Data Catalog: II/246 [Google Scholar]
- Davis, P. J., Siess, L., & Deschamps, R. 2013, A&A, 556, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Deca, J., Marsh, T. R., Østensen, R. H., et al. 2012, MNRAS, 421, 2798 [NASA ADS] [CrossRef] [Google Scholar]
- Deschamps, R., Siess, L., Davis, P. J., & Jorissen, A. 2013, A&A, 557, A40 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Drake, N. A., & Pereira, C. B. 2008, AJ, 135, 1070 [NASA ADS] [CrossRef] [Google Scholar]
- Driebe, T., Schoenberner, D., Bloecker, T., & Herwig, F. 1998, A&A, 339, 123 [Google Scholar]
- Edvardsson, B., Andersen, J., Gustafsson, B., et al. 1993, A&A, 275, 101 [NASA ADS] [Google Scholar]
- Famaey, B., Pourbaix, D., Frankowski, A., et al. 2009, A&A, 498, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gilmore, G., Randich, S., Asplund, M., et al. 2012, The Messenger, 147, 25 [NASA ADS] [Google Scholar]
- Gorlova, N., Van Winckel, H., Vos, J., et al. 2014, Setting a new standard in the analysis of binary stars, ed. A. Tkachenko, EAS Conf. Ser. [arXiv:1403.2287] [Google Scholar]
- Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Sci. Rev., 130, 105 [Google Scholar]
- Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Heber, U. 2009, ARA&A, 47, 211 [NASA ADS] [CrossRef] [Google Scholar]
- Hurley, J. R., Pols, O. R., & Tout, C. A. 2000, MNRAS, 315, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Hurley, J. R., Tout, C. A., & Pols, O. R. 2002, MNRAS, 329, 897 [NASA ADS] [CrossRef] [Google Scholar]
- Jofre, P., Heiter, U., Soubiran, C., et al. 2014, A&A, 564, A133 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Jorissen, A., Van Eck, S., Mayor, M., & Udry, S. 1998, A&A, 332, 877 [NASA ADS] [Google Scholar]
- Jorissen, A., Frankowski, A., Famaey, B., & van Eck, S. 2009, A&A, 498, 489 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kupka, F. G., Ryabchikova, T. A., Piskunov, N. E., Stempels, H. C., & Weiss, W. W. 2000, Baltic Astron., 9, 590 [Google Scholar]
- Landsman, W., Simon, T., & Bergeron, P. 1993, PASP, 105, 841 [NASA ADS] [CrossRef] [Google Scholar]
- Lind, K., Bergemann, M., & Asplund, M. 2012, MNRAS, 427, 50 [NASA ADS] [CrossRef] [Google Scholar]
- Malina, R. F., Marshall, H. L., Antia, B., et al. 1994, AJ, 107, 751 [NASA ADS] [CrossRef] [Google Scholar]
- Mashonkina, L., Korn, A. J., & Przybilla, N. 2007, A&A, 461, 261 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Masseron, T. 2006, Ph.D. Thesis, Observatoire de Paris [Google Scholar]
- Masseron, T., Plez, B., Van Eck, S., et al. 2014, A&A, submitted [Google Scholar]
- Merle, T., Thévenin, F., Pichon, B., & Bigot, L. 2011, MNRAS, 418, 863 [NASA ADS] [CrossRef] [Google Scholar]
- Mermilliod, J.-C., Andersen, J., Latham, D. W., & Mayor, M. 2007, A&A, 473, 829 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Østensen, R. H., & Van Winckel, H. 2011, in Evolution of Compact Binaries, eds. L. Schmidtobreick, M. R. Schreiber, & C. Tappert, ASP Conf. Ser., 447, 171 [Google Scholar]
- Østensen, R. H., & Van Winckel, H. 2012, in Fifth Meeting on Hot Subdwarf Stars and Related Objects, eds. D. Kilkenny, C. S. Jeffery, & C. Koen, ASP Conf. Ser., 452, 163 [Google Scholar]
- Plez, B. 2012, Astrophysics Source Code Library, 5004 [Google Scholar]
- Pojmanski, G. 1997, Acta Astron., 47, 467 [NASA ADS] [Google Scholar]
- Pompéia, L., Masseron, T., Famaey, B., et al. 2011, MNRAS, 415, 1138 [NASA ADS] [CrossRef] [Google Scholar]
- Pounds, K. A., Allan, D. J., Barber, C., et al. 1993, MNRAS, 260, 77 [NASA ADS] [Google Scholar]
- Pye, J. P., McGale, P. A., Allan, D. J., et al. 1995, MNRAS, 274, 1165 [NASA ADS] [Google Scholar]
- Raskin, G., van Winckel, H., Hensberge, H., et al. 2011, A&A, 526, A69 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Reddy, B. E., Tomkin, J., Lambert, D. L., & Allen de Prieto, C. 2003, MNRAS, 340, 304 [NASA ADS] [CrossRef] [Google Scholar]
- Rojas, M., Drake, N. A., Pereira, C. B., & Kholtygin, A. F. 2013, Astrophysics, 56, 57 [NASA ADS] [CrossRef] [Google Scholar]
- Siess, L., Izzard, R. G., Davis, P. J., & Deschamps, R. 2013, A&A, 550, A100 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Siess, L., Davis, P. J., & Jorissen, A. 2014, A&A, 565, A57 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Soberman, G. E., Phinney, E. S., & van den Heuvel, E. P. J. 1997, A&A, 327, 620 [NASA ADS] [Google Scholar]
- Soker, N. 2000, A&A, 357, 557 [NASA ADS] [Google Scholar]
- Strassmeier, K., Washuettl, A., Granzer, T., Scheck, M., & Weber, M. 2000, A&AS, 142, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Tout, C. A., & Eggleton, P. P. 1988, MNRAS, 231, 823 [NASA ADS] [CrossRef] [Google Scholar]
- Udry, S., Mayor, M., Maurice, E., et al. 1999a, in Precise Stellar Radial Velocities, eds. J. B. Hearnshaw, & C. D. Scarfe, IAU Colloq., 170, ASP Conf. Ser., 185, 383 [Google Scholar]
- Udry, S., Mayor, M., & Queloz, D. 1999b, in Precise Stellar Radial Velocities, eds. J. B. Hearnshaw, & C. D. Scarfe, IAU Colloq., 170, ASP Conf. Ser., 185, 367 [Google Scholar]
- van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vennes, S., Mathioudakis, M., Doyle, J. G., Thorstensen, J. R., & Byrne, P. B. 1995, A&A, 299, L29 [NASA ADS] [Google Scholar]
- Vennes, S., Christian, D. J., Mathioudakis, M., & Doyle, J. G. 1997, A&A, 318, L9 [NASA ADS] [Google Scholar]
- Vennes, S., Christian, D. J., & Thorstensen, J. R. 1998, ApJ, 502, 763 [NASA ADS] [CrossRef] [Google Scholar]
- Vos, J., Østensen, R. H., Degroote, P., et al. 2012, A&A, 548, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vos, J., Østensen, R. H., Németh, P., et al. 2013, A&A, 559, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
All Tables
Line list used for determining atmospheric parameters and chemical composition of atomic species.
All Figures
![]() |
Fig. 1 The Ca ii H and K lines, from the HERMES spectrum on HJD 2 455 066.748 (2009, August 23). The emission present in the cores is somewhat more intense than the one presented in Fig. 5 of Vennes et al. (1997). |
In the text |
![]() |
Fig. 2 The orbital solution for IP Eri. The lower panel shows the observed (O) minus the calculated (C) orbits. |
In the text |
![]() |
Fig. 3 Line by line abundance analysis for Fe i and Fe ii lines as a function of the excitation potential χ and the reduced equivalent width W/λ. The grey areas represent the standard deviations around the mean abundances of Fe i and Fe ii. |
In the text |
![]() |
Fig. 4 Example of CN-line fitting. |
In the text |
![]() |
Fig. 5 Example of Ba ii line fitting (Ba ii λ585.37 nm). |
In the text |
![]() |
Fig. 6 Example of Ce ii line fitting (Ce ii λ533.01 nm). |
In the text |
![]() |
Fig. 7 Comparison of light-element abundances in IP Eri (black dot), in field stars (Edvardsson et al. 1993, grey plusses), and in mild barium stars (Rojas et al. 2013, black crosses). |
In the text |
![]() |
Fig. 8 Comparison of s-process element abundances of IP Eri with field stars (Edvardsson et al. 1993, grey plusses; Reddy et al. 2003, grey squares; Pompéia et al. 2011, grey asterisks) and mild barium stars (Rojas et al. 2013, black crosses). |
In the text |
![]() |
Fig. 9 Evolutionary channels for the formation of a He WD. The dashed lines refer to channels where the eccentricity can be preserved (see text for details). |
In the text |
![]() |
Fig. 10 The positions of IP Eri and its twin system HR 1608 (black circled crosses; Beavers & Eitter 1988, quoted by Landsman et al. 1993) in the eccentricity – period diagram are compared with those of the long-period sdB binaries, with (mostly) pre-mass-transfer binaries (normal K and M giants), and with post-mass-transfer binaries (Barium and S stars). Symbols are as indicated in the figure label. |
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.