The Subdwarf B + White Dwarf Binary KPD 1930 + 2752, a Supernova Type Ia Progenitor Candidate
Total Page:16
File Type:pdf, Size:1020Kb
15th European Workshop on White Dwarfs ASP Conference Series, Vol. 372, 2007 R. Napiwotzki and M. R. Burleigh The Subdwarf B + White Dwarf Binary KPD 1930 + 2752, a Supernova Type Ia Progenitor Candidate S. Geier,1 S. Nesslinger,1 U. Heber,1 N. Przybilla,1 R. Napiwotzki,2 and R.-P. Kudritzki3 1Dr.–Remeis–Sternwarte, Institute for Astronomy, University Erlangen-Nuremberg, Sternwartstr. 7, 96049 Bamberg, Germany 2Centre of Astrophysics Research, University of Hertfordshire, College Lane, Hatfield AL10 9AB, UK 3Institute for Astronomy, University of Hawaii, 2680 Woodlawn Drive, Honolulu, HI 96822, USA Abstract. The nature of the progenitors of type Ia supernovae is still under controversial debate. KPD 1930+2752 is one of the best SN Ia progenitor can- didates known today. The object is a double degenerate system consisting of a subluminous B (sdB) star and a massive white dwarf. Maxted, Marsh, & North (2000) conclude that the system mass exceeds the Chandrasekhar mass. This conclusion rests on the assumption that the sdB mass is 0.5 M⊙. However, recent binary population synthesis calculations suggest that the mass of an sdB star may range from 0.3 M⊙ to more than 0.7 M⊙. It is therefore important to mea- sure the mass of the sdB star simultaneously with that of the white dwarf. Since the rotation of the sdB star is tidally locked to the orbit the inclination of the system can be constrained. An analysis of the ellipsoidal variations in the light curve allows to tighten the constraints derived from spectroscopy. We derive the mass–radius relation for the sdB star from a quantitative spectral analysis. The projected rotational velocity is determined for the first time from high-resolution spectra. In addition a reanalysis of the published light curve is performed. The atmospheric and orbital parameters are measured with unprecedented accuracy. −1 In particular the projected rotational velocity vrot sin i = 92.3 ± 1.5 kms is determined. The mass of the sdB is limited between 0.45 M⊙ and 0.52 M⊙. The total mass of the system ranges from 1.36 M⊙ to 1.48 M⊙ and hence is likely to exceed the Chandrasekhar mass. So KPD 1930 + 2752 qualifies as an excellent double degenerate supernova Ia progenitor candidate. 1. Introduction Supernovae of type Ia (SN Ia) play an important role in the study of cosmic evolution. They are regarded as the best standard candles for the determina- tion of the cosmological parameters. The nature of their progenitors is still under debate. There is general consensus that only the thermonuclear explosion of a white dwarf (WD) is compatible with the observed features of SN Ia. A white dwarf has to accrete mass from a close companion to reach the Chan- drasekhar limit of 1.4 M⊙. Two main scenarios of mass transfer are currently under discussion. In the so-called single degenerate (SD) scenario (Whelan & Iben 1973), the mass donating component is a red giant/subgiant, which fills 393 394 Geier et al. its Roche lobe and is continually transferring mass to the white dwarf. Accord- ing to the so-called double degenerate (DD) scenario (Iben & Tutukov 1984) the mass donating companion is a white dwarf, which eventually merges with the primary due to orbital shrinkage caused by gravitational wave radiation. KPD 1930 + 2752 was identified as a pulsating subdwarf B star (sdBV). Multi- periodic variations with short periods and low amplitudes were detected. In addition to these pulsations a strong variation at a much longer period of about 4100 s was found. This variation could be identified as an ellipsoidal deformation of the sdB star most likely caused by a massive companion. Bill`eres et al. (2000) predicted the period of the binary to be two times the period of this brightness variation. This was proven by Maxted et al. (2000), who measured a radial ve- locity curve of KPD 1930 + 2752 which matched the proper period. The radial velocity amplitude combined with an assumption of the canonical mass for sdB stars MsdB = 0.5 M⊙ led to a lower limit for the mass of the system derived from the mass function. This lower limit M ≥ 1.47 M⊙ exceeded the Chandrasekhar mass of 1.4 M⊙. The system is considered to become double degenerate when the subdwarf eventually evolves to a white dwarf. Orbital shrinkage caused by gravitational wave radiation will lead to a merger of the binary. Maxted et al. concluded that KPD 1930 + 2752 could be a good candidate for the progenitor of a type Ia supernova. This conclusion rests on the assumption that the mass of the primary is 0.5 M⊙. SdB stars have been identified as core helium burning stars at the blue end of the horizontal branch. According to canonical stellar evolution the mass is fixed by the onset of the core helium flash to about half a solar mass (Heber 1986). However, binary population synthesis calculations (Han et al. 2003) suggest, that the mass range for sdB stars is much wider, ranging from 0.3 M⊙ to more than 0.7 M⊙. Therefore we aim at measuring the mass of the sdB star simultaneously with that of the white dwarf. Since the rotation of the sdB star is tidally locked to the orbit the inclination of the system can be constrained, if the sdB radius and the projected rotational velocity can be measured with high precision. An analysis of the ellipsoidal variations in the light curve allows to tighten the constraints derived from spectroscopy. We derive the mass–radius relation for the sdB star from a quantitative spectral analysis. The projected rotational velocity is determined for the first time from high-resolution spectra. In addition a reanalysis of the published light curve is performed. 2. Observations and Data Analysis With the 10 m Keck I Telescope at the Mauna Kea Observatory two hundred high-resolution spectra were obtained using the High Resolution Echelle Spec- trometer (HIRES). Two spectra were taken with the ESO Very Large Telescope UT2 (Kueyen) equipped with the UV-Visual Echelle Spectrograph (UVES). Ad- ditional 150 low-resolution spectra were obtained with the 2.2 m Telescope at the Calar Alto Observatory using the Calar Alto Faint Object Spectrograph (CAFOS). The data was reduced using the ESO–MIDAS package. The radial velocities were obtained by cross correlation with a model spectrum at rest wave- length. Combining with data from literature we derive: K = 341 ± 1kms−1, P = 0.0950933 ± 0.0000015 d. The sdB+WD Binary KPD 1930+2752 395 For the analysis of the spectra of KPD 1930+2752 we used LTE model grids for 10 times solar metallicity. As the surface gravity is of utmost importance for our analysis we also calculated new grids of models and synthetic spectra to account for NLTE effects and metal line blanketing simultaneously. The mean parameter values from the model atmosphere fits are: Teff = 35 200 ± 500 K and log g = 5.61 ± 0.06 dex. The primary aim of the high-resolution time-series spectroscopy was to measure the projected rotational velocity of KPD 1930 + 2752 as accurately as possible. For this purpose the spectra obtained with HIRES and UVES were used. The projected rotational velocity was measured by convolving a synthetic spectrum calculated from the best fit model atmosphere with a rotational broad- ening ellipse for appropriate vrot sin i (see Figure 1). The result is very accurate: −1 vrot sin i = 92.3 ± 1.5kms . 3. Mass and Inclination KPD 1930 + 2752 is obviously affected by the gravitational forces of the com- panion, demonstrated by its ellipsoidal deformation. Since the period of the photometric variations is exactly half the period of the radial velocity variations the rotation of the sdB star is tidally locked to the orbit. Having determined the gravity and projected rotational velocity, we have three equations at hand that constrain the system, with the sdB mass MsdB being the only free parameter. Besides the mass function 3 3 2 Mcomp(sin i) KsdBP f(MsdB,Mcomp)= 2 = (1) (MsdB + Mcomp) 2πG these are v P sin i = rotsini (2) 2πR M G R = sdB (3) s g With log g obtained from the model atmosphere analysis, the radius of the star R was calculated using the standard mass–radius relation (Eq. 3). Together with the orbital period of the system P and the projected rotational velocity vrotsini = vrot sin i the inclination of the system sin i was derived for different values of the sdB mass. Because the rotation is tidally locked to the orbit, the rotational period of the sdB equals the orbital period of the system. Therefore the absolute value of the rotational velocity could be calculated. From the measured projected rotational velocity the inclination of the system follows (Eq. 2). With the sdB mass as a free parameter, the measured radial velocity semi-amplitude KsdB and orbital period P the mass function was solved numerically (Eq. 1) to derive the mass of the companion Mcomp and calculate the total mass of the binary. The fact that sin i cannot exceed unity gave a lower limit for the mass of the sdB MsdB ≥ 0.45 M⊙. The total mass of the system exceeds the Chandrasekhar limit for almost all assumptions of MsdB. The inclination angle of the system is close to 90◦ implying that KPD 1930 + 2752 may be an eclipsing binary. 396 Geier et al. 1.00 1.80 113 113 112 1.70 112 0.98 eclipse limit 111 eclipse limit 111 110 sun 110 109 109 108 ) / M 1.60 108 0.96 107 comp 107 sin(i) 106 106 1.50 105 + M 105 χ2 χ2 sdB 0.94 (M 1.40 0.92 1.30 0.48 0.52 0.56 0.60 0.48 0.52 0.56 0.60 MsdB / Msun MsdB / Msun Figure 1.