Issue |
A&A
Volume 614, June 2018
|
|
---|---|---|
Article Number | A77 | |
Number of page(s) | 12 | |
Section | Stellar structure and evolution | |
DOI | https://doi.org/10.1051/0004-6361/201629789 | |
Published online | 18 June 2018 |
The MUCHFUSS photometric campaign
1
Institute for Astronomy and Astrophysics, Kepler Center for Astro and Particle Physics, Eberhard Karls University,
Sand 1,
72076
Tübingen, Germany
e-mail: schaffenroth@astro.uni-tuebingen.de
2
Institut für Astro- und Teilchenphysik, Universität Innsbruck,
Technikerstr. 25/8,
6020
Innsbruck, Austria
3
Dr. Remeis-Observatory & ECAP, Astronomical Institute, Friedrich-Alexander University Erlangen-Nürnberg,
Sternwartstr. 7,
96049
Bamberg, Germany
4
Department of Physics, University of Warwick,
Coventry
CV4 7AL, UK
5
Argelander Institute for Astronomy,
Auf dem Hügel 71,
53121
Bonn, Germany
Received:
26
September
2016
Accepted:
9
February
2018
Hot subdwarfs (sdO/Bs) are the helium-burning cores of red giants, which have lost almost all of their hydrogen envelope. This mass loss is often triggered by common envelope interactions with close stellar or even substellar companions. Cool companions like late-type stars or brown dwarfs are detectable via characteristic light-curve variations like reflection effects and often also eclipses. To search for such objects, we obtained multi-band light curves of 26 close sdO/B binary candidates from the MUCHFUSS project with the BUSCA instrument. We discovered a new eclipsing reflection effect system (P = 0.168938 d) with a low-mass M dwarf companion (0.116 M⊙). Three more reflection effect binaries found in the course of the campaign have already been published; two of them are eclipsing systems, and in one system only showing the reflection effect but no eclipses, the sdB primary is found to be pulsating. Amongst the targets without reflection effect a new long-period sdB pulsator was discovered and irregular light variations were found in two sdO stars. The found light variations allowed us to constrain the fraction of reflection effect binaries and the substellar companion fraction around sdB stars. The minimum fraction of reflection effect systems amongst the close sdB binaries might be greater than 15% and the fraction of close substellar companions in sdB binaries may be as high as 8.0%. This would result in a close substellar companion fraction to sdB stars of about 3%. This fraction is much higher than the fraction of brown dwarfs around possible progenitor systems, which are solar-type stars with substellar companions around 1 AU, as well as close binary white dwarfs with brown dwarf companions. This might suggest that common envelope interactions with substellar objects are preferentially followed by a hot subdwarf phase.
Key words: subdwarfs / binaries: eclipsing / binaries: spectroscopic / brown dwarfs / stars: fundamental parameters
© ESO, 2018
1 Introduction
Stars do not form in isolation and therefore most of them are members of binary or multiple systems (see Duchêne & Kraus 2013, for a review). The separations of a significant fraction of those binary pairs are so small that the stars will interact with each other during their lifetimes. Several studies have determined the frequency of sufficiently close companions around solar-type stars (see Winn & Fabrycky 2015; Sahlmann et al. 2011, and references therein). The fraction of stellar companions closer than 10 AU is about 13%. While interactions in the main sequence phase are possible for the closest systems, it is more likely that the stars start to interact in the post-main sequence phase as soon as the more massive star expands and evolves to become a red giant increasing its radius by two to three orders of magnitude. Depending mostly on the mass ratio, the resulting interaction will lead to mass-transfer, which can be either stable or unstable. In the most extreme cases, this interaction can change the evolution of the stars dramatically and form objects that can hardly be explained otherwise (see González Martínez-País et al. 2014, Chapter 2, for a review).
An important class of such objects are the hot subdwarf stars (sdO/Bs), which are evolved, core helium-burning objects with thin hydrogen-dominated envelopes and masses around 0.5 M⊙ (see Heber 2009, 2016, for a review). The progenitor star has to lose almost its entire hydrogen envelope in the red-giant phase to form such an object. About half of the sdB stars reside in close binaries with short periods from as little as 1 h to a few days (Maxted et al. 2001). Such close binaries are believed to be formed by a common-envelope and spiral-in phase (Han et al. 2002, 2003), as the separation in these systems is much smaller than the size of the red-giant progenitor star. In this scenario of unstable mass-transfer, a low-mass companion is swallowed by the close-by red giant. Subsequently, the companion spirals towards the core of the red giant while removing the hydrogen envelope. Although the common-envelope ejection channel is not properly understood in detail (see Ivanova et al. 2013, for a review), it provides a reasonable explanation for the strong mass loss required to form sdB stars. Such close post-CE binaries (PCEBs) are also found among the white dwarf population. However, since the majority of stars end their lives as white dwarfs, the relative fraction of PCEBs is much smaller than for the sdOB stars.
It has been known for decades that the lowest-mass stars close to the hydrogen burning limit can expel a common envelope (e.g., Kilkenny et al. 1981). Substellar objects with masses below that limit have been found in close orbits around stars as well. About 10% of the solar-type stars have giant planets with masses greater than 0.1 MJup in orbits smaller than 10 AU. The more massive brown dwarf (BD) companions on the other hand are found only around <1% of the solar-type stars. This paucity of BDs is also known as the brown dwarf desert. All those close substellar companions will eventually interact with their host stars. The discovery of three sdB systems with close brown dwarf companions (Geier et al. 2011c; Schaffenroth et al. 2014, 2015) provides the best evidence that substellar objects are able to interact with a star in a common-envelope phase and form an sdB star.
Han et al. (2002, 2003) performed a binary population study to predict the fractions of hot subdwarf stars coming from different formation channels using Monte-Carlo simulations with different model parameters. Their favored model predicts a fraction of 18% of the hot subdwarf stars coming from the common envelope channel, which produces hot subdwarf stars with low-mass stellar or substellar companions, most of the companions being M dwarfs. Acounting for selection effects, they predict that the measurable fraction of low-mass stellar or substellar companions to hot subdwarf stars should be about 40%, when the selection is sensitive to short-period binaries like in Maxted et al. (2001), for example.
Constraining the fraction of low-mass stellar and substellar companions to hot subdwarf stars can hence help for a better understanding of the formation of hot subdwarf stars and the common envelope phase. The current sample of hot subdwarf binaries consists of 142 binaries, 26 of them with low-mass stellar or substellar companions (Kupfer et al. 2015), which relates to a fraction of 15%. However, the selection is quite inhomogeneous and is therefore not suited to derive the fraction of low-mass stellar or substellar companions to hot subdwarf stars.
Several different typical light variations are found in sdB binaries (more details in Sect. 2.1) allowing us to restrict the nature of their companions. Here we present a photometric follow-up of targets from the MUCHFUSS project. In Sect. 2 we describe the observations and present the already published results of our photometric follow-up campaign. In Sect. 3.1, we present the new discovery of an eclipsing sdB binary with a low-mass M dwarf companion and the discovery of a new sdB pulsator in a binary. Moreover, we show two He sdOs with light variations. In Sect. 4, we try to confine the fraction of substellar companions. Finally, in Sect. 5 we discuss our results.
2 The MUCHFUSS photometric follow-up campaign
The Massive Unseen Companions to Hot Faint Stars from SDSS (MUCHFUSS) project was initially designed to find massive compact companions like massive white dwarfs, neutron stars or even black holes as companions to hot subdwarf stars. To discover rare objects, a huge initial dataset is necessary. The enormous SDSS database was therefore the starting point for the survey. Hot subdwarfs are found most easily by applying a color cut to Sloan photometry. The spectra of the color-selected stars were subsequently classified by visual inspection.
Since the spectra of hot subdwarfs in close binary systems do not show any lines of the companions, we selected all single-lined sdO/Bs with spectra that are sufficiently bright (V < 18.5 mag) to have spectra with reasonable signal-to-noise ratio (S∕N > 5). The SDSS spectra are coadded from three individual spectra with an exposure time of 15 min taken consecutively. Measuring the radial velocities of those individual spectra and performing a large spectroscopic follow-up campaign, we selected close binary sdB stars, which show significant RV-variability (Geier et al. 2011a, 2015, 2017) for the photometric follow-up. In addition to this thoroughly selected sample (see Table A.1) we also observed several hot subdwarf stars as backup targets (see Tables B.1 and B.2). This selection turned out to be not only well suited to find compact companions, but also low-mass main-sequence or substellar companions with very short orbital periods (e.g., Geier et al. 2011b; Schaffenroth et al. 2014).
2.1 Light variations found in sdB binaries
In our spectroscopic follow-up campaign we measured the radial velocity curves and orbital parameters of the binaries. However, in such single-lined spectroscopic binaries only minimum masses can be derived for the companions, as the inclination is unknown. To unravel the nature of the companions, photometric follow-up is helpful, because the diverse classes of companions can cause characteristic signatures in the binary light curves. The so-called reflection effect is observed in close binary systems consisting of a hot primary and a cool companion. As the secondaries in these systems are supposed to orbit synchronously, the hemisphere of the cool companion facing the hot primary is constantly irradiated and heated up, which leads to an increased flux around the secondary eclipse (see Wilson 1990; Budaj 2011, for a detailed discussion of this effect). The amplitude of the reflection effect scales with the temperature ratio and the radii of the primary and secondary star, as well as the inverse orbital separation (Wilson 1990; Budaj 2011). Hence, the reflection effect is strong when both components of a close binary system have a very small separation, similar radii, and a high temperature difference. In systems consisting of a low-mass main sequence star or a brown dwarf and a hot, compact star like a hot subdwarf, these conditions are fulfilled. The reflection effect is visible as a sinusoidal variation with a period equal to the orbital period (see Fig. 1).
Due to the short periods of such sdB binaries, these systems also have a high probability of being eclipsing. Eclipsing systems are of high value because they allow us to determine the masses and radii of both components with a combined photometricand spectroscopic analysis. Moreover, the separation of the system can be determined. Eclipsing binaries with low-mass stellar or substellar companions are called HW Virginis (HW Vir) systems after their prototype. Only 19 of these systems are published so far, which makes this type of system quite rare. Due to the unique shape of thelight curve, almost all of them have been found in photometric surveys or by chance while looking for pulsations.
In Fig. 1, the dependency of the reflection effect amplitude on the inclination angle is illustrated with synthetic light curves adopting the parameters of the prototype HW Vir, which represents a typical HW Vir system (see Table C.1 and Fig. A.1). Even for an inclination angle as small as 10°, variations ofalmost 2% are predicted, which are detectable from the ground. Multi-band light curves are very useful to confirm the reflection effect. Since the temperature of the irradiated hemisphere of the companion is heated up to a temperature of about 10 000–20 000 K (Vučković et al. 2016), which is lower than the temperature of the sdB, the observed variation is wavelength dependent and becomes stronger at longer wavelengths, where the contrast in brightness between the irradiated companion and the hot subdwarf becomes stronger. Due to the much smaller radius of the companion, the reflection effect is not visible in sdB + WD systems (see Fig. 1).
However, in the shortest period sdB binary systems with WD companions, the tidal influence of the rather massive and compact companion can lead to an ellipsoidal deformation of the sdB. This results in a sinusoidal variation with half the orbital period (e.g., Wilson & Vaccaro 1997). The amplitude of the ellipsoidal variation scales with the mass of the companion and the inverse orbital period. Since it is caused by the distorted shape of the primary, the wavelength dependence of the variation is much lessthan for the reflection effect. Multi-band photometry, therefore, allows us to distinguish between variations due to the reflection effect and those due to ellipsoidal deformation. Moreover, Fig. 1 shows that the reflection effect in sdB binaries is much stronger than the ellipsoidal deformation for a given orbital period.
In some cases multiperiodic light variations due to pulsations of the sdB primary are detected in addition to the binary signatures. Two main types of sdB pulsators are known. Kilkenny et al. (1997) discovered the sdBV r stars, also called V361 Hya stars, which show nonradial p-mode pulsations with short periods of ~100–300 s and amplitudes from a few parts per million to about 5%. The second class was found by Green et al. (2003). These stars, called sdBV s stars or V1093 Her stars, show high-radial-order g-mode pulsations with periodsof ~2000–8000 s and very small amplitudes below 0.5%. Since the pulsation periods of the V1093 Her stars are very similar to the orbital periods of the closest sdB binaries and can show a complex wavelength dependency, such pulsations can be misclassified as shallow reflection effects or ellipsoidal variations.
Recently, apparently irregular light variations have been discovered in some sdO stars, the hotter siblings of the sdBs. It is not yet clear where these strange variations come from, but they look quite similar to variations found in cataclysmic variables (Green et al. 2014).
![]() |
Fig. 1 Different light-curve models showing the typical light variations of hot subdwarf binaries superimposed by a typical non-detection light curve of the photometric Muchfuss follow-up. All light-curvemodels were calculated with the MORO code (Drechsel et al. 1995) in photometric band RB . We calculated models of a sdB+dM system with the parameters of HW Vir (see Table C.1) representing a typical HW Virginis system at different inclinations showing the reflection effect and of a sdB binary system with a low-mass WD with the same period and RV semi-amplitude at 90° and 70°. Variations in the light-curve of the sdB+WD systems are due to ellipsoidal deformation of the sdB. |
![]() |
Fig. 2 Transmission curves of the BUSCA passbands taken from the BUSCA homepage (http://www.caha.es/CAHA/Instruments/BUSCA/bands.txt) weighted by the quantum efficiency of the detector(https://www.caha.es/guijarro/BUSCA/QE_thick_thin.txt). |
2.2 Observations
The number of known hot subdwarf stars with cool low-mass companions (reflection effect systems) is still small and most of them have been discovered due to their characteristic light curves in surveys looking for pulsating stars or from the data archives of planetary transit and transient surveys (e.g., Schaffenroth et al. 2013; Almeida et al. 2012).
However, the fraction of low-mass stellar and substellar companions in hot subdwarf binaries cannot be derived based on this heterogeneous sample. This motivated us to conduct a photometric follow-up of 26 spectroscopically selected sdB binary candidates from the MUCHFUSS project (Geier et al. 2015, 2017). In addition to the selection criteria described in Sect. 2, we picked systems brighter than g ≃ 18 mag to obtain light curves with an S/N sufficient to also find small light variations. Moreover, we preferentially chose systems which already had additional spectroscopic follow-up and therefore highly significant RV variations. If no such systems were observable at a given time, backup targets were chosen based on their visibility and magnitude depending on the observing conditions.
This follow-up was done mostly with Bonn University Simultaneous CAmera (BUSCA; see Reif et al. 1999) on the 2.2 m-telescope located at the Calar Alto Observatory in Spain. This instrument turned out to be perfect for our purposes, as it is possible to observe in four bands simultaneously. We did not use any filters but the intrinsic transmission curve given by the beam splitters (UB , BB , RB , IB ). In this way all visible light is used most efficiently. The transmission curves of the BUSCA bandpasses are shown in Fig. 2.
The data were taken during several observing runs listed in Table 1. In total, 30 nights were granted for this programme with BUSCA. For each star time-resolved photometry with a duration between 1.5 and 2.5 h and exposure times from 30 to 180 s per image was taken. The strategy was tailored to find light variations with periods of 2–3 h, which matches the typical orbital periods of the HW Vir systems.
Defining smaller windows around the stars within the full images, the read-out time was reduced from 2 min for the full-frame image to about 15–20 s. We observed our target and four comparison stars in 60 × 60 pixel windows to perform differential photometry. The comparison stars were selected using the SDSS Navigate Tool1. We chose comparison stars with similar magnitudes (Δm < 2 mag) in all SDSS bands from u to z. The data reduction was done using IRAF2. The standard CCD reduction was done using the tools for bias-, and flatfield-correction and the light curves were extracted using the aperturephotometry package daophot. To correct for changing airmass and conditions, the light curve of the target was divided by the light curves of the comparison stars.
Log of observations with BUSCA.
2.3 Summary of published light variations
In the course of our photometric follow-up campaign we already discovered and published several objects. The spectra of SDSS J082053.53+000843.4 (J0820) and SDSS J162256.66+473051.1 (J1622) showed RV variations with short periods and the photometric follow-up confirmed them to be eclipsing with a period of 0.096 d for J0820 (Geier et al. 2011c)3 and 0.069 dfor J1622 (Schaffenroth et al. 2014). The companions of both systems have masses of 0.045–0.068 M⊙ and 0.064 M⊙, respectively.Those systems were the first confirmed sdB+BD systems.
Finally, we discovered the sdB binary FBS 0117+396 showing the reflection effect as well as pulsations. In addition to two short-periodpulsation modes, we detected several modes with long periods. FBS 0117+396 is therefore the first hybrid pulsator in a reflection effect binary with an M-dwarf companion (Østensen et al. 2013).
3 New discoveries
3.1 J192059+372220 – a new HW Vir system
During the course of our photometric follow-up with BUSCA we discovered the new HW Vir system J192059+372220 (J1920), which was selected as an RV-variable sdB. The observations were taken and reduced as described in Sect. 2.2 during the two runs in Sept/Oct 2010/2011. Its light curve, which is shown in Fig. 3, displays a prominent reflection effect and grazing eclipses. We obtained 39 medium-resolution (grism T13, R ~ 4000, wavelength coverage 3450–5750 Å) spectra with the TWIN spectrograph on 26/27/29/30 May 2012 at the 3.5 m telescope at the Calar Alto observatory, which were reduced with the MIDAS package4.
Those spectra were used to determine the radial velocity curve (Fig. 4). The RV was measured using the SPAS routine (Hirsch 2009) by fitting Gaussians and Lorentzians to the Balmer and helium lines. A circular orbit was assumed due to the short period, and therefore a sine curve was fitted to the measured RVs to determine the RV semi-amplitude K = 57.2 ± 2.3 km s−1, the system velocity γ = 18.0 ± 1.7 km s−1, and the orbital period 0.168938 ± 0.00011 d. The results are summarized in Table 2.
The atmospheric parameters of the sdB primary were determined by fitting synthetic spectra, which were calculated using LTE model atmospheres with solar metallicity and metal line blanketing (Heber et al. 2000), to the Balmer and helium lines. Each spectrum was fitted separately, because in several other HW Vir systems an apparent change of the atmospheric parameters with the orbital phase has been detected (Schaffenroth et al. 2013; Heber et al. 2004), which is caused by the time-variable contributionof the irradiated companion. This effect is only visible when the atmospheric parameters can be derived with sufficient accuracy. In our analysis this effect was not visible and thus we co-added all spectra, which were shifted to RV zero before, to improve the S/N. The results are summarized in Table 2. The effective temperature Teff = 27500 ± 1000 K and the surface gravity of logg = 5.4 ± 0.1 are typical for an sdB situated on the extreme horizontal branch.
From the BUSCA light curves we could determine the orbital period to be 0.168965(43) d, which is consistent with the period determined from the RV curve. Unfortunately, the S/N of our data is not sufficent for a proper light curve analysis. A preliminary analysis using MORO (MOdified ROche program, see Drechsel et al. 1995) as described in, for example, Schaffenroth et al. (2014) indicates that the inclination angle of the system is about 67°. Adopting this value and assuming a canonical mass of 0.47 M⊙ for the sdB, we derive a companion mass of 0.116 M⊙ from the binary mass function. This mass estimate is consistent with a late M dwarf making J1920 a quite typical HW Vir system.
![]() |
Fig. 3 BUSCA light curve of the eclipsing sdB+dM system J1920 (UB, BB, RB, IB, from top to bottom). |
![]() |
Fig. 4 Radial velocity curve of J1920 with the best fit superimposed. The residuals are shown in the lower panel. |
Atmospheric and fundamental parameters of J1920.
3.2 GALEX J050735+034815 – a newly discovered pulsator
For the bright backup target5 GALEX J050735+034815 (J0507) (Németh et al. 2012) we detected a small sinusoidal variation with a period of about 2.7 h, which is shown in Fig. 5. Since the amplitude of the variation is similar in the V B and RB bands, a reflection effect can be excluded. Initially, we suspected the variation to be caused by an ellipsoidal deformationand the sdB to be in a close binary with a WD companion similar to the sdB+WD binary CD−30°11223, which was also found in the MUCHFUSS project (Geier et al. 2013).
However, most recently Kawka et al. (2015) determined the binary parameters of this system based on spectroscopy. The orbital period of ~0.528 d is much longer than the period of the light variation we detected, which can therefore not be an ellipsoidal variation. Németh et al. (2012) determined the atmospheric parameters of the sdB star. With an effective temperature of Teff = 23 990 K and a surface gravity of logg = 5.42 the star is situated in the instability region of the long-period sdB pulsators. Since period and amplitude of the observed variation are consistent with a pulsation mode, we conclude that J0507 is most likely an V1093 Her pulsator.
![]() |
Fig. 5 BUSCA light curve of J0507 in BB and RB showing most probable long-period g-mode pulsations. |
![]() |
Fig. 6 BUSCA light curve of He-sdO J1415 in RB showing strong irregular variations. |
3.3 Variations of He-sdOs
As already shown in Green et al. (2014) and Geier et al. (2015), several helium-rich sdO stars have been found to show irregular RV and light variations. To investigate this further we observed several He-sdOs. J141549+111213 (J1415) shows irregular variations (see Fig. 6 and Geier et al. 2015). The origin of these variations is still unclear. Moreover, the He-sdO J112414+402637 (J1124) shows a low-amplitude sinusoidal variation (see Fig. 7). The cause of this variation is also unclear and further RV and light-curve observations are necessary to investigate these systems. No significant variations have been detected in seven other sdO and He-sdO stars (see Table B.2).
![]() |
Fig. 7 BUSCA light curve of the He-sdO J1124 showing low-amplitude light variations (BB and RB ). |
4 The fraction of low-mass stellar and substellar companions
As shown before, 4 out of the 26 observed sdB or sdOB binaries with significant RV shifts show a reflection effect. Three of them show eclipses and two of them have substellar companions. For 21 of the observed sdB or sdOB systems, it is possible to exclude typical sdB+dM systems with parameters similar to HW Vir and also sdB+BD systems with parameters similar to J1622 applying the method outlined in Appendix A. One of the light curves is too noisy and/or the measured RV shifts are too small to be conclusive. In this case, a low-mass stellar and substellar companion cannot be firmly excluded.
Therefore, we conclude that 81% of our sample are not HW Vir-type systems, 15% are HW Vir or reflection effect systems, 8% of the sdB binaries of our sample have substellar companions and for another 4% no conclusions can be drawn. Since we can only firmly exclude systems with short orbital periods, the fractions of low-mass stellar and substellar companions have to be regarded as lower limits only.
Although HW Vir systems have been studied for decades, no other sample allows us to constrain their fraction of reflection effect and HW Vir binaries in a similar way. The discovery history of those systems (see Table C.1) is very inhomogeneous. Their conspicuousness is due to their characteristic light curves. Some have been known for decades (e.g., AA Dor, HW Vir), while others have been found in surveys for pulsating sdB stars (e.g., NY Vir, HS 2231+2441, HS 0705+6700, EC 10246−2707). The most recent discoveries are from large light-curve archives of transit (e.g., ASAS 10232, NSVS 14256825), microlensing (e.g., BUL−SC16 335) or transient (e.g., PTF1 J072455.75+125300.3) surveys. This sample can therefore not be used to draw any quantitative conclusions.
The most comprehensive dataset of sdO/B light curves to date is the one collected by Østensen et al. searching for pulsations (e.g., Østensen et al. 2010b). The most recent catalog contains 776 sdO/Bs with light-curve data (Østensen priv. comm.). Only 12 systems of the MUCHFUSS target list are included, because the SDSS sample is fainter than the rest. Four of them were also observed by us with BUSCA. But more importantly, those light curves have an average duration of just 30 min (the pulsations periods are just a few min) and have often been detrended to filter out long-period variations (some binaries have been missed because of that). This dataset is therefore not suited for our purposes.
Light curves of unprecedented quality and duration have been obtained by the Kepler space mission. However, due to the pre-selection of potential planet host stars, only 29 single-lined sdBs and sdOBs have been observed; among them one HW Vir (2M 1938+4603) and three reflection effect systems (Østensen et al. 2010c). This fraction of 14% seems to be consistent with our result. However, the Kepler targets have not been selected based on RV variations and the number of spectroscopically confirmed close binaries is significantly smaller. The Kepler fraction of reflection effect and HW Vir systems is therefore even higher than the fraction derived here.
It must be pointed out that our sample and the Kepler sample are still small and the derived number fractions are therefore affected by small number statistics. Larger and more homogeneous samples like the ones from the Kepler K2 and the proposed TESS mission will shed more light on this important question.
5 Conclusions
We found another eclipsing sdB binary J1920 with a low-mass M dwarf companion increasing the number of known systems to 20. Moreover, a known sdB binary was found to contain a long-period pulsating sdB star.
Additionally, we tried to constrain the nature of the companions for the other observed sdB binaries and the fractions of cool, low-mass and brown dwarf companions. More than 15% of the close sdB binaries have cool, low-mass companions and more than 8% of the sdB binaries seem to have substellar companions. This fraction is much lower than predict by Han et al. (2003); however, we are only sensitive to short periods and unseen companions such as very low-mass stars or substellar objects and white dwarfs, which means our fraction is only a lower limit.
Several studies have shown that about one third of the sdBs have composite spectra and are therefore members of wide binary systems with solar-type main sequence companions. About half of the single-lined sdBs show significant RV variations on short time scales (see Heber 2016, for a review). Taking those fractions into account, we derive a close brown dwarf fraction of about 3% around sdBs.
The frequency of close brown dwarfs around sun-like stars, the likely progenitors of most sdBs, is <1% and therefore likely lower (Sahlmann et al. 2011). This higher fraction makes it more likely that substellar companions, which were destroyed during the common-envelope phase, could also be responsible for the formation of single sdB stars (Geier et al. 2011c; Schaffenroth et al. 2014).
Only eight close PCEB WD+BD systems are known in double-lined binaries, as the WD is much smaller and hence fainter than an sdB star and the BD particulary in the infrared of comparable brightness. Although no reliable number fractions are known yet, the fraction of close substellar companions to the more than 40 000 spectroscopically identified WDs is likely to be significantly smaller than the one derived for hot subdwarfs. This alone is a strong indication that common envelope interactions with substellar objects are preferentially followed by a hot subdwarf phase. Taking a closer look at the WD masses of the eight known WD+BD systems WD 0137−349B (0.39 M⊙, Maxted et al. 2006), GD 1400 (0.67 M⊙, Farihi & Christopher 2004), WD 0837+185 (0.798 M⊙, Casewell et al. 2012), NLTT 5306 (0.44 M⊙ Steele et al. 2013), CSS 21055 (0.53 M⊙, Littlefair et al. 2014), SDSS J155720.77+091624.6 (0.447 M⊙, Farihi et al. 2017), SDSS J1205−0242 (0.39 M⊙, Parsons et al. 2017) and SDSS J1231+0041 (0.56 M⊙, Parsons et al. 2017), this conclusion is strengthened even further. Six of the known close WD+BD binaries (75%) have masses in the range 0.39–0.56 M⊙, which is very close to the typical mass of core helium-burning hot subdwarfs. The WD mass distribution shows a strong peak at about 0.6 M⊙, while WDs with masses around 0.5 M⊙ are very rare (e.g., Tremblay & Bergeron 2009). It is therefore possible that those six systems have previously evolved through an sdB phase.
To end up on the extreme horizontal branch, the hydrogen envelope of the progenitor star must be removed at the same time as the core of the red giant ignites helium burning. This happens right at the tip of the red giant branch, where the star reaches its maximum radius and the envelope is most loosely bound. Perhaps those conditions are favorable for very low-mass companions to eject the envelope. We conclude that sdB+BD binaries are a likely outcome of the interaction between stars and substellar objects with a suitable initial separation and therefore it is very important to study such interactions.
Acknowledgements
Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA), operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofisica de Andalucia (CSIC) with BUSCA at 2.2 m telescope and TWIN at the 3.5 m telescope. Also based on observations with the ISIS spectrograph at the William Herschel Telescope located at the Roques de los Muchachos observatory at La Palma, Canary Islands. V.S. acknowledges funding by the Deutsches Zentrum für Luft- und Raumfahrt (grant 50 OR 1110), the Erika-Giehrl-Stiftung, and the Deutsche Forschungsgemeinschaft under grant GE 2506/9-1. S.G. acknowledges funding by the Deutsche Forschungsgemeinschaft under grant GE 2506/8-1.
Appendix A Estimation of the fraction of low-mass stellar and substellar companions
To estimate the true fraction of sdB binaries with cool, low-mass companions in our sample, we not only have to count the three HW Virs and one pulsating, reflection effect binary we found (see Sect. 2.2), but also have to quantify how many similar systems (if any) remained undetected in the other 22 binary candidates whose light curves do not showany significant variations. Since the targets were selected spectroscopically we not only have the light curves but also RV shifts measured for all of them. Here we use both to put constraints on the companion type.
In Fig. A.1 we show the relation between the companion mass and the orbital period for the HW Vir systems. From this plot it can be seen that HW Vir represents a typical HW Vir system. Furthermore, J1622 is a good representative for an HW Vir system with brown dwarf companion. Hence, we choose this systems as prototype systems for the further investigation. We calculated synthetic light curves adopting the parameters of HW Vir and J1622 (Table C.1) varying the inclinations with the MORO code (Drechsel et al. 1995). For any given inclination angle we then derived the amplitude of the light variation. In addition, we calculated the RV amplitudes of those systems for different inclination angles. Figure A.2 shows the expected semi-amplitudes of the light variations in RB as well as the expected semi-amplitudes of the radial velocity variations for different inclinations. By fitting parabolas to the points we derived the following relations between the inclination i and the relative semi-amplitude of the light variation in the RB band (650 nm) A as well as the semi-amplitude of the radial velocity K, thereby i is given in degrees. K and A are higher in the case of HW Vir despite the longer period of this binary with respect to J1622 due to the higher mass and larger radius of the companion:
Comparing the two panels of Fig. A.2 it becomes clear that a close sdB+dM/BD binary seen under a small inclination angle will show small light variations and at the same time small RV variations. If a binary candidate with small or undetected light variations should show high RV variations, it cannot be a cool low-mass object, but must be a compact companion like a WD. In this case we can firmly exclude a close sdB+dM/BD binary. We will now apply this criterion to our sample.
To derive upper limits for possible, but undetected light variations we normalized the light curves in the RB -band and calculated mean values and standard deviations. These standard deviations are also a measure of the typical photometric uncertainties for each light curve. Some of the light curves show a linear trend with changing airmass. This is due to the fact that the target sdB star is usually much bluer than the available comparison stars and that the atmospheric extinction is wavelength dependent. We refrained from correcting for this trend, which leads to an overestimation of the standard deviation and, hence, a more conservative uncertainty estimate.
![]() |
Fig. A.1 Period-companion mass relation of the known HW Vir systems. Our prototype systems (HW Vir for an HW Vir system with low-mass main sequence companion and J1622 for an HW Vir with brown dwarf companion) are marked in red. The parameters of the known HW Vir systems are summarized in Table C.1. |
The derived standard deviations are now adopted as upper limits for the semi-amplitudes A of hidden reflection effects (A ≦ σ). Using Eqs. (A.1b) and (A.2b) we then derive the maximum possible inclinations for hidden reflection effect systems similar to HW Vir and J1622 (using Amax(imax) = σ). These inclinations imax are then converted to upper limits for the corresponding RV semi-amplitudes Kmax(imax) and RV shifts 2 ⋅ Kmax using Eq. (A.1a) and (A.2a). If the measured radial velocity shift is higher than the expected one (Δ vrad,max > 2 ⋅ Kmax), we can exclude systems with parameters similar to HW Vir or J1622.
The results for all non-detection light curves are displayed in Table A.1. The table is sorted by the false detection probability for RV variability lnp (Geier et al. 2015). We divided the table into solved systems with known orbital parameters in bold face, RV variable systems with ln p< −9.2 which corresponds to a false-alarm probability <0.01% and RV-variable candidates with lnp < −3.0 (<5%). Table B.1 shows the results for all RV constant and back-up targets, which have been observed during the MUCHFUSS photometric follow-up. Moreover, we observed nine sdOs or He-sdOs selected in the MUCHFUSS project. They are summarized in Table B.2. Two of them showed variations as presented in Sect. 3.3.
Due to bad weather, the light curves of some of our targets do not have the pursued length of 1.7 (the period of J1622) to 2.8 h (the period of HW Vir). However, if the light curve covers only parts of the orbit, the size of the measured amplitude of a possible reflection effect depends on the phase in which the system was observed. For example if we assume that only half of the orbit was observed, we would only measure half of the amplitude in the worst case. This effect has to be considered to draw firm conclusions from the light curves. In the case of J204613−045418, we took a light curve of only 0.67 h. However, Geier et al. (2011b) determined the period of the system to be 0.24311 d, which is much longer than the period of HW Vir or J1622. They also derived a semi-amplitude K1 of 134.3 ± 7.8 km s−1, which is much higher than the K1 of 87.9± 7.8 km s−1. Therefore, we can exclude a typical HW Vir system and a substellar companion. The same is valid for the other systems with known periods (marked in bold). In the case of J153411+543345 only 70% of the orbital period of J1622 was covered. Hence, the amplitude of a hidden reflection effect could be as high as 20.0 o for a system with the parameters of J1622 or 23.7 for the parameters of HW Vir. This translates to a KJ1622,max of about 19 km s−1 and KHW Vir,max of about 25 km s−1, still much smaller than the measured RV shift. The same was done for the other systems, which were not observed for an adequate length of time. The results of this analysis can be found in Table A.1.
Results of the analysis of the observed light curves without variations in the MUCHFUSS follow-up.
Appendix B Back-up targets
![]() |
Fig. A.2 Expected relative semi-amplitude of the reflection effect in RB (left panel) and expected RV semi-amplitude of the hot subdwarf (right panel) for different inclinations. They were calculated with Eqs. (A.1) and (A.2) for the parameters of a typical HW Vir system like the prototype HW Vir (+) and for an sdB+BD system like J1622 (x). |
Observations and light variations of the observed back-up targets.
Observations and light variations of He sdO and sdO from the MUCHFUSS project.
Appendix C Parameters of the known HW Vir systems
Parameters of the known HW Vir systems.
References
- Almeida, L. A., Jablonski, F., Tello, J., & Rodrigues, C. V. 2012, MNRAS, 423, 478 [NASA ADS] [CrossRef] [Google Scholar]
- Barlow, B. N., Kilkenny, D., Drechsel, H., et al. 2013, MNRAS, 430, 22 [NASA ADS] [CrossRef] [Google Scholar]
- Beuermann, K., Dreizler, S., Hessman, F. V., & Deller, J. 2012, A&A, 543, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Budaj, J. 2011, AJ, 141, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Casewell, S. L., Burleigh, M. R., Wynn, G. A., et al. 2012, ApJ, 759, L34 [NASA ADS] [CrossRef] [Google Scholar]
- Derekas, A., Németh, P., Southworth, J., et al. 2015, ApJ, 808, 179 [NASA ADS] [CrossRef] [Google Scholar]
- Drechsel, H., Haas, S., Lorenz, R., & Gayler, S. 1995, A&A, 294, 723 [NASA ADS] [Google Scholar]
- Drechsel, H., Heber, U., Napiwotzki, R., et al. 2001, A&A, 379, 893 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Duchêne, G. & Kraus, A. 2013, ARA&A, 51, 269 [Google Scholar]
- Farihi, J. & Christopher, M. 2004, AJ, 128, 1868 [NASA ADS] [CrossRef] [Google Scholar]
- Farihi, J., Parsons, S. G., & Gänsicke, B. T. 2017, Nat. Astron., 1, 0032 [NASA ADS] [CrossRef] [Google Scholar]
- For, B.-Q., Green, E. M., Fontaine, G., et al. 2010, ApJ, 708, 253 [NASA ADS] [CrossRef] [Google Scholar]
- Geier, S., Hirsch, H., Tillich, A., et al. 2011a, A&A, 530, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geier, S., Maxted, P. F. L., Napiwotzki, R., et al. 2011b, A&A, 526, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geier, S., Schaffenroth, V., Drechsel, H., et al. 2011c, APJ, 731, L22 [NASA ADS] [CrossRef] [Google Scholar]
- Geier, S., Marsh, T. R., Wang, B., et al. 2013, A&A, 554, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geier, S., Østensen, R. H., Heber, U., et al. 2014, A&A, 562, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geier, S., Kupfer, T., Heber, U., et al. 2015, A&A, 577, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geier, S., Kupfer, T., Heber, U., et al. 2017, A&A, 602, C2 [CrossRef] [EDP Sciences] [Google Scholar]
- González Martínez-País, I., Shahbaz, T., & Casares Velázquez J. 2014, Accretion Processes in Astrophysics (Cambridge: Cambridge University Press) [Google Scholar]
- Green, E., Johnson, C., Wallace, S., et al. 2014, eds. V. van Grootel, E. Green, G. Fontaine, & S. Charpinet, in ASP Conf. Ser., 481, 161 [NASA ADS] [Google Scholar]
- Green, E. M., Fontaine, G., Reed, M. D., et al. 2003, ApJ, 583, L31 [NASA ADS] [CrossRef] [Google Scholar]
- Han, Z., Podsiadlowski, P., Maxted, P. F. L., & Marsh, T. R. 2003, MNRAS, 341, 669 [NASA ADS] [CrossRef] [Google Scholar]
- Han, Z., Podsiadlowski, P., Maxted, P. F. L., Marsh, T. R., & Ivanova, N. 2002, MNRAS, 336, 449 [NASA ADS] [CrossRef] [Google Scholar]
- Heber, U. 2009, ARA&A, 47, 211 [NASA ADS] [CrossRef] [Google Scholar]
- Heber, U. 2016, PASP, 128, 082001 [NASA ADS] [CrossRef] [Google Scholar]
- Heber, U., Reid, I. N., & Werner, K. 2000, A&A, 363, 198 [NASA ADS] [Google Scholar]
- Heber, U., Drechsel, H., Østensen, R., et al. 2004, A&A, 420, 251 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hirsch, H. 2009, Phd thesis, Friedrich Alexander Universität Erlangen Nürnberg [Google Scholar]
- Hong, K., Lee,J. W., Lee, D.-J., et al. 2017, PASP, 129, 014202 [NASA ADS] [CrossRef] [Google Scholar]
- Hoyer, D., Rauch, T., Werner, K., Hauschildt, P. H., & Kruk, J. W. 2015, A&A, 578, A125 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ivanova, N., Justham, S., Chen, X., et al. 2013, A&ARv., 21, 59 [Google Scholar]
- Kawka, A., Vennes, S., O’Toole, S., et al. 2015, MNRAS, 450, 3514 [NASA ADS] [CrossRef] [Google Scholar]
- Kilkenny, D., Hill, P. W., & Penfold, J. E. 1981, MNRAS, 194, 429 [NASA ADS] [CrossRef] [Google Scholar]
- Kilkenny, D., Koen, C., O’Donoghue, D., & Stobie, R. S. 1997, MNRAS, 285, 640 [Google Scholar]
- Kupfer, T., Geier, S., Heber, U., et al. 2015, A&A, 576, A44 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lisker, T., Heber, U., Napiwotzki, R., et al. 2005, A&A, 430, 223 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Littlefair, S. P., Casewell, S. L., Parsons, S. G., et al. 2014, MNRAS, 445, 2106 [NASA ADS] [CrossRef] [Google Scholar]
- Maxted, P. f. L., Heber, U., Marsh, T. R., & North, R. C. 2001, MNRAS, 326, 1391 [NASA ADS] [CrossRef] [Google Scholar]
- Maxted, P. F. L., Napiwotzki, R., Dobbie, P. D., & Burleigh, M. R. 2006, Nature, 442, 543 [NASA ADS] [CrossRef] [Google Scholar]
- Németh, P., Kawka, A., & Vennes, S. 2012, MNRAS, 427, 2180 [NASA ADS] [CrossRef] [Google Scholar]
- Németh, P., Ziegerer, E., Irrgang, A., et al. 2016, ApJ, 821, L13 [NASA ADS] [CrossRef] [Google Scholar]
- Østensen, R., Oreiro, R., Drechsel, H., et al. 2007, in 15th European Workshop on White Dwarfs, eds. R. Napiwotzki, & M. R. Burleigh, ASP Conf. Ser., 372, 483 [NASA ADS] [Google Scholar]
- Østensen, R. H., Oreiro, R., Hu, H., Drechsel, H., & Heber, U. 2008, Hot Subdwarf Stars and Related Objects, eds. U. Heber, C. S. Jeffery, & R. Napiwotzki, ASP Conf. Ser., 392, 221 [NASA ADS] [Google Scholar]
- Østensen, R. H., Green, E. M., Bloemen, S., et al. 2010a, MNRAS, 408, L51 [NASA ADS] [CrossRef] [Google Scholar]
- Østensen, R. H., Oreiro, R., Solheim, J.-E., et al. 2010b, A&A, 513, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Østensen, R. H., Silvotti, R., Charpinet, S., et al. 2010c, MNRAS, 409, 1470 [NASA ADS] [CrossRef] [Google Scholar]
- Østensen, R. H., Geier, S., Schaffenroth, V., et al. 2013, A&A, 559, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Parsons, S. G., Hermes, J. J., Marsh, T. R., et al. 2017, MNRAS, 471, 976 [NASA ADS] [CrossRef] [Google Scholar]
- Polubek, G., Pigulski, A., Baran, A., & Udalski, A. 2007, in 15th European Workshop on White Dwarfs, eds. R. Napiwotzki, & M. R. Burleigh, ASP Conf. Ser., 372, 487 [NASA ADS] [Google Scholar]
- Pribulla, T., Dimitrov, D., Kjurkchieva, D., et al. 2013, Information Bulletin on Variable Stars, 6067, 1 [NASA ADS] [Google Scholar]
- Reif, K., Bagschik, K., de Boer, K. S., et al. 1999, eds. M. M. Blouke, & G. M. Williams, SPIE Conf. Ser., 3649, 109 [NASA ADS] [CrossRef] [Google Scholar]
- Sahlmann, J., Ségransan, D., Queloz, D., et al. 2011, A&A, 525, A95 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schaffenroth, V. 2015, Phd thesis, Friedrich Alexander Universität Erlangen Nürnberg [Google Scholar]
- Schaffenroth, V., Geier, S., Drechsel, H., et al. 2013, A&A, 553, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schaffenroth, V., Geier, S., Heber, U., et al. 2014, A&A, 564, A98 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schaffenroth, V., Barlow, B. N., Drechsel, H., & Dunlap, B. H. 2015, A&A, 576, A123 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schindewolf, M., Levitan, D., Heber, U., et al. 2015, A&A, 580, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Steele, P. R., Saglia, R. P., Burleigh, M. R., et al. 2013, MNRAS, 429, 3492 [NASA ADS] [CrossRef] [Google Scholar]
- Tremblay, P.-E., & Bergeron, P. 2009, ApJ, 696, 1755 [NASA ADS] [CrossRef] [Google Scholar]
- Van Grootel, V., Charpinet, S., Brassard, P., Fontaine, G., & Green, E. M. 2013, A&A, 553, A97 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vučković, M., Aerts, C., Østensen, R., et al. 2007, A&A, 471, 605 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Vučković, M., Bloemen, S., & Østensen, R. 2014, eds. V. van Grootel, E. Green, G. Fontaine, & S. Charpinet, ASP Conf. Ser., 481, 259 [NASA ADS] [Google Scholar]
- Vučković, M., Østensen, R. H., Németh, P., Bloemen, S., & Pápics, P. I. 2016, A&A, 586, A146 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wilson, R. E. 1990, ApJ, 356, 613 [NASA ADS] [CrossRef] [Google Scholar]
- Wilson, R. E., & Vaccaro, T. R. 1997, MNRAS, 291, 54 [NASA ADS] [CrossRef] [Google Scholar]
- Winn, J. N., & Fabrycky, D. C. 2015, ARA&A, 53, 409 [NASA ADS] [CrossRef] [Google Scholar]
- Wolz, M. 2016, Teacher’s thesis, Friedrich Alexander Universität Erlangen Nürnberg [Google Scholar]
http://www.eso.org/sci/software/esomidas//. Furthermore, 11 spectra with grism R600 giving similar resolution (coverage from 3500 to 5300 Å) from the ISIS spectrograph at the William Herschel Telescope at the Roque de los Muchachos Observatory on La Palma were taken on the 9 and 11 July 2012, which were reduced using IRAF2. The exposure time of the TWIN and ISIS spectra was 20 min.
see Table B.1 for a list of observed backup targets.
All Tables
Results of the analysis of the observed light curves without variations in the MUCHFUSS follow-up.
All Figures
![]() |
Fig. 1 Different light-curve models showing the typical light variations of hot subdwarf binaries superimposed by a typical non-detection light curve of the photometric Muchfuss follow-up. All light-curvemodels were calculated with the MORO code (Drechsel et al. 1995) in photometric band RB . We calculated models of a sdB+dM system with the parameters of HW Vir (see Table C.1) representing a typical HW Virginis system at different inclinations showing the reflection effect and of a sdB binary system with a low-mass WD with the same period and RV semi-amplitude at 90° and 70°. Variations in the light-curve of the sdB+WD systems are due to ellipsoidal deformation of the sdB. |
In the text |
![]() |
Fig. 2 Transmission curves of the BUSCA passbands taken from the BUSCA homepage (http://www.caha.es/CAHA/Instruments/BUSCA/bands.txt) weighted by the quantum efficiency of the detector(https://www.caha.es/guijarro/BUSCA/QE_thick_thin.txt). |
In the text |
![]() |
Fig. 3 BUSCA light curve of the eclipsing sdB+dM system J1920 (UB, BB, RB, IB, from top to bottom). |
In the text |
![]() |
Fig. 4 Radial velocity curve of J1920 with the best fit superimposed. The residuals are shown in the lower panel. |
In the text |
![]() |
Fig. 5 BUSCA light curve of J0507 in BB and RB showing most probable long-period g-mode pulsations. |
In the text |
![]() |
Fig. 6 BUSCA light curve of He-sdO J1415 in RB showing strong irregular variations. |
In the text |
![]() |
Fig. 7 BUSCA light curve of the He-sdO J1124 showing low-amplitude light variations (BB and RB ). |
In the text |
![]() |
Fig. A.1 Period-companion mass relation of the known HW Vir systems. Our prototype systems (HW Vir for an HW Vir system with low-mass main sequence companion and J1622 for an HW Vir with brown dwarf companion) are marked in red. The parameters of the known HW Vir systems are summarized in Table C.1. |
In the text |
![]() |
Fig. A.2 Expected relative semi-amplitude of the reflection effect in RB (left panel) and expected RV semi-amplitude of the hot subdwarf (right panel) for different inclinations. They were calculated with Eqs. (A.1) and (A.2) for the parameters of a typical HW Vir system like the prototype HW Vir (+) and for an sdB+BD system like J1622 (x). |
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.