A modern study of HD 166734: a massive supergiant system Based on observations collected at the European Southern Observatory (La Silla, Chile) with FEROS and TAROT and on data collected at the San Pedro Mártir observatory (Mexico).

A modern study of HD 166734: a massive supergiant system thanks: Based on observations collected at the European Southern Observatory (La Silla, Chile) with FEROS and TAROT and on data collected at the San Pedro Mártir observatory (Mexico).

L. Mahy Space sciences, Technologies, and Astrophysics Research (STAR) Institute, Université de Liège, Quartier Agora, Bât B5c, Allée du 6 août, 19c, B-4000 Liège, Belgium
mahy@astro.ulg.ac.be Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, Bus 2401, B-3001 Leuven, Belgium F.R.S.-FNRS Postdoctoral researcher
   Y. Damerdji Centre de Recherche en Astronomie, Astrophysique et Géophysique, route de l’Observatoire BP 63 Bouzareah, 16340 Algiers, Algeria Space sciences, Technologies, and Astrophysics Research (STAR) Institute, Université de Liège, Quartier Agora, Bât B5c, Allée du 6 août, 19c, B-4000 Liège, Belgium
   E. Gosset Space sciences, Technologies, and Astrophysics Research (STAR) Institute, Université de Liège, Quartier Agora, Bât B5c, Allée du 6 août, 19c, B-4000 Liège, Belgium
mahy@astro.ulg.ac.be F.R.S.-FNRS Senior Research Associate
   C. Nitschelm Unidad de Astronomía, Facultad de Ciencias Básicas, Universidad de Antofagasta, Antofagasta, Chile    P. Eenens Departamento de Astronomía, Universidad de Guanajuato, Apartado 144, 36000 Guanajuato, GTO, Mexico    H. Sana Instituut voor Sterrenkunde, KU Leuven, Celestijnenlaan 200D, Bus 2401, B-3001 Leuven, Belgium    A. Klotz Université de Toulouse, UPS-OMP, IRAP, Toulouse, France
Received September 15, 1996; accepted March 16, 1997
Key Words.:
Stars: early-type - Stars: binaries: general - Stars: fundamental parameters - Stars: individual: HD 166734
offprints: L. Mahy


Aims:HD 166734 is an eccentric eclipsing binary system composed of two supergiant O-type stars, orbiting with a 34.5-day period. In this rare configuration for such stars, the two objects mainly evolve independently, following single-star evolution so far. This system provides a chance to study the individual parameters of two supergiant massive stars and to derive their real masses.

Methods:An intensive monitoring was dedicated to HD 166734. We analyzed mid- and high-resolution optical spectra to constrain the orbital parameters of this system. We also studied its light curve for the first time, obtained in the VRI filters. Finally, we disentangled the spectra of the two stars and modeled them with the CMFGEN atmosphere code in order to determine the individual physical parameters.

Results:HD 166734 is a O7.5IfO9I(f) binary. We confirm its orbital period but we revise the other orbital parameters. In comparison to what we found in the literature, the system is more eccentric and, now, the hottest and the most luminous component is also the most massive one. The light curve exhibits only one eclipse and its analysis indicates an inclination of . The photometric analysis provides us with a good estimation of the luminosities of the stars, and therefore their exact positions in the Hertzsprung-Russell diagram. The evolutionary and the spectroscopic masses show good agreement with the dynamical masses of 39.5  for the primary and 33.5  for the secondary, within the uncertainties. The two components are both enriched in helium and in nitrogen and depleted in carbon. In addition, the primary also shows a depletion in oxygen. Their surface abundances are however not different from those derived from single supergiant stars, yielding, for both components, an evolution similar to that of single stars.


1 Introduction

Massive stars constitute a small part of the stellar population but have a considerable influence on their environment. About 70% of them (Sana et al. 2012), and even 90% if we consider orbital periods longer than 3500 days (Sana et al. 2014), are moreover expected to be binaries or multiple systems. This multiplicity can affect the way that they evolve through tidally-induced-rotational mixing, exchange of matter and angular momentum through Roche lobe overflow, or even the coalescence of both components. In the case of young systems with large separations, the different components evolve independently, following single-star evolution (Pavlovski & Hensberge 2005) whilst, when the components interact with each other, their evolution can completely differ from that of single stars, affecting their size, mass, surface chemistry and/or even their rotational rate.

To better understand how massive stars evolve and improve evolutionary models, comparisons between the latter and spectroscopic observations are essential. Nevertheless, for single stars, the modeling by atmosphere codes does not provide all the information about their fundamental properties, especially on the real masses of the stars. Herrero et al. (1992) raised the attention to the mass discrepancy problem. This term refers to a disagreement between the mass determined by comparing the position of the star in the Hertzsprung-Russell (HR) diagram according to the predictions of stellar evolution calculations and the mass obtained from the spectroscopically derived values of the surface gravity and of the radius of the star. In this context, detached eclipsing binary systems are of crucial importance to provide direct measurements of the physical properties of stars such as their radius or their dynamical mass, allowing us to test the predictive abilities of evolutionary models for single stars.

HD 166734 (V411 Ser) is particularly promising for such an analysis. Unfortunately, this system has neither a measured inclination nor a cluster membership. It was located by Conti et al. (1980) at a distance of 2.3 kpc from Earth. It is composed of two supergiant O-type stars in a wide eccentric orbit. This configuration, rather rare for such stars, suggests few interactions between each other, making them perfect objects for testing the theoretical schemes of stellar evolution. HD 166734 has not been investigated since the seventies when Conti et al. (1980) provided for the first time the parameters of its orbit. These authors described this system as an O7.5If primary and an O9I secondary moving on an eccentric orbit () with a period of 34.54 days. These two objects had similar luminosities but the secondary presents a slightly larger minimum mass than the primary (31  vs 29 , respectively). On the basis of these minimum masses Conti et al. (1980) predicted eclipses in the light curve of HD 166734. Even though Otero & Wils (2005) mentioned HD 166734 as photometrically variable, the light curve has never been subject to a dedicated analysis.

The present paper aims at showing the advances that have been achieved in the analysis of this interesting system and is organized as follows. The spectroscopic and photometric data are described in Section 2, together with their reduction procedure. The radial velocity (RV) determination and the orbital solution of HD 166734 are presented in Section 3. Section 4 reports the analysis of the VRI light curves. In Section 5, we disentangle the spectra of each component and we model them with an atmosphere code to constrain the individual parameters and the surface abundances of each star. Section 6 is devoted to the winds of both components and to their interactions. Section 7 provides the evolutionary statuses of both components and discusses the distance to HD 166734. Finally our conclusions are given in Section 8.

2 Observations and data reduction

2.1 Spectroscopy

HD 166734 was observed through a multi-epoch campaign with two different instruments. We first acquired high-resolution optical spectra with FEROS (Fiber-fed Extended Range Optical Spectrograph) successively mounted on the ESO-1.52m (for the observations taken before 2002) and on the ESO/MPG-2.2m telescopes at La Silla (Chile). This instrument has a resolving power of about 48 000 and the detector was a 2k 4k EEV CCD with a pixel size of 15m 15m. It provides échelle spectra composed of 39 orders, allowing us to cover the [3800–9200]Å wavelength domain. The data reduction was performed with an improved version of the FEROS pipeline as described by Sana et al. (2006). The data normalization was performed by fitting polynomials of degree 4 or 5 to carefully chosen continuum windows.

We completed this dataset by acquiring medium-resolution échelle spectra with the ESPRESSO spectrograph. This instrument is mounted on the 2.12m telescope at Observatorio Astronómico Nacional of San Pedro Mártir (SPM) in Mexico. The spectra are slitted over 27 orders to cover the [3800–6950]Å region. The spectral resolving power is of about 18 000. To increase the signal-to-noise ratio of our data, we took several consecutive observations and we combined them. The data were reduced by using the échelle package available within the ESO-MIDAS software.

The journal of the observations is provided in Table 4 in the Appendix.

2.2 Photometry

We obtained the light curve of HD 166734 in the VRI filters with the TAROT (Télescope à Action Rapide pour l’Observation des phénomènes Transitoires) robotic observatory located at La Silla. It is a fully autonomous 25 cm aperture telescope () covering a field of view thanks to an Andor CCD camera. Its spatial sampling is /pixel. The reduction method of the TAROT observations is described in Damerdji et al. (2007). The object was observed between March 2011 and September 2016 in the V filter, between February 2010 and September 2016 in the R filter and between February 2011 and September 2016 in the I filter. These three light curves are composed of 5021, 5178, and 4928 points, respectively.

3 Radial velocity measurements and orbital solution

Figure 1 shows FEROS spectra as a function of the orbital phase (computed from orbital parameters determined below and listed in Table 1). It illustrates the variability of the line profiles for two particular regions featuring the He i 4471 absorption line as well as the C iii 5696 emission line.

Figure 1: FEROS spectra in the region of He i 4471 (left) and C iii 5696 (right) shifted vertically as a function of phase (times 3 for clarity) from 0.0 at the bottom to 1.0 at the top.

To construct the orbital solution of HD 166734, we first measure the RVs of both components. For this purpose, we fit their spectral lines with Gaussian profiles. The computations of the RVs are done by adopting the rest wavelengths from Conti et al. (1977) for the lines with a rest wavelength shorter than 5000 Å, and those of Underhill (1995) for the others. Then we refine those values (notably for the spectra where the lines of both components are blended) by using a new method designed by one of the authors (YD). The latter is based on a minimum distance method considering a comparison between the observed spectrum and a composite (two-components) theoretical template. We perform this measurement for several lines individually, no matter their ionization stage. The final set of RVs (provided in Table 4 in the Appendix) is obtained by computing the mean values of the radial velocities measured on all the various spectral lines.

SB2 solution
Primary Secondary
 [day] 34.537723 0.001330
 [HJD  2 450 000] 2195.064 0.036
1.179 0.016
 [km s] 7.12 0.96
 [km s] 142.01 1.53 167.39 1.79
 [] 76.23 0.89 89.85 1.05
 [] 27.92 0.81 23.69 0.68
rms [km s]
111 The orbital solution is computed from 129 RV measurements (listed in the Appendix). refers to the periastron passage. , and denote the apparent systemic velocity, the semi-amplitude of the radial velocity curve and the projected separation between the center of the star and the center of mass of the binary system.
Table 1: Orbital solution of HD 166734. The given errors correspond to 1 - .

The determination of the orbital period is performed from the Heck-Manfroid-Mersch technique (hereafter HMM, Heck et al. 1985) that was revised by Gosset et al. (2001, see their appendix). This technique is basically a Fourier method generalized to address the unevenly sampled time series. The HMM power spectrum is determined independently of the mean of the data and thus does not suffer from the biases of the Scargle (1982) method. As input of this technique, we first use the RVs measured on the N iv 4058 line that is only visible in the primary’s spectrum. The periodogram exhibits two outstanding peaks at  d and at  d with a strong neighbor at  d. However, given their large width, the two last-quoted frequencies overlap. We then use HMM on our final set of RVs (provided in Table 4), on the basis of the difference between the secondary RVs and the primary RVs (RV–RV), and detect a frequency at  d, giving the period of  days and a neighbor at  d, corresponding to a period of 34.534560 days. We determine an error of  d (corresponding to an error on the period of  days) calculated on our dataset from the expression given by Lucy & Sweeney (1971, see also or , for further explanations). We compare these values to the result obtained by Conti et al. (1980). We stress that these authors determined an orbital frequency of  d from the difference between the RVs of the two stars, but that the computed periodogram also shows a neighbor at  d that overlaps with the other peak. To decrease the widths of the peaks, we thus need a dataset covering a longer timescale. We merge the two datasets and we determine in the periodogram once again two outstanding peaks at  d, corresponding to a period of 34.627837 days and at  d, corresponding to a period of 34.537542 days. By coupling the dataset of Conti et al. (1980) with our dataset, the error on the frequency becomes  d which corresponds to an error on the period of  days. The periodogram is exhibited in Fig. 11 in the Appendix. We then use the Liège Orbital Solution Package (LOSP222LOSP is developed and maintained by H. Sana and is available at http://www.stsci.edu/hsana/losp.html. The algorithm is based on the generalization of the SB1 method of Wolfe et al. (1967) to the SB2 case along the lines described in Rauw et al. (2000) and Sana et al. (2006).) to determine the SB2 orbital solution of the system. The refined value of the orbital period, provided by this program (in order to improve the of the fit) is days and is adopted in the following. We thus confirm the period computed by Conti et al. (1980) on the basis of the difference between their secondary RVs and their primary RVs and we emphasize that our refined value is similar to the period associated to the alias found in the periodogram computed from the two coupled datasets. The parameters are listed in Table 1 whilst the RV curves are shown in Fig. 2. From these orbital parameters, the inferior conjunction is expected at and the superior conjunction is at .

Figure 2: Radial velocity curves of HD 166734 computed with the RVs from the present paper (provided electronically). The primary is represented by filled symbols whilst the open ones indicate the RVs of the secondary. Circles represent the RVs measured on the FEROS spectra and diamonds indicate the RVs measured on the ESPRESSO data. corresponds to the periastron passage.

Our orbital solution presents an orbit more eccentric (0.62 vs. 0.46, respectively) than that derived by Conti et al. (1980). Moreover, our monitoring (more intensive near the periastron passage) allows us to better constrain the semi-amplitudes, and thus the minimal mass of each component. In this new orbital solution, the primary is the most massive component of the system, which was not the case in the analysis of Conti et al. (1980). This result seems thus to suggest that the two components of HD 166734 show a ‘more classical’ evolution than was first suspected (Conti et al. 1980).

4 Photometry

Before studying the light curve to constrain the parameters of both objects, we submit the three photometric datasets of HD 166734, obtained in the VRI filters, to a detailed Fourier analysis by HMM in order to search for the periods. Surprisingly, no outstanding peak is detected around 0.0289 d as would have been expected from the orbital frequency. Instead, we detect in each of the three datasets a series of peaks spaced with the same frequency of 0.0289 d (similar to the orbital frequency), starting at  d until  d (Fig. 3). The dominant member varies from one filter to another, going from  d in the V and R filters to  d in the I filter (see Figs. 12, 13 in the Appendix). This series of frequencies corresponds in fact to harmonics of the orbital frequency, even though the fundamental and the first-order harmonics are not detected. This large number of harmonics is explained by the small size of the eclipse with respect to the orbital cycle. Moreover, these harmonics have aliases between 0.5 and 1.0 d and at higher frequencies. We also apply HMM on the publicly available ASAS data and still find most of these harmonics there.

Figure 3: Top: Power spectrum of the V light curve of HD 166734. The diamond indicates the location of the orbital frequency while the dots represent its harmonics. Bottom: Spectral window computed from the V light curve.

Because of the inclination of the orbit and the orientation of its apses with respect to our line of sight, the light curve of HD 166734 exhibits only one eclipse during its orbital cycle, immediately after the periastron passage () with a total duration in orbital phase of . The phase of this eclipse is in agreement with the estimation of Otero & Wils (2005) but differs from the phase of the conjunction (see Section 3). This difference can come from the errors on the orbital parameters. We emphasize that if we fix the position of the eclipse to the conjunction it does not change the values of the orbital parameters within the uncertainties. This eclipse shows a difference of magnitude of about 0.22 mag and is preceeded by an increase of the total brightness of the system by about 0.03 mag at the periastron passage. We display in Fig. 4 the light curve of HD 166734 in the three different filters sampled over the orbital period. We stress that these light curves contain outliers (with an unknown origin) that do not change their general shape.

To limit the possible effects of the degeneracies that can be produced because of the presence of only one eclipse we use the PHOEBE (PHysics Of Eclipsing BinariEs, v0.31a, Prša & Zwitter 2005) software while keeping the parameters of the orbit fixed at their values obtained through the spectroscopic analysis (Section 3, Table 1). PHOEBE allows one to model the light curve and the RV curves of the object at the same time. This software is based on Wilson & Devinney’s code (Wilson & Devinney 1971) and uses Nelder & Mead’s Simplex fitting method to adjust all the input parameters and to find the best fit to the light curve. Furthermore, it allows one to model multiple bandpasses at once. We include the reflection effects between the two stars in the computation of the light curve and we also fix the effective temperatures to 32000 K and 30500 K, for the primary and the secondary components, respectively, as determined in Section 5. The photometric parameters derived from the modeling of the light curve are given in Table 2. The error bars are computed by exploring the parameter space. For this purpose, we fix one parameter and allow the others to vary in order to reach the minimum of the corresponding to a 68.3% confidence level (1-). Fig. 4 displays the best-fit of the light curve obtained by PHOEBE in the three different bandpasses.

Figure 4: Light curves of HD 166734 obtained with TAROT in the V (top panel), R (second panel) and I (third panel) filters. The red lines represent the best fit computed with PHOEBE. The magnitude outside the eclipse has been set to 0 in the three filters. The representation shown at the bottom is a zoom-in on the eclipse observed in the V light curve (plotted at the top of this Figure).

The system is seen under an inclination of . The eclipse occurs when the primary is in front of the secondary (secondary eclipse) along the line of sight. The two components do not fill their Roche lobe over the entire orbit. They fill about 90%, within an error of 4%, at the periastron displaying a small deformation because of tidal interactions. With this inclination the masses obtained from the orbital solution are estimated to  for the primary and  for the secondary. Their radii are computed to  and , respectively.

The secondary component thus appears slightly more evolved (cooler and about the same size) than the primary. The modeling of the light curve also reports a primary component brighter than the secondary. From their radii and their effective temperatures, we determine stellar luminosities of for the primary and of for the secondary, infering a brightness ratio between the two stars of .

Primary Secondary
3.499 3.026
1.0 1.0
1.0 1.0
333 is the radius toward the pole, is the radius measured toward the other component, is the radius measured toward the direction perpendicular to the two above-mentioned axes, toward the Lagrangian point and is the geometrical mean between , and . refers to the value of the Roche-model potential used in the WD program. are the bolometric albedos (coefficients of reprocessing of the emission of a companion by ”reflection”) and represents the gravity darkening coefficients.
Table 2: Photometric parameters of HD 166734 obtained with PHOEBE. The uncertainties correspond to 1–.

5 Individual parameters

5.1 Disentangling

We use the orbital parameters derived in Section 3 as input of the disentangling program. We apply the Fourier approach of Hadrava (1995) to separate the spectral contributions of both components. This method uses Nelder & Mead’s Downhill Simplex on the multidimensional parameter space to reach the best fit between the recombined component spectra and the observed data. We emphasize that we only apply this technique on the FEROS spectra because they have a better signal-to-noise ratio and a higher resolution.

HD 166734 is a detached eclipsing binary. During the eclipse, the secondary component is however not totally eclipsed. There exists therefore an ambiguity in the determination of the continuum level. The method of Fourier presented by Hadrava (1995) induces oscillations in the continuum of the disentangled spectrum of each component when one of them is not completely eclipsed, as is the case for HD 166734. A renormalization of the disentangled spectra must thus be done with the light ratio derived from the light curve in Section 4 and provided in Table 2. A discussion of the spectral disentangling can be found in Pavlovski & Hensberge (2010) and in Pavlovski & Southworth (2012).

Figure 5: Best CMFGEN fit (red) of the disentangled optical spectrum (black) of the primary component of HD 166734.

To better constrain the continuum level of both disentangled spectra, we decided, however, to compare them with those obtained from the method of Marchenko et al. (1998), revised by González & Levato (2006). As already stressed by Mahy et al. (2015), the technique of González & Levato (2006) gives reliable results outside broad lines such as the Balmer lines. In the wings of those lines, the disentangled spectra exhibit intrinsic excesses or deficiencies of the flux, making them difficult to normalize, and infering larger uncertainties on the shape of their wings. The two techniques have thus their respective advantages, and their comparison allows us to remove the artefacts linked to the disentangling method.

Figure 6: Same as Fig. 5 but for the secondary component of HD 166734.

As mentioned above, the resulting spectra are normalized to remove any oscillation of the continuum linked to the Fourier disentangling technique and then corrected for the brightness ratio of . The final disentangled primary and secondary spectra are displayed in Fig. 5 and in Fig. 6, respectively.

5.2 Spectral-type classification

From the final disentangled spectra, we determine the equivalent widths of several lines to apply Conti’s criteria (Conti & Alschuler 1971; Conti 1973; Mathys 1988, 1989) and derive the spectral classifications of both components. We measure ratios of and , which gives spectral classifications of O7.5 for the primary and of O9 for the secondary. We can then determine the luminosity class of each component using . We find that this value is equal to 0.93 and to 0.65 for the primary and the secondary components, respectively, indicating that both stars are supergiant. We observe the lines of He ii 4686 and N iii 4634–41 in strong emission in the spectrum of the primary. According to Walborn (1971), we add the ‘f’ suffix to the spectral type of the primary. For the secondary, we observe in its final disentangled spectrum that the absorption of the He ii 4686 line is compensated by the emission, whilst the N iii 4634–41 lines are moderately in emission, infering a ‘(f)’ suffix to its spectral classification.

We also compare the component spectra of both stars of the system to the spectral atlas of Walborn & Fitzpatrick (1990), confirming their classification.

In summary, we attribute to the primary an O7.5 If spectral type whilst the secondary is classified as an O9 I(f) star, with an uncertainty of one subtype for each component. These classifications are in agreement with those obtained by Conti et al. (1980), within the uncertainties.

5.3 Modeling

To determine the physical parameters of both components, we use the CMFGEN atmosphere code (Hillier & Miller 1998) on the disentangled spectra corrected for the brightness ratio. This code provides non-LTE atmosphere models including winds and line-blanketing. CMFGEN requires input from an estimate of the hydro-dynamical structure connected to a velocity law of the form , where is the terminal velocity of the wind. Our final models include the following chemical elements: H i, He i-ii, C ii-iv, N ii-v, O ii-v, Al iii, Ar iii-iv, Mg ii, Ne ii-iii, S iii-iv, Si ii-iv, Fe ii-vi, and Ni ii-v with the solar composition of Grevesse et al. (2007) unless otherwise stated. CMFGEN allows a super-level approach to reduce the memory requirements. Therefore, we include about 1771 super levels for a total of about 8075 levels. The formal solution of the radiative transfer equation allows us to compute an emergent spectrum. To do so, we use a depth-independent microturbulent velocity varying linearly from 20 km s at the photosphere to at the top of the atmosphere. The choice of 20 km s for the microturbulence velocity at the photosphere is coherent with other studies made on supergiant massive star populations (see Bouret et al. 2012 or Crowther & Evans 2009). Once the emergent spectrum is calculated we convolve it twice, once with a rotational profile to include the projected rotational velocity () and once with a Gaussian profile to include macroturbulence () in the line profiles.  is determined from the Fourier transform method given by Simón-Díaz & Herrero (2007) whilst  is fixed by eye from the He i 4713 and O iii 5592 lines. We compute values of  km s and  km s for the  and the  of the primary whilst values of  km s and  km s are obtained for the secondary, respectively.

The effective temperature of each component is estimated from the ratio between He i and He ii, mainly from the He i 4471 and He ii 4542 lines, even though the other helium lines are also considered. The typical uncertainty on our determination is 1000 K. The surface gravities of both components are determined from the wings of the Balmer lines H and H. The typical errors on the surface gravities of both components are 0.1 dex. H and H are not taken into account because they are more affected by the stellar winds, or maybe even by the interaction between the winds of the two components. Furthermore, the shapes of these two lines can vary with respect to the orbital phase. In the final disentangled spectra, we obtain mean profiles computed over the entire orbital period that are not necessarily meaningful. The stellar luminosity is determined from the radii derived through the photometry (see Sect. 4) and from the effective temperature determined for each component through the CMFGEN best-fit (in an iterative way). The mass-loss rates and the clumping filling factors are adjusted to reproduce the H and H lines as well as the He ii 4686 line whilst the shape of their emission is sensitive to the exponent of the wind velocity law. The terminal velocity has been computed from the Newtonian escape velocity (where is the gravitational constant, the Eddington factor computed from Vink et al. (2001) by assuming an electron scattering cross-section per unit mass of 0.32 cm g (Lamers & Leitherer 1993), the mass of the star and its radius) by assuming .

The best-fit model has been obtained by minimizing the calculated by varying different parameters in the parameter space. To this end, we generate a non-uniform grid composed of several dozen models for each star. Once all the fundamental parameters are constrained, we run models with different surface abundances (for He, C, N, and O). Specific lines of each element are selected from the disentangled spectra and we quantitatively compare these lines to the synthetic spectra by means of a analysis from which we derive the surface abundance and their uncertainties. This method is similar to that proposed by Martins et al. (2015), except that, rather than renormalizing the function so that the minimum has a value of 1.0, we subtract an appropriate quantity to set the minimum to 1.0. The uncertainty is approximately set to the abundances for which . The lines used to constrain the abundances are C iii 4068–4070, and C iii 4153–56–63 for carbon, N iv 4058 (only for the primary), and N iii 4511–4515–4524 for nitrogen and O iii 5592 for oxygen. We exclude the C iii 4647–50–51, C iii 5696, and C iv 5801–5812 lines since their formation process depends on fine details of atomic physics and on the modeling (Martins & Hillier 2012). Moreover, special care is taken for the helium abundance. Some of the helium lines are also affected by the microturbulent velocity taken to compute the emerging spectrum (McErlean et al. 1998). A compromise is thus made to ensure the best CMFGEN fit. The individual parameters of both components are given in Table 3 and we display, in Figs. 5 and 6, the CMFGEN best-fit model for the primary and the secondary, respectively.

Primary Secondary
[M yr]
 [km s]
0.07 0.07
[km s] 30 30
1.0 1.0
 [km s]
 [km s]
[km s]
444(He/H), (C/H), (N/H) and (O/H) in number. indicates how rapidly the clumping should be switched on.
Table 3: Individual parameters of each component of HD 166734

The individual parameters agree with two supergiant objects according to the tables of Martins et al. (2005), and Bouret et al. (2012). The best-fit is achieved with  K and 30500 K and and 3.10 for the primary and the secondary, respectively. We also derive an overabundance of helium and nitrogen as well as a depletion of carbon in both objects, as expected from their luminosity class. Furthermore, we emphasize that the oxygen is depleted for the primary whilst it is solar for the secondary.

The wind parameters are reminiscent of those of Galactic supergiant O-type stars, according to Bouret et al. (2012) even though stars as cool as the two components of HD 166734 were not analyzed by these authors. These parameters are thus in agreement with those of single massive supergiant objects.

6 Wind-wind interactions

The spectral disentangling has yielded the mean spectrum of each component. From these spectra, the emission observed in the He ii 4686 and the H lines seems to be linked to the primary component over the entire orbital cycle. Our CMFGEN analysis shows that this component has a larger mass-loss rate than the secondary, whilst its terminal velocity is assumed to be similar.

Figure 7: Top: Positions of the maximum of H emission as a function of the orbital phase. The red dots represent the positions of the maximum, whilst the blue line is the radial velocity curve of the primary component. Bottom: Spectral variation of the H line as a function of the orbital phase. The emission vanishes at .

Upon inspection of the H profile in the observed spectra, it is clear that the maximum of emission of this line is associated to the primary on the whole orbit (see the upper panel of Fig. 7). A closer look at this profile at different phases (see the lower panel of Fig. 7) shows that at the periastron () no emission is observed in the spectra, but this emission is present at and even at . During this short time interval, when the two components are closest to each other, it appears that the emission generated by the two winds vanishes.

The configuration of the system indicates that at the periastron the distance between the centers of both components is 71.2 . Therefore, the distance that remains between the two photospheres (of about 16.9 ) is not large enough for the winds to strongly collide.

7 Discussion

7.1 Stellar evolution

The stellar luminosities determined from photometry as well as the effective temperatures derived from the CMFGEN modeling provide the exact locations of both stars in the Hertzsprung-Russell diagram. They are displayed in Fig 8. Through our analysis, the orbital parameters of HD 166734 as well as the stellar parameters and the derived surface abundances of its two components suggest that they have evolved mainly as single stars (at least without Roche lobe overflow event). They are thus interesting objects to be compared with the current evolutionary tracks for massive stars.

The input physics and its implementation in six different evolutionary codes as well as their comparison with observed data were deeply discussed by Martins & Palacios (2013). In the present discussion, we compare the observational parameters of the two components of HD 166734 with the tracks of Geneva (Ekström et al. 2012) and with those of Brott et al. (2011, referred as STERN models by and hereafter in this paper).

7.1.1 Geneva models

The evolutionary tracks of Ekström et al. (2012) are computed with a rotational rate of 0.4, a metallicity of , an overshoot parameter of 0.1 and do not account for magnetism. The upper panel of Fig. 8 displays the evolutionary tracks and the isochrones of Geneva computed with rotation. Given the error bars on the luminosity and on the effective temperature, the two components are located, according to these tracks, close to the characteristic hook at the end of the main sequence. These positions therefore infer two different values for the initial mass of the two components. Indeed, within the error bars, the stars can be located either at the end of the main sequence, immediately before the stage where their temperature starts to increase, or before the yellow supergiant phase. Depending on their exact locations, the initial masses will be expected to be 44 or 47  for the primary and 37 or 39  for the secondary. According to the theoretical predictions, these stars would have ages between 5.2 and 5.4 Myrs for the primary and between 5.8 and 6.0 Myrs for the secondary (depending on the locations of both stars on the hook at the end of the main sequence). Furthermore, these tracks yield a radius of 27  for the primary, and of 26  for the secondary.

We also determine the present evolutionary masses of each component. For the primary we compute a mass of  whilst for the secondary we obtain a mass of . These masses can be compared to the spectroscopic masses computed from the luminosity, the effective temperature and the gravity of each star to be  for the primary and  for the secondary as well as to the dynamical masses of  and , for the primary and the secondary, respectively. The large rotational rate () is however not realistic for massive stars in the Milky Way (Simón-Díaz & Herrero 2014). This infers an overestimation of the ages of the stars, and an underestimation of their evolutionary masses. From the tracks of Geneva, the two components of HD 166734 are thus located close to the terminal age main sequence (TAMS) but given the fast evolution at the end of the main sequence, this part of the HR diagram is expected to be empty (McEvoy et al. 2015). The inclusion of extra mixing (either by tidal effects at periastron or by increasing the overshoot parameter) in the evolutionary models would allow to increase the lifetime of the stars on the main sequence. Even though the initial rotational rate is large for the Geneva models, the spin down is quick. Indeed, for our individual stars taken as single, the actual equatorial velocities are predicted to be 15 km s for the primary and 18 km s for the secondary, which is too small with respect to the rotational velocities of our stars. The spin down is notably ruled by the mass-loss rate included in the models, suggesting that the predictions concerning the mass-loss in the Geneva models seems to be overestimated or that another phenomenon is at work.

The initial rotational velocity of about 300 km s seems therefore too large for both components. We thus quickly compare the positions of the stars to the Geneva evolutionary tracks computed without rotation. From these models, the stars are on the main sequence, with estimated ages of about 3.0 Myrs for the primary and 3.4 Myrs for the secondary. The primary has an initial mass of 59  and an actual mass of 51.4  whilst the secondary has an initial mass of 50  and an actual mass of 45.3 ; thus larger than what we found for the two components from evolutionary tracks with rotation as expected from the theory.

7.1.2 STERN models

We decided to compare our observational values with the STERN evolutionary models (Brott et al. 2011). These models constitute a grid with a sufficiently fine mesh, attenuating the errors due to the interpolation. This contrasts with the Geneva models for which one rotational rate is taken into account. The STERN models are computed with a metallicity of , an overshoot parameter of 0.335 and the rotation is taken as a free parameter. The BONNSAI facility555The BONNSAI web-service is available at http://www.astro.uni-bonn.de/stars/bonnsai (Schneider et al. 2014) was used as a Bayesian approach to derive the distribution of stellar parameters for each component of HD 166734. The best fits provide an initial rotational velocity of 150 km s for both objects, giving an equatorial rotational velocity of 110 km s at their current stage (in adequation with our observational parameters). The lower panel of Fig. 8 displays the locations of the stars on the HR diagram build from the STERN evolutionary tracks and the isochrones computed with an initial rotational velocity of 160 km s. For the primary, the initial mass is expected to be 56.1, giving an evolutionary mass of 47.8. For the secondary, the initial mass is 47.4and the evolutionary mass is 41.2. The STERN models provide an age of about 3.0 Myrs for the primary and of about 3.5 Myrs for the secondary; in agreement with the non-rotating tracks of Geneva. As emphasized by Martins & Palacios (2013), the characteristic hook at the end of the main sequence occurs at lower temperatures because of the large amount of overshooting. These authors also mentioned that an overshooting parameter of 0.335, as used for the computation of the STERN models, is too large for stars with masses above 20 . Therefore, the ages estimated from these models are probably underestimated and their evolutionary masses overestimated. However, the STERN tracks better reproduce the rotational velocities of the supergiants.

Figure 8: Positions in the HR diagram of the two components of HD 166734. Top: Evolutionary tracks (black lines) and isochrones (gray lines) from Ekström et al. (2012), computed with a rotational rate of 300 km s. Bottom: Evolutionary tracks (black lines) and isochrones (gray lines) from Brott et al. (2011), computed with a rotational rate of 160 km s. The error bars represent the 1-.

7.1.3 Surface abundances

Whatever the evolutionary tracks we use, the secondary component appears more evolved than the primary. An enhanced mixing through tides at periastron could explain this difference in age (see Sect. 7.2). As we already mentioned in Section 5, the primary and the secondary exhibit an overabundance in helium and in nitrogen as well as a depletion in carbon. Furthermore, the primary is also depleted in oxygen. The abundances thus confirm the evolved status of each object. However, they do not differ from abundances of single evolved massive stars. When we compare the surface abundances of the two components of HD 166734 to the values of only the supergiant stars determined by Martins et al. (2015) in the framework of the MiMeS project and to the values derived by Bouret et al. (2012) in their analysis of a sample of Galactic supergiant stars, we see a clear correlation between N/O and N/C (see the upper panel of Fig. 9) that does not differ from the results of Bouret et al. (2012) or of Martins et al. (2015), represented in the upper panel of Fig. 9 by black dots. This agrees with the expectations of nucleosynthesis through the CNO cycle. The ratios N/C and N/O follow those of single stars according to the evolutionary models, within the uncertainties. The lower panel of Fig. 9 shows that the ratio N/C increases as surface gravity decreases, as was observed for single supergiant Galactic massive stars (Bouret et al. 2012; Martins et al. 2015). The ratios of N/C and N/O determined from the two components of HD 166734 are thus correlated and in the same range as those measured for presumably single supergiant objects, meaning that no or at least few interactions occur between these two components.

Figure 9: Top: (N/C) (by number) as a function of (N/O). The expected trends for the case of CN and CNO equilibrium are shown by the dashed lines. Bottom: (N/C) (by number) as a function of . The two components of HD 166734 are represented by red dots, whilst the sample of Galactic supergiant stars analyzed by Bouret et al. (2012) and Martins et al. (2015) is represented by black dots. The evolutionary tracks are from Ekström et al. (2012) and are computed for a rotational rate (initial over critical) of 0.4. We use evolutionary tracks with initial masses larger than 10 .

7.2 Synchronization of the system

HD 166734 is a highly-eccentric binary system where the secondary appears to be more evolved than the primary. This problem linked to the ages can be due to tidal mixing at periastron. This tidal mixing would affect the secondary more than the primary because of its size. In addition, some matter could be transferred to the primary with the consequence of increasing its spin, inducing more mixing and therefore a higher temperature; this increase of the temperature would modify the apparent age of the primary by becoming younger. This effect should be however limited because of the configuration of the system.

We then may wonder whether these tidal interactions (that would be responsible for this tidal mixing) can affect the rotations of the components and can possibly synchronize the system (tidal interactions can indeed be responsible for a long timescale evolution of the binary orbital elements and the stellar rotation). To characterize these interactions in HD 166734, we focus on the theory of Hut (1981, 1982). If we compute the ratio between the rotational and the orbital angular momentum (Eq. 14 in Hut 1981), we find for the present stars of HD 166734 that a large part of the angular momentum is contained in the orbit. An important classification parameter () representing the ratio of the orbital angular momentum to the rotational angular momentum at equilibrium was introduced by Hut (1981). According to the value of this parameter, the circularization will take much more time than the synchronization. In the case of HD 166734, for which the orbital separation between the two components is relatively large, the parameter is extremely large, meaning that the circularization time is expected to be longer than the synchronization time. Always according to Hut’s theory, the eccentricity of the orbit should decrease very slowly as well as the semimajor axis ( and ) whilst the rotational angular velocities of both stars should increase as a function of time () but all this evolution is expected on timescales longer than the current age of the system. In eccentric binaries, the term of synchronization between the rotation and the revolution is rather replaced by the pseudo-synchronization at the moment of the periastron. This pseudo-synchronization of the system is reached when the rotational angular velocity of a star is at times the instantaneous orbital angular velocity at periastron (see Hut 1981). Indeed, assuming that the rotational axes of both components are perpendicular to the orbital plane, we use the combination of the inclination () of the system, the equatorial rotational velocities ( and  km s) and the radii ( and ) of both components, obtained in Section 4, to compute rotational periods of 13.0 days for the primary and of 12.3 days for the secondary. Both stars appear to have their rotations almost synchronized with each other but these values are smaller than the orbital period. When we compare the rotational angular velocities with the instantaneous orbital angular velocity at periastron, we compute ratios of about 0.6 for the primary and secondary components. We therefore conclude that pseudo-synchronization (Hut 1981) has not yet been achieved for HD 166734.

7.3 Distance to HD 166734

In Section 4, we determined the bolometric magnitudes of for the primary and for the secondary, giving stellar luminosities of and , respectively. Given their spectral types and their effective temperatures, we can compute, from Martins & Plez (2006), bolometric corrections of and for the primary and the secondary, respectively. From these values, we compute absolute magnitudes of for the primary and for the secondary.

From the observed magnitude of HD 166734 of 8.42 and its of 1.07 (Ducati 2002), accounting for (Martins & Plez 2006), and a brightness ratio of 1.28 between the primary and the secondary components, we estimate, for HD 166734, a distance modulus of 11.85 and thus a distance of  kpc, confirming the distance estimated by Conti et al. (1980).

Figure 10: Section of the FEROS spectrum of HD 166734 showing the Na i interstellar lines, covering a velocity range of to  km s. The dotted lines indicate, from left to right, RVs of , , ,  km s for both panels.

Next, we compare this distance with the kinematic distance estimated from the Galactic rotation curve. In Fig. 10, we show the Na i absorption line profiles. The Na i lines are relatively saturated on a width ranging from to  km s and with a shallower absorption profile between and  km s. From the Galactic rotation model of Fich et al. (1989), this velocity range gives a mean distance of  kpc. This distance is in agreement with the distances estimated from the stellar luminosities or derived by Conti et al. (1980).

8 Conclusions

We present the results of the spectroscopic and photometric campaigns devoted to HD 166734. The spectra taken close to the periastron passage have allowed us to refine the orbital solution of the system, removing the existing doubts on the masses of both components. Our analysis has indeed shown that HD 166734 is a binary system with an orbital period of 34.537723 days and an eccentricity of 0.618. From the photometry, we determine an inclination of for the system and actual masses of  and  for the primary and the secondary.

The individual spectra of each component, obtained by spectral disentangling, indicate a primary with a spectral classification of O7If, an effective temperature of  K, a  of and a luminosity of . We found that the secondary has the spectral type O9I(f) and its stellar parameters are  K, , and . We thus conclude, unlike Conti et al. (1980), that the hottest and the most luminous star is also the most massive one. Both components show modified helium, carbon, nitrogen and oxygen surface abundances. They are both enriched in helium and nitrogen and are depleted in carbon. In addition, the primary is also depleted in oxygen. Their surface abundances do not however differ from the predictions of the ‘classical’ evolutionary paths for single stars. This means that no or, at most, few interactions occur between the two stars, making them perfectly suitable as test objects for evolutionary models.

Systems containing a pair of massive evolved O-type stars are relatively rare. Despite their number, they represent excellent tests for atmospheric and evolutionary models. Analyses of the kind we did for HD 166734 are thus crucial to validate these models. In this context, we may quote other systems that can help with this validation. An interesting system is R139 (Taylor et al. 2011) composed of an O6Iafc primary and of an O6Iaf secondary on a wide eccentric orbit. This system is the most massive binary system known to host two evolved Of supergiants but the physical parameters of each component are not complete enough to allow a deep comparison with the evolutionary models. Finally, WR21a (Tramper et al. 2016; Gosset & Nazé 2016) is a very massive system composed of two younger O3/WN5 and O3 stars. Although this system is less evolved, it has orbital period ( days) and eccentricity () comparable to HD 166734.

In the present paper, we compare the fundamental parameters of each component with evolutionary tracks of Geneva (Ekström et al. 2012) and of STERN (Brott et al. 2011). We conclude from our observations that the tracks of Geneva have an overly large initial rotational rate, that shortens the lifetime of the star on the main sequence and it seems that the mass-loss used to compute the tracks is also too large to reproduce the actual rotational velocities of the stars. At the same time, the STERN tracks have an overly large overshoot parameter, but the spin down of the stars is sufficient to agree with the actual rotational velocities of the stars. Finally, the ages and the evolutionary masses for both components are different as a function of the tracks used in the analysis. Even though other systems (such as those quoted above) must still be investigated, HD 166734 allowed us to show that the evolutionary models are not reliable enough to reproduce the exact characteristics of supergiant stars.

The authors thank the anonymous referee for his/her constructive remarks on the paper. This research was supported by the Fonds National de la Recherche Scientifique (F.R.S.-F.N.R.S.), by the PRODEX XMM contract (Belspo), and through the ARC grant for Concerted Research Actions, financed by the French Community of Belgium (Wallonia-Brussels Federation). L.M. thanks Pavel Mayer for his advice in the light curve modeling and John Hillier for making CMFGEN available. We are grateful to the staff of La Silla ESO Observatory for their technical support.


  • Bouret et al. (2012) Bouret, J.-C., Hillier, D. J., Lanz, T., & Fullerton, A. W. 2012, A&A, 544, A67
  • Brott et al. (2011) Brott, I., de Mink, S. E., Cantiello, M., et al. 2011, A&A, 530, A115
  • Conti (1973) Conti, P. S. 1973, ApJ, 179, 161
  • Conti & Alschuler (1971) Conti, P. S. & Alschuler, W. R. 1971, ApJ, 170, 325
  • Conti et al. (1977) Conti, P. S., Leep, E. M., & Lorre, J. J. 1977, ApJ, 214, 759
  • Conti et al. (1980) Conti, P. S., Massey, P., Ebbets, D., & Niemela, V. S. 1980, ApJ, 238, 184
  • Crowther & Evans (2009) Crowther, P. A. & Evans, C. J. 2009, A&A, 503, 985
  • Damerdji et al. (2007) Damerdji, Y., Klotz, A., & Boër, M. 2007, AJ, 133, 1470
  • Ducati (2002) Ducati, J. R. 2002, VizieR Online Data Catalog, 2237
  • Ekström et al. (2012) Ekström, S., Georgy, C., Eggenberger, P., et al. 2012, A&A, 537, A146
  • Fich et al. (1989) Fich, M., Blitz, L., & Stark, A. A. 1989, ApJ, 342, 272
  • González & Levato (2006) González, J. F. & Levato, H. 2006, A&A, 448, 283
  • Gosset & Nazé (2016) Gosset, E. & Nazé, Y. 2016, A&A, 590, A113
  • Gosset et al. (2001) Gosset, E., Royer, P., Rauw, G., Manfroid, J., & Vreux, J.-M. 2001, MNRAS, 327, 435
  • Grevesse et al. (2007) Grevesse, N., Asplund, M., & Sauval, A. J. 2007, Space Sci. Rev., 130, 105
  • Hadrava (1995) Hadrava, P. 1995, A&AS, 114, 393
  • Heck et al. (1985) Heck, A., Manfroid, J., & Mersch, G. 1985, A&AS, 59, 63
  • Herrero et al. (1992) Herrero, A., Kudritzki, R.-P., Vilchez, J. M., et al. 1992, in Lecture Notes in Physics, Berlin Springer Verlag, Vol. 401, The Atmospheres of Early-Type Stars, ed. U. Heber & C. S. Jeffery, 21
  • Hillier & Miller (1998) Hillier, D. J. & Miller, D. L. 1998, ApJ, 496, 407
  • Hut (1981) Hut, P. 1981, A&A, 99, 126
  • Hut (1982) Hut, P. 1982, A&A, 110, 37
  • Lamers & Leitherer (1993) Lamers, H. J. G. L. M. & Leitherer, C. 1993, ApJ, 412, 771
  • Lucy & Sweeney (1971) Lucy, L. B. & Sweeney, M. A. 1971, AJ, 76, 544
  • Mahy et al. (2011) Mahy, L., Gosset, E., Baudin, F., et al. 2011, A&A, 525, A101
  • Mahy et al. (2015) Mahy, L., Rauw, G., De Becker, M., Eenens, P., & Flores, C. A. 2015, A&A, 577, A23
  • Marchenko et al. (1998) Marchenko, S. V., Moffat, A. F. J., & Eenens, P. R. J. 1998, PASP, 110, 1416
  • Martins et al. (2015) Martins, F., Hervé, A., Bouret, J.-C., et al. 2015, A&A, 575, A34
  • Martins & Hillier (2012) Martins, F. & Hillier, D. J. 2012, A&A, 545, A95
  • Martins & Palacios (2013) Martins, F. & Palacios, A. 2013, A&A, 560, A16
  • Martins & Plez (2006) Martins, F. & Plez, B. 2006, A&A, 457, 637
  • Martins et al. (2005) Martins, F., Schaerer, D., & Hillier, D. J. 2005, A&A, 436, 1049
  • Mathys (1988) Mathys, G. 1988, A&AS, 76, 427
  • Mathys (1989) Mathys, G. 1989, A&AS, 81, 237
  • McErlean et al. (1998) McErlean, N. D., Lennon, D. J., & Dufton, P. L. 1998, A&A, 329, 613
  • McEvoy et al. (2015) McEvoy, C. M., Dufton, P. L., Evans, C. J., et al. 2015, A&A, 575, A70
  • Montgomery & O’Donoghue (1999) Montgomery, M. H. & O’Donoghue, D. 1999, Delta Scuti Star Newsletter, 13, 28
  • Otero & Wils (2005) Otero, S. A. & Wils, P. 2005, Information Bulletin on Variable Stars, 5644
  • Pavlovski & Hensberge (2005) Pavlovski, K. & Hensberge, H. 2005, A&A, 439, 309
  • Pavlovski & Hensberge (2010) Pavlovski, K. & Hensberge, H. 2010, in Astronomical Society of the Pacific Conference Series, Vol. 435, Binaries - Key to Comprehension of the Universe, ed. A. Prša & M. Zejda, 207
  • Pavlovski & Southworth (2012) Pavlovski, K. & Southworth, J. 2012, in IAU Symposium, Vol. 282, From Interacting Binaries to Exoplanets: Essential Modeling Tools, ed. M. T. Richards & I. Hubeny, 359–364
  • Prša & Zwitter (2005) Prša, A. & Zwitter, T. 2005, ApJ, 628, 426
  • Rauw et al. (2000) Rauw, G., Sana, H., Gosset, E., et al. 2000, A&A, 360, 1003
  • Sana et al. (2012) Sana, H., de Mink, S. E., de Koter, A., et al. 2012, Science, 337, 444
  • Sana et al. (2006) Sana, H., Gosset, E., & Rauw, G. 2006, MNRAS, 371, 67
  • Sana et al. (2014) Sana, H., Le Bouquin, J.-B., Lacour, S., et al. 2014, ApJS, 215, 15
  • Scargle (1982) Scargle, J. D. 1982, ApJ, 263, 835
  • Schneider et al. (2014) Schneider, F. R. N., Langer, N., de Koter, A., et al. 2014, A&A, 570, A66
  • Simón-Díaz & Herrero (2007) Simón-Díaz, S. & Herrero, A. 2007, A&A, 468, 1063
  • Simón-Díaz & Herrero (2014) Simón-Díaz, S. & Herrero, A. 2014, A&A, 562, A135
  • Taylor et al. (2011) Taylor, W. D., Evans, C. J., Sana, H., et al. 2011, A&A, 530, L10
  • Tramper et al. (2016) Tramper, F., Sana, H., Fitzsimons, N. E., et al. 2016, MNRAS, 455, 1275
  • Underhill (1995) Underhill, A. B. 1995, ApJS, 100, 433
  • Vink et al. (2001) Vink, J. S., de Koter, A., & Lamers, H. J. G. L. M. 2001, A&A, 369, 574
  • Walborn (1971) Walborn, N. R. 1971, ApJS, 23, 257
  • Walborn & Fitzpatrick (1990) Walborn, N. R. & Fitzpatrick, E. L. 1990, PASP, 102, 379
  • Wilson & Devinney (1971) Wilson, R. E. & Devinney, E. J. 1971, ApJ, 166, 605
  • Wolfe et al. (1967) Wolfe, Jr., R. H., Horak, H. G., & Storer, N. W. 1967, The machine computation of spectroscopic binary elements, ed. O. Struve & M. Hack, 251

Appendix A Journal of the spectroscopic observations of HD 166734

HJD–2 450 000 Instrument RV RV HJD–2 450 000 Instrument RV RV
2201.6084 ESPRESSO 0.189 94.7 -102.5 3643.5890 FEROS 0.940 -186.1 199.6
2201.6311 ESPRESSO 0.190 91.0 -103.0 3643.6064 FEROS 0.941 -190.1 198.6
2201.6636 ESPRESSO 0.191 93.3 -101.2 3797.9001 FEROS 0.408 52.4 -44.3
2202.5888 ESPRESSO 0.218 85.9 -88.3 3800.9080 FEROS 0.495 37.3 -29.4
2202.6127 ESPRESSO 0.219 84.4 -86.2 3862.8646 FEROS 0.289 57.2 -85.5
2202.6560 ESPRESSO 0.220 79.5 -85.4 3864.9075 FEROS 0.348 43.5 -83.1
2381.9031 FEROS 0.410 41.1 -63.3 3867.9057 FEROS 0.435 37.7 -30.5
2381.9212 FEROS 0.410 35.3 -60.3 3867.9232 FEROS 0.436 37.3 -24.2
2382.9269 FEROS 0.439 37.8 -59.8 3874.8779 FEROS 0.637 12.1 -7.2
2782.8462 FEROS 0.019 -67.2 77.3 3874.8953 FEROS 0.638 7.2 9.1
2783.7920 FEROS 0.046 47.8 -76.0 3880.7990 FEROS 0.809 -75.1 95.1
2784.9020 FEROS 0.078 63.7 -87.5 3880.8165 FEROS 0.809 -73.3 95.6
3130.7969 FEROS 0.093 88.9 -94.8 3899.7138 FEROS 0.356 42.7 -79.7
3131.8345 FEROS 0.123 93.4 -118.2 3899.7312 FEROS 0.357 45.4 -82.3
3132.8529 FEROS 0.153 89.7 -103.5 3912.8945 FEROS 0.738 -35.5 49.4
3133.8712 FEROS 0.182 89.1 -108.5 3913.8940 FEROS 0.767 -56.9 62.5
3134.8491 FEROS 0.210 82.8 -79.7 3914.6723 FEROS 0.789 -59.0 67.9
3288.6241 ESPRESSO 0.663 -5.1 6.6 3918.5369 FEROS 0.901 -149.4 157.2
3288.6393 ESPRESSO 0.663 -4.8 9.6 3918.5543 FEROS 0.902 -152.6 157.1
3288.6539 ESPRESSO 0.664 2.1 -3.2 3919.5103 FEROS 0.929 -179.5 191.4
3289.6186 ESPRESSO 0.692 -31.3 72.1 3919.5277 FEROS 0.930 -179.9 193.2
3289.6364 ESPRESSO 0.692 -37.8 58.5 3950.6415 FEROS 0.831 -88.7 95.8
3289.6513 ESPRESSO 0.693 -26.6 54.6 3950.6589 FEROS 0.831 -87.4 97.7
3291.6744 ESPRESSO 0.751 -43.7 71.7 3951.5085 FEROS 0.856 -107.5 117.9
3291.6890 ESPRESSO 0.752 -49.0 63.1 3951.5304 FEROS 0.856 -106.5 117.4
3291.7036 ESPRESSO 0.752 -45.9 63.3 3952.7421 FEROS 0.892 -137.8 157.9
3509.8746 FEROS 0.069 49.8 -82.9 3952.7595 FEROS 0.892 -140.9 159.6
3511.9189 FEROS 0.128 85.5 -113.3 3953.5369 FEROS 0.915 -160.9 182.8
3512.8354 ESPRESSO 0.155 81.2 -90.9 3953.5543 FEROS 0.915 -161.4 181.6
3512.8482 ESPRESSO 0.155 84.4 -91.8 3954.5829 FEROS 0.945 -181.0 209.6
3512.8653 ESPRESSO 0.155 85.3 -90.6 3954.6004 FEROS 0.945 -181.1 210.9
3512.8804 ESPRESSO 0.156 86.2 -94.8 4193.9278 ESPRESSO 0.875 -113.6 154.3
3512.9230 FEROS 0.157 99.8 -99.3 4193.9416 ESPRESSO 0.875 -104.1 162.8
3513.8269 ESPRESSO 0.183 76.2 -102.5 4194.9639 ESPRESSO 0.905 -138.2 177.0
3513.8443 ESPRESSO 0.184 73.3 -93.2 4194.9788 ESPRESSO 0.905 -130.8 183.8
3513.8590 ESPRESSO 0.184 74.7 -85.4 4194.9951 ESPRESSO 0.906 -134.6 181.1
3515.8300 ESPRESSO 0.241 84.9 -86.9 4196.9348 ESPRESSO 0.962 -195.0 233.4
3515.8472 ESPRESSO 0.242 62.6 -84.9 4196.9570 ESPRESSO 0.963 -195.7 232.2
3515.8619 ESPRESSO 0.242 66.9 -87.0 4196.9760 ESPRESSO 0.963 -198.2 226.2
3516.8296 ESPRESSO 0.270 61.4 -81.5 4197.8627 ESPRESSO 0.989 -180.8 198.8
3516.8445 ESPRESSO 0.271 51.1 -80.2 4197.8852 ESPRESSO 0.989 -178.5 195.9
3516.8592 ESPRESSO 0.271 50.9 -82.1 4197.9099 ESPRESSO 0.990 -176.1 193.4
3517.8129 ESPRESSO 0.299 60.6 -67.5 4199.8838 ESPRESSO 0.047 42.1 -47.5
3517.8275 ESPRESSO 0.299 60.6 -73.6 4199.9010 ESPRESSO 0.048 38.2 -49.8
3517.8443 ESPRESSO 0.300 60.0 -74.7 4199.9157 ESPRESSO 0.048 38.2 -43.5
3550.6987 FEROS 0.251 75.5 -99.5 4249.9649 ESPRESSO 0.497 26.2 -17.9
3550.7161 FEROS 0.251 76.9 -99.7 4249.9768 ESPRESSO 0.498 26.9 -19.1
3552.8692 FEROS 0.314 56.7 -70.2 4250.8564 ESPRESSO 0.523 24.5 -30.7
3552.8852 FEROS 0.314 53.3 -73.5 4250.8863 ESPRESSO 0.524 23.0 -24.5
3588.6468 FEROS 0.350 51.6 -69.6 4251.7829 ESPRESSO 0.550 7.9 0.1
3588.6642 FEROS 0.350 52.8 -70.9 4251.8100 ESPRESSO 0.551 3.8 0.2
3597.6673 FEROS 0.611 4.9 -7.3 4251.8433 ESPRESSO 0.552 6.1 0.2
3597.6847 FEROS 0.611 6.1 -8.7 4982.8318 ESPRESSO 0.717 -36.9 40.9
3601.6977 FEROS 0.727 -24.2 32.2 4983.7593 ESPRESSO 0.743 -11.0 81.1
3601.7151 FEROS 0.728 -24.7 31.2 4984.7515 ESPRESSO 0.772 -42.1 76.8
3617.5842 FEROS 0.187 87.8 -96.2 4985.8245 ESPRESSO 0.803 -59.3 97.0
3617.6017 FEROS 0.188 88.3 -95.5 4986.8138 ESPRESSO 0.832 -68.3 106.2
3620.5870 FEROS 0.274 57.8 -88.1 4987.8062 ESPRESSO 0.861 -109.4 141.2
3620.6044 FEROS 0.275 57.9 -87.4 4988.7544 ESPRESSO 0.888 -133.0 163.0
3625.6001 FEROS 0.420 41.9 -70.1 4989.7606 ESPRESSO 0.917 -160.6 186.2
3625.6175 FEROS 0.420 35.5 -70.0 4992.7391 ESPRESSO 0.003 -134.9 113.5
3628.6180 FEROS 0.507 25.0 -12.0 5021.7712 FEROS 0.844 -88.0 110.3
3628.6355 FEROS 0.507 22.7 -14.6 5025.6253 FEROS 0.956 -187.1 217.4
3640.5021 FEROS 0.851 -105.4 123.6 5028.6212 FEROS 0.042 15.9 3.5
3640.5196 FEROS 0.852 -105.4 123.7
Table 4: Journal of the observations of HD 166734.

Appendix B Periodograms

Figure 11: Top: Power spectrum computed on the basis of RV–RV. Bottom: Spectral window computed from the spectroscopic observations.
Figure 12: Top: Power spectrum of the R light curve of HD 166734. The diamond indicates the location of the orbital frequency while the dots represent its harmonics. Bottom: Spectral window computed from the R light curve.
Figure 13: Top: Power spectrum of the I light curve of HD 166734. The diamond indicates the location of the orbital frequency while the dots represent its harmonics. Bottom: Spectral window computed from the I light curve.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description