Draft version November 19, 2020 Typeset using LATEX twocolumn style in AASTeX63 DY Pegasi: An SX Phoenicis Star in a Binary System with an Evolved Companion Hui-Fang Xue1 and Jia-Shu Niu2, 3, 4 1Department of Physics, Taiyuan Normal University, Jinzhong, 030619, China 2Institute of Theoretical Physics, Shanxi University, Taiyuan 030006, China 3State Key Laboratory of Quantum Optics and Quantum Optics Devices, Shanxi University, Taiyuan 030006, China 4Collaborative Innovation Center of Extreme Optics, Shanxi University, Taiyuan, Shanxi 030006, China ABSTRACT In this work, the photometric data from the American Association of Variable Star Observers are collected and analyzed on the SX Phoenicis star DY Pegasi (DY Peg). From the frequency analysis, −1 −1 we get three independent frequencies: f0 = 13:71249 c days , f1 = 17:7000 c days , and f2 = −1 18:138 c days , in which f0 and f1 are the radial fundamental and first overtone mode, respectively, while f2 is detected for the first time and should belong to a nonradial mode. The O − C diagram of the times of maximum light shows that DY Peg has a period change rate (1=P0)( dP0= dt) = −(5:87 ± 0:03) × 10−8 yr−1 for its fundamental pulsation mode, and should belong to a binary system that has an orbital period Porb = 15425:0 ± 205:7 days. Based on the spectroscopic information, single star evolutionary models are constructed to fit the observed frequencies. However, some important parameters of the fitted models are not consistent with that from observations. Combing with the information from observation and theoretical calculation, we conclude that DY Peg should be an SX Phoenicis star in a binary system and accreting mass from a dust disk, which was the residue of its evolved companion (most probably a hot white dwarf at the present stage) produced in the asymptotic giant branch phase. Further observations are needed to confirm this inference, and it might be potentially a universal formation mechanism and evolutionary history for SX Phoenicis stars. 1. INTRODUCTION tudes (≥ 0m:1; see Fu et al.(2008)). There are several SX Phoenicis (SX Phe) stars, a subgroup of the high- scenarios proposed to illustrate the formation mecha- amplitude δ Scuti stars (HADS), are old Population II nism and evolutionary history of SX Phe stars (see, e.g., stars. They always pulsate in single or double radial Rodr´ıguez& L´opez-Gonz´alez(2000)), but the origin of modes (such as SW Ser, AE UMa, etc.), but some also them is still unknown up to now. show nonradial modes coupling with the radial modes. DY Peg is an SX Phe star with a low metallicity Because of the insufficient amount and the generally ([Fe=H] = −0:8, Burki & Meylan(1986) and Hintz et al. poor photometric precision of the observation data, (2004); [Fe=H] = −0:56, Pe~naet al.(1999)). The vari- whether any low-amplitude pulsations exist besides the ability of DY Peg was discovered by Morgenroth(1934), dominant radial modes in most SX Phe stars is still un- whereafter a good amount of photometric monitoring known. Although most SX Phe stars, which are char- was done to record and analyze its behavior of light- acterized by high amplitudes of pulsation, low metal- ness variation (see, e.g., Iriarte(1952); Meylan et al. licity, and large spatial motion, are found to be mem- (1986); Percy et al.(2007)). Based on the secular ob- arXiv:2008.02542v3 [astro-ph.SR] 18 Nov 2020 bers of globular clusters (Rodr´ıguez& L´opez-Gonz´alez servations, the period change of DY Peg was continu- 2000), some of them have been discovered in the gen- ously studied in history (see, e.g., Quigley & Africano eral star fields (Rodr´ıguez& Breger 2001). In particu- (1979); Mahdy & Szeidl(1980); Pena & Peniche(1986); lar, pulsations in the majority of the field SX Phe vari- Derekas et al.(2003); Hintz et al.(2004); Derekas et al. ables display very simple frequency spectra with short (2009); Fu et al.(2009)). Li & Qian(2010) did a more periods (≤ 0d:08) and large visual peak-to-peak ampli- detailed research on the period change of this star, in which they reported the variation of the period can be described by a secular decrease of the period at a rate Corresponding author: Jia-Shu Niu of −6:59 × 10−13 days cycle−1, and a perturbation from [email protected] a companion star in an eccentric orbit with a period of 15414.5 days. Unlike the period change that has 2 Xue & Niu been studied adequately, the pulsation frequency of DY The first pulsation mode with frequency f0 = −1 Peg was not detected accurately beside the fundamen- 13:71249 c days and amplitude a0 = 240:3 mmag tal mode. Garrido & Rodriguez(1996) and Pop et al. dominates the light curves of DY Peg. The secondary −1 (2003) reported that DY Peg should be a double-mode pulsation mode with frequency f1 = 17:7000 c days pulsator, while it was not confirmed in subsequent works and a small amplitude a1 = 5:2 mmag is obvious in (Fu et al. 2009; Barcza & Benk}o 2014). the present work, which was not confirmed in previous In the following sections, we extract some impor- works because of the low signal-to-noise (see Garrido & tant information from observations, construct theoret- Rodriguez(1996), Pop et al.(2003), Fu et al.(2009), ical models and present some discrepancies between ob- and Barcza & Benk}o(2014)). servation and theoretical calculation. Then, we propose The ratio of f0=f1 = 0:775 agrees well with the the- some inferences to relieve the discrepancies and back- oretical calculation on the fundamental and first over- track the evolutionary history of DY Peg. This paper is tone radial modes (∼ 0:77, see Petersen & Christensen- organized as follows: Section 2 presents the data reduc- Dalsgaard(1996); Poretti et al.(2005)), illustrating DY tion procedures from observations; theoretical models Peg does pulsate in the two radial modes. are constructed and the calculation results are shown in What is interesting is that a third independent fre- −1 Section 3; Section 4 gives the discussion and conclusions. quency solution with frequency f2 = 18:138 c days and amplitude a2 = 2:8 mmag is detected in this work, 2. OBSERVATIONS AND ANALYSIS for the first time. The frequency f2 is close to f1 with a smaller amplitude, therefore, we suggest that f2 should 2.1. Observations be a nonradial mode.1 For a definite mode identification The time-series photometric data in the V band on of f2, multicolour photometry or time resolved high res- DY Peg is downloaded from the American Associa- olution spectroscopy is needed. tion of Variable Star Observers (AAVSO) International Database (Kafka 2020), which cover from 2003 August 2.3. The O − C Diagram 2 to 2019 December. After the heliocentric corrections of the Julian date and magnitude shifts elimination be- Based on the observations between 2003 and 2019 (see tween different nights, the light curves are used to ex- Figure1), the light curves around the maxima were fit- tract the times of maximum light (see Figure1). A ted by a fourth polynomial. We have obtained 139 times portion of the light curves covering a period of 32 days of maximum light from these light curves, and estimated (from 2011 October 30 to 2011 December 1) are used to their uncertainties via Monte Carlo simulations. The make frequency analysis. Table1 lists the detailed in- newly determined times of maximum light and the un- formation of the observations for the frequency analysis, certainties are listed in Table C1. In Li & Qian(2010), and Figure2 shows the relevant light curves. 412 times of light maximum of DY Peg obtained from photoelectric or CCD data had been collected, which are also used in our O − C analysis. In addition, 138 times 2.2. Frequency Analysis of maximum light in the V band are collected from the The software Period04 (Lenz & Breger 2005) is used literature, which are listed in Table C2. In total, 689 to perform Fourier transformations and frequency pre- times of maximum light spanning 70 years are used to whitenning process for the light curves of DY Peg. Fig- perform the O − C analysis in this work. 3 ure3 shows the spectral window and Fourier amplitude As it has been shown by Li & Qian(2010), a lin- spectra of the pre-whitenning process. The statistical ear or quadratic fit cannot reproduce the times of light criterion of an amplitude signal-to-noise ratio is set to maximum precisely. Consequently, we fit the times of be 4.0 for judging the reality of a newly discovered peak light maximum with a quadratic plus a function of sines, in the Fourier spectra. The noises are determined as the mean amplitudes around each peak with a box of 6 1 c days−1. The ratio of f0=f2 = 0:756 can rule out the assumption that f2 is a radial mode. In total, 14 statistically significant frequencies have 2 Because the amplitude of f0 is about 46 times larger than that been detected, including 3 independent frequencies of f1 and the times of maximum light are dominated by the −1 −1 f mode, O − C method can be used effectively to analyze the (f0 = 13:71249 c days , f1 = 17:7000 c days , and 0 pulsating and orbital parameters in this case. f = 18:138 c days−1), together with 11 harmonics or 2 3 In the following analysis, we give typical uncertainties of 0.0006 linear combinations of them. The solid curves in Fig- days and 0.0005 days to the times of maximum light in the liter- ure2 show the fits with the multi-frequency solution ature detected by photoelectric photometers and CCD cameras, which is listed in Table2.