This research paper presents a comprehensive analysis of RR Lep and BF Vel two southern shortperiod semidetached oscillating Algols oEA stars that have been shown to be triple systems. We determined the spectral types of the primary components and calculated the radial velocities from spectra observed with the Australian National Universitys 2.3 m telescope and Wide Field Spectrograph. The spectra of the Na I D doublet confirmed the presence of tertiary components that were apparent in the broadening function analyses. During primary eclipses chromospherical activity in their secondary components was apparent in their HSUBSUB spectra it was also apparent in the Na I D spectra of BF Vel with its almost complete annular eclipse. Groundbased telescopes were used for observations in several pass bands for photometric analyses. These observations were complemented by data from the Transiting Exoplanet Survey Satellite TESS mission allowing us to model the light curves followed by a detailed analysis of pulsations. Eclipsetiming variation ETV analyses of both systems were used to determine the most likely mechanisms modulating the orbital period. We found mass values of MSUB1SUB  2.9 MSUBSUB and MSUB2SUB  0.75 MSUBSUB for the components of RR Lep along with MSUB1SUB  1.93 MSUBSUB and MSUB2SUB  0.97 MSUBSUB for those of BF Vel. By integrating information from photometry spectroscopy and ETV analysis we found that tertiary components revolve around both systems. The primary star of RR Lep pulsates in 36 frequencies of these five were identified as independent modes with the dominant one being 32.28 dSUP1SUP. The pulsating component of BF Vel oscillates in 37 frequencies with the frequency 46.73 dSUP1SUP revealed as the only independent mode. For both systems a number of frequencies were found to be related to the orbital frequency. Their physical properties were compared with other oEA stars in massradius and HertzsprungRussell diagrams while the pulsational properties of their  Sct components were compared with currently known systems of this type within the orbitalpulsation period and loggpulsation period diagrams.
['binaries: close', 'binaries: eclipsing', 'stars: fundamental parameters', 'stars: individual: BF Vel', 'stars: individual: RR Lep', 'stars: variables: δ Scuti', 'Astrophysics - Solar and Stellar Astrophysics']
Modelling of eclipsing binary systems with pulsating components and tertiary companions BF Vel and RR Lep
{'Modelling of eclipsing binary systems with pulsating components and tertiary companions: BF Vel and RR Lep': 'A. Liakos 1 , D. J. W. Moriarty 2 , 3 , A. Erdem 4 , 5 , J. F. West 3 , and P. Evans 6  \n- 1 Institute for Astronomy, Astrophysics, Space Applications and Remote Sensing, National Observatory of Athens, Metaxa & Vas. Pavlou St., GR-15236, Penteli, Athens, Greece  \ne-mail:  \nalliakos@noa.gr  \n- 2 School of Mathematics and Physics, The University of Queensland, Queensland 4072, Australia\n- 3 Astronomical Association of Queensland, St. Lucia, Queensland, 4067, Australia\n- 4 Astrophysics Research Center and Ulupınar Observatory, Çanakkale Onsekiz Mart University, TR-17100, Çanakkale, Türkiye\n- 5 Department of Physics, Faculty of Arts and Sciences, Çanakkale Onsekiz Mart University, Terzio˘glu Kampüsü, TR-17100, Çanakkale, Türkiye\n- 6 El Sauce Observatory, Coquimbo Province, Chile  \nReceived XX January 2024; accepted XX January 2024', 'ABSTRACT': "This research paper presents a comprehensive analysis of RR Lep and BF Vel, two southern short-period semi-detached oscillating Algols (oEA stars), which are shown to be triple systems. Spectral types of the primary components were determined and radial velocities calculated from spectra observed with the Australian National University's 2.3 m telescope and Wide Field Spectrograph. Spectra of the Na I D doublet confirmed the presence of tertiary components which were apparent in the broadening function analyses and, with H α spectra during primary eclipses, indicated chromospherical activity in their secondary components. Ground-based telescopes were used for observations in several pass bands for photometric analyses. These observations were complemented by data from the Transiting Exoplanet Survey Satellite (TESS) mission to enable the modelling of the light curves, followed by a detailed analysis of pulsations. Eclipse-timing variation (ETV) analyses of both systems were used to determine the most likely mechanisms modulating the orbital period. We found mass values M 1 = 2 . 9 M ⊙ and M 2 = 0 . 75 M ⊙ for the components of RR Lep, and M 1 = 1 . 93 M ⊙ and M 2 = 0 . 97 M ⊙ for those of BF Vel. By integrating information from photometry, spectroscopy and ETV analysis, we found that tertiary components revolve around both systems. The primary star of RR Lep pulsates in 36 frequencies, of which five were identified as independent modes, with the dominant one being 32.28 d -1 . The pulsating component of BF Vel oscillates in 37 frequencies, with the frequency 46.73 d -1 revealed as the only independent mode. For both systems, many frequencies were found to be related to the orbital frequency. Their physical properties were compared with other oEA stars in Mass-Radius and Hertzsprung-Russell diagrams, and the pulsational properties of their δ Sct components were compared with currently known systems of this type within the orbital-pulsation period and log g -pulsation period diagrams.  \nKey words. binaries: eclipsing - stars: fundamental parameters - binaries: close - stars: oscillations - stars: variables: delta Scuti stars: individual (BF Vel) - stars: individual (RR Lep)", '1. Introduction': "The δ Sct stars exhibit rapid and multi-periodic pulsations occurring within a frequency range of 4-80 d -1 . Their pulsation behaviour is predominantly driven by radial and low-order nonradial p -modes, triggered by the κ -mechanism (cf. Eddington 1926; Baker & Kippenhahn 1962; Zhevakin 1963; Breger 2000; Aerts et al. 2010; Balona et al. 2015, and references therein). Furthermore, they also pulsate in high-order non-radial modes, influenced by turbulent pressure within the hydrogen convective zone (Antoci et al. 2014). These stars typically have masses ranging from 1.5-2.5 M ⊙ and belong to spectral classes between A and early F. They are primarily situated within the classical instability strip, spanning from the main-sequence to the giants branch (Aerts et al. 2010).  \nBinary stars are stellar systems comprising two gravitationally bound components orbiting a common center of mass. They serve as crucial astrophysical tools for determining fundamental stellar characteristics, especially mass, radius, and temperature. Close binaries, characterized by orbital periods ranging from less  \nthan one day up to tens of days, have properties that potentially allow for mass transfer between the components (cf. Kopal 1959; De Loore & Doom 1992; Kallrath & Milone 1999; Hilditch 2001, and references therein). Eclipsing binaries (EBs) constitute a specific subset, where the orbital plane aligns with the observer's line of sight, resulting in mutual eclipses of the components. These eclipses induce periodic brightness variations, with their depth, shape, and duration o ff ering insights into various orbital and physical properties such as inclination, temperature, radii, and mass function (cf. Prša 2018). Furthermore, longterm monitoring of eclipses using the Eclipse-Timing Variation (ETV) technique reveals orbital period modulations attributable to physical mechanisms like mass transfer or loss. Spectroscopic binaries (SBs) exhibit detectable periodic Doppler shifts in their spectra due to the components' orbital motion around the common center of mass. These spectral line shifts enable determination of the components' radial velocities and hence their mass ratio. SBs are classified as single-lined (SB1) or double-lined (SB2) respectively, depending on whether absorption or emission lines from one or both components are detected. Notably,  \nSB1(2) + EB systems are particularly significant as they provide direct measurements of both orbital and absolute physical parameters of the binary components.  \nMkrtichian et al. (2002) coined the term 'oEA stars' (oscillating EBs of Algol type) to describe eclipsing pairs with a δ Sct mass accretor component of spectral type (B)A-F. The impact of mass accretion on pulsations remains a key unresolved question in this branch of asteroseismology (cf. Tkachenko et al. 2009; Bowman et al. 2019; Mkrtichian et al. 2022). The association between orbital ( P orb) and dominant pulsation ( P pul) periods for binaries, including those with a δ Sct member, was first noted by Soydugan et al. (2006a). Following a comprehensive six-year survey of candidate systems, Liakos et al. (2012) refined this relationship with an expanded dataset. Zhang et al. (2013) presented the initial theoretical endeavor to establish this connection. Further progress came with the work of Liakos & Niarchos (2017), who compiled an updated catalogue containing approximately 200 known cases, from which only 20 semi-detached and 15 detached systems had precise absolute parameters. They also proposed the existence of subgroups exhibiting distinct properties based on the Roche geometry of their host systems and a boundary of 13-days in the orbital period, beyond that the pulsations are not a ff ected by the companion's presence. Kahraman Aliçavu¸s et al. (2017), focusing solely on EBs, suggested an approximately doubled threshold for P orb.  \nThe exceptional precision (in the order of 10 -4 mag) and continuous time coverage of photometric data from the Transiting Exoplanet Survey Satellite mission (TESS; Ricker et al. 2015) o ff er a valuable resource for in-depth investigations of stellar oscillations. This capability is particularly potent when coupled with multi-band photometry and high-resolution spectroscopy. The temporal resolutions of data collected by the TESS missions range from 20 seconds to 30 minutes, proving highly e ff ective for identifying short-period frequencies, a characteristic feature of δ Sct stars. Notably, owing to their precision, they also enable the detection of low-amplitude pulsations at the level of a few µ mag (cf. Kurtz et al. 2020; Kim et al. 2021; Liakos et al. 2022). Furthermore, these missions a ff ord extensive temporal coverage, sometimes spanning years, which proves invaluable for determining the timing of minima in EBs. This, in turn, facilitates a more thorough investigation of their orbital period variations.  \nAccording to the catalogue of Liakos & Niarchos (2017) and the updated list hosted by AL (which will be published in a future work), there are currently more than 1000 known cases of binaries with δ Sct components. The vast majority (more than 90%) of these systems lack information about the physical and (or) detailed pulsation parameters of their oscillating members. At present, there are approximately only 100 cases with well determined parameters (i.e. they are of the SB1(2) + EB types) and they are distributed almost equally into detached and semidetached systems, with less than 50 cases per group. From these systems, only 33 oEA stars are SB2 + EB systems. Therefore, new studies providing accurate results regarding the absolute parameters of the pulsating stars accompanied with high photometric data from space missions are extremely important for Asteroseismology.  \nThis paper is a continuation of the detailed study on individual southern oEA stars with a δ Sct component. We report analyses of two EBs, namely RR Leporis and BF Velorum. These systems are previously known EBs exhibiting pulsations (see next paragraphs for historical elements), but the absolute properties of their components have not been determined accurately to date. Therefore, since both of them were found to be  \nSB2 + EB systems, this work enriches the sample of oEAs with well-determined absolute stellar properties by 6%.  \nRR Leporis had been observed photoelectrically by Bookmyer et al. (1986), Abhyankar & Vyas (1989), Vyas & Abhyankar (1989), and Samec et al. (1989). In the latter work there was a slight evidence for pulsational behaviour, which was confirmed later by the CCD observations of Dvorak (2009). Liakos & Niarchos (2013) using systematic multi-filter CCD observations, obtained full light curves from which they derived an orbital model solution and determined the dominant pulsation frequency as 33.3 d -1 . The light curves of Vyas & Abhyankar (1989) were analysed by Erdem & Öztürk (2016), who found that orbital period changes of the system resulted in a downward parabolic behaviour superimposed by a periodic curve. The latter was attributed to the possible existence of a low-mass tertiary component. Kahraman Aliçavu¸s et al. (2024), using Échelle spectra, reported the system to be a single lined binary and calculated the radial velocity curve of the primary component; they calculated a photometric model of the system using TESS data. They derived mass values of 1.79 M ⊙ and 0.5 M ⊙ for the primary and secondary components of the system, respectively. Moreover, they found 14 pulsation frequencies, from which six were identified as independent modes.  \nThe system BF Velorum is a rather neglected system. The only available light curves come from Manimanis et al. (2009). In that work, the first photometric model of this system was presented and the pulsational behaviour of the primary component was described.  \nFor these two systems we have undertaken extensive data collection, including medium resolution spectroscopic data (Sect. 2), as well as multi-band ground-based and space-borne photometric observations (Sect. 3). We aimed to: a) model the orbital properties of the systems (Sect. 4); b) Determine stellar parameters with high precision (Sect. 5); c) Perform EclipseTiming Variation (ETV) analysis to identify the orbital period modulating mechanisms (Sect. 6); d) Detect pulsation frequencies and estimate oscillation modes (Sect. 7). The findings and conclusions derived from these analyses are discussed in Sect. 8.", '2. Spectroscopy': "Spectra were observed with the stellar aperture on the wide field spectrograph (WiFeS) mounted at the Nasmyth focal position on the Australian National University's 2.3 m telescope at Siding Spring Observatory (SSO), Australia. The detailed observations log is given in Table 1. A brief overview of the spectrograph is given here; for full details see Dopita et al. (2007, 2010). The 560 nm dichroic beam splitter was used. Fairchild Imaging CCD cameras with a resolution of 4096 × 4096 pixels and a pixel scale of 0.5 arcsec were used; one was optimised for the blue channel and the other for the red channel. The Ne-Ar arc lamp was used for wavelength calibration; the precision of the wavelength calibration was improved by observing arc spectra before and after sets of target spectra. Spectra were reduced with P ywifes , a data reduction pipeline that was written for the WiFeS (Childress et al. 2014).  \nOne dimensional spectra were extracted from the 3D data cubes produced by the P ywifes routines using a zero pixel radius around the spatial pixel where the target was centred and visualised in Q fitsview 1 . This procedure was used also to avoid pixels with noise from probable cosmic ray hits that  \nTable 1: Log of spectroscopic observations with the wide field spectrograph on the ANU 2.3 m telescope.  \n4050 4100 4150 4200 4250 4300 4350 Wavelength (A)  \n<!-- image -->  \nFig. 1: Rectified spectra of the primary components of BF Vel and RR Lep and for comparison, the A5V star HD 23194.  \nwere not fully removed in the reduction routines. The sky background flux values were generally about two orders of magnitude lower than those of the target spectra. Flux units are ergs cm -2 s -1 arcsec -2 Å -1 . See Moriarty et al. (2019) and Liakos et al. (2022) for more details of the equipment and procedures used.", '2.1. Spectral classification': "Spectra for classification were observed with the B3000 grating ( R = 3000) which has a wavelength range between 3200 and 5900 Å that includes the important Ca II K line. The spectra of the primary components of BF Vel and RR Lep were observed with relatively short exposures during the secondary eclipse (phase 0.5) to minimise interference by their secondary and tertiary components.  \nSpectral types were determined using the X class programme and the M kclass 2 reference libraries (Gray & Corbally 2014). X class provides a direct visual comparison of the target spectrum with reference spectra from the M kclass libraries, which encompass integer spectral temperature types from O6 to M5, and a wide range of luminosity classes between V and Ia. The target spectra were normalised to unity at 4503 Å to match the library spectra. The spectra were classified initially by visual comparisons with selected reference spectra across the wave-  \nength range of 3800-4600 Å. In addition, careful reference was made to the structure of the Balmer lines, and to several individual metallic absorption lines, including Ca II K, Ca I 4226, and Mg II 4481. They were also compared to other A5V reference spectra, such as that of HD 23194, drawn from the MILES Library, and the examples provided in the work of Gray & Corbally (2009).  \nThe best fit for the primary components of BF Vel and RR Lep was A5V. The spectra were rectified to deliver a linear continuum baseline using the 'Autorectify' function in the XMK25 display routine in X class (see Fig. 1).", '2.2. Radial velocities': "The B7000 and R7000 gratings were used to capture spectra for radial velocity determinations and for indications of chromospherical activity during the primary eclipse of BF Vel. Exposure times with the B7000 grating were long enough to allow detection of the faint secondary stars' metal lines. At the same time as long exposure times were used with the B7000 grating, shorter exposure times were used with the R7000 grating to minimise blending and enhance the separation of Na I D lines of the primary star from those of the tertiary component in each system.  \nRadial velocities (RVs) of the systems and the binary components were calculated with the broadening function method implemented in the R avespan code written by B. Pilecki (Ruci'nski 1992, 2002; Pilecki et al. 2017) and are listed in Table A.1.  \nFig. 2: Examples of broadening function curves. The orbital velocity of the primary components is marked with a blue line and that of the secondary components with a magenta line. The systemic velocity is indicated by the dashed line. (a) BF Vel phase 0.28, 2020 Nov. 1, exp. time = 240 s. (b) BF Vel phase 0.73, 2020 Nov. 2, exp. time = 240 s. (c) RR Lep phase 0.25, 2020 Dec. 1, exp. time = 240 s. (d) RR Lep phase 0.76, 2020 Nov. 3, exp. time = 420 s.  \n<!-- image -->  \nTemplate spectra were taken from the synthetic spectra library of Coelho et al. (2005). The procedures are described in general by Moriarty et al. (2019) and Liakos et al. (2022), with R avespan settings chosen specifically for the systems in this paper as follows. The spectra with a wavelength range 4600-4840 Å and 4880-5540 Å were normalised in R avespan ; a preference setting of velocity range 1 . 6 + 60 was chosen as that gave the best fit of broadening function peak velocities to velocities observed in the Na I D spectra; other preference settings were the gravity coe ffi cient 3.0, metallicity 0 and resolution 2 km s -1 . The temperature settings of the synthetic spectrum templates were 5250 K for BF Vel and 5000 K for RR Lep. These settings excluded the Balmer lines and far blue spectral region to emphasise the lines of the late spectral type secondary component of each system.  \nThe presence of a tertiary star in each of these systems was evident in the broadening function curve of the primary component, either as a double peak with one peak having a low or zero velocity or as an asymmetrical curve (see Fig. 2). The RVs of the secondary and tertiary components of BF Vel were shown clearly at the quadrature phases in Fig. 2a, b. The RVs of the secondary component of RR Lep were clearly evident in the broadening function analyses, whereas broadening function peaks of the primary and tertiary components of RR Lep were merged (Fig. 2c, d). As described below, spectra of the Na I D doublet and H α were also used to determine RVs and provide initial values of the primary components' velocities as inputs in R avespan .", '2.3. H α spectra': "The systemic velocity of each system determined with radial velocity analyses was subject to uncertainty due to blending of the primary components' spectral lines with those of the tertiary components. As the spectral class of the primary components is A5V, their Balmer lines are much stronger than those of the secondary and tertiary components which have later spectral  \n<!-- image -->  \nFig. 3: H α spectra at phase 0.5 (secondary eclipse). (a) BF Vel, exp. time = 240 s; (b) RR Lep, exp. time = 60 s.  \n<!-- image -->  \ntypes, indicated by their sodium spectra (see Sect. 2.4). Therefore, spectra during secondary eclipses were observed to provide independent values for the systemic velocities (see Fig. 3). With relatively short exposures, the Balmer lines at phase 0.5 were dominated by those of the bright A5V primary stars, whereas during primary eclipses, the Balmer lines are a composite of the three components in each system. The H α line centre of BF Vel was at 6563.1 Å, that is at 14 km s -1 , which is a little above the range of 4 . 7 ± 1 km s -1 determined by the modelling (Sect. 4). The H α line is broadened by the combination of the primary component's rapid rotational velocity and blending with that of the tertiary star. The H α line centre of RR Lep was at 6561.7 Å, indicating a systemic velocity of ∼ -50 km s -1 . The uncertainty of the velocities determined from the H α spectra is about ± 5 km s -1 due to the low-to-medium resolution of R = 7000 of the wide field spectrograph and the broadening of the Balmer lines by the rapid rotational velocities of the primary components of the systems (Sect. 4). The H α line of RR Lep is broadened to the red due to the tertiary component.", '2.4. Sodium D doublet spectra': "Spectra of the Na I D doublet were used to check the velocities calculated in R avespan . At quadrature phases, the lines of the binary pairs are separated from those of the tertiary components to an extent dependent on their orbital velocities and length of the observation exposures. The sodium D lines of the tertiary star in the BF Vel system were apparent, although partially blended with those of the primary component, at phase 0.75 with an exposure time of 300 s (Fig. 4a). Note at that phase, the D2 line of the secondary component of BF Vel was not distinguishable from noise. One night later, at the opposite quadrature and with a shorter exposure of 180 s, there was a distinct separation of the tertiary star's D2 line from that of the primary star; the tertiary D1 line was also noticeable, although partially blended with the lines of the primary and secondary components (Fig. 4b). Similarly, at phase 0.69 with a 300 s exposure, the narrow lines of the tertiary star were blended, but visible, and at phase 0.25 with a shorter exposure, were clearly separated from the primary component's sodium lines (Fig. 4c, d). The D2 line of the secondary component was apparent at phase 0.69. However, its D1 line was not apparent at phase 0.25 on 2020 Mar. 16; see Sect. 2.5 for an explanation of the variability in the secondary component's sodium doublet spectra.  \nWith a maximum radial velocity of -116 ± 10 km s -1 (i.e. RV1 , total = K 1 + V 0; see Tables 3 and A.1) at phase 0.25, the RR Lep primary star's sodium D lines were distinct from those of the tertiary star (Fig. 5a, b). In contrast, at phase 0.75 the  \nFig. 4: Spectra of the BF Vel Na I D doublet at orbital phases ( ϕ ) around quadrature; R = 7000. Downward arrows mark the calculated positions of the line centres for the primary star (S1) and secondary (S2) that give the best fit to the RAVESPAN analyses. Lines of the tertiary component are indicated with magenta arrows angled up. Exposure times: (a) 300 s; (b) 180 s; (c) 300 s; (d) 105 s.  \n<!-- image -->  \npositive radial velocity of the primary star was close to that of the tertiary star, resulting in blending of its spectral lines with those of the tertiary star. The lines of the primary star around phase 0.75 were broadened to the red of the sharp lines of the tertiary star (Fig. 5c, d). The Na I D lines of the tertiary component were strong relative to those of the primary component at phase 0.17 (Fig. 5e). The tertiary H α spectrum at phase 0.17 was blended with that of the primary component, although noticeable with broadening to the red (Fig. 5f). The secondary component's flux was too weak for its H α line to be apparent.", '2.5. Chromospherical activity': "The spectrum of BF Vel during a primary eclipse is dominated by the secondary component, which has almost the same diameter as the primary star. With an orbital inclination of close to 90 · , the eclipse is annular. Emission in the H α line centre of BF Vel was observed during the primary eclipse on 2017 Feb. 14 (Fig. 6a); the H α line centre was infilled on 2020 Mar. 15 (Fig. 6b). Emission was also observed on the same dates in the Na I D line centres of BF Vel (Fig. 6c, d). These spectra indicate the secondary component is chromospherically active. Evidence of chromospherical activity in the secondary component of BF Vel was also apparent in the variability in its sodium D spectral lines at other phases (Fig. 4). Its Balmer lines are weak compared to those of the primary star, whereas its sodium D lines are stronger. The spectra also include blending of the tertiary component's lines.  \nIn contrast to BF Vel, the spectral lines of the primary star of RR Lep contribute to the spectrum during a primary eclipse as the orbital inclination is ∼ 81 · . Spectra in the range 5700-  \nFig. 5: Spectra of the Na I D doublet of RR Lep during quadrature phases (a-e) and a corresponding H α spectrum at phase 0.17 (f); R = 7000. Exposure times: (a, b) 240 s; (c) 360 s; (d) 120 s; (e, f) 240 s.  \n<!-- image -->  \n-6700 Å were reduced with a wavelength interval of 0.42 Å to enhance resolution in the Na I D and H α regions. The Na I D lines of the tertiary component in RR Lep are strong and distinct from those of the binary components (Fig. 7a, c). The infilling in the H α line of RR Lep at phases 0 and 0.5 indicates that both the primary and secondary components are chromospherically active (Fig. 7b, d). The strength of these lines indicates that the spectral type of the tertiary star is about G8.", '3.1. Ground-based photometry': "The system RR Lep was observed with the 40 cm Cassegrain telescope of the University of Athens Observatory, Greece which is equipped with the ST-10XME CCD. The observations were made in the B and V pass bands (Bessell specification) in nine nights between January and March 2012. The light curves (LCs) have been already published in Liakos & Niarchos (2013) but we  \nFig. 6: Spectra of BF Vel during primary eclipses. Indications of chromospherical activity in the secondary component of BF Vel during primary eclipse with emission in the H α line centre (a, b) and emission in the Na I D lines on the same dates (c, d). Exposure times were 60 s.  \n<!-- image -->  \nFig. 7: Spectra of RR Lep during eclipses. Downward arrows mark the calculated position of the combined binary components' lines. (a) The magenta arrows indicate the line centres of the tertiary star's Na I D doublet. (b) The H α line is partly infilled broadened to the red; it is marked with a magenta arrow. (c) Secondary eclipse; the Na I D lines of the primary star broadened to the red by the tertiary star (magenta). (d) The H α line broadened to the red. Exposure times were 60 s.  \n<!-- image -->  \nincluded them in the present work as additional information for the modelling.  \nThe Evans 0.51 m telescope at El Sauce Observatory in Chile was used to observe BF Vel in BVRI Johnson-Cousins filters  \non 26 nights, commencing 13 Dec. 2022. The telescope was equipped with a Moravian C3-26000 CMOS camera. The first five data sets were observed at 2-min cadence in alternating B and V filters. Subsequent observations for B , V , and R c filters were at 1-min cadence and I c filter observations at 90-sec cadence. Standard calibration di ff erential photometry was then performed on each data set using 'AstroImageJ V5' (Collins et al. 2017) resulting complete LCs for B , V , and R c filters and an almost complete LC for the I c filter.", '3.2. TESS photometry': 'RRLep (TIC 169532543) was observed with a 2-min cadence in Sector 5 (Nov. 15 to Dec. 11, 2018) and with a 10-min cadence in Sector 32 (Nov. 20 to Dec. 16, 2020). BF Vel (TIC 190894967) was observed with a 30-min cadence in Sector 8 (Feb. 2 to 27, 2019), with a 30-min cadence in Sector 9 (Feb. 28 to Mar. 25, 2019), with a 2-min cadence in Sector 35 (Feb. 9 to Mar. 6, 2021), and with a 200-s cadence in Sector 62 (Feb. 12 to Mar. 10, 2023). These TESS data were downloaded from the Mikulski Archive for Space Telescopes (MAST) (cf. Jenkins et al. 2016). The straight Simple Aperture Photometry (SAP) data were used, since the Pre-search Data Conditioning Simple Aperture Photometry (PDCSAP) is optimised for planet transits, and it appears that the PDCSAP detrending sometimes leads to additional e ff ects. The LCs from Sector 5 for RR Lep and Sector 35 for BF Vel were used to ensure that a reasonably high Nyquist frequency could be sustained for the pulsation analysis. The Sector 35 includes the data in two orbits, as in other sectors. When comparing the data taken in these two orbits, the primary minimum depths are slightly di ff erent. This a ff ects the residual LCs (and therefore the pulsation analysis) when the entire data set is resolved. Therefore, these two data obtained in two orbits were solved separately. Photometry of these two binary systems is publicly accessible (Ricker et al. 2015).', '4. Modelling radial velocity and light curves': "The numerical integration method of Wilson & Devinney (1971, WD) was applied to solve the LCs and RV curves of BF Vel and RR Lep. This method models the RV curve and (or) the LC of a given binary star, considering mainly the ellipticity and proximity e ff ects, and the equipotential surfaces of the components. Radial velocities were not determined with R avespan close to the orbital phases 0.0 and 0.5, where the components occult each other and proximity e ff ects are dominant.  \nIn the WD code (and most LC modelling programmes) the e ff ective temperature ( T 1) of the primary star is usually taken as a fixed parameter when solving the LCs of a given binary system. Therefore, reliable determination of T 1 is important. We examined three methods to determine the e ff ective temperature of the primary components of BF Vel and RR Lep.  \n(i) To calculate the intrinsic (unreddened) colours of both systems, the colour indices were first derived from the study of Høg et al. (2000) to be B -V = 0 . 13 ± 0 . 08 mag for BF Vel and B -V = 0 . 15 ± 0 . 04 mag for RR Lep. Then, the colour excesses of E ( B -V ) = 0 . 095 ± 0 . 015 mag for BF Vel and E ( B -V ) = 0 . 078 ± 0 . 010 mag for RR Lep were computed using the stellar absorption values of AV = 0 . 295 ± 0 . 048 mag for BF Vel and AV = 0 . 242 ± 0 . 030 mag for RR Lep in the Gaia DR3 catalogue (Gaia Collaboration et al. 2022). Finally, the intrinsic colours were calculated to be ( B -V )0 = 0 . 035 ± 0 . 081 mag for BF Vel and ( B -V )0 = 0 . 072 ± 0 . 041 mag for RR Lep. When compared with the calibration data sets of Pecaut & Mamajek  \nA. Liakos et al.: Modelling of eclipsing binary systems with pulsating components and tertiary companions: BF Vel and RR Lep  \nTable 2: Results of optimal curve-fitting to the radial velocity and light curves of BF Vel using WD modelling.  \n(2013) and Eker et al. (2018), these intrinsic colours indicate that the spectral types of these systems would be close to A1-2V, and therefore are not compatible with our spectral results.  \n(ii) Kahraman Aliçavu¸s et al. (2024) reported the temperature of the primary component of RR Lep as 7800 ± 150 K from the analysis of H β lines in the HERMES spectra, which they observed at phases close to the secondary eclipse.  \n(iii) In Sect. 2.1 of the present study, the spectral type of their primary components was assigned as A5V using the line matching method. From the calibration data sets of Pecaut & Mamajek (2013) and Eker et al. (2018), the temperature corresponding to the A5V spectral type is 8100 K and 8200 K, respectively. Therefore, the TESS LCs of BF Vel were solved by using the temperatures of 8100 K and 8200 K for T 1, while the TESS LCs of RR Lep were solved with temperatures of 7800 K, 8100 K, and 8200 K for T 1. Then, the temperature with the smallest chisquare value (i.e. 8100 K) was adopted and kept fixed for the primary components. The temperature of the secondary components was adjusted in the WD iterations.  \nOther parameters that were adjusted in the WD iterations included: the semi-major axis ( a ), mass ratio ( q ), and orbital inclination ( i ) of the binary system, phase shift ( ∆ ϕ ; which allows the WD code to adjust for a zero point error in the ephemeris used to compute the phases; the unit is the orbital period), systemic velocity of the binary ( V 0), non-dimensional surface potential parameter of each component ( Ω 1, Ω 2), fractional luminosity of the primary component ( L 1), and third light contribution to the total light of the system ( L 3).  \nA quadratic limb-darkening law was assumed; the limbdarkening coe ffi cients ( x i and y i) of the components were taken from Claret (2017) for the TESS LCs. For the ground-based LCs, the limb darkening coe ffi cients were taken from the tables of Claret & Bloemen (2011), according to the e ff ective temperature of the components and the filters used. The bolometric gravity darkening exponents were adopted from Lucy (1967) to be 0.32 for convective atmospheres of the components ( T < 7200 K) and from von Zeipel (1924) to be 1.0 for radiative atmospheres of the components ( T ≥ 7200 K). The bolometric albedos of the components were set to 0.5 for convective atmospheres and to 1.0 for radiative atmospheres (Ruci'nski 1969). These parameters were kept constant during all iterations.  \nThe iterations were started in MODE 2 in the WD code, which corresponds to a detached configuration; however, the final simultaneous RV + LC solution was switched to MODE 5, which applies to semi-detached binaries where the secondary component fills its Roche lobe. Therefore, the cited Ω 2 retains its lobe-filling value and was not adjusted in this model.  \nThe parameters of the final WD models for the RVs and ground-based and TESS LCs of BF Vel and RR Lep are listed in Tables 2 and 3, respectively. The first two rows of these tables contain the ephemeris elements used for phasing both the LC and RV data. Since the amplitudes of the RVs, K 1 and K 2, are not directly given by the WD method, they were calculated using the relevant parameters in the RV + LCs fitting. The value of χ 2 was obtained from P i ( y i , o -y i , c ) 2 /∆ y 2 i (Bevington 1969), where y i , o and y i , c are the observed and calculated light levels (or  \nTable 3: Results of optimal curve-fitting to the RV curves and light curves of RR Lep using WD modelling.  \nRVs) at a given phase, respectively, and ∆ y i is an error estimate for the measured values of y i , o. The reduced χ -squared value was obtained from χ 2 red = χ 2 /ν , where ν is the number of degrees of freedom of the data set, that is ν = N obs -N param, where N obs is the number of observations and N param is the number of parameters in the model. The comparison of the RV measurements and photometric observations with the WD models are shown in Figs. 8, 9, 10, and 11 for BF Vel and in Figs. 13, 14, and 15 for RR Lep. The Roche geometries of BF Vel and RR Lep, acquired using the program B inary M aker (Bradstreet & Steelman 2002), are plotted in Figs. 12 and 16, respectively. In these configurations, which show that BF Vel and RR Lep are semi-detached binary systems, their primary components fill 84% and 85% of their Roche lobe, respectively, while the secondary components fill their Roche lobe completely.  \nIn the WD method, the photometric data observed in various bands for a given binary star can be resolved either separately or all together simultaneously (with the RV data, if any). Due to the systemic uncertainties caused by di ff erent permeabilities of the filters, di ff erent pixel sensitivity and other instruments used, the model obtained from simultaneous solutions of various band photometric data cannot represent the observations better than the model derived from separate band solutions. Since we will perform a pulsation analysis using the residuals obtained after subtracting the eclipsing model from the observed LCs in this study, we solved each band's observed LC separately.  \nAs can be seen in Table 2 for BF Vel, there is no significant di ff erence or discrepancy within the error limits between the model parameters in single-band LC + RVsolutions. However for RR Lep, values for the temperature of the secondary compo-  \nnent ( T 2) resulting from the solution of ground-based photometric observations are higher than those obtained from TESS solutions (according to Table 3). In particular, in the solution of the ground-based B -band light curve, T 2 is ∼ 600 K higher than temperatures given in the solution of TESS data. The CROWDSAP parameter, which is the ratio of the target flux to the total flux in the photometric aperture, is 0.985 for TESS Sector 5 and and 0.986 for Sector 32. These values show that only approximately 1.5% of the observed flux in the SAP aperture does not come from the target (and the analysis would just add that to the third light), so this does not seem to be the reason. On the other hand, when the third light was neglected in the solution of the groundbased B light curve, T 2 decreases to 4900 K, but, unfortunately, the di ff erence was still present. When all light and RV curves were solved simultaneously, a value of 4650 K was obtained for T 2. However, the theoretical curve of this solution did not match the B -band observations particularly well (i.e. the fitting was not good). Therefore it seems that, the higher temperature of the secondary component in the ground-based photometric model may be due to systematic errors in the ground-based data.  \nIn the TESS LC of Sector 62 of BF Vel there is a noticeable asymmetry at its maximum values (see Fig. 11). The light flux at Maximum II ( Φ = 0 . 75) is slightly lower than at MaximumI ( Φ = 0 . 25). Two di ff erent approaches were used to model this asymmetry: (i) The assumption of a cool spot on the secondary component due to the chromospherical activity detected in Sect. 2.5. (ii) The assumption of a hot spot that may form on the primary component due to mass transfer from the secondary to the primary, which causes the orbital period of the system BF Vel to increase in the modelled parabolic change in  \nFig. 8: Radial velocities of BF Vel with the WD model fitting. RVs of the primary and the secondary components are marked as filled and hollow symbols, respectively. Residuals to the model are plotted in the bottom panel.  \n<!-- image -->  \nFig. 9: Ground-based BVI light curves of BF Vel with the WD model fitting. Residuals to the model are plotted on the bottom panel. The V and I light curves and their residuals are shifted downward to enhance visibility.  \n<!-- image -->  \nthe ETV analysis in Sect. 6.2 and fits the semi-detached configuration. The WD code combined with the Monte Carlo (MC) optimization procedure (Zola et al. 2004) were used to solve the TESS Sector 62 LC under these two di ff erent spot approaches. The most important contribution of the MC procedure in the  \nFig. 10: Two light curves of TESS Sector 35 (with a 2-min cadence) obtained in two di ff erent consecutive orbits of BF Vel and the WD model fitting. Residuals to the model are plotted on the bottom panel.  \n<!-- image -->  \nFig. 11: Light curves of TESS Sector 62 (with a 200-s cadence) for BF Vel and the WD model fitting. A solid gray line was added to the light flux value 1.0 to display the asymmetry between maxima in the light curve. Residuals to the model are plotted on the bottom panel.  \n<!-- image -->  \nFig. 12: Roche geometry of BF Vel at orbital phase 0.75. Crosses denote the barycenters of the stars and the system.  \n<!-- image -->  \nFig. 13: Radial velocities of RR Lep with the WD model fitting. RVs of the primary and the secondary components are marked as filled and hollow symbols, respectively. Residuals to the model are plotted in the bottom panel.  \n<!-- image -->  \nFig. 14: Ground-based BV light curves of RR Lep with the WD model fitting. Residuals to the model are plotted on the bottom panel. The V LC and its residuals have been shifted downward to enhance visibility.  \n<!-- image -->  \nWD + MCmethod is the use of a range of input values for the free parameters, rather than fixed input values. Thus, the WD + MC method searches for the model that provides the best fit of the observed LC in the solution space, by conducting hundreds of thousands of iterations in the selected input ranges. A Chi-square test was also performed using the model and observed data, and  \nFig. 15: TESS Sector 5 (with a 2-min cadence; upper panel) and Sector 32 (with a 10-min cadence; lower panel) light curves of RR Lep with the WD model fitting. Residuals to the model are plotted in the bottom figures.  \n<!-- image -->  \nFig. 16: Roche geometry of RR Lep at orbital phase 0.75. Crosses denote the barycenters of the stars and the system.  \n<!-- image -->  \nit was found that the cool spot model reflects the observations better. The results are given in Table 2, and a comparison of the observed LC and the cool spot model synthetic LC is presented in Fig. 11. The spot parameters given in Table 2 are: The spot's latitude ( β ), longitude ( λ ), angular radius ( γ ), and temperature factor ( κ ). As a LC solution only was developed with WD + MC for the Sector 62 data, there are no RV curve solution parameters in the last column of Table 2. When the LCs of TESS Sectors 35 and 62 are compared with each other, it is apparent that the depths of minima are quite di ff erent. The depth of Minimum I ( Φ = 0 . 0) was determined as 0.95 mag for Sector 35 and 1.20 mag for Sector 62, while the depth of Minimum II ( Φ = 0 . 5) was determined as 0.76 mag for Sector 35 and 0.98 mag for Sector 62. The reason for the di ff erent eclipse depths in the TESS data might be caused by contamination within the photometric aperture. The CROWDSAP parameter for the Sector 35 is 0.916, which indicates that is 8.4% of the flux in the aperture does not come from the target. Unfortunately, this parameter is not available for the Sector 62 data. This di ff erence in the WD + MC iterations is attributed to the third light ( L 3) and not to the spot e ff ect for Sector 62 data, although the third light contribution is at the level of 0.11 for Sector 35 data, it reduced L 3 to zero for Sector 62 data. On the other hand, while the fractional luminosity of the primary component ( L 1) was 0.76 for Sector 35 data, L 1 increased to 0.87 for Sector 62 data. Thus, the third light contribution of 0.11, which cannot be seen in Sector 62, was included with L 1 in WD + MCresults. Apart from this, the fact that the other parameters representing the LC for the Sector 62 data were similar to the parameters for Sector 35 data and other band data within error limits, does not change the results of the general solution for BF Vel.", '5. Absolute Parameters and Evolutionary status': "The physical parameters of the components of both systems were calculated based on simultaneous modelling of the V -band LC and RVs (see Sect. 4) in order to simplify calculation of the bolometric magnitude ( M bol = M V + BC V), and then derive the total luminosity of the system ( M bol = M bol , ⊙ -2 . 5 log L / L ⊙ ) and the distance to the system from the distance module ( M V = V + 5 log π + 5 -A V). The physical parameters thus obtained are given in Table 4 along with their errors. The nominal solar values, which were adopted by IAU 2015 Resolutions B2 and B3, were used in our calculations. Bolometric corrections for the components were taken from the study of Eker et al. (2018), according to their e ff ective temperatures.  \nThe results of the RV + LC analyses show that the total luminosity of RR Lep is 28.17 L ⊙ and that of BF Vel is 12.43 L ⊙ (see Tables 2, 3, and 4). In the Gaia DR3 catalogue, the system luminosities produced from the FLAME module are given as 37.32 L ⊙ for RR Lep and 12.74 L ⊙ for BF Vel.  \nAn absolute magnitude of M G = 1 . 97 mag for BF Vel was calculated using the G -band mean magnitude m G = 10 . 77 mag, parallax π = 1 . 9776 mas from the Gaia DR3 catalogue and the interstellar absorption (or extinction) A G = 0 . 275 mag, produced from the ESP-HS module (Creevey et al. 2023). Based on these data, the BF Vel system luminosity is 12.78 L ⊙ , assuming zero bolometric correction. The system's luminosity and the distance value calculated from the V -band LC + RV solution for BF Vel (12.43 L ⊙ and 514 ± 32 pc) agree with the values given in the Gaia DR3 catalog (12.74 or 12.78 L ⊙ and 506 ± 7 pc).  \nIn a similar manner for RR Lep, using m G = 9 . 94 mag, π = 2 . 4453 mas, and A G = 0 . 223 mag, we derived M G = 1 . 66 mag, and a system luminosity of 17.09 L ⊙ , assuming zero bolomet-  \nTable 4: Absolute parameters of RR Lep and BF Vel.  \ncorrection. However, the system's luminosity and distance value obtained from the V -band LC + RV solution (28.17 L ⊙ and 537 ± 42 pc) do not match the values given in the Gaia DR3 catalogue (37.32 or 17.09 L ⊙ and 409 ± 3 pc). It should be noted that the interstellar absorption values from the FLAME module and the ESP-HS module in the Gaia DR3 catalogue are quite di ff erent (0.789 and 0.223 mag, respectively), and similarly, the RR Lep system's luminosities also show a large di ff erence. Furthermore, there is still a di ff erence of approximately 100 pc between the distance calculated from our V -band LC + RV solution and that from the Gaia DR3 parallax. If we substitute the interstellar absorption used in our calculations with the value of 0.8 mag, as calculated from the FLAME module, this di ff erence is eliminated. We conclude that the aforementioned discrepancies arise due to the strength of the tertiary component's spectral lines, which were blended with those of the primary component at most phases (Sect. 2). Around the first quadrature, at phases 0.17 and 0.25, spectra of the tertiary component were clearly evident as its sodium D lines were considerably stronger than those of the primary star.  \nThe physical properties of the primary components of these systems and their pulsational behaviour (Sect. 7) indicate that they are δ Sct stars. In order to compare their properties with similar cases, we used mostly the sample of Liakos & Niarchos (2017) for the oEA stars and various works for individual cases from the literature. The gathered sample contains 96 oEA systems (including BF Vel and RR Lep), however there is information about the absolute parameters of only 72 of these. In particular, 33 of them are double-lined spectroscopic and eclipsing systems (SB2 + EB), which are the most important tools for calculating absolute parameters of stars. Nine of them are SB1 + EB systems and one is SB1 + ellipsoidal system (SB1 + el); the latter systems can be considered as the second most informative cases. The masses of the pulsating stars in the other 29 oEA systems were assumed, based on their spectral type. The locations of the components of RR Lep and BF Vel in the Mass-Radius (M-R) and Hertzsprung-Russell (H-R) evolutionary diagrams are given  \nFig. 17: Mass-Radius (top panel) and Hertzsprung-Russell (bottom panel) evolutionary diagrams. In both diagrams the black symbols denote the δ Scuti members of oEA stars that belong either to spectroscopic single- or double-lined (SB1 or SB2) and eclipsing (EB) or ellipsoidal (el) systems (dots) or only to eclipsing systems (squares). The sample was mainly gathered from Liakos & Niarchos (2017). Magenta and cyan symbols refer to the primary (filled) and secondary (open) components of RR Lep and BF Vel, respectively. Solid black lines indicate the ZAMS and TAMS limits and the coloured solid lines (B = blue, R = red) the boundaries of the instability strip (IS; Soydugan et al. 2006b).  \n<!-- image -->  \nin Fig. 17. The purpose of these diagrams is to a) check the evolutionary stages of both components of each studied system and b) compare the evolutionary stages of their primaries with other δ Sct stars in oEA systems. For the latter reason, we included only the pulsating members of the oEA stars and not their nonpulsating companions.  \nThe position of the primary component of RR Lep in the M-R diagram indicates that it is well inside the main-sequence boundaries and is identified as the most massive δ Sct star of the current sample. The primary of BF Vel is located on the ZAMS boundary in the same diagram. The positions of the secondary components in both diagrams show that they have evolved out of the main-sequence. The primary components lie well inside the classical instability strip as shown in the H-R diagram. As both BF Vel and RR Lep are semi-detached systems, their primary stars are accreting mass from their companions, which have evolved to the sub-giant stage of evolution.", '6. Eclipse Timing Variation analysis': "For the analysis of orbital period variations we employed the 'LITE' code of Zasche et al. (2009). The code uses statistical weights ( w ) on the times of minima (ToM) based on the method followed for the observations of the eclipses. Older ToM derived from visual or photographic data were assigned w = 1, the photoelectric w = 7, and the CCD w = 10. Past ToM of each studied system were obtained from online databases 3 , 4 , 5 , 6 , 7 and were enriched with those derived from our ground-based data, TESS data and data from other ground-based surveys. The ToM used in the following analyses are given in Table B.1.  \nThe code is able to fit the data points with Light-Time E ff ect (LITE; Irwin 1959) and parabolic curves. The LITE curve has seven free parameters, namely the ephemeris of the EB ( T 0 and P ), the amplitude of the variation A , and the orbital period P 3, the time of the periastron passage JD 0, the argument of periastron Ω , and the orbital eccentricity e 3 of the tertiary component. The parabolic fitting has three free parameters; the ephemeris of the binary and the parabolic coe ffi cient C 2. We used ephemerides from Kreiner (2004) for both systems; they were left free to adjust during the iterations. The selection of the final model is based also on the Bayesian information criterion (BIC). The results of the Eclipse-Timing Variation (ETV) analyses for both systems were applied as inputs in the 'InPeVeb' software (Liakos 2015) to calculate the parameters of the most likely orbital modulating mechanisms.", '6.1. RR Lep': "The ETV analysis of RR Lep is based on 140 ToM spanning from the mid 1930's up to 2023. A periodic variation due to a tertiary companion is not particularly obvious in the ETV diagram. However, as a noticeable third light was detected in both the spectroscopic and photometric observations, a LITE curve was selected for fitting the ToM. Given that the system has a semi-detached configuration, a parabola was also tested as a possible solution.  \nThe LITE solution is plotted in the upper panel of Fig. 18 and its parameters are listed in Table 5. This solution suggests that a tertiary component in eccentric orbit with a period of ∼ 74 yr orbits the EB. The mass function of the third body suggests a minimum mass of ∼ 0 . 33 M ⊙ (i.e. in coplanar orbit with the EB). However, following the formalism of Liakos & Niarchos (2012), we found that if the tertiary member is a main-sequence star, then, in order to contribute ∼ 6%of the total luminosity (see Table 3), its mass should be ∼ 1 . 14 M ⊙ . The latter is possible if its inclination is ∼ 19 . 5 · with respect to the EB's orbital plane. Alternatively, a higher inclination angle of the third body's orbit, which would yield a lower mass value, is also possible under the assumption that the tertiary component is an evolved star (i.e. subgiant or giant).  \nThe parabolic fitting resulted in a downward parabola (lower panel of Fig. 18) that suggests an orbital period decrease. This contradicts the conventional semi-detached configuration (i.e. the less massive star transfers mass to the more massive component). However, the derived curvature of the parabola may indicate non-conservative mass transfer ( ˙ M tr) and mass loss ( ˙ M loss) due to magnetic braking (Hilditch 2001; Erdem et al. 2005).  \nFig. 18: Fitting of a LITE curve (upper plot) and of a downward parabola (lower plot) on the ToM of RR Lep. The bigger the symbol the bigger the statistical weight. Open and filled symbols denote the primary and secondary ToM, respectively. Lower panels show the residuals of each solution.  \n<!-- image -->  \nFig. 19: Fitting of a LITE curve (upper plot) and of an upward parabola (lower plot) on the ToM of BF Vel. Symbols have the same meaning as in Fig. 18.  \n<!-- image -->  \nTable 5: LITE parameters of the ETV analysis of RR Lep and BF Vel.Table 6: Mass transfer / loss parameters for both systems as alternative explanations for their ETVs.  \na  \nassumed  \nTherefore using the method of Erdem et al. (2005), we assumed that the secondary component is the mass donor star of the system and transfers mass with a typical rate of 10 -8 M ⊙ yr -1 to the primary, while simultaneously losing mass due to magnetic braking. The gyration constant of this component was assumed as k = 0 . 32, while its Alfvén radius as 10 R 2 (cf. Soydugan et al. 2011; Liakos & Niarchos 2013). Results regarding this solution are given in Table 6.  \nThe BIC values for the proposed solutions resulted as -1075 for the LITE and -1073 for the parabola, indicating that LITE is more appropriate solution. Moreover, as our photometric and spectroscopic data demonstrate the presence of a tertiary component, and a parabolic fitting in the ToM of RR Lep needs many assumptions to explain the ETV, we also conclude that the correct solution for the observed orbital period changes of this system is the LITE.", '6.2. BF Vel': "There are 61 ToM currently available for BF Vel. All of them, except for two, cover the time period after the year 2000; the other two are from the late 1930s and were based on photographic plates. Taking these old ToM into account, we fitted a LITE curve to the data points. The selection of this curve is based on the detection of third light in both spectroscopic and photometric observations. The fitted curve is illustrated in the upper plot of Fig. 19 and its parameters are listed in Table 5. According to this solution, a third body orbits the EB with a period of ∼ 25 yr. The f ( m 3) of this solution derives a minimum mass of ∼ 0 . 36 M ⊙ . Therefore, similar to that discussed above for RR Lep, its inclination with respect to the EB's orbit should be ∼ 24 · in order for its mass as a dwarf to be ∼ 1 M ⊙ ,  \nwhich is su ffi cient to explain the observed ∼ 8% contribution in the system's total luminosity (see Table 2). Alternatively, a less massive, but more evolved star (i.e. more luminous) at a higher inclination angle, is also possible. It should be mentioned that a P 3 of ∼ 36 . 75 yr is also possible. The fit on the data points using this value yielded slightly larger P res 2 (i.e. 0.0022) but also a much smaller f ( m 3) value (i.e. 0.0021 M ⊙ ) that cannot easily explain the observed light contribution from the tertiary component. Therefore we neglected this solution.  \nAs a second solution, the data points could be fitted by an upward parabola, which agrees with the semi-detached status of the system. For this case, the fitting is given in the lower plot of Fig. 19 and the corresponding parameters in Table 6. However, that does not accommodate the presence of the tertiary component. The BIC value of LITE ( -604) is lower than that of the parabolic fitting ( -573). Therefore, the LITE curve is the better solution, not only statistically, but also because it explains the presence of the tertiary component.", '7. Pulsation Analyses': 'The frequency search for both systems relied only on the TESS LC residuals because of their photometric accuracy, time resolution and continuity in time (i.e. no time gaps). Based on the spectral types and the physical properties of the components of both systems, as described in the previous sections, only the primary stars have properties of δ Scuti stars. Therefore, the data points between orbital phases 0.9-0.1, when the pulsating components are eclipsed, were excluded from the frequency analyses. The software PERIOD04 (classical Fourier analysis; Lenz & Breger 2005) was used to extract the pulsation frequencies. Although the typical oscillation range for δ Scuti stars is 4-80 d -1 (Breger 2000; Bowman & Kurtz 2018), we extended the search from 0 to 80 d -1 in order to detect possible δ Scutiγ Doradus hybrid behaviour, or low-frequency modes that are connected to g -mode pulsations or the orbital frequency ( f orb) or other frequency combinations. The noise level of each detected frequency was computed around it with a spacing of 2 d -1 and a box size of 2 (Balona 2014). After the computation of a frequency, the residuals were pre-whitened for the next one. The search was stopped when the signal-to-noise ratio of the detected frequency was approximately 5 (S / N ∼ 5; Baran et al. 2015; Bowman & Michielsen 2021). The TESS data sets used for the analyses were observed in short-cadence mode (i.e. 2-min time resolution) and span over approximately 25 continuous days. The Nyquist frequency for these data sets is ∼ 360 d -1 and the frequency resolution (1 . 5 / T , where T is the observations time range in days; cf. Loumos & Deeming 1978) is ∼ 0 . 06 d -1 .  \nThe results of the analyses are divided into two main parts; the independent frequencies and the combination frequencies. In addition to the Fourier modelling of the pulsations, we also estimated the most likely l -degrees of the independent frequencies. For this, we used the method of Breger (2000) for calculating the pulsation constant Q of each frequency (see also Liakos et al. 2022) by using the advantage of the known absolute properties of the pulsating star (Sect. 5). Then, we compared the Q values with the models of Fitch (1981) in order to derive the most likely oscillation modes. Table 7 includes the results only for the independent frequencies detected in each system. In particular, we list: the increasing number of the frequency ( n ), its value ( f n), amplitude ( A ), phase ( Φ ), S / N, Q , and l degrees. Table C.1 contains the same information as Table 7, except the Q and the l -degrees, but for the combination frequencies.', '7.1. RR Lep': 'Two TESS time series are currently available for RR Lep, namely Sector 5 and 32, and they span over 25.6 d. The data set of Sector 5 has a time resolution of 2 minutes and that of Sector 32, 10 minutes. Furthermore, these data sets were acquired with a separation of 1.94 years. Therefore only the data of Sector 5 were used (app. 13800 data points) in the analysis to minimise introduction of alias frequencies (Breger 2000).  \nThe pulsating member of RR Lep was found to oscillate in more than 30 frequencies. However, only five, namely f 1, f 2, f 5, f 6, and f 7, are identified as independent. We detected many harmonics of the f orb ( = 1.0924 d -1 ). The other frequencies are combinations of others and many of them result from the rotational splitting of the independent frequencies by the f orb and its harmonics due to the orbital motion of the pulsator around the common centre of mass. The rotationally split modes would not be spaced exactly equally by the orbital frequency due to the influence of the Coriolis force on the pulsations (Ledoux 1951). According to the Q values of the independent frequencies, all of them are identified as non-radial pressure modes with their l -degrees to range between 1 and 3. The ratio P pul / P orb for these frequencies is lower than 0.046, which agrees with the conclusion of Zhang et al. (2013, i.e. P pul / P orb < 0 . 07) regarding the non-radial nature of these frequencies. On the other hand, we note that f 5 / f 1 ∼ f 5 / f 7 ∼ 0 . 73, f 5 / f 2 ∼ 0 . 74, and f 5 / f 6 ∼ 0 . 78. These frequency ratios are typical for the radial fundamental and the first overtone modes (Stellingwerf 1979; Oaster et al. 2006; Aerts et al. 2010). Therefore, f 5 could be a radial fundamental mode. However due to the very close values of the other independent frequencies, we cannot be certain which of them could be the first overtone. The periodogram of RR Lep, including some of the strongest detected frequencies, is shown in two plots with di ff erent scales in Fig. 20, and a representative Fourier fit on the data points is illustrated in Fig. 22. Detailed results for the frequency search are given in Tables 7 and C.1.  \nThe comparison of the results of Liakos & Niarchos (2013) with the present ones shows that the authors of that work had found as the most dominant frequency the + 1 alias of the present f 1. Hence, for reasons of completeness, we also re-analysed the ground-based B and V residual data of Liakos & Niarchos (2013) based on the LC + RV model presented herein (Sect. 4). The scope of this re-analysis was the determination of the oscillation modes using the amplitude ratio and the phase di ff erence of a given frequency value from two di ff erent photometric filters. The new analysis of these data yielded the same alias frequency. Therefore, we fixed it at the value found from the analysis of TESS data and we tried to detect other frequencies as well. Unfortunately, given that the sample of data points and the photometric accuracy of these data sets are rather inferior in comparison to that of TESS, in addition to the dominant frequency only a second one, which di ff ers between the filters, with S / N > 4 (Table 8) was detected. The mode identification for the dominant frequency was made using the FAMIAS software (Zima 2008) and the MADδ Sct models (Montalbán & Dupret 2007). Particularly, the closest model to the stellar and pulsational parameters of the primary of RR Lep was found for a star with a mass of 2.2 M ⊙ (i.e. the maximum value in the MAD models), a temperature of 8470 K, and log g = 3 . 94 cm s -2 . This model suggests an l -degree equal to 3, similarly to that found for f 1 of the TESS data (Table 7).  \nTable 7: Independent oscillation frequencies of the primary components of RR Lep and BF Vel based on the TESS data sets.  \nFig. 20: Periodogram of the pulsating component of RR Lep. The strongest frequencies are indicated.  \n<!-- image -->', '7.2. BF Vel': 'The frequency search of BF Vel was based only on the shortcadence (2-min time resolution) TESS data set of Sector 35. This data set spans over approximately 24 days, but there is a 6-day time gap. The total number of points used for the analysis, after the exclusion of the primary eclipses, is 10589. Unfortunately, our ground-based photometric data cover much shorter time spans (i.e. ∼ 2.5 d in B and ∼ 4 d in V and I filters) and have much lower photometric accuracy, therefore, they were excluded from the analysis. The TESS data set of Sector 62, although it has the same time resolution and covers about the same time span as that of Sector 35, it includes ∼ 30% fewer data points and has a time di ff erence of approximately two years with that of Sector 35. Therefore, as the combination of these data sets into one common frequency analysis is subject to alias, we excluded this data set from our analyses.  \nThe periodogram of the oscillating member of BF Vel (Fig. 21) is dominated by the orbital frequency and its harmon-  \nFig. 21: The same as Fig. 20 but for BF Vel.  \n<!-- image -->  \nTable 8: Oscillation frequencies of the primary component of RR Lep using the B and V ground-based data of Liakos & Niarchos (2013).  \nics. We detected most of the harmonics of f orb up to 38 f orb. Only one independent frequency was detected, that is f 3 = 46 . 73 d -1 , and also a few frequencies that indicate rotational splitting by the f orb and its harmonics. In total, 37 oscillation frequencies were detected, but we cannot be certain of the intrinsic origin of all of them, since many of them were found to be connected to an extremely low-frequency value (i.e. f 6) or to high-order multiples of the f orb. The Q value of f 3, based on the absolute parameters of the primary component and the same method as in the case of RR Lep, was found to be 0 . 014 ± 0 . 002 d. According to models of Fitch (1981), f 3 is either a p -mode or a radial mode and particularly the 4 th overtone of the radial fundamental mode. The pulsations model is given in Tables 7 and C.1, while a sample of the Fourier fitting is presented in Fig. 22.  \nFig. 22: Fourier fittings on a small fraction of the available TESS data for each system studied.  \n<!-- image -->', '8. Summary, discussion, and conclusions': "Multi-epoch, multi-filtered ground-based and high-accuracy space-borne photometry of the systems RR Lep and BF Vel have been combined with medium resolution spectroscopy to derive the physical properties of their components and, furthermore, to study their pulsational behaviour. Both systems are identified as oEA stars with their primary components displaying δ Sct type pulsations. Radial velocities of both components of the binary pairs were determined with broadening function analyses of the blue spectra. These analyses also indicated the presence of tertiary companions, which were confirmed with spectra of H α and Na I D lines. Orbital period modulations were analysed to determine probable orbital period modulating mechanisms. The secondary stars of both systems have evolved to the subgiant stage, and are chromospherically active.  \nAlthough for both systems a significant contribution of third light was detected in the photometric and spectroscopic observations and analyses, they lack clear periodic variations in their orbital periods that could be attributed to the existence of their tertiary companions. However, we conclude that the LITE is the most appropriate orbital period modulating mechanism.  \nThe primary star of RR Lep pulsates in five independent frequencies, with the dominant being 32.275 d -1 . All these frequencies are identified as p -modes, but f 5 could also be a radial fundamental mode. In addition, 31 combination frequencies were detected, many of which are connected to the orbital frequency of the system. This star is the most massive (2.9 M ⊙ ) among the δ Sct stars in the sample of the oEA systems. Our results regarding the pulsation modelling partially agree with those of Kahraman Aliçavu¸s et al. (2024). They found seven independent frequencies with the most dominant one to be 31.866 d -1 followed by the 32.276 d -1 as the second strongest one. We also detected these two frequencies ( f 1 and f 2) but in reverse order (i.e. the amplitude of f 1 was found larger than that of f 2). The frequencies 23.665 d -1 and 32.492 d -1 were detected in both studies ( f 5 and f 7 in our work). The main discrepancies between the two works concern the current f 6, which was not detected by Kahraman Aliçavu¸s et al. (2024), and their classification for the 25.782 d -1 , 26.449 d -1 , and 31.986 d -1 . The latter frequencies were also detected in our work (i.e. f 8, f 9, and f 21), but  \nwe found that they are just combinations of other frequencies. Larger discrepancies occur between our work and that of Liakos & Niarchos (2013), who detected only two frequencies, namely 33.280 d -1 and 24.318 d -1 . The first one in their study is probably the alias of the present f 1 (i.e. f 1 = 33 . 280 -1 d -1 ), which could be caused by the time gaps in the ground-based LCs that they used (alias frequencies, cf. Breger 2000). Although we reanalysed the data of Liakos & Niarchos (2013) by setting the dominant frequency f 1 to that found from the TESS data, we did not detect any other frequencies in the B and V ground-based data sets similar to those detected in the TESS data. However, for f 1, both methods of mode identification (i.e. models of Fitch (1981) and the MAD models) yielded the same l -degree value ( l = 3). The TESS data that we used for modelling is far more precise than that used in previous studies with ground-based data; however, the latter are useful for checking the main results, not for pulsation modelling.  \nOne independent frequency mode ( = 46.731 d -1 ) was found for the δ Sct star of BF Vel, and more than 35 other combination frequencies were detected. The majority of these frequencies are splittings of the orbital frequency due to rotation of the star around the systemic centre of mass, but many of them are probably leftovers from other stronger frequencies. Our results have large discrepancies compared with those of Manimanis et al. (2009). In particular, they detected seven frequencies with the dominant one being 44.94 d -1 , which is close to our f 3. The other frequencies reported in that work are either harmonics of others or they were not characterized at all. With the exception of their 7.48 d -1 , they did not detect any other slow frequencies, in contrast to our results. However, it should be noted that the results of Manimanis et al. (2009) were based on a much smaller data sample (i.e. 350 points in B filter spanning over three consecutive nights) in comparison with the TESS data that we used. The short time span of their data series could also be the reason they did not detect slow frequencies connected to the f orb.  \nIn order to compare the absolute and pulsation properties of the oscillating components of BF Vel and RR Lep with similar systems, we placed them among other δ Scuti stars of other oEA systems in two diagrams. The first diagram (Fig. 23) is the orbital-pulsation periods ( P orb -P puls) correlation and the second (Fig. 24) is the pulsation period-evolutionary stage ( P puls -log g )  \nFig. 23: P orb -P puls correlation for δ Scuti stars of oEA systems (star symbols) and the locations of the pulsating components of RRLep (magenta) and BF Vel (cyan). Black solid line represents the linear fitting of Liakos & Niarchos (2017).  \n<!-- image -->  \nFig. 24: P puls -log g correlation for δ Scuti stars that belong to SB1(2) + EB(el) (black dots) and EB (black squares) oEA systems and the locations of the pulsating components of RR Lep (magenta) and BF Vel (cyan). Black solid line represents the linear fitting of Liakos & Niarchos (2017).  \n<!-- image -->  \ncorrelation (cf. Soydugan et al. 2006a; Liakos et al. 2012; Liakos & Niarchos 2015, 2017). For both diagrams, the majority of the sample and the empirical linear fittings were taken from Liakos & Niarchos (2017). The sample of the P orb -P puls diagram includes 96 oEA systems, while that of P puls -log g includes 75 δ Scuti stars categorized according to the type of their host systems (i.e. SB1(2) + EB(el) or EB; for details see Sect. 5). The reason for the di ff erent number of stars in these samples is that the absolute parameters have not been determined for every oEA system. The pulsating components of BF Vel and RR Lep follow the empirical trends in these plots. The primary of BF Vel is among the ten fastest pulsators and the ten youngest stars of these samples, while its host system is among the five systems with the shortest orbital periods. The δ Scuti component of RR Lep is, similarly to that of BF Vel, a relatively young and fast pulsating star in comparison with the rest of these samples. These findings are in very good agreement with the main conclusions of Liakos (2017) regarding a) the connection of the dominant frequency with the age of the star (i.e. log g ; the younger the star the faster its pulsation frequency) and b) the connection between the orbital and pulsations periods (i.e. the shorter the systemic orbital period the faster the pulsation frequency of the oscillating component).  \nFurther investigations of similar systems are strongly encouraged. Expanding the study sample of oEA stars with precisely determined absolute properties will enhance our understanding of how interactions between components in close binary systems impact pulsations.  \nAcknowledgements. It is a pleasure to express our appreciation of the high quality and ready availability, via the Mikulski Archive for Space Telescopes (MAST), of data collected by the TESS mission. Funding for the TESS mission is provided by the NASA Explorer Program. This research partly made use of the SIMBAD database, operated at CDS, Strasbourg, France, and of NASA's Astrophysics Data System Bibliographic Services. A. L. acknowledges financial support from the NOA's internal fellowship 'SPECIES' (No. 5094). D. J. W. M. and J. F. W. acknowledge grants for time on the ANU 2.3 m telescope from the Edward Corbould Research Fund of the Astronomical Association of Queensland. D. J. W. M. is grateful to Prof. Michael Drinkwater for his extensive assistance with the planning and data analysis of the ANU 2.3 m telescope observations. D. J. W. M. thanks C. M. Moriarty, A. Mohit, Y. Rist, L. Coetzee, and S. Sweet for their assistance with the spectroscopic observations. We thank the anonymous reviewer for the valuable comments that improved the quality of this work.", 'Data Availability': 'The majority of data included in this article are available as listed in the paper or from the online supplementary material it refers to. The TESS data are available online from the MAST repository ( https://mast.stsci.edu/portal/ Mashup/Clients/Mast/Portal.html . All the times of minima used for the O-C analysis in this work will be made available at the CDS ( http://cdsweb.u-strasbg.fr/ ).', 'References': "Abhyankar, K. D. & Vyas, M. L. 1989, Space Sci. Rev., 50, 376  \nAerts, C., Christensen-Dalsgaard, J., & Kurtz, D. W. 2010, Asteroseismology (Springer Netherlands) Antoci, V., Cunha, M., Houdek, G., et al. 2014, ApJ, 796, 118 Baker, N. & Kippenhahn, R. 1962, ZAp, 54, 114 Balona, L. A. 2014, MNRAS, 439, 3453 Balona, L. A., Daszy'nska-Daszkiewicz, J., & Pamyatnykh, A. A. 2015, MNRAS, 452, 3073 Baran, A. S., Koen, C., & Pokrzywka, B. 2015, MNRAS, 448, L16 Bevington, P. R. 1969, Data reduction and error analysis for the physical sciences (McGraw-Hill Education) Bookmyer, B. B., Faulkner, D. R., & Samec, R. G. 1986, Information Bulletin on Variable Stars, 2873, 1 Bowman, D. M., Johnston, C., Tkachenko, A., et al. 2019, ApJ, 883, L26 Bowman, D. M. & Kurtz, D. W. 2018, MNRAS, 476, 3169 Bowman, D. M. & Michielsen, M. 2021, A&A, 656, A158 Bradstreet, D. H. & Steelman, D. P. 2002, in American Astronomical Society Meeting Abstracts, Vol. 201, American Astronomical Society Meeting Abstracts, 75.02 Breger, M. 2000, in Astronomical Society of the Pacific Conference Series, Vol. 210, Delta Scuti and Related Stars, ed. M. Breger & M. Montgomery, 3 Childress, M. J., Vogt, F. P. A., Nielsen, J., & Sharp, R. G. 2014, Ap&SS, 349, 617 Claret, A. 2017, A&A, 600, A30 Claret, A. & Bloemen, S. 2011, A&A, 529, A75 Coelho, P., Barbuy, B., Meléndez, J., Schiavon, R. P., & Castilho, B. V. 2005, A&A, 443, 735 Collins, K. A., Kielkopf, J. F., Stassun, K. G., & Hessman, F. V. 2017, AJ, 153, 77 Creevey, O. L., Sordo, R., Pailler, F., et al. 2023, A&A, 674, A26 De Loore, C. W. H. & Doom, C. 1992, Structure and evolution of single and binary stars, Vol. 179 (Springer Science + Business Media Dordricht) Dopita, M., Hart, J., McGregor, P., et al. 2007, Ap&SS, 310, 255 Dopita, M., Rhee, J., Farage, C., et al. 2010, Ap&SS, 327, 245 Dvorak, S. 2009, Communications in Asteroseismology, 160, 64 Eddington, A. S. 1926, The Internal Constitution of the Stars (Cambridge: Cambridge University Press) Eker, Z., Bakı¸s, V., Bilir, S., et al. 2018, MNRAS, 479, 5491 Erdem, A., Budding, E., Demircan, O., et al. 2005, Astronomische Nachrichten, 326, 332 Erdem, A. & Öztürk, O. 2016, New A, 48, 33 Fabricius, C., Høg, E., Makarov, V. V., et al. 2002, A&A, 384, 180 Fitch, W. S. 1981, ApJ, 249, 218  \n- Gaia Collaboration, Klioner, S. A., Lindegren, L., et al. 2022, A&A, 667, A148 Gray, R. O. & Corbally, C. J. 2014, AJ, 147, 80\n- Gray, R. O. & Corbally, Christopher, J. 2009, Stellar Spectral Classification (Princeton University Press)  \nHilditch, R. W. 2001, An Introduction to Close Binary Stars (Cambridge University Press)  \nHøg, E., Fabricius, C., Makarov, V. V., et al. 2000, A&A, 355, L27 Irwin, J. B. 1959, AJ, 64, 149  \nJenkins, J. M., Twicken, J. D., McCauli ff , S., et al. 2016, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 9913, Software and Cyberinfrastructure for Astronomy IV, ed. G. Chiozzi & J. C. Guzman, 99133E  \nKahraman Aliçavu¸s, F., Soydugan, E., Smalley, B., & Kubát, J. 2017, MNRAS, 470, 915  \nKahraman Aliçavu¸s, F., Aliçavu¸s, F., Çoban, Ç. G., Handler, G., & De Cat, P. 2024, MNRAS, 527, 4076  \nKallrath, J. & Milone, E. F. 1999, Eclipsing binary stars : modeling and analysis (Springer-Verlag New York)  \nKim, S.-L., Lee, J. W., Lee, C.-U., et al. 2021, AJ, 162, 212  \nKopal, Z. 1959, Close binary systems (John Wiley & Sons)  \nKreiner, J. M. 2004, Acta Astron., 54, 207  \nKurtz, D. W., Handler, G., Rappaport, S. A., et al. 2020, MNRAS, 494, 5118 Ledoux, P. 1951, ApJ, 114, 373  \nLenz, P. & Breger, M. 2005, Communications in Asteroseismology, 146, 53 Liakos, A. 2015, in Astronomical Society of the Pacific Conference Series, Vol. 496, Living Together: Planets, Host Stars and Binaries, ed. S. M. Rucinski, G. Torres, & M. Zejda, 286  \nLiakos, A. 2017, A&A, 607, A85  \nLiakos, A., Moriarty, D. J. W., Blackford, M. G., et al. 2022, A&A, 663, A137  \nLiakos, A. & Niarchos, P. 2012, New A, 17, 634  \nLiakos, A. & Niarchos, P. 2013, Ap&SS, 343, 123  \nLiakos, A. & Niarchos, P. 2015, in Astronomical Society of the Pacific Conference Series, Vol. 496, Living Together: Planets, Host Stars and Binaries, ed. S. M. Rucinski, G. Torres, & M. Zejda, 195  \nLiakos, A. & Niarchos, P. 2017, MNRAS, 465, 1181  \nLiakos, A., Niarchos, P., Soydugan, E., & Zasche, P. 2012, MNRAS, 422, 1250  \nLoumos, G. L. & Deeming, T. J. 1978, Ap&SS, 56, 285  \nLucy, L. B. 1967, ZAp, 65, 89  \nManimanis, V. N., Vamvatira-Nakou, C., & Niarchos, P. G. 2009, Ap&SS, 323, 115  \nMkrtichian, D., Gunsriviwat, K., Lehmann, H., et al. 2022, Galaxies, 10, 97  \nMkrtichian, D. E., Kusakin, A. V., Gamarova, A. Y., & Nazarenko, V. 2002, in  \nAstronomical Society of the Pacific Conference Series, Vol. 259, IAU Col-  \nloq. 185: Radial and Nonradial Pulsationsn as Probes of Stellar Physics, ed.  \nC. Aerts, T. R. Bedding, & J. Christensen-Dalsgaard, 96  \nMontalbán, J. & Dupret, M. A. 2007, A&A, 470, 991  \nMoriarty, D. J. W., Liakos, A., Drinkwater, M. J., et al. 2019, Ap&SS, 365, 1  \nOaster, L., Smith, H. A., & Kinemuchi, K. 2006, PASP, 118, 405  \nPecaut, M. J. & Mamajek, E. E. 2013, ApJS, 208, 9  \nPilecki, B., Gieren, W., Smolec, R., et al. 2017, ApJ, 842, 110  \nPrša, A. 2018, Modeling and Analysis of Eclipsing Binary Stars; The theory and design principles of PHOEBE (IOP Publishing Ltd)  \nRicker, G. R., Winn, J. N., Vanderspek, R., et al. 2015, Journal of Astronomical Telescopes, Instruments, and Systems, 1, 014003  \nRuci'nski, S. M. 1969, Acta Astron., 19, 245  \nRuci'nski, S. M. 1992, AJ, 104, 1968  \nRuci'nski, S. M. 2002, AJ, 124, 1746  \nSamec, R. G., Fuller, R. E., Bookmyer, B. B., & Faulkner, D. R. 1989, PASP, 101, 180  \nSoydugan, E., ˙ Ibanoˇglu, C., Soydugan, F., Akan, M. C., & Demircan, O. 2006a, MNRAS, 366, 1289  \nSoydugan, E., Soydugan, F., ¸Senyüz, T., Püsküllü, Ç., & Demircan, O. 2011, New A, 16, 72  \nSoydugan, E., Soydugan, F., Demircan, O., & ˙ Ibanoˇglu, C. 2006b, MNRAS, 370, 2013  \nStellingwerf, R. F. 1979, ApJ, 227, 935  \nTkachenko, A., Lehmann, H., & Mkrtichian, D. E. 2009, A&A, 504, 991  \nvon Zeipel, H. 1924, MNRAS, 84, 665  \nVyas, M. L. & Abhyankar, K. D. 1989, A&AS, 81, 81  \nWilson, R. E. & Devinney, E. J. 1971, ApJ, 166, 605  \nZasche, P., Liakos, A., Niarchos, P., et al. 2009, New A, 14, 121  \nZhang, X.-B., Luo, C.-Q., & Fu, J.-N. 2013, ApJ, 777, 77  \nZhevakin, S. A. 1963, ARA&A, 1, 367  \nZima, W. 2008, Communications in Asteroseismology, 157, 387  \nZola, S., Rucinski, S. M., Baran, A., et al. 2004, Acta Astron., 54, 299", 'Appendix A: Radial velocities': 'Table A.1 contains the RV measurements derived from our spectroscopic observations. The values in the O -C columns represent the deviations of the individual measurements from the fitted RV curves. Orbital phases were calculated using the ephemeris given in Tables 2 and 3.  \nTable A.1: Radial velocity measurements of the components of RR Lep and BF Vel.', 'Appendix B: Times of minima': 'Table B.1 includes the ToM used for the ETV analyses of RR Lep and BF Vel (Sect. 6). Particularly, we list: the HJD of the ToM, the type of observation ( Obs . type ) based on which the respective data for its calculation were obtained (i.e. vis = visual, pg = photographic, pe = photoelectric, CCD; TESS; ASAS = CCD), and the literature source (e.g. name(s) of observer(s), campaign, mission).  \nTable B.1: Times of minima (ToM) of RR Lep and BF Vel used in the ETV analyses.  \nTable B.1 (continued)', 'Appendix C: Combination pulsation frequencies': 'Table C.1 contains the increasing number of the frequency ( n ), its value ( f n), amplitude ( A ), phase ( Φ ), S / N, and the most possible frequency combination ( Combo ). This table is complementary to Table 7 and their integration is the complete model of the pulsational behaviour of each studied system (Sect. 7).', 'T able C.1: Combination pulsation frequencies of the primary components of RR Lep and BF V el.': 'y  \nfrequenc  \nintrinsic  \nan  \nnot  \nand  \ner  \nv  \nlefto  \npossible  \n<!-- image -->'}