The eclipsing Algol OGLE-LMC-DPV-065

On the long cycle variability of the Algol OGLE-LMC-DPV-065 and its stellar, orbital and disk parameters

R.E. Mennickent, M. Cabezas, G. Djurašević, T. Rivinius, P. Hadrava, R. Poleski, I. Soszyński, L. Celedón, N. Astudillo-Defru, A. Raj, J. G. Fernández-Trincado, L. Schmidtobreick, C. Tappert, V. Neustroev, I. Porritt
Universidad de Concepción, Departamento de Astronomía, Casilla 160-C, Concepción, Chile
Astronomical Institute of the Academy of Sciences of the Czech Republic, Boční II 1401/1, Prague, 141 00, Czech Republic
Astronomical Observatory, Volgina 7, 11060 Belgrade 38, Serbia
Isaac Newton Institute of Chile, Yugoslavia Branch
European Southern Observatory
Department of Astronomy, Ohio State University, 140 W. 18th Ave., Columbus, OH 43210, USA
Warsaw University Observatory, Al. Ujazdowskie 4, 00-478 Warszawa, Poland
Departamento de Matemática y Física Aplicadas, Universidad Católica de la Santísima Concepción, Alonso de Rivera 2850, Concepción, Chile
Indian Institute of Astrophysics, II Block Koramangala, Bangalore 560034, India
Institut Utinam, CNRS UMR6213, Univ. Bourgogne Franche-Comté, OSU THETA, Observatorie de Besançon, BP 1615, 25010 Besançon Cedex, France
Instituto de Astronomía y Ciencias Planetarias, Universidad de Atacama, Copayapu 485, Copiapó, Chile
Instituto de Física y Astronomía, Universidad de Valparaíso, Chile
Astronomy research unit, P.O. Box 3000 90014 University of Oulu, Finland
Turitea Observatory, Palmerston North, New Zealand
E-mail: Based on ESO proposal 092.D-0385(A) and CNTAC 2014B-13.

OGLE-LMC-DPV-065 is an interacting binary whose double-hump long photometric cycle remains hitherto unexplained. We analyze photometric time series available in archive datasets spanning 124 years and present the analysis of new high-resolution spectra. A refined orbital period is found of 100316267 00000056 without any evidence of variability. In spite of this constancy, small but significant changes in timings of the secondary eclipse are detected. We show that the long period continuously decreases from 350 to 218 days during 13 years, then remains almost constant for about 10 years. Our study of radial velocities indicates a circular orbit for the binary and yields a mass ratio of 0.203 0.001. From the analysis of the orbital light curve we find that the system contains 13.8 and 2.81 M  stars of radii 8.8 and 12.6 R and absolute bolometric magnitudes -6.4 and -3.0, respectively. The orbit semi-major axis is 49.9 R and the stellar temperatures are 25460 K and 9825 K. We find evidence for an optically and geometrically thick disk around the hotter star. According to our model, the disk has a radius of 25 R, central and outer vertical thickness of 1.6 R and 3.5 R, and temperature of 9380 K at its outer edge. Two shock regions located at roughly opposite parts of the outer disk rim can explain the light curves asymmetries. The system is a member of the double periodic variables and its relatively high-mass and long photometric cycle make it similar in some aspects to Lyrae.

stars: binaries: eclipsing, close, spectroscopic, stars: activity, circumstellar matter, fundamental parameters

1 Introduction

Stellar magnetic dynamos are relatively common in solar type stars, and magnetic activity in binaries containing GK dwarfs is well documented in the RS CVn systems (Hall, 1989). The situation in Algol-type variables is less clear. Algols are semi-detached binaries with intermediate mass components, where the less massive star, dubbed secondary or donor, is more evolved than the more massive star, named gainer or primary. Some authors have proposed that orbital period changes observed in some close binaries might be driven by magnetic cycles through the Applegate (1992) mechanism; the angular momentum of the star and the binary is redistributed during the magnetic cycle producing the observed orbital period changes (Lanza, Rodonò, & Rosner, 1998; Lanza & Rodonò, 1999). Further studies indicate that the presence of a dynamo may modulate the mass transfer rate in Algol systems, leading to a characteristic impact of the dynamo cycle on the system luminosity (Bolton, 1989; Meintjes, 2004). In this context the existence of a group of hot Algols showing a long photometric cycle lasting on average about 33 times the orbital period might be relevant, since this variability has been recently interpreted in terms of a magnetic dynamo (Schleicher & Mennickent, 2017). If this hypothesis turns to be correct, one may deduce that the stellar dynamo is also active in the hot, rapidly rotating (orbitally synchronized) A-type giants in some semidetached Algols. In fact, for the Algol binary V393 Sco indirect evidence for magnetism in the secondary star has been deduced from the presence of chromospheric emission lines (Mennickent, Schleicher, & San Martin-Perez, 2018). These authors note that the spin-up of the donor during mass-transfer stage increases its dynamo number, likely enhancing the probability of occurrence of a magnetic dynamo at the semi-detached stage.

The aforementioned group of hot Algols showing long photometric cycles additional to their orbital variability is named Double Periodic Variables (DPVs, Mennickent et al., 2003; Mennickent, Otero, & Kołaczkowski, 2016; Poleski et al., 2010; Pawlak et al., 2013; Mennickent, 2017). DPVs are semidetached binaries typically consisting of a A/F/G giant star filling its Roche lobe and transferring mass onto a B-type primary surrounded by a circumprimary disk. Among Galactic DPVs, one famous example is Lyrae (Guinan, 1989; Harmanec et al., 1996; Harmanec, 2002).

Few extragalactic DPVs have been studied at some detail. Among them, the case of OGLE-LMC-DPV-065 (OGLE05200407-6936391; R.A. = 05:20:04.07, Dec. = -69:36:39.1) is notable, since it is one of the brightest DPVs in the LMC ( = 14.74, = -0.07), and shows a remarkable change in the long period from 350 to 210 days in 15 years that clearly stands out among the rest of the DPVs. In addition, the system is eclipsing, with a 1.4 mag deep main eclipse and a comparatively large amplitude of the long cycle of 0.3 mag, in the band. To date, there is no other DPV with such a remarkable change in the long cycle. The orbital period has been reported as = 10031645 0.000033 (Poleski et al., 2010).

The above credentials make OGLE-LMC-DPV-065 an ideal target for a deeper study. If the long cycle is related to changes in the mass transfer driven by a magnetic dynamo, it might show up in spectroscopic and photometric signatures during the long cycle. For this reason we conducted a long-term spectroscopic monitoring of this target with UVES at the VLT (Sec. 2.3). For the sake of order and clarity we have divided our work in two parts. In this first paper we analyze the available photometric time series making use of archive data, present our new high-resolution spectroscopic observations, calculate the system and orbital parameters and provide a solution for the stellar radius, mass, luminosity, surface gravity along with a characterization of the accretion disk. In a second forthcoming paper we will provide an analysis of the spectroscopic changes during the long cycle and a study of the evolutionary stage of the binary. We notice that a short and preliminary spectroscopic study of this object based on the data presented in this paper has been presented in a recent conference (Cabezas et al., 2019).

2 Observations and methods

2.1 Photometric observations

We included OGLE-II (Szymanski, 2005)111 and OGLE-III/IV data222 OGLE-III/IV data kindly provided by the OGLE team.. The OGLE-IV project is described by Udalski, Szymański, & Szymański (2015). Poleski et al. (2010) published the OGLE-II and OGLE-III -band333 and -band444 data of this star. We also considered 460 magnitudes from the Digitalized Harvard plates (DASCH project)555 covering 59.4 years, since August 1893 to January 1953. In addition, we obtained new photometry with the CTIO 1.3m telescope operated by the SMARTS consortium in service mode between November 2014 and October 2015, with the ANDICAM camera and filters and . Another data set was collected by Ian Porritt in Turitea Observatory, New Zealand, with the 0.36-meter Meade telescope and a yellow filter. These new data were reduced in the usual way, removing bias and performing flat field corrections in the images and calculating differential magnitudes with respect to constant comparison stars. Finally, 664 -band magnitudes were included from the ASAS-SN catalogue. The photometric observations analyzed in this paper amount to 3099 data points, cover 124 years and are summarized in Table 1.

Source N HJD HJD mag std. band
DASCH 460 12697.8482 34399.4995 14.996 0.219
OGLE-II 915 50455.6744 51873.7744 14.898 0.218
OGLE-III 504 52123.9345 54953.5268 14.907 0.246
OGLE-IV 73 55326.4931 57710.7482 14.901 0.302
CTIO 97 56964.7927 57327.7354 14.901 0.317
OGLE-II 95 50467.7237 51631.5633 14.908 0.244
OGLE-III 90 52990.6851 54948.4703 14.929 0.245
Turitea 106 56342.9193 56467.8381 14.918 0.368
ASAS-SN 664 56789.4535 57974.8870 14.918 0.116
CTIO 95 56964.7954 57327.7381 14.918 0.330
Table 1: Summary of photometric observations considered in this paper. The number of measurements, starting and ending times for the series, average magnitude and standard deviation (in magnitudes) are given. Single point uncertainties in the -band and -band for OGLE data are between 4 and 6 mmag. The zero point for HJD is 2 400 000. See comment on the average magnitudes in the text.

2.2 Light curve disentangling

We separated the light curve into long- and a short-period components. For that we used an algorithm especially designed to disentangle multi-periodic light curves through the analysis of their Fourier component amplitudes. The method is described in Mennickent et al. (2012) and a short summary is given here. The main frequency is found with a period searching algorithm, this is usually the orbital frequency. A least-square fit is then applied to the light curve with a fitting function consisting of a sum of sine functions representing the main frequency and their additional significant harmonics. Afterwards the residuals are inspected for a new periodicity . This new periodicity (the long cycle in the case of DPVs) and their harmonics are included in the new fitting procedure. Finally, we obtain the light curve represented by a sum of Fourier components of frequency and and their harmonics. Data residuals with respect to the second and first theoretical light curves are the photometric series representing the orbital and long cycle, respectively.

2.3 Spectroscopic observations

We were granted 25 hours for spectroscopic observations of the target with the ESO Ultraviolet and Visual Echelle Spectrograph UVES at the Kueyen telescope in the Paranal Observatory in service mode. This is a cross-dispersed echelle spectrograph designed to operate with high efficiency from the atmospheric cut-off at 300 nm to the long- wavelength limit of the CCD detectors (about 1100 nm). The light beam is split into two arms, UV-Blue and Visual-Red, within the instrument. The two arms can be operated separately or in parallel with a dichroic beam splitter. The instrument provides accurate calibration of the wavelength scale down to an accuracy of at least 50 m/s.

With the aim of covering both the orbital as well as the long cycle 27 spectra were secured between October 1, 2013 and February 1, 2015 with the dichroic#2 mode in the ranges 37604985, 57007520 and 76659445 Å. The slit widths of 09 at the blue channel and 08 at the red channels provided a resolving power of 50 000 and 55 000, respectively. The object was observed at typical airmass 1.4 and with 3000 s exposure time per single exposure. A typical S/N of 65 was achieved at 480 nm. A summary of the spectroscopic observations is given in Table 7.

Figure 1: Light curve for OGLE II (black), III (red), and IV (green) datasets along with CTIO (blue) data. The upper envelope is produced by the long-period variability. The data points fainter than about 15.2 mag are taken during primary eclipse. The data intervals given in Table 4 are also shown.
Figure 2: (up): Long period versus long cycle number. (down): Observed minus calculated times of the long cycle maximum versus long cycle number. Data are based on a test period of 240 days and are given in Table 3. Vertical dashed lines indicate spectroscopic observation times.

3 Data analysis

3.1 The orbital period

Zero points can be an issue when combining photometry obtained at different sites with different detectors, filters and sky conditions. We have shifted OGLE IV, ASAS-SN, Turitea, DASCH and CTIO magnitudes to the OGLE-II and OGLE-III averages before performing the analysis described in this paper. In Table 1 we provide only the original magnitude average for OGLE-II and OGLE-III data and also for the DASCH -band magnitude. The light curve in the -band is shown in Fig. 1.

We conducted a search for the orbital period using standard methodologies: eclipse times were measured interactively in the light curve with the computer cursor and a straight line fit was performed with the measured (epoch, time) pairs; the resulting slope gave the orbital period and their error. We also used the Period04 program, that calculates errors based on a Monte Carlo technique (Lenz & Breger 2005). The -band residuals were obtained after removing the long-term cycle (see next section). The periodicities found in different datasets are given in Table 2. We can see that the data are consistent with a constant orbital period; we find the following ephemerides for the main eclipse:

Method note
100316450 00000330 - Poleski et al. (2010)
100317800 00000780 Eclipse timings band
100316259 0.0000230 Period04 band
100316173 0.0000128 Period04 residuals
100316257 0.0000066 Period04 and bands
100316267 00000056 Weighted mean
Table 2: Summary of the search for the orbital period .

3.2 The orbital and long cycle light curves

Maxima of the long cycle were measured directly from the light curves and compared with the ephemerides for a 240 day test period, as reported in Table 3. The observed minus calculated () diagram, constructed with the observed times of maxima () and the predicted ones (), shows that the long cycle decreased at the beginning of the observations then remained more or less constant during about 14 cycles (Fig. 2). We notice that considering the MACHO data analyzed by Mennickent et al. (2005, HJD: 2448900-2451500, not included in this paper), which is previous to the OGLE data reported here, the long period has decreased from about 350 to 218 days continuously during about 13 years, before entering in a phase of almost constant period, that lasted for slightly more than 10 yr.

At every epoch we defined a local long cycle period , subtracting the observed maximum timing from the previous one and dividing by the number of elapsed times, as given in Table 3. After inspection of Fig. 2 we choose six data intervals characterized by a more or less constant long cycle and large number of observations (Table 4). This procedure allowed to apply the disentangling to every data block considering the variability of the long cycle. The resulting disentangled light curves are shown in Fig. 3, they reveal that the long cycle is double-humped and that it shape remains relatively constant. In addition, the orbital light curve shows a small but significant variability (Figs. 4 and 5): (i) on the 5th interval between 2 455 804 and 2 456 405 the system is brighter at quadratures, and produces larger scatter in the long cycle light curve, (ii) on the first interval the main eclipse seems to be shallower, (iii) significant variability is observed during secondary eclipse; the secondary eclipse seems to occur earlier in interval 1 than in interval 2, and (iv) the shape of the eclipses vary minimally during the maximum, the minimum and the secondary maximum of the long cycle, perhaps the egress of the main eclipse around phases 0.1-0.2 is shallower during the low stage. The changes in timing of minima during the secondary eclipse might indicate changes in the photo-center of the eclipsed or eclipsing source, or changes in circumstellar matter or the donor hemisphere facing the gainer. An unseen/undetected body that dynamically affects the photo-center is another possibility.

We did the same exercise with the -band but we had to use a smaller number of intervals due to the smaller number of observations in this band. The intervals are documented in Table 5. The long cycle usually has a smaller amplitude than in the -band and the orbital light curve shows subtle variability. These changes are better visualized in the combined light curve (Fig. 5).

Figure 3: (up): Orbital phase curves at the intervals 1 to 6 defined in Table 4. (Down): Long cycle phase curves with different periods. Intervals 1 to 6 defined in Table 4 are illustrated from the top-left to the down-right panels. In both panels the magnitude is differential -band and the color is used to indicate time strings of nearby data-points.
Figure 4: Up: A zoom into the secondary (left) and primary (right) eclipses for data intervals I1 and I2, as defined in Table 4. Vertical axis shows the differential -band magnitude. Down: Eclipses during main and secondary long cycle maxima (0.9 1.1 and 0.4 0.6, respectively) and minima ( 0.2 0.4 and 0.6 0.8). Vertical axis shows the differential -band magnitude.
Figure 5: Up: Orbital (left) and long cycle (right) -band differential light curves at the intervals defined in Table 5. The color is used to indicate time strings of nearby data-points. Down: Combined data -band light curve. Colors indicate time strings of subsequent points.
Figure 6: Disentangled donor spectrum and the best fit model at two different spectral regions.
N (d) (d) band/source
-8 571.68 4.49 697.89 -126.21 OGLE
-7 859.75 7.05 937.89 -78.14 288.06 OGLE
-6 1125.58 9.48 1177.89 -52.31 265.83 OGLE
-4 1645.23 8.04 1657.89 -12.66 259.83 OGLE
-2 2131.23 10.49 2137.89 -6.66 243.00 OGLE
0 2617.89 15.47 2617.89 0.00 243.33 OGLE
2 3077.87 7.04 3097.89 -20.02 229.99 OGLE
3 3311.28 4.98 3337.89 -26.61 233.41 OGLE
5 3747.12 7.00 3817.89 -70.77 217.92 CTIO
5 3749.41 6.93 3817.89 -68.48 219.06 OGLE
6 3973.60 2.52 4057.89 -84.29 224.19 OGLE
8 4403.01 10.98 4537.89 -134.88 214.71 OGLE
8 4408.67 7.00 4537.89 -129.22 217.54 CTIO
10 4835.37 9.55 5017.89 -182.52 213.35 OGLE
13 5477.80 5737.89 -260.09 214.14 OGLE
17 6343.85 6697.89 -354.04 216.51 Turitea
17 6347.25 6697.89 -350.64 217.36 OGLE
17 6353.84 8.00 6697.89 -344.05 219.01 CTIO
20 7011.15 4.00 7417.89 -406.74 219.10 OGLE
20 7011.20 2.00 7417.89 -406.69 219.12 CTIO
20 7011.23 10.00 7417.89 -406.66 219.13 ASAS-SN
22 7448.75 7897.89 -449.14 218.76 OGLE
22 7453.28 10.00 7897.89 -444.61 221.03 ASAS-SN
Table 3: Observed and calculated times of long cycle maxima. is the cycle number. HJD are referred to 2 450 000. stands for and for . Errors are given when available. We have used the linear ephemerides with zero point 2617.89 and period 240 days. An estimate of the local period is given subtracting the observed timing from the previous one and dividing by the number of elapsed cycles.
label range (d)
I1 538 455 1000 282 855.66
I2 377 1000 1900 258 1128.82
I3 109 1900 2800 240 2618.78
I4 125 2800 3500 231 3316.75
I5 330 3500 6500 216 3753.70
I6 104 6500 7500 219 7009.62
Table 4: Data intervals used for the long-period analysis. HJD are referred to 2 450 000. is the number of -band data points. Times for long cycle maxima are given.
label range (d)
V1 95 468 1632 265 580.73
V2 90 2991 4949 219 2858.65
V3 95 6965 7328 220 6789.63
Table 5: Data intervals used for disentangling the -band light curve. HJD are referred to 2 450 000. is the number of data-points. Times for long cycle maxima are given.

3.3 Spectra components and orbital/system parameters

In order to obtain the radial velocities and orbital parameters we used the KOREL code (Hadrava, 1995, 1997) based on the method of Fourier disentangling, yielding directly the orbital parameters together with the decomposed spectra of the multiple stellar system under study. In addition, we also used the code FOTEL (Hadrava, 1990) to estimate the errors of the orbital parameters.

We notice that the system can be classified as SB2, i.e. both stellar components are detected in the spectrum, in particular in helium and hydrogen lines. The detected components correspond to an early B-type (primary or gainer), and an early A-type (secondary or donor). The method of spectra disentangling does not use any template or another information about the laboratory wavelengths of the spectral lines, therefore the systemic velocity is set to zero. For this reason we adopted an average of systemic velocities calculated by Gaussian adjustments for different spectral lines of each component. For the gainer we obtain km s and for the donor km s. The lines used in this calculation are shown in Table 6.

We notice that S II/III measurements systematically differ from other lines, suggesting a different formation place. For this reason they are not included in the above calculation. The disentangled spectra are shown in Figs. 6 and 7.

We performed the calculation of radial velocity in seven regions of every spectrum. These regions were chosen because they include several narrow, unblended and well identified metallic lines. All our spectra were prepared with a routine written in IRAF6 and the sampling auxiliary code PREKOR (Hadrava, 2004) was used. To diminish the numerical errors of the disentangling we sampled each spectral region in the maximum number of bins allowed by the code, viz. 4096. This results in the average resolution 0.726 km s per bin, which is higher than the original resolution on the spectrograph detector.

Radial velocities obtained with KOREL for the cases of circular orbit are given in Table 7 and their best fit is shown in Fig. 8. The operation of the KOREL code is described in Hadrava (2004).

The radial velocity for each components is given by


where the sum is realized on the orbits that influence the movement of the star. The true anomaly is calculated according to


where is obtained from the solution of Kepler’s equation.

The orbital parameters obtained by disentangling of the seven spectral regions are summarized in Table 8. The Solution I, which we accept for our modeling of photometry, has been obtained using an independent disentangling of each region separately and then calculating mean solutions and standard deviations of each parameter. Solution II is the simultaneous (”multi-region”) disentangling of all the regions together. The errors of the parameters were obtained using the Bayesian estimate, i.e. from the moments of the Bayesian probability distribution (Hadrava, 2016) We have also solved the radial-velocity curve using the FOTEL code with the input radial velocities obtained from the disentangling. The resulting values of parameters were within the error-bars of the Solution I, but their errors were for about one order underestimated, so we skipped this solution. Finally, the multi-region Solution III is to verify that a possible eccentricity of the orbit can be neglected.

Spectral line
ion (Å) (km s) (km s)
SiII 4128.054 - 281.136
SiII 4130.894 - 278.431
HeI 4143.76 279.219 -
SII 4153.068 290.429 -
NII 4227.74 276.360 -
FeII 4233.172 - 280.476
NII 4236.91 277.24 -
NII 4241.78 272.029 -
CrII 4242.364 - 275.769
ScII 4246.822 - 275.291
SIII 4253.589 289.218 -
OII 4414.884 275.352 -
OII 4416.974 274.255 -
HeI 4437.551 276.146 -
TiII 4443.794 - 277.682
NII 4447.04 273.909 -
TiII 4533.960 - 281.188
SII 4552.410 289.434 -
SiIII 4567.840 276.484 -
SiIII 4574.757 274.978 -
TiII 4549.617 - 284.985
CrII 4558.650 - 277.747
TiII 4563.757 - 281.192
TiII 4571.968 - 282.472
OII 4590.971 274.638 -
OII 4596.174 273.442 -
NII 4607.153 272.240 -
SiII 4621.418 271.289 -
FeII 4629.336 - 278.677
OII 4638.854 276.561 -
OII 4641.811 276.864 -
OII 4649.138 276.155 -
OII 4699.21 270.561 -
OII 4705.355 276.678
HeI 4713.143 278.230 -
FeII 4731.453 - 280.792
Mean (no SII/III) 275.1322.325 279.7632.832
Mean (all) 277.0305.464 -
Table 6: RV zero points derived from different lines.
Date-ut HJD
-2450000.0 km km
2013-10-02 6567.7689 0.547 0.982 12.082        0.789 -59.521        3.382
2013-10-04 6569.8218 0.751 0.992 43.182        1.020 -209.288    5.953
2013-10-06 6571.8038 0.949 0.001 14.536        0.830 -67.676        2.214
2013-10-07 6572.8075 0.049 0.005 -13.402        0.618 61.623        3.562
2013-10-19 6584.7162 0.236 0.060 -43.399        1.163 207.660        7.547
2013-10-22 6587.8479 0.548 0.074 12.647        0.685 -60.792        3.658
2013-12-22 6648.5833 0.603 0.351 25.914        0.940 -124.089    5.008
2013-12-24 6650.7561 0.819 0.361 39.975        0.886 -190.796    4.926
2013-12-31 6657.5987 0.501 0.393 0.115            0.649 0.015        2.346
2014-01-04 6661.6647 0.907 0.411 24.468        0.615 -117.341    2.902
2014-01-18 6675.6735 0.303 0.475 -41.331        0.952 198.184        6.020
2014-01-19 6676.6249 0.398 0.480 -26.363        0.566 126.834        2.904
2014-02-11 6699.5976 0.688 0.584 40.177        0.934 -193.185    5.941
2014-02-15 6703.5569 0.083 0.603 -21.199        0.864 103.465        3.704
2014-02-16 6704.5526 0.182 0.607 -39.697        0.839 189.726        6.078
2014-09-01 6901.8820 0.853 0.508 34.533        0.659 -169.289    4.214
2014-11-03 6964.8370 0.128 0.796 -31.257        0.812 149.933        5.166
2014-11-20 6981.7855 0.818 0.873 39.944        0.924 -191.460    5.220
2014-11-22 6983.7327 0.012 0.882 -3.961        1.684 13.814        2.907
2014-11-25 6986.7551 0.313 0.896 -40.347        0.907 193.721        5.611
2014-11-26 6987.7523 0.413 0.900 -23.026        0.504 111.558        3.302
2014-11-27 6988.7666 0.514 0.905 3.562            0.391 -16.488        3.285
2014-12-08 6999.6148 0.595 0.954 24.437        1.068 -116.546    5.316
2014-12-09 7000.7687 0.710 0.960 42.190        1.088 -202.936    5.340
2014-12-14 7005.6695 0.199 0.982 -41.140        0.994 198.112        6.064
2015-01-20 7042.5709 0.877 0.151 30.758        0.632 -147.288    3.797
2015-02-01 7054.5922 0.075 0.205 -19.228        0.538 94.275        3.972
Table 7: RVs for the primary and secondary components from KOREL solutions in a circular orbit. The radial velocity is the average from each spectral region and we considered the standard deviation as error. Orbital phases are given for the orbital ephemerides given by Eq. (1) and long cycle phases for a long period of 219 days and = 57009.62.
Parameter I II III
[d] 10.0316267 (fixed)
* 92.310.02 92.3050.004 94.79
[km s] 42.600.97 42.440.33 42.45
[km s] 210.56.4 214.11.8 213.5
0.2030.008 0.1980.002 0.199
[deg] 90 90 178.7
0 0 0.021
Table 8: Orbital parameters obtained in different solutions
Figure 7: The component spectra in our seven selected regions disentangled using  KOREL. For each panel, the upper spectrum corresponds to the primary  component and the lower spectrum to the secondary component. The  wavelength scale is shifted for the mean -velocity to correspond to the  laboratory wavelengths of the lines.
Figure 8: Fit of theoretical RVs to the average velocities from Table 7.

Once disentangled the donor spectrum, we compared it with a grid of solar-metallicity synthetic models constructed with SPECTRUM666 grayro/spectrum/spectrum.html and search for the synthetic spectrum minimizing residuals. We find the best fit with a stellar spectrum of = 9825 75 K, = 53 3 km s and log g = 3.2 0.2. Comparisons of the donor disentangled spectrum with the best fit model are shown in Fig. 6. Similarly, from the region 41204199 Å  we obtained a model with T = 22000 K and = 70.6 km s for the gainer.

4 Models for the system

4.1 Model for an optically thick disk around the gainer

Part of the phenomenology of DPVs has been associated with the presence of an optically thick disk around the gainer, probably feed by a Roche-lobe filling donor star (e.g. Garcés L. et al., 2018). Consistently, we model the orbital light curve of OGLE-LMC-DPV-065 considering the stellar fluxes of the two stars, the contribution of an accretion disk around the primary and eventually the light contribution of hot/bright spots located in the outer disk rim. The basic elements of the binary system model with a plane-parallel disk and the corresponding light curve synthesis procedure are described by Djurašević (1992a, b, 1996). The code has been successfully applied to several close binaries including the well-studied binary system Lyrae (e.g. Djurašević et al., 2010, 2012; Garrido et al., 2013; Mennickent et al., 2012; Mennickent & Djurašević, 2013).

We assume that the disk is optically and geometrically thick and that its outer edge is approximated by a cylindrical surface. The vertical thickness of the disk can change linearly with radial distance, allowing different disk’s conical shapes: plane- parallel, concave or convex. The geometrical parameters of the disk are its radius (), its vertical thickness at the outer edge (d) and the vertical thickness at the inner boundary (d). The cylindrical edge of the disk is characterized by its temperature, T, and the conical surface of the disk by a radial temperature profile inspired in the temperature distribution proposed by Zola (1991):


We further assume that the disk is in physical contact and thermal equilibrium with the gainer, so its inner radius and corresponding temperature are equal to the radius and temperature of the star (, T). The temperature of the disk at the edge (T) and the temperature exponent (a), as well as the radii of the gainer () and of the disk () are free parameters, determined by solving the inverse problem (see Section 4.2).

Motivated by previous research on DPVs (Mennickent et al. 2016), our model includes active regions on the edge of the disk. These active regions are usually revealed in Doppler tomography maps of disks in Algol binaries (Richards, 2004). These regions have higher local temperatures than the disk, and produce a non-uniform distribution of radiation. We consider two active regions: a hot spot (h) and a bright spot (b), characterized by their temperatures T, angular dimensions (radius) and longitudes . The longitude is measured from the line joining the gainer and donor centers in the direction opposite to the orbital motion in the orbital plane. These parameters are also determined by solving the inverse problem. We also consider a possible departure of symmetry of light emerging from the hot spot due, for instance, to the impact of the gas stream in the disk. This deviation is described by the angle between the line perpendicular to the local disk edge surface and the direction of the hot spot maximum radiation in the orbital plane. The second spot in the model (here named bright spot), simulates the spiral structure of an accretion disk, observed in hydrodynamical calculations (Heemskerk, 1994). The tidal forces exerted by the donor star produce a spiral shock, observed as one or two extended spiral arms in the outer disk regions. This bright spot can also be interpreted as a region where the disk significantly deviates from the circular shape.

Two potential limitations of the code need to be briefly mentioned: the lack of treatment of the donor irradiation by the hot spot, and the lack of inclusion of a possible not eclipsed additional third light, considering that the long-cycle light was already removed with the process of disentangling. However, the very good fit obtained (based on minimization) suggests that these additional light sources, if present, are much fainter than the stars and the disk/spots. In addition, while the donor irradiation by the hot spot is not included, the much larger effect of the donor irradiation by the gainer is implemented in our code.

4.2 The fit to the orbital light curve

The fit to the orbital light curve was performed using the inverse-problem solving method based on the simplex algorithm, and the model for the binary system described in the previous section. The inverse-problem method is the process of finding the set of parameters that will optimally fit the synthetic light curve to the observed data. We used the Nelder-Mead simplex algorithm (Press et al., 1992) with the optimization described by Dennis & Torczon (1991). For details see Djurašević (1992b).

Based on results of the previous sections we fixed the spectroscopic mass ratio to = 0.203 and the donor temperature to T = 9825 K. In addition, we set the gravity darkening coefficient and the albedo of the gainer and the donor to = 0.25 and A = 1.0, following von Zeipel’s law for radiative shells and complete re-radiation (von Zeipel, 1924). The limb darkening for the components was calculated as described by Djurašević et al. (2010).

We assume that the donor is rotating synchronously, i.e. the non-synchronous rotation coefficient, defined as the ratio between the actual and the Keplerian angular velocity is f = 1.0. This is justified since it is assumed that the donor has filled its Roche lobe (i.e. the filling factor of the donor was set to F = 1.0), then the synchronization is expected as consequence of the tidal forces.

The case for the gainer is different, since the accreted material from the disk is expected to transfer angular momentum increasing the rotational speed of the gainer up to the critical velocity as soon as even a small fraction of the mass has been transferred (Packet, 1981; de Mink, Pols, & Glebbeek, 2007; Deschamps et al., 2013). For this reason we assumed critical rotation for the gainer, and estimated a non-synchronous rotation factor f = 8.9 in the critical rotation regime.

The best fit along with the residuals, individual donor, disk and gainer flux contributions and the view of the optimal model at orbital phases 0.25, 0.50 and 0.75, are shown in Fig. 9. We note that the residuals show no dependence on orbital or long-cycle phases, except a larger random scatter around main eclipse. Parameters are given in Table 9 and the sensitivity of with some parameters is illustrated in Fig. 10. We find that at quadrature and -band, the gainer contributes 27% more flux than the donor and the disk only 48% of the donor to the total flux.

We find that the system contains a 13.8 M  star and a 2.81 M  star with absolute magnitude = -6.4 and -3.0 respectively, separated by 49.9 R. The stellar temperatures are T = 25460 K and T (fixed) = 9825 K. The best-fitting model contains an optically and geometrically thick disk around the hotter, more massive gainer star. With a radius of 25 R, the disk is 2.8 times larger than the central star ( 8.8 R). The disk has a convex shape, with central thickness d 1.6 R and edge thickness d 3.5 R. The temperature of the disk decreases from T = 9380 K at its edge, to T = 25460 K at the inner radius.

We notice that the hot spot has 205 angular radius and covers 12% of the disk outer rim, and it is situated at longitude = 3124, roughly between the components of the system, at the place where the gas stream falls on to the disk (Lubow & Shu, 1975). The temperature of the hot spot is approximately 18 per cent higher than the disk edge temperature, i.e. T = 11068 K. Although including the hot spot region into the model improves the fit, it cannot explain the light curve asymmetry completely. By introducing one additional bright spot, larger than the hot spot and located on the disk edge at = 1149, the fit becomes much better. This bright spot is 26% hotter than the disk at its edge, i.e. T = 11819 K and has an angular radius of 279, covering 16% of the disk outer rim.

5 Discussion

Only a few DPVs have been studied spectroscopically in detail and therefore few of them posses relatively well-determined orbital and stellar parameters; 9 Galactic DPVs and the LMC DPV OGLE05155332-6925581 are documented by Mennickent, Otero, & Kołaczkowski (2016) and recently stellar and orbital parameters were provided for V 495 Cen by Rosales Guzmán et al. (2018). Our study of OGLE-LMC-DPV-065 presented in this paper is the second spectroscopic study of an LMC DPV.

In Fig. 11 we compare OGLE-LMC-DPV-065 data with those of other DPVs and classical Algols, these later taken as reference. It is clear that DPVs are hotter and more massive than ordinary Algols, a fact already noticed in previous studies. In addition, it is clear that OGLE-LMC-DPV-065 is a comparatively massive and hot DPV, in many aspects similar to Lyrae. In Table 10 we provide a comparison between these systems based on the results of Mennickent & Djurašević (2013), although see also the recent research on the Lyrae disk by Mourard et al. (2018) for a complementary approach confirming the existence of a hot spot and obtaining roughly the same disk size but from an interferometric study.

The similarity is especially significant in inclination angle, stellar masses, surface gravities and time scale of the long-cycle length. Both systems are found in a mass transfer stage, harbor a comparatively hot accretion disk and massive B A type stars for the DPV standard (Fig. 9). The radial extension of the disk is also similar along with the location of the hot and bright spots. As a jet has been detected in Lyrae (Harmanec et al., 1996; Ak et al., 2007; Lomax & Hoffman, 2011), it is then possible that the same structure exists in OGLE-LMC-DPV-065 and could be related to the long-cycle through a magnetic dynamo as suggested by Schleicher & Mennickent (2017). On the other hand, an important difference is the remarkable long-cycle change observed in OGLE-LMC-DPV-065 which is not observed in Lyrae. The large amplitude of the long-cycle in OGLE-LMC-DPV-065 is also remarkable. In comparison, the long-cycle in Lyrae is of low amplitude and relatively constant in period. Orbital period changes can be explained in terms of conservative mass transfer in a binary system. Hence it is possible that both systems are in different stages of the mass transfer episode. A much larger mass transfer in Lyrae might explain why this binary shows a variable orbital period, whereas OGLE-LMC-DPV-065, eventually with a smaller mass transfer rate, does not. In addition, Lyrae has a larger and brighter secondary star, which might also play a role in the observed differences between both systems. These issues will be investigated in a forthcoming paper.

If a magnetic dynamo is the cause for the long-cycle, then these two systems with similar parameters but different long cycle light curve morphology, constitute constrains to be satisfied by any competent detailed physical model of the long variability. Our next study will explore this point, establishing the evolutionary stage of OGLE-LMC-DPV-065 and analyzing the spectroscopic changes during the long cycle. We will also present numerical calculations aimed to test the hypothesis of variable mass transfer driven by a magnetic dynamo as proposed by Schleicher & Mennickent (2017).

6 Conclusions

We have analyzed the variability of the eclipsing Algol OGLE-LMC-DPV-065 considering new and published photometric data spanning 124 years. The orbital and long-cycle light curves have been disentangled and characterized. We also presented the first spectroscopic study of this binary system obtaining the mass ratio and temperature of the cooler stellar component. These quantities served as fixed input parameters in our model of the light curve, that was done following an inverse-problem methodology. The best solution shows a reasonable fit to the light curve providing additional parameters for the binary, the stellar components and the circumprimary accretion disk. The main results of our research can be summarized as follows:

  • We find a refined orbital period of 100316267 00000056 without any evidence of variability.

  • Small but significant changes in timings of the secondary eclipse are detected. They might be caused by circumstellar material.

  • The long-cycle is characterized by a double hump light-curve at and bands, of amplitude about 0.3 and 0.2 mag, respectively, whose general shape is more or less constant, with only minor variability.

  • We find that after a continuous decrease of the long-period during about 13 years, from 350 to 218 days, it remained almost constant by about 10 years.

  • The study of radial velocities indicates a binary in a circular orbit with mass ratio of 0.203 0.001.

  • We find that the system consists of a pair of stars of 13.8 and 2.81 M  of radii 8.8 and 12.6 R and absolute bolometric magnitudes -6.4 and -3.0, respectively.

  • We find stellar temperatures of 25460 K and 9825 K for the gainer and donor, respectively.

  • We find an orbital semi-major axis of 49.9 R.

  • We find evidence of an accretion disk with a radius of 25 R, central thickness 1.6 R and edge thickness 3.5 R. The temperature of the disk decreases from 25460 K at the inner radius to 9380 K at its outer edge.

  • As happens in other DPVs, two hot shock regions located at roughly opposite parts of the outer disk rim can explain the light curves asymmetries.

  • OGLE-LMC-DPV-065 resembles in some aspects to the well-studied binary Lyrae. However, its orbital period does not change, this could indicate a smaller mass transfer rate.

FIXED PARAMETERS: - mass ratio of the components, - temperature of the less-massive (cooler) donor, - filling factor for the critical Roche lobe of the donor, - non-synchronous rotation coefficients of the gainer and donor respectively, - gravity-darkening coefficients of the gainer and donor, - albedo coefficients of the gainer, donor and disk.

Quantities: - number of observations, - final sum of squares of residuals between observed (LCO) and synthetic (LCC) light-curves, - root-mean-square of the residuals, - orbit inclination (in arc degrees), - disk dimension factor (ratio of the disk radius to the critical Roche lobe radius along y-axis), - disk-edge temperature, , , - disk thicknesses (at the edge and at the center of the disk, respectively) in the units of the distance between the components, - disk temperature distribution coefficient, - filling factor for the critical Roche lobe of the hotter, more-massive gainer (ratio of the stellar polar radius to the critical Roche lobe radius along z-axis for a star in critical rotation regime), - temperature of the more-massive (hotter) gainerr, - hot and bright spots’ temperature coefficients, and - spots’ angular dimensions and longitudes (in arc degrees), - angle between the line perpendicular to the local disk edge surface and the direction of the hot-spot maximum radiation, - non-synchronous rotation coefficients of the gainer in critical rotation regime, - dimensionless surface potentials of the hotter gainer and cooler donor, , - stellar masses and mean radii of stars in solar units, - logarithm (base 10) of the system components effective gravity, - absolute stellar bolometric magnitudes, , , , - orbital semi-major axis, disk radius and disk thicknesses at its edge and center, respectively, given in the solar radius units.

Table 9: Results of the analysis of DPV065 light-curves obtained by solving the inverse problem for the Roche model with an large accretion disk partially obscuring the more-massive (hotter) gainer in critical non-synchronous rotation regime.
Figure 9: Observed (LCO), synthetic (LCC) light-curves and the final O-C residuals between the observed and synthetic light curves of OGLE-LMC-DPV-065; fluxes of donor, gainer and of the accretion disk, normalized to the donor flux at phase 0.25; the views of the model at orbital phases 0.18, 0.55, 0.70 and 0.98, obtained with parameters estimated by the light curve analysis. This is for a gainer in critical non-synchronous rotation regime.
Figure 10: Plots showing the dependence of .
Figure 11: Comparison of physical data for semi-detached Algols from Dervişoǧlu, Tout, & Ibanoǧlu (2010, primaries open blue circles and secondaries open red circles) and DPVs (Mennickent, Otero, & Kołaczkowski, 2016, (primaries blue crosses and secondaries red squares)). The zero-age main sequence for Z = 0.02 is plotted with a solid black line and evolutionary tracks for single stars with initial masses (in solar masses) labelled at the track footprints are also shown (Pols et al., 1998). The positions of the gainer and donor of OGLE-LMC-DPV-065 are indicated by color arrows, whereas Lyrae is labeled with black arrows. The best evolutionary tracks for the primary and secondary of the DPV HD 170582 are also plotted by two solid lines in the left-hand panel (Mennickent, Otero, & Kołaczkowski, 2016).


We thanks the referee, Denis Mourard, who helped to improve the first version of this manuscript. This paper uses photometric data acquired under CNTAC proposal CN2014B-13. This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. R.E.M. gratefully acknowledges support by VRID-Enlace 218.016.004-1.0, FONDECYT 1190621, and the Chilean Centro de Excelencia en Astrofísica y Tecnologías Afines (CATA) BASAL grant AFB-170002. The OGLE project has received funding from the Polish National Science Centre grant MAESTRO no. 2014/14/A/ST9/00121. G. D. gratefully acknowledges the financial support of the Ministry of Education and Science of the Republic of Serbia through the project 176004, Stellar physics. J.G.F.-T. is supported by FONDECYT N. 3180210. N.A-D. acknowledges support from FONDECYT #3180063. We thanks Shelby Owens for reducing data of the Turitea Observatory.

Table 10: Comparison between the OGLE-LMC-DPV-065 data obtained in this paper and those of Lyrae obtained by Mennickent & Djurašević (2013) and references therein.


  • Applegate (1992) Applegate J. H., 1992, ApJ, 385, 621
  • Ak et al. (2007) Ak H., et al., 2007, A&A, 463, 233
  • Barría et al. (2013) Barría D., Mennickent R. E., Schmidtobreick L., Djurašević G., Kołaczkowski Z., Michalska G., Vučković M., Niemczura E., 2013, A&A, 552, A63
  • Bolton (1989) Bolton C. T., 1989, SSRv, 50, 311
  • Cabezas et al. (2019) Cabezas, M., Hadrava, P., Mennickent, R. E., Rivinius, T., 2019, KOREL disentangling of the eclipsing Algol OGLE-LMC-DPV-065, Contrib. Obs. Skalnaté Pleso 49, in press.
  • de Mink, Pols, & Glebbeek (2007) de Mink S. E., Pols O. R., Glebbeek E., 2007, AIPC, 948, 321
  • Dennis and Torczon (1991) Dennis J.E., Torczon V., 1991, Direct search methods on parallel machines, in SIAM Journal on Optimization, Vol. 1, Issue 4, pages 448-474
  • Dervişoǧlu, Tout, & Ibanoǧlu (2010) Dervişoǧlu A., Tout C. A., Ibanoǧlu C., 2010, MNRAS, 406, 1071
  • Deschamps et al. (2013) Deschamps R., Siess L., Davis P. J., Jorissen A., 2013, A&A, 557, A40
  • Djurašević (1992a) Djurašević G., 1992a, Ap&SS, 196, 267
  • Djurašević (1992b) Djurašević G., 1992b, Ap&SS, 197, 17
  • Djurašević (1996) Djurašević G., 1996, Ap&SS, 240, 317
  • Djurašević et al. (2010) Djurašević G., Latković O., Vince I., Cséki A., 2010, MNRAS, 409, 329
  • Djurašević et al. (2012) Djurašević G., Vince I., Antokhin I. I., Shatsky N. I., Cséki A., Zakirov M., Eshankulova M., 2012, MNRAS, 420, 3081
  • Garcés L. et al. (2018) Garcés L J., Mennickent R. E., Djurašević G., Poleski R., Soszyński I., 2018, MNRAS, 477, L11
  • Garrido et al. (2013) Garrido H. E., Mennickent R. E., Djurašević G., Kołaczkowski Z., Niemczura E., Mennekens N., 2013, MNRAS, 428, 1594
  • Guinan (1989) Guinan E. F., 1989, SSRv, 50, 35
  • Hadrava (1990) Hadrava P., 1990, Contributions of the Astronomical Observatory Skalnate Pleso, 20, 23
  • Hadrava (1995) Hadrava P., 1995, A&AS, 114, 393
  • Hadrava (1997) Hadrava P., 1997, A&AS, 122, 581
  • Hadrava (2004) Hadrava P., 2004, PAICz, 92, 15
  • Hadrava (2016) Hadrava P., 2016, Astrophysics and Space Science Library, Vol. 439, p. 113  
  • Hall (1989) Hall D. S., 1989, SSRv, 50, 219
  • Harmanec et al. (1996) Harmanec P., et al., 1996, A&A, 312, 879
  • Harmanec (2002) Harmanec P., 2002, AN, 323, 87
  • Harmanec et al. (1996) Harmanec P., et al., 1996, A&A, 312, 879
  • Heemskerk (1994) Heemskerk M. H. M., 1994, A&A, 288, 807
  • Lanza, Rodonò, & Rosner (1998) Lanza A. F., Rodonò M., Rosner R., 1998, MNRAS, 296, 893
  • Lanza & Rodonò (1999) Lanza A. F., Rodonò M., 1999, A&A, 349, 887
  • Lomax & Hoffman (2011) Lomax J. R., Hoffman J. L., 2011, BSRSL, 80, 689
  • Lubow & Shu (1975) Lubow S. H., Shu F. H., 1975, ApJ, 198, 383
  • Meintjes (2004) Meintjes P. J., 2004, MNRAS, 352, 416
  • Mennickent (2017) Mennickent R. E., 2017, SerAJ, 194, 1
  • Mennickent & Djurašević (2013) Mennickent R. E., Djurašević G., 2013, MNRAS, 432, 799
  • Mennickent et al. (2003) Mennickent R. E., Pietrzyński G., Diaz M., Gieren W., 2003, A&A, 399, L47
  • Mennickent et al. (2005) Mennickent R. E., Cidale L., Díaz M., Pietrzyński G., Gieren W., Sabogal B., 2005, MNRAS, 357, 1219
  • Mennickent et al. (2012) Mennickent R. E., Djurašević G., Kołaczkowski Z., Michalska G., 2012, MNRAS, 421, 862
  • Mennickent, Otero, & Kołaczkowski (2016) Mennickent R. E., Otero S., Kołaczkowski Z., 2016, MNRAS, 455, 1728
  • Mennickent, Schleicher, & San Martin-Perez (2018) Mennickent R. E., Schleicher D. R. G., San Martin-Perez R., 2018, PASP, 130, 094203
  • Mourard et al. (2018) Mourard D., et al., 2018, A&A, 618, A112
  • Packet (1981) Packet W., 1981, A&A, 102, 17
  • Pawlak et al. (2013) Pawlak M., et al., 2013, AcA, 63, 323
  • Poleski et al. (2010) Poleski R., Soszyński I., Udalski A., Szymański M. K., Kubiak M., Pietrzyński G., Wyrzykowski Ł., Ulaczyk K., 2010, AcA, 60, 179
  • Pols et al. (1998) Pols O. R., Schröder K.-P., Hurley J. R., Tout C. A., Eggleton P. P., 1998, MNRAS, 298, 525
  • Press et al. (1992) Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 1992, Numerical Recipes in Fortran, second ed., in: The Art of Scientific Computing, vol. 120, Cambridge University Press
  • Richards (2004) Richards M. T., 2004, AN, 325, 229
  • Rosales Guzmán et al. (2018) Rosales Guzmán J. A., Mennickent R. E., Djurašević G., Araya I., Curé M., 2018, MNRAS, 476, 3039
  • Schleicher & Mennickent (2017) Schleicher D. R. G., Mennickent R. E., 2017, A&A, 602, A109
  • Szymanski (2005) Szymanski M. K. 2005, AcA, 55, 43
  • Udalski, Szymański, & Szymański (2015) Udalski A., Szymański M. K. & Szymański G.2015, AcA, 65, 1
  • von Zeipel (1924) von Zeipel H., 1924, MNRAS, 84, 702
  • Zola (1991) Zola S., 1991, AcA, 41, 213
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description