Free Access
Issue
A&A
Volume 570, October 2014
Article Number A3
Number of page(s) 9
Section Stellar structure and evolution
DOI https://doi.org/10.1051/0004-6361/201322888
Published online 06 October 2014

© ESO, 2014

1. Introduction

Recent detection of yellow supergiant star as the possible progenitor of supernovae have raised serious questions about our understanding of the evolution of massive stars. Bersten et al. (2012) and Benvenuto et al. (2013) present recent evidence that the yellow supergiant star found at the location of supernova SN 2011dh in the famous nearby galaxy M51 was the SN progenitor. The authors of these papers have produced a self-consistent model to explain how the star exploded. In their model, the exploding massive star was a member of a close binary system. A comprehensive analysis of similar binary systems will advance our knowledge of binary systems as possible progenitors of supernova explosions.

The spectrum of HD 50975 was classified as a composite spectrum of two stars, F8 Ib supergiant and an early type star, using the objective prism with dispersion of 108 Å mm-1 at Hγ by MacConnell & Bidelman (1976). The same spectral type was found for a primary star by Barbier et al. (1978) on the basis of spectra with dispersion of 200 Å mm-1. Barbier et al. (1978) analyzed photometric data gathered using the ultraviolet satellite TD1 and noted that HD 50975 is an “unexplained case”. HD 50975 is located in the region of G4 subdwarfs according to uvby photometry (Olsen 1979), and low-amplitude light variations are suspected by Olsen (1983). The period (P) and amplitude (A) of light variations were determined by Koen & Eyer (2002), P = 95.1 days and A = 0.013 mag. Radial velocity variations were detected for HD 50975 by Bartkevicius & Sperauskas (1994).

In the present paper a comprehensive analysis and modeling of the collected long-term radial velocity measurements, photometric data, spectral energy distribution (SED), cross-correlation function (CCF), and high-resolution spectra was performed for HD 50975 to calculate orbital elements, atmospheric parameters, abundances, luminosities, and masses. A spectrum in the ultraviolet and blue wavelength region was analyzed in detail to clarify the nature of unseen companion.

2. Observations

Radial-velocity monitoring of HD 50975 was started in 1989 using the CORAVEL spectrometer of the Vilnius University installed on the 1.65 m and 0.63 m telescopes at the Moletai Observatory (Lithuania). The CORAVEL spectrometer (Upgren et al. 2002) is based on the principles of the photoelectric radial-velocity scanner developed by Griffin (1967) and operates by scanning a spectrum of star across the mask and obtaining an online cross-correlation velocity. The standard deviation of a single observation for late-type stars brighter than about the 11th magnitude is usually better than 0.8 km s-1. The velocities have been standardized using the observations of IAU radial velocity standards (Udry et al. 1999) and are close to the system of velocities published by Nidever et al. (2002). The difference in zero point was found to be 0.14 km s-1 for F-G-K type stars and exhibit an rms scatter of 0.5 km s-1.

High-resolution spectra for HD 50975 were obtained with the coudé échelle spectrometer MAESTRO on the 2 m telescope at the Observatory on the Terskol Peak in Northern Caucasus (altitude of 3100 m) equipped with a Wright Instruments CCD detector. Three spectra were obtained between November 2012 and January 2013. In addition, spectra of several hot, rapidly rotating stars were observed to identify the telluric absorption lines. All the spectra covered region from about 3600 to 9600 Å and have a resolving power R ~ 45 000. Unfortunately, our efforts to observe a high-resolution spectrum shortward of 3600 Å was unsuccessful. All spectra were bias-subtracted, flat-field-corrected, and converted to one-dimensional spectra using the standard DECH20T package1. The equivalent widths of absorption lines were measured by Gaussian fitting to the observed profiles.

Table 1

Radial velocity measurements for HD 50975A, along with the uncertainties, phases, and residuals for the orbital curve.

3. Analysis and results

3.1. Orbital solution and light curve

Orbital parameters for HD 50975 were calculated using 53 radial velocities observed in the time span from February 1989 to April 2008. The measured heliocentric radial velocities and their errors are presented in Table 1, together with the time of observations, phases, and residuals (O–C). Orbital elements are listed in Table 2. All of the radial velocities upon which the orbit is based have been attributed the same (unit) weight. Radial velocity curve for HD 50975A is shown in Fig. 1 (upper panel). In addition, two light curves are plotted in the middle and bottom panels. Published photometric data are phased with the orbital period of P = 190.22 d and Tmax of radial velocity as an epoch for phase = 0. The photometric data were retrieved from the Hipparcos and ASAS databases: 74 measurements taken from the Hipparcos epoch photometry (Perryman & ESA 1997) and 410 measurements from the ASAS-3 V band photometry (Pojmanski 2002). We used ASAS photometry observed by a two-pixel aperture because of the lowest reduction errors and scatter of measurements. In addition, the measurements of D (uncertain) quality and crude errors (3σ threshold) were rejected. The remaining 389 measurements were phased with the orbital period and are plotted in Fig. 1 (middle panel). Double-peaked light curve synchronized with the maximum and minimum of radial velocity variations gives evidence about the ellipsoidal shape of the primary star.

thumbnail Fig. 1

Observed heliocentric radial velocities for HD 50975A (filled circles) as a function of phase, with the velocity curve corresponding to the calculated orbital elements (upper panel). ASAS (middle panel) and Hipparcos (lower panel) photometry phased using the same orbital period gives evidence of the ellipsoidal shape of the primary star.

Table 2

The orbital elements along with standard deviations for HD 50975A.

3.2. IUE spectra

The spectrum of the hot secondary star dominates the composite spectrum of HD 50975 in the far ultraviolet spectral region. Low-dispersion spectra observed by SWP and LWR cameras (SWP 04082 and LWR 03613) are available for HD 50975 in the IUE archive. We used IUE NEWSIPS data recalibrated by Massa & Fitzpatrick (2000) in order to minimize systematic errors. To avoid the 2200 Å interstellar extinction bump and the possible influence of primary star HD 50975A we employed only spectra observed by SWP camera. In addition, spectra of four comparison stars of spectral types between B1 and B4 were extracted from the database. Normalized spectra are plotted in Fig. 2. Criteria of spectral classification developed by Rountree & Sonneborn (1991) in the far ultraviolet region are used for spectral classification of the secondary star HD 50975B. On the basis of temperature sensitive features of Si ii at 1260 Å, Si iii at 1300 Å, C ii at 1340 Å, and Si iv at 1400 Å, the temperature subtype was found to be B2 ± 0.5 for HD 50975B.

The profile of hydrogen Lyman α line was synthesized using code STARSP and Kurucz (1993) models (see Sect. 3.6.3) for the set of effective temperatures between 15 000 and 25 000 K. The best fit between the observed and calculated profiles (see Fig. 2) was found for the atmospheric model with an effective temperature, Teff = 17 000 K, and gravity, log (g) = 4.5 (cgs), in support of spectral type ~B3 for HD 50975B.

thumbnail Fig. 2

Spectrum of HD 50975 retrieved from IUE archive in the far ultraviolet region where radiation of the hot secondary star dominates. Four standard stars of spectral types between B0.5 and B3 are presented for comparison purposes. The profiles of Lyman α line at 1215 Å calculated using atmospheric models with Teff = 17 000 (solid line) and 25 000 K (dashed line) are shown.

3.3. The uvby photometry

A set of uvby colors for the primary component HD 50975A was calculated using unreddened colors of HD 50975 (Hauck & Mermilliod 1998) for a range of color excesses E(BV) from 0.05 to 0.25 mag. For (by), m1, and c1 of HD 50975A, the folowing limits are accepted:

  • 1)

    The range of (by) values from 0.34–0.44 mag according mean colors for F7-F9 supergiants (Gray et al. 2001), which correspond the uncertainty in spectral subtype ±1. This range was increased by 0.02 mag because of possible errors in observations and reduction.

  • 2)

    The values of c1 and m1 are within ±0.1 mag for the standard line of supergiants (Gray 1991). The appropriate contribution of the secondary star was subtracted, and the difference in absolute magnitudes between primary and secondary components was varied in steps of 0.01 mag. Normal ubvy colors for the secondary star were adopted from Crawford (1978). In addition, the relation (by)0 = 0.0867(c1 + m1) + 0.125 given by Shulov (1986) was used to calculate the normal (by) index for B star because Crawford (1978) provided only m1 and c1 indices. The results of calculations are plotted in Fig. 3 with the sequence of yellow supergiants given by Gray (1991) and the main sequence F-type stars given by Crawford (1975). In addition, a sample of supergiants of spectral types from F6 to G0 classified by Gray et al. (2001) are plotted for comparison purposes. Their colors are adopted from the photometric catalog compiled by Hauck & Mermilliod (1998) and are unreddened using E(by) values given by Gray et al. (2001) and reddening ratios given by Crawford & Mandwewala (1976). The range of parameters for the secondary star was found to be as follows: spectral type B0.5 V–B1.5 V, difference in magnitudes between primary and secondary star, ΔV = 1.8 − 2.4 mag, and extinction, E(BV) = 0.1 − 0.18 mag.

thumbnail Fig. 3

Modeled position of HD 50975A (crosses) in the c1 − (by) diagram calculated for the set of parameters (spectral type, magnitude, extinction) of the secondary star. The solid line represents main sequence F stars; dashed line – sequence of yellow supergiants; filled circles – a sample of single supergiants of spectral types from F6 to G0.

3.4. Spectral energy distribution

A set of the composite SEDs of two stars was calculated in the spectral region from 1250 to 5450 Å using the SEDs of MK spectral types given in Sviderskiene (1988) and was compared with the observed SED of HD 50975, which is based on the IUE spectrum and uvby photometry. The IUE spectrum was averaged over 50 Å intervals centered on the same wavelength as for MK-type SED. The uvby magnitudes were converted into absolute fluxes using the Gray (1998) calibration. The SED of MK types expressed in terms of intensities per unit wavelength interval were transformed into absolute fluxes by equating intensities to fluxes in the y passband (λ0 = 5456 Å) of the uvby system. Assuming a spectral type of F8  I for the primary star HD 50975A, the best fit between observed and calculated SED was achieved for the secondary of spectral type B1 V, a ratio of fluxes of 0.12 at 5500 Å for the companion stars, and E(BV) = 0.12 mag (see Fig. 4). We also estimated extinction using IUE data and the relationship of Vieira & Ponz (1995) applicable for O3–G5 single stars:

where F1,F2,Fc are fluxes estimated at λ1 = 1600 Å,λ2 = 2400 Å,λc = 2175 Å, respectively. The contribution of F8 Ib star to the flux in the ultraviolet (up to 1% at 1600 Å, 3–5% at 2175 Å, and 5–8% at 2400 Å) was extracted. The extinction was found to be E(BV) = 0.14 ± 0.04 using this method for both small and large apertures and for the averaging windows of 50 Å.

thumbnail Fig. 4

Comparison of the observed (IUE observations: filled circles, uvby photometry unreddened by E(bv) = 0.12: open circles), and modeled SED for HD 50975. The solid line represents a composite SED of two stars of spectral type F8 Iab and B1 V with a ratio of fluxes of 0.12 at 5500 Å.

3.5. Coravel cross-correlation function

The shape of CCF measured with the CORAVEL spectrometer is a convolution of averaged profiles of absorption lines and slots of mask on which the spectrum of star was focused. Therefore, the area defined by CCF is an analogue for an equivalent width (EW) of spectral line. The continuum contribution from the hot secondary partially fill the photospheric lines of the primary, making them less deep than they should be. This extra contribution (so-called veiling) can be measured by adding an extra continuum to a template star in order to fit the observed lines. Veiling of the spectrum of primary star reduces the equivalent widths. The decrease in EW can be used to estimate a difference in magnitudes between primary and unseen secondary components. Because of the different SEDs of the companion stars, the amount of veiling of the primary star increases toward short wavelengths.

The CCF of HD 50975 are plotted in Fig. 5, along with those for one main sequence star of spectral type G0V and single supergiant of spectral type F8 Ib (HD 194093), i.e. the same spectral type as for HD 50975. The spectral type of a supergiant was estimated by Gray et al. (2001). The EW of CCF for HD 50975 is about 15% lower than measured for the F8 Ib comparison star. Numerical simulations of the cross-corelation between the spectrum of F8 Ib, the B2 star, and the spectral transmission function of the mask was carried out. The decrease in EW by 15% corresponds to the difference of about 2.5 mag in magnitude between the primary and secondary stars in the case of similar metallicity of HD 50975 and the comparison star.

thumbnail Fig. 5

Cross-corelation function measured using CORAVEL spectrometer for HD 50975 (filled circles), along with those for two comparison stars: G0 V main sequence star (solid line) and single F8 I supergiant (dashed line). The dotted line represents CCF of the single F8 I star decreased in the equivalent width by 15%.

3.6. High-resolution spectroscopy

Absorption lines of iron are selected and measured in the observed high-resolution spectrum of HD 50975 with the goal of estimating the atmospheric parameters and metallicity of the primary star HD 50975A. The red and near infrared spectral region from about 6100 to 9600 Å was exploited to avoid a significant contribution of the secondary star. The modeled and scaled flux, Fν(B2) /Fν(F8) = 0.1 at 5500 Å, of the companion stars limits the contribution of the flux of B2 star below 8% in the red and near-infrared wavelength region (6.5 ± 1.5% from 6000 to 8000 Å). Such veiling should cause a minor decrease in EWs and a calculated metallicity (<0.1 dex). In addition, unblended iron lines are measured in the visual and blue spectral regions with the goal of quantiting the amount of veiling in different spectral regions.

A search for absorption lines for the secondary star was performed in the blue wavelength region shortward of 4500 Å. Unfortunately, in the blue region the spectrum of primary star is very crowded by absorption lines and detection of weak lines of the secondary star was plagued. A high-resolution spectrum shortward of 3600 Å is needed to detect photospheric lines of B2 star and to calculate the mass ratio of components.

3.6.1. Effective temperature of primary star

The excitation balance of iron was analyzed to obtain a color-independent estimate of the temperature, and Teff was chosen so that the derived iron abundance vs. excitation potential has a zero slope (Fig. 6). This excitation analysis confirms a temperature Teff = 5900 ± 150 K for HD 50975A, which is close to the mean effective temperature of F8 I supergiants, Teff = 5800 K (see Lyubimkov et al. 2010).

thumbnail Fig. 6

Abundance of iron derived from the Fe i lines in red wavelength region vs. the excitation potential of the lower energy level of the line, illustrating the estimation of the excitation temperature for HD 50975A.

The intensity of the sodium D12 doublet is a good indicator of stellar effective temperature; however, the sodium photospheric lines are asymmetric probably because of unresolved circumstellar absorption, hence plaguing examination of temperature using this criterion. In addition, the interstellar component is visible on the red wings of D1 and D2. The results of spectral synthesis for two effective temperatures and for the final calculated gravity, microturbulence, and abundances are presented in Fig. 7. As can be seen, the synthesized profiles of photospheric lines confirm an effective temperature of HD 50975A higher than 5500 K. The EW was measured to be 250 mÅ for interstellar D1 component, and an independent estimate of reddening was made for HD 50975 using relation between equivalent width and reddening (see Munari & Zwitter 1997), E(BV) ≃ 0.09 ± 0.01 mag.

thumbnail Fig. 7

Observed spectrum of HD 50975 (thick line) around the Na i D12 doublet, along with the calculated photospheric spectra for two different effective temperatures, Teff = 6000 K (thin line) and 5500 K (dashed line). The wavelengths are corrected for stellar velocity. Sharp interstellar components are visible on the red wings of the photospheric lines.

3.6.2. Gravity and microturbulence

The available spectroscopic methods for estimating the surface gravity are based on the study of the ionization equilibrium and damping wings of strong lines. Gravity for the primary star HD 50975A was determined by requiring that Fe i and Fe ii lines provide the same iron abundance (see Fig. 8). From the ionization balance, log g = 1.4 ± 0.3 (cgs) was found. The wings of the Hα line synthesized using Kurucz (1993) model in general confirm the gravity calculated using ionization balance (see Fig. 9).

thumbnail Fig. 8

Abundances derived from the Fe i and Fe ii (neutral lines – open circles; ionized lines – filled circles) lines in red wavelength region plotted vs. the EW of the line, illustrating the estimation of the microturbulent velocity for HD 50975A.

Table 3

Final atmospheric parameters and relative (to the Sun) abundances for primary star HD 50975A calculated using lines in the red and near-infrared wavelength region.

thumbnail Fig. 9

Observed spectrum of HD 50975 (dots) around Hα line, along with the one calculated using atmospheric model with Teff = 6000 K, log (g) = 1.5 (cgs), and [Fe/H] = −0.3.

The procedure for estimating the microturbulent velocity is illustrated in Fig. 8 for the primary star of HD 50975. As seen, a good agreement between the iron abundances derived from individual lines was found using ξt = 4.1 ± 0.3 km s-1. The final atmospheric parameters for HD 50975A are given in Table 3.

3.6.3. Abundances

The abundances were computed using the standard LTE line analysis program WIDTH6, written by Kurucz. This code employs an input model atmosphere to compute the strength of a given atomic line formed in such an atmosphere. The model atmospheres were taken from the grid of Kurucz (1993). The main goal was to estimate the iron abundance in the atmosphere of HD 50975. In addition, we calculated abundances for chemical elements that are represented by a set of unblended lines redward of 6100 Å. Spectral lines stronger than 300 mÅ were not used for abundance calculations. Oscillator strengths for the lines are compiled from a variety of sources by Luck (priv. comm.). A majority of the values are from high-precision laboratory measurements. The VALD database (Piskunov et al. 1995) was used to examine some problematic lines. In Table 4 the EWs are given for measured lines, along with the oscillator strengths and excitation potentials. The resulting mean relative abundances normalized with the solar-system abundances of Grevesse & Sauval (1998) are listed in Table 3. The abundance [X/H] is determined as the difference between the X abundance of the star and the Sun, i.e. [X/H] = log ϵ(X) − log ϵ(X), where the abundance log ϵ corresponds to the standard logarithmic scale with the hydrogen abundance log ϵ(H) = 12.00. Some wavelength regions around important atomic lines were calculated using the spectral synthesis code STARSP (Tsymbal 1996) based on the Kurucz’s code SYNTHE. The code for the spectrum synthesis adopts atomic and molecular partition functions and the continuous opacity package from ATLAS9 (Kurucz 1992). The systematic errors in abundances produced by the uncertainties in Teff (±150 K), log g (±0.3 dex), and microturbulent velocity (±0.3 km s-1) would lead to errors, less than 0.15 dex for the calculated elements.

The amount of veiling of EWs of the primary star because the continuum of the secondary star was examined on the basis of iron abundance calculated using Fe i lines over a broad wavelength region from 4000 to 9600 Å. A significant trend in the iron abundance from the wavelength of the line used was not detected in the region from 6100 to 9600 Å, [Fe/H] = −0.26 ± 0.06 dex, which agrees with the modeled fluxes of the companion stars. Indeed, it is difficult to detect minor ~3% variations in the flux in the observed spectrum because of uncertainties in atomic data, setting of continuum, and measurements of EWs. However, the mean iron abundance calculated using lines in the wavelength region from 5400 to 5600 Å (around λeff of V band) was found to be lower, [Fe/H] = −0.36 ± 0.04 dex (9 lines), which confirms the increasing flux from the secondary star and reduction of EWs of HD 50975. The continuum definition is perfect in the visual spectral region, and the difference of 0.1 dex in the mean abundance is because of lower EWs by about 5% on average. Two other samples of iron lines in the wavelength regions from 5030 to 5230 Å and from 4150 to 4260 Å give the mean iron abundance of [Fe/H] = −0.48 ± 0.04 dex (12 lines) and [Fe/H] ≃ −1.0 ± 0.14 dex (3 lines), respectively, in support of increasing veiling of EWs for HD 50975A in the blue spectral region. Unfortunately, in the blue region shortward of 4300 Å, the continuum definition is uncertain and a number of unblended lines are low, so the calculated mean iron abundance is uncertain.

Table 4

Wavelengths, excitation potentials, gf-values, and equivalent widths of lines measured in the spectrum of HD 50975A.

4. Mass and luminosity

The mass and bolometric magnitude of the secondary star HD 50975B was accepted to be a mean value for B2 main-sequence stars according Harmanec (1988), MB = 8.6 M and Mbol = −4.5 mag. The adopted mass for the secondary is close to what was derived from calibration given by Schmidt-Kaller (1982). The maximum (upper) mass of the primary star HD 50975A was estimated using the mass of secondary star, the mass function, and orbital inclination i ≤ 80 deg, MA ≤ 12.7 M. The inclination angle of 80 deg was adopted because the light curve does not show any indication of eclipses, and it takes the theoretical dimensions of both stars and the orbit into account. On the other hand, the minimal mass of the primary HD 50975A should be higher than the mass of the secondary star because of its more advanced stage of evolution.

An independent determination of the bolometric magnitude for HD 50975A can be made with the well known relation log (M/M) = log (g) − 4log (Teff) − 0.4 Mbol + 12.51, if the gravity and temperature are known from spectroscopy. When using the effective temperature Teff = 5900 K and adopting a spectroscopic gravity of log g = 1.4, the relation leads to the bolometric magnitude, Mbol = −5.5 ± 0.3 mag, for the primary star in the case of mass 8.6 MMA ≤ 12.7 M.

Extremely strong oxygen lines of an infrared triplet at 7771.94, 7774.17, and 7775.39 Å in the spectrum of HD 50975A confirm high luminosity for the primary star. The EW of the line at 7771.94 Å , which is the least blended member of the triplet, was measured to be about 600 mÅ, and the combined EW of all three components EW = 1.6 Å agrees with luminosity around log L/L ≃ 4 (Kovtyukh et al. 2008).

The final parameters for the both components of HD 50975 are provided in Table 5. The mean efective temperature of B2 main sequence stars according to De Jager & Nieuwenhuijzen (1987) was accepted for HD 50975B. Uncertainties are evaluated in accordance with the uncertainty of ±0.5 in the spectral subtype for secondary component, ±0.5 mag in the difference between magnitudes ΔV, and ±0.05 mag in E(BV). In Fig. 10 the positions of both binary components in the Hertzsprung-Russell (HR) diagram are compared to the evolutionary tracks calculated by Padova group (see Bertelli et al. 2009). Two tracks are given for 8 M and 10 M with a solar and a depleted metallicity of Z = 0.008 ([Fe/H] = −0.32 dex) to illustrate the uncertainty. In addition, a sample of F-G type stars adopted from Kovtyukh et al. (2010) is shown to display the area occupied by single F supergiants of different masses in the HR diagram: open circles are their original MV, crosses – MV derived from trigonometric parallaxes. As can be seen, the initial mass of the primary star HD 50975A is about 10 M according the determined stellar parameters and evolutionary tracks. The isochrone of 25 Myr fits both components of HD 50975 well in the log (Teff) vs. log (L/L) plane (see Fig. 10).

Table 5

Fundamental parameters for the primary and secondary components of HD 50975A.

thumbnail Fig. 10

Position of the primary and secondary components of HD 50975 (filled circles) in the log (Teff) vs. log (L/L) plane. Evolutionary tracks retrieved from Bertelli et al. (2009) labeled by their initial mass in solar units are given for solar composition (Z = 0.019, thin solid line) and depleted metallicity (Z = 0.008, dotted line). The thick gray line presents a segment of the 25 Myr isochrone. A sample of single F-G type stars are plotted for comparison.

5. Discussion and conclusions

A comprehensive analysis of the long-term radial velocity measurements collected for the primary star of a spectroscopic binary HD 50975A was performed with the goal of calculating orbital elements. The orbital period was found to be 190.22 ± 0.01 days, and the mass function was calculated to be high, f(m) = 1.402 ± 0.005 M. The absence of eclipses in the light curve gives evidence about an orbital inclination of i ≤ 80 deg. Unfortunately, high-resolution spectroscopy shortward of 3600 Å and a search for absorption lines of the secondary star were unsuccessful. Therefore, calculating the mass ratio was postponed to the future.

Analysis and modeling of the low-resolution IUE spectrum of HD 50975 in the far ultraviolet region, along with those for early type standard stars, revealed the nature of the hot secondary. HD 50975B is a star of spectral type B2 ±0.5 near ZAMS with an effective temperature Teff ≃ 21 000 K. Comparison of the SED for HD 50975 with those for F8 supergiant and hot secondary confirm the spectral type about B1 for the secondary star HD 50975B. The mass and bolometric magnitude was adopted to be M ≃ 8.6 M and Mbol ≃ − 4.5 mag, according to the mean parameters for B2 main sequence stars (Harmanec 1988). Using Mbol = −4.5 mag and BC = 1.5 mag (see Flower 1996), MV ≃ − 3.0 mag for HD 50975B.

The high-resolution spectrum of HD 50975 in the large spectral region was studied in detail with the goal of estimating the atmospheric and fundamental parameters of both stellar components. The red and near infrared spectral region from 6100 to 9600 Å was analyzed and modeled using the method of atmospheric models with the goal of estimating atmospheric parameters, abundances, luminosity, and distance for the primary star HD 50975A. The effective temperature and surface gravity were found to be typical of F8 supergiants: Teff = 5900 ± 150 K, and log (g) = 1.4 ± 0.3 (cgs). The spectroscopic temperature and gravity lead to a bolometric magnitude for the primary star, Mbol = −5.5 ± 0.3 mag for mass 8.6 MMA ≤ 12.7 M. Strong lines of the oxygen infrared triplet confirm high luminosity. The spectroscopic parameters compared with stellar evolutionary tracks comfirm the initial mass of about 10 M for the primary. The distance modulus, (VMV) = 13.0 mag, and total absorption, AV = 0.3 mag, places HD 50975 at the distance of 3.5 kpc. Using the estimated distance and galactic longitude, b = + 5 deg, the distance was estimated to be about 300 pc above the Galactic plane. The iron abundance calculated using equivalent widths and atmospheric models confirms slight metal deficiency in the atmosphere of primary, [Fe/H] = −0.26 dex. Modeling of the flux of the both companion stars and of the measurements of veiled EWs of iron lines give evidence for a minor contribution of secondary HD 50975B in the red and near infrared spectral region. The calculated abundances for Si and Ca agree with the mild iron deficiency on average. The iron abundances for single A, F, and G-supergiants with the masses between 4 and 15 M are between log ϵ (Fe) = 7.28 and 7.67 in comparison with 7.45 for the Sun (Lyubimkov et al. 2010). The abundance of europium was found to be enhanced in the atmosphere of HD 50975A, [Eu/H] = + 0.61 ± 0.07 dex. Thus, the site where HD 50975 formed seems to be enriched by r-process elements in the past.

The absolute magnitude MV = −3.0 mag for HD 50975B and the distance modulus estimated for HD 50975A lead to the apparent magnitude V ≃ 10.0 mag for the secondary star HD 50975B in the V band, therefore the difference in magnitudes between the secondary and primary stars could be about 2.5 mag. Detected veiling of EW by about 5% around λeff = 5450 Å of V band confirms a bit larger difference in the apparent magnitudes. Some discrepancy between the estimations on the basis of CCF and photometry seems to come from the different wavelength range for the CORAVEL mask and from the neglected slight metal deficiency for HD 50975A.

High luminosity, log (L/L) = 4.1 ± 0.1, places HD 50975A among the brightest yellow supergiants. Unfortunately, the parallax of HD 50975, π = 2.68 ± 0.54 mas (van Leeuwen 2007), leads to the distance of kpc. At such a distance, log (L/L) = 2.5 for the primary star HD 50975A and 1.1 for the secondary star HD 50975B, which are discrepant with the results gathered from spectroscopy. In addition, the minimum mass of secondary, which matches the calculated value of f(m) = 1.40 M, is about 6 M in the case of MA>MB and i< 80 deg. On the other hand, the position of HD 50975 at the distance of 0.37 kpc should be variable because of the orbital motion and a change in the position during a year should exceed the parallax. Curiously, the old Hipparcos parallax, π = 1.38 ± 0.98 mas (Perryman & ESA 1997), does not contradict the high luminosity estimated from spectroscopy. We concluded that the parallax and proper motions are uncertain for HD 50975A.

The progenitor of the SN 2011dh was detected in archival pre-explosion Hubble Space Telescope images (see Maund et al. 2011). It is consistent with an F8 supergiant (log (L/L) = 4.92 ± 0.2, Teff = 6000 ± 280 K). Through comparison with stellar evolution tracks, this corresponds to a single star at the end of core C-burning with an initial mass of 13 ± 3 M. Early photometric and spectroscopic observations show fast evolution similar to the transitional type IIb SN 2008 ax and suggest that a large amount of the progenitor’s hydrogen envelope was removed before explosion. There are two main mechanisms proposed for stars to lose their outer layers: strong winds and mass transfer in interacting binary systems. Bersten et al. (2012) carried out stellar evolution calculations for a system of two massive stars in a close orbit, in which phases of mass transfer occur. By assuming a system of stars with initial masses of 16 and 10 M, and an initial period of 125 days, they were able to obtain a configuration for the mass-donor star that closely matches the observations in the pre-supernova images and has a core mass that is consistent with their hydrodynamic model. The parameters of the binary system HD 50975 are not very far from those simulated for the SN 2011dh progenitor. HD 50975 is near a semi-detached configuration with a radius of primary of RA = 107 R, estimated using the fundamental relation among R, L, and Teff, and a distance between both components of 370 R, calculated using the period and masses given in Tables 2 and 5. The radius of the Roche lobe for HD 50975A was calculated to be 120 R according to a formula by Eggleton (1983). Double-peaked light curve of HD 50975A, which is synhronized with the maximum and minimum of radial velocity variations, confirms a gravitational distortion of the primary star.


Acknowledgments

L.Z. acknowledges the support from the Research Council of the Lithuania under the grant No.VIZIT-5-TYR-052. We thank S. Velichko and A. Bondar for observations on the 2 m telescope at the observatory on the Terskol Peak. The referee, T. Zwitter, is thanked for his constructive comments on the paper. This research made use of the SIMBAD database, operated at the CDS, Strasbourg, France. Some of the data presented in this paper were obtained from the Multimission Archive at the Space Telescope Science Institute (MAST). STScI is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS5-26555. Support for MAST for non-HST data is provided by the NASA Office of Space Science via grant NAG5-7584 and by other grants and contracts.

References

  1. Barbier, R., Dossin, F., Jaschek, M., et al. 1978, A&A, 66, L9 [NASA ADS] [Google Scholar]
  2. Bartkevicius, A., & Sperauskas, J. 1994, Balt. Astron., 3, 49 [NASA ADS] [Google Scholar]
  3. Bersten, M. C., Benvenuto, O. G., Nomoto, K. E., et al. 2012, ApJ, 757, 31 [NASA ADS] [CrossRef] [Google Scholar]
  4. Benvenuto, O. G., Bersten, M. C., & Nomoto, K. 2013, ApJ, 762, 74 [NASA ADS] [CrossRef] [Google Scholar]
  5. Bertelli, G., Nasi, E., Girardi, L., & Marigo, P. 2009, A&A, 508, 355 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Crawford. D. L. 1975, AJ, 80, 955 [NASA ADS] [CrossRef] [Google Scholar]
  7. Crawford. D. L. 1978, AJ, 83, 48 [NASA ADS] [CrossRef] [Google Scholar]
  8. Crawford, D. L., & Mandwewala, N. 1976, PASP, 86, 91 [Google Scholar]
  9. De Jager, C., & Nieuwenhuijzen, H. 1987, A&A, 177, 217 [NASA ADS] [Google Scholar]
  10. Eggleton, P. P. 1983, ApJ, 268, 368 [NASA ADS] [CrossRef] [Google Scholar]
  11. Flower, J. P. 1996, ApJ, 469, 355 [NASA ADS] [CrossRef] [Google Scholar]
  12. Gray, R. O. 1991, A&A, 252, 237 [NASA ADS] [Google Scholar]
  13. Gray, R. O. 1998, AJ, 116, 482 [NASA ADS] [CrossRef] [Google Scholar]
  14. Gray, R. O., Napier, M. G., & Winkler, L. I. 2001, AJ, 121, 2148 [NASA ADS] [CrossRef] [Google Scholar]
  15. Grevesse, N., & Souval, A. J. 1998, Space Sci. Rev., 85, 161 [NASA ADS] [CrossRef] [Google Scholar]
  16. Griffin, R. 1967, ApJ, 148, 465 [NASA ADS] [CrossRef] [Google Scholar]
  17. Harmanec, P. 1988, Bull. Astron. Inst. Czechoslovakia, 39, 329 [Google Scholar]
  18. Hauck, B., & Mermilliod, M. 1998, A&AS, 129, 431 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Koen, C., & Eyer, L. 2002, MNRAS, 331, 45 [NASA ADS] [CrossRef] [Google Scholar]
  20. Kovtyukh, V. V., Chekhonadskikh, F. A., Luck, R. E., et al. 2010, MNRAS, 408, 1568 [Google Scholar]
  21. Kovtyukh, V. V., Gorlova, N. I., & Belik, S. I. 2012, MNRAS, 423, 3268 [NASA ADS] [CrossRef] [Google Scholar]
  22. Kurucz, R. L.1992, in Stellar Populations of Galaxies, eds. B. Barbuy, & A. Renzini (Dordrecht: Kluwer), IAU Symp., 149, 225 [Google Scholar]
  23. Kurucz, R. L. 1993, CD-ROM 13 & 18, Smithsonian Astrophys. Obs. [Google Scholar]
  24. Lyubimkov, L. S., Lambert, D. L., Rostopchin, S. I., Rachkovskaya, T. M., & Poklad, D. B. 2010, MNRAS, 402, 1369 [NASA ADS] [CrossRef] [Google Scholar]
  25. MacConnell, D. J., & Bidelman, W. P. 1976, AJ, 81, 225 [NASA ADS] [CrossRef] [Google Scholar]
  26. Massa, D., & Fitzpatrick, E. L. 2000, ApJS, 126, 517 [NASA ADS] [CrossRef] [Google Scholar]
  27. Maund, J. R., Fraser, M., Ergon, M., et al. 2011, ApJ, 739, L37 [NASA ADS] [CrossRef] [Google Scholar]
  28. Munari, U., & Zwitter, T. 1997, A&A, 318, 269 [NASA ADS] [Google Scholar]
  29. Nidever, D. L., Marcy, G. W., Butler, R. P., Fische, D. A., & Vogt, S. S. 2002, ApJS, 141, 503 [NASA ADS] [CrossRef] [Google Scholar]
  30. Olsen, E. H. 1979, A&AS, 37, 367 [NASA ADS] [Google Scholar]
  31. Olsen, E. H. 1983, A&AS, 54, 55 [NASA ADS] [Google Scholar]
  32. Perryman, M. A. C., & ESA 1997, The Hipparcos and Tycho catalogues (Noordwijk, Netherlands: ESA Publications Division), ESA SP, 1200 [Google Scholar]
  33. Piskunov, N. E., Kupka, F., Ryabchikova, T. A., Weiss, W. W., & Jeffery, C. S. 1995, A&AS, 112, 525 [NASA ADS] [Google Scholar]
  34. Pojmanski, G. 2002, Acta Astron., 52, 397 [NASA ADS] [Google Scholar]
  35. Rountree, J., & Sonneborn, G. 1991, ApJ, 369, 515 [NASA ADS] [CrossRef] [Google Scholar]
  36. Schmidt-Kaler, T. 1982, in Landolt-Boernstein, Vol. VI/2b, eds. K. Schaifers, & H. H. Voigt (Springer-Verlag), 1 [Google Scholar]
  37. Shulov, O. S. 1986, Sov. Astron., 30, 435 [NASA ADS] [Google Scholar]
  38. Sviderskiene, Z. 1988, Vilnius Obs. Bull., 80, 7 [NASA ADS] [Google Scholar]
  39. Tsymbal, V. V. 1996, ASP Conf. Ser., 108, 198 [Google Scholar]
  40. Udry, S., Mayor, M., Maurice, E., et al. 1999, in Precise Stellar Radial Velocities, eds. J. B. Hearnshaw, & C. D. Scarfe, ASP Conf., 185, 383 [Google Scholar]
  41. Upgren, A., Sperauskas, J., & Boyle, R. P. 2002, Balt. Astron., 11, 91 [Google Scholar]
  42. van Leeuwen, F. 2007, A&A, 474, 653 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  43. Vieira, E. F., & Ponz, J. D. 1995, A&AS, 111, 393 [NASA ADS] [Google Scholar]

All Tables

Table 1

Radial velocity measurements for HD 50975A, along with the uncertainties, phases, and residuals for the orbital curve.

Table 2

The orbital elements along with standard deviations for HD 50975A.

Table 3

Final atmospheric parameters and relative (to the Sun) abundances for primary star HD 50975A calculated using lines in the red and near-infrared wavelength region.

Table 4

Wavelengths, excitation potentials, gf-values, and equivalent widths of lines measured in the spectrum of HD 50975A.

Table 5

Fundamental parameters for the primary and secondary components of HD 50975A.

All Figures

thumbnail Fig. 1

Observed heliocentric radial velocities for HD 50975A (filled circles) as a function of phase, with the velocity curve corresponding to the calculated orbital elements (upper panel). ASAS (middle panel) and Hipparcos (lower panel) photometry phased using the same orbital period gives evidence of the ellipsoidal shape of the primary star.

In the text
thumbnail Fig. 2

Spectrum of HD 50975 retrieved from IUE archive in the far ultraviolet region where radiation of the hot secondary star dominates. Four standard stars of spectral types between B0.5 and B3 are presented for comparison purposes. The profiles of Lyman α line at 1215 Å calculated using atmospheric models with Teff = 17 000 (solid line) and 25 000 K (dashed line) are shown.

In the text
thumbnail Fig. 3

Modeled position of HD 50975A (crosses) in the c1 − (by) diagram calculated for the set of parameters (spectral type, magnitude, extinction) of the secondary star. The solid line represents main sequence F stars; dashed line – sequence of yellow supergiants; filled circles – a sample of single supergiants of spectral types from F6 to G0.

In the text
thumbnail Fig. 4

Comparison of the observed (IUE observations: filled circles, uvby photometry unreddened by E(bv) = 0.12: open circles), and modeled SED for HD 50975. The solid line represents a composite SED of two stars of spectral type F8 Iab and B1 V with a ratio of fluxes of 0.12 at 5500 Å.

In the text
thumbnail Fig. 5

Cross-corelation function measured using CORAVEL spectrometer for HD 50975 (filled circles), along with those for two comparison stars: G0 V main sequence star (solid line) and single F8 I supergiant (dashed line). The dotted line represents CCF of the single F8 I star decreased in the equivalent width by 15%.

In the text
thumbnail Fig. 6

Abundance of iron derived from the Fe i lines in red wavelength region vs. the excitation potential of the lower energy level of the line, illustrating the estimation of the excitation temperature for HD 50975A.

In the text
thumbnail Fig. 7

Observed spectrum of HD 50975 (thick line) around the Na i D12 doublet, along with the calculated photospheric spectra for two different effective temperatures, Teff = 6000 K (thin line) and 5500 K (dashed line). The wavelengths are corrected for stellar velocity. Sharp interstellar components are visible on the red wings of the photospheric lines.

In the text
thumbnail Fig. 8

Abundances derived from the Fe i and Fe ii (neutral lines – open circles; ionized lines – filled circles) lines in red wavelength region plotted vs. the EW of the line, illustrating the estimation of the microturbulent velocity for HD 50975A.

In the text
thumbnail Fig. 9

Observed spectrum of HD 50975 (dots) around Hα line, along with the one calculated using atmospheric model with Teff = 6000 K, log (g) = 1.5 (cgs), and [Fe/H] = −0.3.

In the text
thumbnail Fig. 10

Position of the primary and secondary components of HD 50975 (filled circles) in the log (Teff) vs. log (L/L) plane. Evolutionary tracks retrieved from Bertelli et al. (2009) labeled by their initial mass in solar units are given for solar composition (Z = 0.019, thin solid line) and depleted metallicity (Z = 0.008, dotted line). The thick gray line presents a segment of the 25 Myr isochrone. A sample of single F-G type stars are plotted for comparison.

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.