Abstract
Introduction:
The first multi-color photometric observations and low-resolution spectra of the Algol-type short-period eclipsing binary systems IWUMa and V367 Gem are presented here. The stellar atmospheric parameters of the primary stars were derived through spectral analysis. Ground-based and TESS-based light curves were analyzed.
Methods:
ULySS software was employed to obtain atmospheric parameters by fitting full spectra with model spectra. Wilson–Devinney code was used to analyze the light curves.
Results:
The photometric solutions suggest that IW UMa and V367 Gem are both typical Algol systems in a semidetached configuration, with secondary stars filling their Roche lobes. The mass ratios of IW UMa and V367 Gem are 0.145 (3) and 0.478 (15), respectively. A third light was detected in the light curve analysis of V367 Gem with a luminosity contribution of . Using OC analysis, we found that the period of IW UMa is likely increasing over the long term, superposed with cyclic variation. The period of V367 Gem is decreasing at a rate of on a long time scale.
Discussion:
The period decrease of V367 Gem may be due to the mass and angular momentum loss from the system. Additionally, a probable cyclic variation has been identified for the period of V367 Gem, which may be caused by the third body.
1 Introduction
Algol-type binaries are a unique subset of semi-detached eclipsing binaries (EB systems) characterized by mass transfer between their components, which significantly influences their evolution. In these systems, the more luminous and massive component, referred to as the “primary”, is typically composed of B-, A-, or F-type stars. The less luminous and less massive “secondary” component may be an evolved giant or subgiant star. In these systems, the primary component is a main sequence star, while the secondary component fills its Roche lobe, leading to mass transfer from the secondary to the primary. These systems began as binary systems with unequal masses. The more massive component (the future secondary) evolved from the main sequence more quickly than the less massive component (the future primary). When the former primary fills its Roche lobe, a significant mass-transfer episode occurs during which a substantial amount of mass from the former primary is transferred to the less massive component, causing it to become more massive than the Roche-lobe-filling star (the former more massive component). Therefore, current Algols are systems undergoing slow mass transfer from a Roche-lobe-filling giant or subgiant low-mass star (the former massive star) to a main sequence, massive, mass-accreting component (the current primary). The mass transfer results in an increase in the orbital period. Secondary components with later spectral types are magnetically active, producing features such as spots, plages, flares, and winds (e.g., Hall, 1989). These systems lose angular momentum due to magnetic stellar winds, suggesting that mass transfer between the components should be non-conservative (e.g., Sarna et al., 1997; 1998). The observed period decrease in the orbital periods may be attributed to mass loss and transfer and angular momentum loss (AML) (Qian, 2000; 2001; Yang and Wei, 2009; Erdem et al., 2010; Soydugan et al., 2011; Ibanoǧlu et al., 2006). In addition, a magnetic braking mechanism is often referenced to explain the observed orbital shrinkage. Chen et al. (2006) suggested that, in addition to magnetic braking, a circumbinary disk may play an important role in AML from Algol-type binaries, which could explain the low-mass ratio systems undergoing rapid mass transfer.
IW UMa (ASASSN-V J064643.50 + 205322.1) is a confirmed Algol-type variable binary with a short orbital period of 0.7747182 (1) d (Kreiner, 2004). IW UMa was identified as eclipsing Am binary candidates by confirming the relationship between different spectral types, based on the low-resolution spectra of the Large Sky Area MultiObject Fiber Spectroscopic Telescope (LAMOST) (Tian et al., 2023). The MK spectral type of primary component is kA5hA5mA9 (Tian et al., 2023). Borovicka and Sarounova (1996) first released the VR band light curves of IW UMa.
V367 Gem (SVS 863, ASASSN-V J092900.54 + 434402.7) is an Algol-type variable binary with a short orbital period of 0.699229 (1) d as revised by Hintz and Brown (2007). V367 Gem was also identified as eclipsing Am binary candidates with an MK spectral type of the primary component classified as kA3hA5mA7 by Tian et al. (2023). Hintz and Brown (2007) also released the light curves and found some minimal color changes with values of (BV) = 0.270.01 and (VR) = 0.180.01. Based on the measured colors and the shape of the light curve, V367 Gem was classified as an Algol-type system by Hintz and Brown (2007). Unfortunately, the two Algol-type eclipsing binaries have not been investigated in more detail.
In this study, we obtain the multi-color high-precision light curves and low-resolution spectra of both targets. Based on the times of the light minimum collected and calculated from the TESS database, the OC diagrams with long time-spanning eclipse timings were obtained for the short-period Algol-type eclipsing binaries IW UMa and V367 Gem. We thus obtained more accurate geometric and physical parameters.
2 Observation and data reduction
2.1 Photometric observations
The new complete multi-color light curves of IW UMa in JohnsonCousin filters (Ažusienis and Straižys, 1969; Bessell, 1983) were observed on 30 January and 24–26 March, 2023 and on 7 February, 2024, with a PI VersArray 1300B LN CCD (size: pixel) photometric system attached to the 80 telescope (Tsinghua-NAOC Telescope) at the Xinglong Station of the National Astronomical Observatories of China (NAOC). New complete JohnsonCousin (Ažusienis and Straižys, 1969; Bessell, 1983) multi-color light curves of V367 Gem were observed on 11–13 January, 2024. The observations were also made with the 80 telescope (Tsinghua-NAOC Telescope) at NAOC. The PHOT package of IRAF was used to process all images. The comparison and check stars were chosen to determine the differential magnitudes. For IW UMa, the non-variable star TYC 2997-1178-1 was selected as the comparison, and TYC 2997-1377-1 and TYC 2997-1501-1 were the check stars. GSC 01342-01119 and UCAC4 555-033241 were selected as comparison stars for V367 Gem, and TYC 1342-417-1 was the check star. Table 1 shows information for the comparison and check stars and for both targets. The comparison stars have a similar brightness to the targets. The observed CCD images of both targets are shown in Figures 1 and 2; the variable star is marked “V”, the comparison star is marked “C”, and the check star is marked “Ch”. The two check stars for IW UMa are marked “Ch1” and “Ch2” in Figure 1, and the two comparison stars for V367 Gem are marked “C1” and “C2” in Figure 2. The light curves of IW UMa and V367 Gem are shown in Figures 3 and 4, respectively. The data observed on different nights are marked with different colors, and the B, V, and band data are represented by a square, circle, triangle, and inverted triangle, respectively. The magnitude difference of the comparison and check stars in all filters are shown at the bottom of the figures. In Figure 3, the mean standard deviations of the magnitude difference in each filter are approximately 0.00678 in the B band, 0.00637 in the V band, 0.00589 in the band, and 0.00527 in the band. In Figure 4, the mean standard deviations of the magnitude difference of each filter are 0.01075 in the B band, approximately 0.01175 in the V band, 0.00695 in the band, and 0.00818 in the band.
TABLE 1
| Targets | Name | Mag | ||
|---|---|---|---|---|
| Variable (V) | IW UMa | V = 12.20 | ||
| Comparison (C) | TYC 2997-1178-1 | V = 12.14 | ||
| Check (Ch1) | TYC 2997-1377-1 | V = 10.95 | ||
| Check (Ch2) | TYC 2997-1501-1 | V = 10.67 | ||
| Variable (V) | V367 Gem | V = 11.35 | ||
| Comparison (1) | GSC 01342-01119 | V = 12.7 | ||
| Comparison (2) | UCAC4 555-033241 | V = 13.80 | ||
| Check (Ch) | TYC 1342-417-1 | V = 11.66 |
Information of IW UMa, V367 Gem, and the corresponding comparison and check stars.
FIGURE 1
FIGURE 2
FIGURE 3
FIGURE 4
2.2 TESS observations
The Transiting Exoplanet Survey Satellite (TESS) is an all-sky transit survey whose principal goal is to detect Earth-sized planets orbiting bright stars that are amenable to follow-up observations to determine planet masses and atmospheric compositions (Ricker et al., 2014). These targets are read out as postage stamps and made available to the community as target pixel files (TPFs) and corrected light curves. These files comprise several time series, including a simple aperture photometry, a corrected light curve, position vectors, and quality flags. The aperture photometry flux series is termed “SAP_FLUX”, while the flux series that has the common instrumental systematics removed using cotrending basis vectors (CBV) files is termed “PDCSAP_FLUX”. SAP_Flux was used to obtain the light curves of IW UMa and V367 Gem in this study.
IW UMa (TESS target identifier: TIC 105439439) and V367 Gem (TESS target identifier: TIC 56897408) have been observed by TESS for many days in different sectors. The data were retrieved using the Python package (Lightkurve Collaboration et al., 2018) from the MAST data archive (http://https://mast.stsci.edu/), which provides three different types of data products. Two kinds of light curves from the authors TESS-SPOC and QLP, and the other is the TESScut cutouts of the calibrated full frame images from the author TESScut. Using the TESScut tool, we downloaded pixel files centered on IW UMa and V367 Gem. Then, the Python procedure dat.to_lightcurve was used to obtain a simple aperture photometry light curve. The flux data are normalized. We thus obtained the light curves of IW UMa from sector 48 with a cadence of 600 s and that of V367 Gem from sectors 71 and 72 with a cadence of 200s. The observation time of the light curves of IW UMa from sector 48 lasted from 29 January to 11 February and from 14 to 25 February 2023. The observation time of light curves of V367 Gem from sectors 71 and 72 lasted 19–28 October, 1–10 and 15–24 November, and 28 November to 7 December 2023. The time-series light curves of each sector and phase-binned light curve of IW UMa and V367 Gem obtained from the TESS database are presented in Figures 5 and 6. The phase-binned light curves of both targets were based on the part of the data in corresponding sectors.
FIGURE 5
FIGURE 6
2.3 Spectroscopic observations
The spectra of IW UMa and V367 Gem were observed using the Beijing Faint Object Spectrograph and Camera (BFOSC), which was mounted on the 2.16 m telescope of the Xinglong station of the National Astronomical Observatories of China (NAOC), Chinese Academy of Sciences (CAS). The observations were made on 14 January 2023, using a low-dispersion spectrometer BFOSC and grism G6. The slit width and line dispersion of grism G6 are 1.8 arcsec and 88 Å , respectively. The observations cover the wavelength range 330–545 nm (Fan et al., 2016). The spectra resolution per pixel is 1.98 Å. IRAF was used to process the observation images and extract the spectra. Normalized flux was obtained, and the atmospheric absorption lines were also corrected. Low-resolution spectra can only show the spectral lines of the primary star, and the lines of the secondary star may be too faint to be detected. The observed spectra of both targets are shown in Figure 7 with black lines.
FIGURE 7
University of Lyon Spectroscopic Analysis Software (ULySS) (Koleva et al., 2009) was employed to obtain atmospheric parameters by fitting full spectra with model spectra generated by an interpolator with the ELODIE library (Prugniel and Soubiran, 2001).
The red lines in Figure 8 represent the fitting spectra. The atmospheric stellar parameters were derived—for IW UMa: = 737870 K, = 4.180.04 cm/s2, [Fe/H] = −0.290.04 dex; for V367 Gem: = 829052 K, = 4.280.04 cm/s2, [Fe/H] = 0.040.05 dex. For the two targets, the spectra were analyzed as single-star spectra, ignoring the secondary (faint) component in the spectra analysis program, the spectra of the primary component, which is casually diluted by the secondary component, will mimic a low metallicity. The atmospheric parameters obtained by ULySS can be used to describe the atmospheric characteristics of the primary star.
FIGURE 8
The physical properties of the primary can be estimated from the atmosphere parameters by matching with stellar isochrones for normal single stars, assuming that the mass-accreting component is in thermal equilibrium (Zhang et al., 2019), using the PARSEC (Padova and Trieste Stellar Evolution Code) isochrone database (Bressan et al., 2012; Chen et al., 2014; 2015; Tang et al., 2014) along with the lognormal form initial mass function by Chabrier (2001). The primary stars’ physical properties can then be determined, including mass and radius—for IW UMa: , , ; for V367 Gem: , , .
The blue-violet spectrum (3800–4600 Å) of IW UMa, V367 Gem, and part of the spectrum of two MK standard spectra with a type of A3 V (Bet Leo) and F0 III are shown in Figure 8. It is clear that the Ca II K lines of IW UMa and V367 Gem are quite weaker in strength than that of the A3 and F0 MK standards, which exhibit typical characteristics similar to Am stars. The low resolution of the spectral data means that many metallic properties are not adequately captured. Indeed, the Am classification of the primary components of Algol undergoing mass accretion and experiencing high atmospheric turbulence dynamics contradicts the theories of Am star formation. The accurate determination of chemical composition should be based on multi-element chemical composition analysis conducted through high-resolution spectroscopy.
3 OC diagram analyses
All available light minimum times of IW UMa and V367 Gem were collected from the literature. Meanwhile, we computed the light minima based on the TESS data for both targets in HJD time—listed in the Appendix of Supplementary Tables A1 and B1 marked “This study*” in the corresponding reference column.
3.1 IW UMa
A total of 95 eclipse times of IW UMa spanning approximately 28 years were used to structure the OC diagram, including two eclipse times observed by us and 62 eclipse times obtained from the TESS data. All eclipse times used for the period analysis are listed in the Supplementary Table A1 Appendix. For the reference items marked “*”, the original publications unavailable online. The corresponding OC were obtained with the following ephemeris:In Equation 2, “Min I” are the observed primary eclipse times and 2460030.1425 (2) is the primary eclipse time of the multi-color light curves.
The OC diagram of IW UMa is shown in Figure 9. In particular, the open square was used to mark the light minimum obtained from the TESS database. The same weight of 1 was taken for all CCD data, including the TESS data in the OC analysis. Using the least squares method, Equation 2 can be used to describe the OC diagram:
FIGURE 9
The error associated with the coefficient of the quadratic term is quite large. The black line in Figure 9 represents the fitting curve. We noticed that this parabolic fit is statistically insignificant. The fitting residuals are shown in the middle panel of Figure 9, which may indicate a cyclic variation. The final residuals are shown at the bottom of Figure 9. The period of IW UMa is likely increasing over the long term superposed with cyclic variation. Considering the large fitting error, this conclusion requires further confirmation.
3.2 V367 Gem
For V367 Gem, 117 eclipse times crossing approximately 25 years were used to calculate the OC diagram, including two new observed and 103 eclipse timings estimated from the TESS data by this study (Supplementary Table B1 Appendix). Equation 3 is used to obtained the corresponding OC and E.
The OC diagram with error bars of V367 Gem based on all the eclipsing times is shown in the topleft panel of Figure 10; the light minimum obtained from the TESS database is marked with green color. In the OC analysis program, the same weight of 1 was adopted for all the CCD data. Using the least square method, the OC of V367 Gem can be described by Equation 4.
FIGURE 10
The quadratic term in the above polymerization indicates that the period is undergoing a long-term decrease at a rate of . The black line in Figure 10 represents the corresponding fitting curve. The fitting residuals are shown in the left-bottom panel of Figure 10. We found that the OC diagram constructed by the light eclipse timings obtained from TESS data may exhibit an incomplete cycle of a sinusoidal variation with a short period. The cyclic variation can be fitted by Equation 5:in which: is the projected semi-major axis and is the velocity of light; , , and are the eccentricity, true anomaly, and longitude of the periastron passage, respectively; is the eccentric anomaly.
The OC diagram based only on the mid-times of eclipses obtained from TESS are shown in the top-right panel in Figure 10, and the residuals are shown in the right-bottom panel. The incomplete cyclic variation fits the OC diagram only based on TESS data well, indicating a cyclic variation with a short period of 0.26 (37) years and an amplitude of A0.002 (24) days. The OC diagram analysis demonstrates that the orbital period of V367 Gem is undergoing a continuous decrease and a probable cyclic variation simultaneously. The cyclic variation shown in the O-C of TESS mid-times requires further confirmation.
4 Light curve analyses
The Wilson–Devinney (WD) program (Wilson and Devinney, 1971; Wilson, 1990; 2012) was used to analyze all the light curves (including these based on the TESS data) of the two targets. The gravity-darkening coefficients and bolometric albedos were taken according to temperature (Ruciński, 1969). The bandpass limb-darkening coefficients were taken from van Hamme (1993), and the logarithmic bolometric coefficients were applied. There are various modes for different binary configurations in the WD program. We acquired the converged solutions of both targets with Mode 5 (the semidetached case with the secondary component filling the critical RL). The corresponding adjustable parameters of Mode 5 are: the monochromatic luminosity of star 1, , , , ; the orbital inclination i; the mean temperature of star 2, ; the dimensionless potential of the primary star, . We also set the third light as one adjustable parameter. In order to confirm a reliable mass ratio, the photometric solutions based on a group of assumed mass ratio values from 0.01 to 1 were obtained with the differential correction program. The search-step length was set to 0.01. The value of q resulted in the lowest sum of weighted squared deviations was chosen as the initial value for q in the WD modeling. Thence, q was treated as an adjustable parameter to achieve the final converged solution.
4.1 IW UMa
The first multi-color light curves of IW UMa were observed. The TESS observation data from sector 48 (covering HJD 2459625 to HJD 2459636) were utilized to generate the light curve for analysis. The light curves are shown at the bottom of Figure 5. In the figures below, the observed multi-color light curve is labeled “LC_”, while the light curve obtained from the TESS database is labeled “LC_TESS”. Different linear ephemerides were used to calculate the phased light curves. For the multi-color light curves, the ephemeris given by Equation 1 was employed. For the TESS-based light curve, the linear ephemeris Equation 6 was used: in which the epoch 2460030.14246 (12) is the new observed eclipse time and 2459625.73832 (168) is the eclipse time based on the TESS database.
The effective temperature of the primary star obtained as = 7378 K from our spectral analysis. The effective temperature of the primary star was fixed at 7378 K in the WD modeling. The gravity-darkening coefficients and bolometric albedos are taken as = 1, = 0.32, = 1, and = 0.5 (Ruciński, 1969). The third light was set as an adjustable parameter, but no converged result was obtained. The process of q search was confirmed. For all assumed values of q, the sums of weighted square deviations based on two sets of light curves are displayed in Figure 11. The minimal values achieved at q = 0.15 and q = 0.20 were based on the observed multi-color light curves and the TESS-based light curves, respectively. Then, q was made as an adjustable parameter, and q = 0.15 and q = 0.20 were taken as the initial value in the WD modeling for the two sets of light curves, respectively. The theoretical light curves for the multi-color light curves and TESS-based light curve are shown in Figures 12, 13 in red lines, respectively. The theoretical light curves fitted the light curves very well. All converged photometric solutions are listed in Table 2. The geometric structure is shown in Figure 14.
FIGURE 11
FIGURE 12
FIGURE 13
TABLE 2
| IW UMa | V367 Gem | |||
|---|---|---|---|---|
| LC_ | LC_TESS | LC_ | LC_TESS | |
| Parameters | ||||
| Mode | 5 | 5 | 5 | 5 |
| 1.0 | 1.0 | 1.0 | 1.0 | |
| 0.32 | 0.32 | 0.32 | 0.32 | |
| 1.0 | 1.0 | 1.0 | 1.0 | |
| 0.5 | 0.5 | 0.5 | 0.5 | |
| q () | 0.1445 (34) | 0.2441 (13) | 0.478 (15) | 0.460 (10) |
| 7378 | 7378 | 8280 | 8280 | |
| 4671 (26) | 4570 (4) | 4579 (32) | 4497 (15) | |
| 80.66 (31) | 74.681 (38) | 80.45 (32) | 80.10 (30) | |
| 2.872 (11) | 3.2616 (58) | 3.601 (40) | 3.371 (22) | |
| 2.0882 | 2.3392 | 2.833 | 2.798 | |
| 0.98019 (2) | — | 0.9764 (5) | — | |
| 0.96153 (6) | — | 0.9470 (11) | — | |
| 0.94494 (10) | — | 0.9152 (19) | — | |
| 0.92734 (16) | — | 0.8785 (32) | — | |
| — | 0.89913 (8) | — | 0.9093 (15) | |
| — | — | 0.3779 (87) | — | |
| — | — | 0.4269 (80) | — | |
| — | — | 0.4664 (81) | — | |
| — | — | 0.5040 (91) | — | |
| — | — | — | 0.9137 (13) | |
| 0.3653 (12) | 0.32998 (52) | 0.3172 (26) | 0.3401 (16) | |
| 0.3764 (13) | 0.33783 (56) | 0.3252 (28) | 0.3506 (18) | |
| 0.3801 (13) | 0.34199 (57) | 0.3325 (29) | 0.3602 (18) | |
| 0.2120 (14) | 0.24665 (38) | 0.2971 (25) | 0.2940 (17) | |
| 0.2204 (15) | 0.25666 (40) | 0.3100 (26) | 0.3068 (18) | |
| 0.2522 (16) | 0.28925 (41) | 0.3426 (26) | 0.3393 (18) | |
| 0.6110 (24) | 0.78662 (95) | 0.9725 (62) | 0.8939 (37) | |
| 0.3178 (38) | 0.2994 (12) | 0.3954 (93) | 0.4833 (75) | |
| 0.0093 | 0.0092 | 0.0077 | 0.0074 | |
Photometric solutions of IW UMa and V367 Gem.
FIGURE 14
4.2 V367 Gem
The first multi-color light curves of V367 Gem were observed. Different linear ephemeris was used to calculate the phased light curves for the observed multi-color light curves and TESS observation. For multi-color light curves, the ephemeris given by Equation 3 was used. For the TESS-based light curves, the linear ephemeris Equation 7 was used:in which the epoch 2460322.22210 (25) is the new observed eclipse times and 2460264.18578 (8) is the eclipse times based on the TESS data.
In the WD modeling, the effective temperature of the primary star was fixed at 8280 K, which was obtained from our spectral analysis. The convergent solutions were confirmed with Mode 5. The third light was set as an adjustable parameter, and a converged solution was obtained. The mass ratio search was also confirmed for V367 Gem based on the two sets of light curves. The sums of the weighted square deviations for all assumed values of q are displayed in Figure 15 based on the two sets of light curves. In order to make the minimal q value more obvious, only the trend of change from q = 0.3 to q = 1 was plotted. The minimal value achieved q = 0.49 and q = 0.46 for observed multi-color light curves and TESS light curves, respectively. We thus make “q” an adjustable parameter and set q = 0.49 and q = 0.46 as the initial value in the WD modeling for the two sets of light curves, respectively. The theoretical light curves with and without a third light for multi-color light curves are shown in Figure 16 in red and green lines, respectively. The residual of the theoretical light curves with a third light is shown at the bottom of Figure 16, where the theoretical light curves with a third light fit the light curves better. For LC_TESS, the converged solution with a third light was obtained. The theoretical light curves with a third light for the TESS light curve are shown in Figure 17 in red lines. All of the converged photometric solutions with a third light are listed in Table 2. The geometric structure of V367 Gem is plotted in Figure 18.
FIGURE 15
FIGURE 16
FIGURE 17
FIGURE 18
5 Discussion and conclusion
The atmospheric parameters of the primary components of the Algol-like short-period eclipsing Am binaries IW UMa and V367 Gem are obtained based on the spectra analyses: for IW UMa, = 737870 K, = 4.180.04 cm/s2, [Fe/H] = -0.290.04 dex; for V367 Gem, = 829052 K, = 4.280.04 cm/s2, [Fe/H] = 0.040.05 dex. The physical properties of the primary can be estimated from the atmospheric parameters by matching with stellar isochrones, using the PARSEC database. For IW UMa: , , ; for V367 Gem: , , . Light curve analysis of both targets used the WD program. The light curves of both targets show Algol-like characteristics. The theoretical light curves with mode 5 fitted the light curves of both targets quite well. The fill-degree factors of the primary components of IW UMa and V367 Gem are 0.3178 (38) and 0.3426 (26) based on the multi-color light curve solutions. The mass ratio of IW UMa is estimated as 0.145 (3), which is quite low. The mass ratio of V367 Gem is 0.478 (15). We detected a third body in the light curve analysis of V367 Gem. The luminosity contribution of the third light is about based on the observed multi-color light curves. The value is higher based on the TESS light curve because of the high background light interference of TESScut files, so such a value is not informative. The orbital inclinations of two targets are for IW UMa and for V367 Gem, which means that the photometric solutions can be reliable. Based on the physical properties of the primary and photometric solutions, the absolute parameters of the two binaries are estimated and are listed in Table 3.
TABLE 3
| IW UMa | V367 Gem | |
|---|---|---|
| Parameters | ||
Absolute parameters of IW UMa and V367 Gem.
All available eclipse times were used to investigate the variations of the orbital period by analyzing the OC diagram. We have attempted a parabolic fit superposed with a cyclic variation to describe the OC diagram of IW UMa but failed to obtain reliable conclusions due to the dispersion of data. It appears that the period of IW UMa is likely increasing over the long term. The period of V367 Gem is decreasing at a rate of on a long time scale. Additionally, the OC diagram derived from the TESS data may show an incomplete cycle of a sinusoidal variation. A cyclic variation was identified with a short period of 0.26 (37) years and an amplitude of A0.002 (24) days. Long-term monitoring is essential to confirm the existence of this cyclic variation. The long-term period increase or decrease change occurs commonly in systems with semi-detached configurations.
IW UMa and V367 Gem are typical Algol systems in a semidetached phase with the secondary star filling the RL. The period increase of IW UMa may be the cause of the mass transfer from the secondary to the primary star (e.g., Qian, 2002; Tian and Zhu, 2019). The period decrease of V367 Gem may mainly be caused by the mass and angular momentum loss (AML) via an enhanced stellar wind of the more evolved secondary star (e.g., Tout and Eggleton, 1988; Liao et al., 2017; 2019). The probable cyclic variation may be caused by the magnetic activity of one or both components (Applegate, 1992) and the light travel time effect (LTTE) through the presence of a tertiary companion. In light curve analysis, the detection of a third light may indicate the existence of a distant companion, which could plausibly explain the cyclic variation observed in V367 Gem due to LTTE. The third body plays an important role in the evolution of the binaries; it will cause the cyclic variation of an orbital period and also will take away the angular momentum of the system, leading to the orbital contraction process and then promoting the evolution of contact binaries (Tokovinin et al., 2006; Qian et al., 2013a; 2014b; 2017). Both targets are in a semidetached configuration, with the secondary star filling the RL, and then will evolve into a marginal-contact phase with poor thermal contact (e.g., Liao et al., 2012; Qian et al., 2013b; 2014a; Liao and Sarotsakulchai 2019), such as GW Gem and V723 Per (Tian and Zhu, 2019).
The components in very-short-period near-contact systems are likely to be rotating synchronously. For a circular orbit, the following equation can be used to obtain the equatorial rotational velocity (Carquillat and Prieur, 2007), where is the radius of the considered component and is the orbital period. The projected rotation velocity of the components of IW UMa and V367 Gem can be estimated. For IW UMa, and ; for V367 Gem, and .
The spectra types and temperature of primary stars in IW UMa and V367 Gem indicate that both systems are AF-type Algols, with primary components classified as A- or F- type stars. We combined our two targets with the AF-type Algols collected by Wang et al. (2022) to investigate the evolutionary stages of these systems. The distributions of a Hertzsprung–Russell (HR) diagram are illustrated in Figure 19. Different symbols are used to differentiate AF-type Algols from our target stars, while various colors indicate the primary and secondary components. The primary components of AF-type Algol systems are predominantly located near the zero-age main sequence, whereas nearly all secondary components have evolved away from the main sequence. The mass-luminosity (M-L) relation, mass-radius (M-R) relation, and mass-temperature (M-T) relation of the components are shown in Figure 20. Different symbols are used to differentiate AF-type Algols and our target stars, and various colors indicate the primary and secondary components, respectively. The blue dash lines indicate the corresponding empirical relationships of main sequence (MS) stars (Eker et al., 2018). The primary components of AF-type Algols closely follow the empirical relationships of main sequence stars. In contrast, the secondary components are brighter, larger, and a little hotter than MS stars of the same mass.
FIGURE 19
FIGURE 20
We noticed that the primary components contribute about and to the luminosity of IW UMa and V367 Gem, respectively. The primary component of IW UMa is approximately 19 times brighter than the secondary, while the radius of the primary component is approximately 1.6 times as large as that of the secondary. Even the primary component of V367 Gem is 13 times brighter than the secondary with nearly the same radius. The secondary components of both targets may have evolved through the red giant phase with an expanded radius.
Statements
Data availability statement
The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.
Author contributions
X-MT: conceptualization, formal analysis, funding acquisition, investigation, methodology, software, writing–original draft, and writing–review and editing. L-QJ: data curation, methodology, resources, visualization, and writing–review and editing. Z-HW: formal analysis, funding acquisition, resources, supervision, and writing–review and editing. J-JW: funding acquisition, investigation, methodology, project administration, validation, and writing–review and editing.
Funding
The authors declare that financial support was received for the research, authorship, and/or publication of this article. This work is supported by the National Natural Science Foundation of China (Grant Nos. 12103030 and 12303038), the Natural Science of Shandong Province (Grant No. ZR2021QA082), the Natural Science Foundation of Xinjiang Uygur Autonomous Region (Grant No. 2022DO1A164), the Joint Research Found (Grant No. U1831109) in Astronomy under cooperative agreement between the National Science Foundation of China (NSFC) and Chinese Academy of Sciences (CAS), and the Young Talent Project of the Yunnan Revitalization Talent Support Program and the CAS Light of West China Program. This work is partly supported by the Research Fund of the Key Laboratory of Aircraft Environment Control and Life Support, MIIT, Nanjing University of Aeronautics and Astronautics (Grant No. KLAECLS–E–202203).
Acknowledgments
We acknowledge the support of the staff of the Xinglong 80 cm and the Xinglong 2.16 m telescopes. This research was partially supported by the Open Project Program of the Key Laboratory of Optical Astronomy, National Astronomical Observatories, Chinese Academy of Sciences. This paper has used data collected by the Transiting Exoplanet Survey Satellite (TESS) mission, for which funding is provided by the NASA Explorer Program. We would like to thank the editor and the referees for the very valuable and useful comments that helped improve this paper.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fspas.2024.1497989/full#supplementary-material
References
1
ApplegateJ. H. (1992). A mechanism for orbital period modulation in close binaries. Astrophysical J.385, 621. 10.1086/170967
2
AžusienisA.StraižysV. (1969). Improved determination of the response curves and parameters of the UBV system. Summary of results. Sov. Astron.13, 316.
3
BessellM. S. (1983). VRI photometry: an addendum. Publ. Astronomical Soc. Pac.95, 480–488. 10.1086/131196
4
BorovickaJ.SarounovaL. (1996). The Period and lightcurve of NSV 4497. Inf. Bull. Var. Stars4402, 1.
5
BressanA.MarigoP.GirardiL.SalasnichB.Dal CeroC.RubeleS.et al (2012). PARSEC: stellar tracks and isochrones with the PAdova and TRieste Stellar Evolution Code. Mon. Notices R. Astronomical Soc.427, 127–145. 10.1111/j.1365-2966.2012.21948.x
6
CarquillatJ. M.PrieurJ. L. (2007). Contribution to the search for binaries among Am stars - VIII. New spectroscopic orbits of eight systems and statistical study of a sample of 91 Am stars. Mon. Notices R. Astronomical Soc.380, 1064–1078. 10.1111/j.1365-2966.2007.12143.x
7
ChabrierG. (2001). The galactic disk mass budget. I. Stellar mass function and density. Astrophysical J.554, 1274–1281. 10.1086/321401
8
ChenW.-C.LiX.-D.QianS.-B. (2006). Orbital evolution of Algol binaries with a circumbinary disk. Astrophysical J.649, 973–978. 10.1086/506433
9
ChenY.BressanA.GirardiL.MarigoP.KongX.LanzaA. (2015). PARSEC evolutionary tracks of massive stars up to 350 M⊙ at metallicities 0.0001 ≤ Z ≤ 0.04. Mon. Notices R. Astronomical Soc.452, 1068–1080. 10.1093/mnras/stv1281
10
ChenY.GirardiL.BressanA.MarigoP.BarbieriM.KongX. (2014). Improving PARSEC models for very low mass stars. Mon. Notices R. Astronomical Soc.444, 2525–2543. 10.1093/mnras/stu1605
11
EkerZ.BakışV.BilirS.SoyduganF.SteerI.SoyduganE.et al (2018). Interrelated main-sequence mass-luminosity, mass-radius, and mass-effective temperature relations. Mon. Notices R. Astronomical Soc.479, 5491–5511. 10.1093/mnras/sty1834
12
ErdemA.DoğruS. S.SoyduganF.ÇiçekC.DemircanO. (2010). Period studies of five neglected algol-type binaries: RW cet, BO Gem, DG lac, SW oph and WY per. New Astron.15, 628–636. 10.1016/j.newast.2010.02.007
13
FanZ.WangH.JiangX.WuH.LiH.HuangY.et al (2016). The Xinglong 2.16-m telescope: current instruments and scientific projects. Publ. Astronomical Soc. Pac.128, 115005. 10.1088/1538-3873/128/969/115005
14
HallD. S. (1989). The relation between RS CVn and Algol. Space Sci. Rev.50, 219–233. 10.1007/BF00215932
15
HintzE. G.BrownP. J. (2007). Revised periods for QS geminorum and V367 geminorum. Publ. Astronomical Soc. Pac.119, 274–283. 10.1086/516602
16
IbanoǧluC.SoyduganF.SoyduganE.DervişoǧluA. (2006). Angular momentum evolution of Algol binaries. Mon. Notices R. Astronomical Soc.373, 435–448. 10.1111/j.1365-2966.2006.11052.x
17
KolevaM.PrugnielP.BouchardA.WuY. (2009). ULySS: a full spectrum fitting package. Astronomy Astrophysics501, 1269–1279. 10.1051/0004-6361/200811467
18
KreinerJ. M. (2004). Up-to-Date linear elements of eclipsing binaries. Acta Astron.54, 207–210.
19
LiaoW. P.QianS. B.LiL. J.ZhouX.ZhaoE. G.ZhangJ. (2017). DI hya: a near-contact binary with a close-in companion. Publ. Astronomical Soc. Pac.129, 034201. 10.1088/1538-3873/aa5869
20
LiaoW. P.QianS. B.LiuN. P. (2012). A CCD photometric study of the contact binary star GSC 03526-01995. Astronomical J.144, 178. 10.1088/0004-6256/144/6/178
21
LiaoW. P.QianS. B.SarotsakulchaiT. (2019). Physical properties of the close-in tertiary in the southern triple-lined system VZ lib. Astronomical J.157, 207. 10.3847/1538-3881/ab17d4
22
LiaoW. P.SarotsakulchaiT. (2019). First photometric investigations of the solar-type binary FV CVn in multiple systems. Publ. Astronomical Soc. Pac.131, 014202. 10.1088/1538-3873/aaeae2
23
Lightkurve CollaborationCardosoJ. V. D. M.HedgesC.Gully-SantiagoM.SaundersN.CodyA. M.et al (2018). Lightkurve: kepler and TESS time series analysis in Python. Astrophys. Source Code Libr. Rec. ascl:1812.013.
24
PrugnielP.SoubiranC. (2001). A database of high and medium-resolution stellar spectra. Astronomy Astrophysics369, 1048–1057. 10.1051/0004-6361:20010163
25
QianS. (2000). Orbital period changes and possible mass and angular momentum loss in two Algol-type binaries: RW Coronae Borealis and TU Herculis. Astronomical J.119, 901–905. 10.1086/301217
26
QianS. (2001). Possible mass and angular momentum loss in algol-type binaries. III. TU cancri, FZ delphini, AY geminorum, VZ leonis, FH orionis, IU persei, XZ persei, and BE vulpeculae. Astronomical J.121, 1614–1622. 10.1086/319404
27
QianS. (2002). Period changes in four algol-type binaries: SX dra, AV del, CU peg and DK peg. Astrophysics Space Sci.282, 399–409. 10.1023/A:1020876723842
28
QianS.-B.HeJ.-J.ZhangJ.ZhuL.-Y.ShiX.-D.ZhaoE.-G.et al (2017). Physical properties and catalog of EW-type eclipsing binaries observed by LAMOST. Res. Astronomy Astrophysics17, 087. 10.1088/1674-4527/17/8/87
29
QianS. B.LiuN. P.LiK.HeJ. J.ZhuL. Y.ZhaoE. G.et al (2013a). BI vulpeculae: a siamese twin with two very similar cool stars in shallow contact. Astrophysical J. Suppl. Ser.209, 13. 10.1088/0067-0049/209/1/13
30
QianS. B.LiuN. P.LiaoW. P.HeJ. J.LiuL.ZhuL. Y.et al (2013b). First photometric investigation of the newly discovered W UMa-type binary star MR com. Astronomical J.146, 38. 10.1088/0004-6256/146/2/38
31
QianS. B.WangJ. J.ZhuL. Y.SnoonthornthumB.WangL. Z.ZhaoE. G.et al (2014a). Optical flares and a long-lived dark spot on a cool shallow contact binary. Astrophysical J. Suppl. Ser.212, 4. 10.1088/0067-0049/212/1/4
32
QianS. B.ZhouX.ZolaS.ZhuL. Y.ZhaoE. G.LiaoW. P.et al (2014b). AL cassiopeiae: an F-type contact binary system with a cool stellar companion. Astronomical J.148, 79. 10.1088/0004-6256/148/5/79
33
RickerG. R.WinnJ. N.VanderspekR.LathamD. W.BakosG. Á.BeanJ. L.et al (2014). “Transiting Exoplanet survey satellite (TESS),” in Space telescopes and instrumentation 2014: optical, infrared, and millimeter wave. Editors OschmannJ.JacobusM.ClampinM.FazioG. G.MacEwenH. A. (Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series), 9143, 914320
34
RucińskiS. M. (1969). The proximity effects in close binary systems. II. The bolometric reflection effect for stars with deep convective envelopes. Acta Astron.19, 245.
35
SarnaM. J.MuslimovA.YerliS. K. (1997). Magnetic activity and evolution of Algol-type stars. Mon. Notices R. Astronomical Soc.286, 209–214. 10.1093/mnras/286.1.209
36
SarnaM. J.YerliS. K.MuslimovA. G. (1998). Magnetic activity and evolution of Algol-type stars - II. Mon. Notices R. Astronomical Soc.297, 760–768. 10.1046/j.1365-8711.1998.01539.x
37
SoyduganF.ErdemA.DoğruS. S.AliçavuşF.SoyduganE.ÇiçekC.et al (2011). Period studies of classical algol-type binaries II: UX Leo, RW Mon, EQ ori, XZ UMa and AX vul. New Astron.16, 253–264. 10.1016/j.newast.2010.11.006
38
TangJ.BressanA.RosenfieldP.SlemerA.MarigoP.GirardiL.et al (2014). New PARSEC evolutionary tracks of massive stars at low metallicity: testing canonical stellar evolution in nearby star-forming dwarf galaxies. Mon. Notices R. Astronomical Soc.445, 4287–4305. 10.1093/mnras/stu2029
39
TianX.-M.WangZ.-H.ZhuL.-Y.YangX.-L. (2023). A new catalog of Am-type chemically peculiar stars based on LAMOST. Astrophysical J. Suppl. Ser.266, 14. 10.3847/1538-4365/acc4b5
40
TianX.-M.ZhuL.-Y. (2019). V723 Persei: a short-period Algol-like near-contact binary. Publ. Astronomical Soc. Jpn.71, 66. 10.1093/pasj/psz044
41
TokovininA.ThomasS.SterzikM.UdryS. (2006). Tertiary companions to close spectroscopic binaries. Astronomy Astrophysics450, 681–693. 10.1051/0004-6361:20054427
42
ToutC. A.EggletonP. P. (1988). Tidal enhancement by a binary companion of stellar winds from cool giants. Mon. Notices R. Astronomical Soc.231, 823–831. 10.1093/mnras/231.4.823
43
van HammeW. (1993). New limb-darkening coefficients for modeling binary star light curves. Astronomical J.106, 2096. 10.1086/116788
44
WangZ. H.ZhuL. Y.YueY. F. (2022). Evolutionary inference and statistical constraints on Algols including SD2-type near contact binaries. Mon. Notices R. Astronomical Soc.511, 488–500. 10.1093/mnras/stac037
45
WilsonR. E. (1990). Accuracy and efficiency in the binary star reflection effect. Astrophysical J.356, 613. 10.1086/168867
46
WilsonR. E. (2012). Spotted star light curves with enhanced precision. Astronomical J.144, 73. 10.1088/0004-6256/144/3/73
47
WilsonR. E.DevinneyE. J. (1971). Realization of accurate close-binary light curves: application to MR cygni. Astrophysical J.166, 605. 10.1086/150986
48
YangY. G.WeiJ. Y. (2009). Evolutionary status of RV trianguli and its related algol-type binaries. Astronomical J.137, 226–235. 10.1088/0004-6256/137/1/226
49
ZhangJ.QianS.-B.WuY.ZhouX. (2019). Unbiased distribution of binary parameters from LAMOST and kepler observations. Astrophysical J. Suppl. Ser.244, 43. 10.3847/1538-4365/ab442b
Summary
Keywords
eclipsing binary, Algol, light curve, period analysis, evolution
Citation
Tian X-M, Jiang L-Q, Wang Z-H and Wang J-J (2025) The first photometric investigation of two short-period Algol-type eclipsing binaries: IW UMa and V367 Gem. Front. Astron. Space Sci. 11:1497989. doi: 10.3389/fspas.2024.1497989
Received
18 September 2024
Accepted
03 December 2024
Published
09 January 2025
Volume
11 - 2024
Edited by
Marcella Marconi, Astronomical Observatory of Capodimonte (INAF), Italy
Reviewed by
Andrzej S. Baran, Missouri State University, United States
Mkrtichian Egishe David, National Astronomical Research Institute of Thailand, Thailand
Updates
Copyright
© 2025 Tian, Jiang, Wang and Wang.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Xiao-man Tian, txmjlx2018@163.com
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.