Absolute Dimensions of the Unevolved F-Type Eclipsing Binary BT Vulpeculae
Total Page:16
File Type:pdf, Size:1020Kb
Accepted for publication in The Astrophysical Journal A Preprint typeset using LTEX style emulateapj v. 12/16/11 ABSOLUTE DIMENSIONS OF THE UNEVOLVED F-TYPE ECLIPSING BINARY BT VULPECULAE Guillermo Torres1, Claud H. Sandberg Lacy2, Francis C. Fekel3, and Matthew W. Muterspaugh4,5 Accepted for publication in The Astrophysical Journal ABSTRACT We report extensive differential V -band photometry and high-resolution spectroscopy for the 1.14 day, detached, double-lined eclipsing binary BT Vul (F0+F7). Our radial-velocity monitoring and light N curve analysis lead to absolute masses and radii of M1 = 1.5439 ± 0.0098 M⊙ and R1 = 1.536 ± N N N 0.018 R⊙ for the primary, and M2 = 1.2196 ± 0.0080 M⊙ and R2 = 1.151 ± 0.029 R⊙ for the secondary. The effective temperatures are 7270±150 K and 6260±180 K, respectively. Both stars are rapid rotators, and the orbit is circular. A comparison with stellar evolution models from the MIST series shows excellent agreement with these determinations, for a composition of [Fe/H] = +0.08 and an age of 350 Myr. The two components of BT Vul are very near the zero-age main sequence. 1. INTRODUCTION Table 1 The discovery of the photometric variability of Times of Minimum Light for BT Vul BT Vul (HD 340072, TYC 2164-161-1, Gaia DR2 1836032243422750464; V = 11.49, SpT F0+F7) was HJD σ O − C made by S. Beljawsky (see Guthnick & Schneller 1939), (2,400,000+) Year Epoch (days) (days) Ecl Type although the original period given (3.988 days) was in- 34740.28 1953.9912 −16421.0 0.0063 +0.004269 1 PG correct. Other than measurements of the times of eclipse 35402.18 1955.8034 −15841.0 0.0063 +0.007879 1 PG since its discovery, no detailed studies of this short-period 43351.786 1977.5682 −8875.0 0.0100 +0.009984 1 V 43779.735 1978.7399 −8500.0 0.0100 +0.008731 1 V system (P =1.14 days) have appeared in the literature. 44022.812 1979.4054 −8287.0 0.0100 +0.009988 1 V In this paper we report the first systematic photometric and spectroscopic monitoring of the binary, leading to a Note. — ‘Ecl’ is 1 for the primary and 2 for the secondary minimum; ‘Type’ is PG, V, and PE for the photographic, visual, full determination of its physical properties. and photoelectric/CCD techniques. The ‘Epoch’ is counted from Section 2 presents our observations, beginning in the reference time of primary minimum in this section. For Section 2.1 with a determination of a highly precise measurements with published uncertainties the errors reported here have been multiplied by scale factors of 3.35 and 2.75 for ephemeris based on more than six decades of timing ob- the photoelectric timings of the primary and secondary minima servations. Section 2.2 then reports our spectroscopic (see text), and scale factors of 4.35 and 1.28 for the visual observations of BT Vul, with the derivation of the radial observations. Measurements with no published errors have been velocities and spectroscopic orbital elements. Our exten- assigned values of 0.0063 days, 0.010 days, and 0.0058 days for the photographic, visual, and photoelectric techniques, respec- sive photometric V -band measurements are described in tively. Sources for the individual measurements are listed at Section 2.3. The analysis of the lightcurves is found in http://var2.astro.cz/EN/brno/eclipsing_binaries.php and Section 3, and is followed in Section 4 by a determina- https://www.bav-astro.eu/index.php/veroeffentlichungen/ tion of the absolute dimensions of the system, and in service-for-scientists/lkdb-engl. (This table is available in Section 5 with a comparison of the physical properties its entirety in machine-readable form.) with current models of stellar evolution. Final remarks tial photometric observations we describe in this paper. appear in Section 6. A linear ephemeris was derived by weighted least squares, in which the published uncertainties, when avail- 2. OBSERVATIONS able, were scaled in order to achieve reduced χ2 values 2.1. Eclipse timings near unity for each type of eclipse and each type of obser- arXiv:2004.00022v1 [astro-ph.SR] 31 Mar 2020 vation. For measurements with no reported uncertainties Times of minimum light have been gathered for BT Vul since 1953 using photographic, visual, and photoelec- we adopted suitable values to the same end (see Table 1). tric/CCD techniques. The available measurements are The linear ephemeris obtained is: collected in Table 1, which contains 115 timings for the Min I(HJD) = 2,453,479.931997(92)+1.141200674(50)E, primary minimum and 40 for the secondary over more than 65 years. Many of them are based on the differen- with the uncertainties indicated in parentheses. Indepen- dent periods determined from the primary and secondary 1 Center for Astrophysics | Harvard & Smithsonian, 60 Garden minima do not differ significantly. A solution allowing for Street, Cambridge, MA 02138, USA; [email protected] separate reference times of primary and secondary min- 2 Physics Department, University of Arkansas, Fayetteville, ima resulted in a phase difference of 0.49987 ± 0.00014, AR 72701, USA 3 Center of Excellence in Information Systems, Tennessee suggesting negligible eccentricity. For the rest of this State University, Nashville, TN 37209, USA work we have therefore assumed the orbit to be circular. 4 Fairborn Observatory, 1327 Duquesne Rd, Patagonia, AZ 85624 2.2. Spectroscopy 5 Current address: Division of Science, Technology, and Math- ematics, Columbia State Community College, 1665 Hampshire Spectroscopic monitoring of BT Vul was carried out Pike, Columbia, TN 38401 at two facilities. Observations at the Center for As- 2 Torres et al. Table 2 Table 3 CfA Radial Velocities for BT Vul Fairborn Radial Velocities for BT Vul HJD RV1 RV2 σ1 σ2 Orbital HJD RV1 RV2 Orbital (2,400,000+) (km s−1) (km s−1) (km s−1) (km s−1) Phase (2,400,000+) (km s−1) (km s−1) Phase 55368.9509 −139.96 136.65 4.56 4.60 0.29074 55867.7441 −115.7 97.0 0.36829 55527.5607 −142.87 135.21 3.54 3.57 0.27576 56010.9947 53.1 −122.1 0.89451 55850.5866 −129.49 121.77 2.50 2.52 0.33369 56023.9539 −151.9 136.7 0.25027 55903.5554 103.22 −177.41 3.99 4.02 0.74866 56058.9383 53.5 −107.2 0.90605 56024.0077 −143.75 133.62 1.86 1.87 0.29741 56060.8952 67.2 −129.2 0.62082 Note. — Orbital phases are counted from the reference time of Note. — Typical uncertainties for the ve- primary eclipse. (This table is available in its entirety in machine- locities are 3.0 km s−1 for the primary and readable form.) 3.4 km s−1 for the secondary. Orbital phases are counted from the reference time of primary trophysics (CfA) began in June of 2010 and contin- eclipse. (This table is available in its entirety ued until October of 2015. They were made with in machine-readable form.) the Tillinghast Reflector Echelle Spectrograph (TRES; Szentgyorgyi & F˝ur´esz 2007; F˝ur´esz 2008), a fiber-fed, the temperatures, and 70/55 kms−1 for the rotational bench-mounted instrument attached to the 1.5m Tilling- broadening. hast reflector at the Fred L. Whipple Observatory on Previous experience at CfA has shown that the raw Mount Hopkins (Arizona, USA). A total of 39 spectra radial velocities can sometimes be affected by systematic were gathered at a resolving power of R ≈ 44,000, and errors that may arise due to lines shifting in and out cover the wavelength region 3800–9100 A˚ in 51 orders. of the spectral window as a function of orbital phase, For the order centered at ∼5187 A,˚ which contains the or due to the barycentric velocity of the Earth (see Mg I b triplet, the signal-to-noise ratios range from 24 to Latham et al. 1996; Torres et al. 1997). Following these 65 per resolution element of 6.8 km s−1. The wavelength authors we applied corrections for this effect based on reference was provided by exposures of a thorium-argon simulations, the details of which may be found in those papers. In the case of the primary the adjustments are lamp before and after each science frame, and the reduc- −1 −1 tions were carried out with a dedicated pipeline. very small (< 0.4 kms ), but they can reach 2 kms For the radial velocity measurements we used the for the secondary. The resulting radial velocities in the two-dimensional cross-correlation algorithm TODCOR heliocentric frame, including corrections, are listed in Ta- (Zucker & Mazeh 1994), with separate templates for the ble 2 along with their uncertainties. The flux ratio we ± primary and secondary selected from a large library of measure using TODCOR is ℓ2/ℓ1 = 0.28 0.02, at the pre-computed synthetic spectra based on model atmo- mean wavelength of our observations (5187 A).˚ spheres by R. L. Kurucz, and a line list tuned to bet- Further spectroscopic observations of BT Vul were ob- ter match the spectra of real stars (see Nordstr¨om et al. tained from 2011 November through 2018 September at 1994; Latham et al. 2002). The synthetic spectra cover Fairborn Observatory, which is situated in southeast Ari- a limited wavelength region of about 300 A˚ centered zona near Washington Camp. The observations were around 5187 A,˚ and the velocity measurements were acquired with the Tennessee State University 2m As- tronomical Spectroscopic Telescope (AST) and a fiber- made using the central 100 A,˚ which contains most of the fed echelle spectrograph (Eaton & Williamson 2007).