The Photometric Investigation of V921 Her Using the Lunar-Based Ultraviolet Telescope of Chang’e-3 Mission
The light curve of V921 Her in ultraviolet band observed by the Lunar-based Ultraviolet Telescope (LUT) is analyzed by the Wilson-Devinney code. Our solutions conclude that V921 Her is an early type marginal contact binary system with an additional close-in component. The binary system is under poor thermal contact with a temperature difference of nearly 700 K between the two components. The close-in component contributes about 19% of the total luminosity in the triple system. Combining the radial velocity study together with our photometric solutions, the mass of the primary star and secondary one is calculated to be , . The evolutionary scenario of V921 Her is discussed. All times of light minimum of V921 Her available in the bibliography are taken into account and the curve is analyzed for the first time. The most probable fitting results are discussed in the paper, which also confirm the existence of a third component ( year) around the binary system. The period of V921 Her is also undergoing a continuously rapid increase at a rate of , which may be due to mass transfer from the less massive component to the more massive one.
In December 2013, the Chinese Lunar Exploration Program carried out the Chang’e-3 mission , which sent a LUT to the surface of the moon. LUT is the first robotic astronomical telescope deployed on the moon surface. It has been working for nearly two years since landing on the moon. Eclipsing binaries are among its observational targets. We have obtained some light curves (LCs) to be analyzed, and V921 Her is a such one.
The light variability of V921 Her (HIP 82344, ) was first discovered by the Hipparcos mission  and the Hipparcos Catalogue gave its spectral type to be A5. Later, Rucinski et al.  carried out the radial velocity study of V921 Her. The results determined that V921 Her was an A-subtype contact binary system with its spectral type to be A7IV. They obtained its mass ratio to be (±0.005). Compared with other W UMa type close binary systems of the same spectral type, V921 Her has a quite long period, which is d. The first photometric analysis of V921 Her was done by Gazeas et al. . They redetermined the mass ratio of the binary system to be and set the effective temperature of the primary star to be . Their results concluded that V921 Her was a shallow contact binary system. However, Karami and Mohebi  and Karami et al.  also obtained the mass ratio of V921 Her, both of which were consistent with those determined by Rucinski et al. .
Orbital period research is a very important part for close binary systems since the period variations contain information about the dynamic interaction occurring between the components of binary systems. However, the period change investigation of V921 Her has been neglected since it was discovered. In the present work, light curve of V921 Her observed by the LUT is analyzed with the Wilson-Devinney (W-D) code, and its period variations are investigated for the first time. The results will give us comprehensive understanding to the physical properties and dynamical evolution of this contact binary system.
2. New Photometric Observations
The LUT is a 150 mm, F/3.75 Ritchey-Chretien telescope working at a Nasmyth focus. It is equipped with a UV-enhanced back-illuminated AIMO 1 k CCD with a field of view to be ; thus the CCD pixel scale is 4.7′′ pixel−1. The passband of filter is about 245–345 nm, with its peak at 250 nm .
The light curve of V921 Her was observed continuously in ultraviolet band by the LUT on February 2, 2015. PHOT subroutine (measured magnitudes for a list of stars) of the IRAF (The Image Reduction and Analysis Facility is hosted by the National Optical Astronomy Observatories in Tucson, Arizona, at URL http://iraf.noao.edu/.) aperture photometry package was used to reduce the observed images .
The phase of our observations is calculated with the following linear ephemeris:
3. Investigation of the Orbital Period
The orbital period variations of a binary star reveal dynamic interactions between its components or the gravitational interactions between the binary system and the close-in component orbiting around it. During this work, all times of light minimum of V921 Her available in the bibliography are taken into account and listed in Table 2. All of the minimum data are observed with CCD camera. The (observed times of light minimum − calculated times of light minimum) values of all light minima are calculated with the linear ephemeris below:The Epoch and corresponding values are listed in the third and fourth columns of Table 2, respectively. The origin in time used in (2) is the first minimum light available in the bibliography for V921 Her.
As shown in Figure 1(a), the linear ephemeris in (2) cannot explain the variations efficiently. We consider long term increase (or decrease) and periodic variations superposed on the linear ephemeris; then where and are corrections to the initial ephemeris and period and is the periodic element given by Irwin  in the case of . Based on the least-squares method, the most probable results of curve fitting are given out.
The resulting ephemeris isWith the quadratic term included in this ephemeris, a continuous period increase, at a rate of , is determined. The sinusoidal term reveals a cyclic period change with a period of 10.2 years and an amplitude of 0.00516 days. The residuals from (3) are showed in Figure 1(c).
It has to be mentioned that the data point of is not used since it is apparently inconsistent with the following two data points. The large group of minima are between and cycles; the rest consists of only 6 very small groups of minima widely dispersed over time. To test the evolution of the diagram, more and more data are necessary in the near future.
The orbital period of binary system is a very important parameter. In the present paper, comparing the period published by other researchers and databases, we chose the initial period to be , since there is only a very small difference among them. The revised period in (4) is 0.8773763(1), which confirms that the period we used is acceptable.
4. Analysis of the Light Curves
To understand the geometrical structure and physical properties of V921 Her, the observed light curve is analyzed using the W-D code (Version 2013) [22–24]. According to the spectral type (A7IV) and mass ratio () derived by Rucinski et al. , the effective temperature of star 1 (star eclipsed at primary minimum light) is set to be , and the mass ratio is fixed to . Radiative and convective outer envelopes are assumed for star 1 and star 2, respectively. Thus, the bolometric albedo ,  and the values of the gravity-darkening coefficients ,  are used. To account for the limb darkening in detail, logarithmic functions are used. The corresponding bolometric and passband-specific limb-darkening coefficients are chosen from Van Hamme’s table .
During the W-D modeling process, the adjustable parameters are as follows: the orbital inclination ; the mean surface temperature of star 2 (); the monochromatic luminosity of star 1 (); the dimensionless potential of star 1 () and star 2 (); and the third light . It is found the solutions converge at both Mode 3 (contact configuration) and Mode 5 (semidetached configuration). The final photometric solutions are listed in Table 3. The theoretical light curves of Mode 3 are displayed in Figure 2 and the theoretical light curves which have not been contaminated by the third light are also plotted with dashed lines. The resulting figure from Mode 5 is very similar to Figure 2. The geometrical structure at phase 0.25 is displayed in Figure 3.
5. Discussion and Conclusions
The light curve solutions of V921 Her show that it converges at both Mode 3 and Mode 5. The solutions of Mode 3 determines a contact configuration with both of its components just filling their critical Roch lobe, and the solutions of Mode 5 determine a semidetached binary system with the secondary component filling the critical Roch lobe and the primary one filling of the critical Roch lobe. However, as displayed in Table 3, the solutions between Mode 3 and Mode 5 have only very small differences in parameters such as orbital inclination ; temperature of star 2 (), potential of star 1 () and star 2 (), luminosity ratio of , and the contribution of the third light . Thus, we can conclude that V921 Her may be a marginal contact or B-type binary system . The temperature difference between its two components is about , which indicates that V921 Her is a poor thermal contact binary system. It is an important target which is exactly in the rare phase predicted by the thermal relaxation oscillation (TRO) theory [30–33].
Compared with the solutions obtained by Gazeas et al. , we determined a much higher orbital inclination and our temperature for star 2 is about lower than theirs, which may be due to the different mass ratio () we used. And also, we add a third light in our solutions, which may also account for these differences. As we conclude that V921 Her is a marginal contact binary system and being on the TRO stage, it is not so strange that a contact configuration was obtained by Gazeas et al. .
Considering the orbital inclination () of ours and the mass function given by Rucinski et al. , , we can easily calculate the mass of the two components to be , . And the orbital semimajor axis is obtained to be . The absolute physical parameters of the two components in V921 Her are listed in Table 4. The absolute physical parameters are calculated according to the solutions of Mode 3. It is almost the same with the values calculated by using the solutions of Mode 5 since they give out nearly equal surface temperature, orbital inclination, and so on.
As discussed by Liao and Qian , the most plausible explanation of the cyclic period changes in close binaries is the light-travel time effect (LTTE) caused by the presence of the tertiary component around the binary system. By assuming a circular orbit (), the projected radius of the orbit that the eclipsing binary rotates around the barycenter of the triple system is calculated with the following equation,where is the amplitude of the oscillation and is the speed of light; that is, . The mass function is computed with the following equation:where and are the gravitational constant and the period of the oscillation. Thus, we obtain the mass function to be . The relationship between the orbital inclination () and the mass () of the tertiary component is plotted in Figure 4(a). The relationship between the orbital inclination () and orbital radius () of the tertiary component is also plotted in Figure 4(b).
D’Angelo et al.  and Pribulla and Ruciński  pointed out that most overcontact binaries existed in multiple systems. A survey of 165 solar-type spectroscopic binaries concluded that the fraction of spectroscopic binaries with additional companions reached for d . As for V921 Her, both light curve solutions and fitting results confirm the existence of a close-in tertiary component orbiting around the binary system. As shown in Figure 2, the third light significantly reduces the depth of the eclipse occultation. According to our results, the luminosity of the third component () is nearly twice the luminosity of the secondary star (). It is a stellar component which contributes about of the total luminosity in the triple system. It is supposed that the presence of a third component has played an important role in the formation and evolution of V921 Her by removing angular momentum from the central binary system during the early dynamical interaction .
The analysis also shows that the orbital period of V921 Her is increasing at a rate of . By considering a conservative mass transfer and combining the estimated mass of V921 Her with the well-known equationthe mass transfer from the less massive component to the more massive one, at a rate of , is determined.
V921 Her is under the key evolutionary stage predicted by the thermal relaxation oscillation (TRO) theory, which may evolve into a contact binary system. It is a very interesting target for testing theoretical scenario of formation and evolution on W UMa type binaries. More photometric and spectroscopic observations are still needed in the future. This object will be monitored for a long time.
The authors declare that they have no competing interests.
The authors appreciate the outstanding work of the LUT team and the support by the team from the ground system of the Chang’e-3 mission. This work is supported by the Chinese Natural Science Foundation (Grant nos. 11133007, 11325315, and 11203066), the Strategic Priority Research Program “The Emergence of Cosmological Structure” of Chinese Academy of Sciences (Grant no. XDB09010202), the Key Research Program of Chinese Academy of Sience (KGED-EW-603), and the Science Foundation of Yunnan Province (Grant no. 2012HC011).
ESA, The Hipparcos and Tycho Catalogues, vol. 1200 of ESA Special Publication, 1997.
K. D. Gazeas, P. G. Niarchos, S. Zola, J. M. Kreiner, and S. M. Rucinski, “Physical parameters of components in close binary systems: VI,” Acta Astronomica, vol. 56, no. 1, pp. 127–143, 2006.View at: Google Scholar
K. Karami, K. Ghaderi, R. Mohebi, R. Sadeghi, and M. M. Soltanzadeh, “Velocity curve analysis of the spectroscopic binary stars RZ Cas, CC Cas, HS Her, HD 93917, V921 Her and Y Cygni by the artificial neural networks,” Astronomische Nachrichten, vol. 330, no. 8, pp. 836–842, 2009.View at: Publisher Site | Google Scholar
K. K. Kwee and H. van Woerden, “A method for computing accurately the epoch of minimum of an eclipsing variable,” Bulletin of the Astronomical Institutes of the Netherlands, vol. 12, p. 327, 1956.View at: Google Scholar
S. Marshall, C. Akerlof, R. Kehoe et al., “The ROTSE project,” Bulletin of the American Astronomical Society, vol. 29, p. 1290, 1997.View at: Google Scholar
T. Pribulla, D. Baludansky, D. Chochol et al., “New minima of selected eclipsing close binaries,” Information Bulletin on Variable Stars, no. 5898, p. 1, 2005.View at: Google Scholar
L. Brát, M. Zejda, and P. Svoboda, “B.R.N.O. Contributions #34,” Open European Journal on Variable Stars, vol. 74, no. 1, p. 1, 2007.View at: Google Scholar
O. Pejcha, “CCD times of minima of several eclipsing binaries,” Information Bulletin on Variable Stars, vol. 5645, p. 1, 2005.View at: Google Scholar
R. H. Nelson, CCD Minima for Selected Eclipsing Binaries in 2004, Information Bulletin on Variable Stars no. 5602, 2005.
T. Krajci, “Photoelectric minima of some eclipsing binary stars,” Information Bulletin on Variable Stars, no. 5806, pp. 1–3, 2007.View at: Google Scholar
S. W. Dvorak, “Times of minima for neglected eclipsing binaries in 2004,” Information Bulletin on Variable Stars, no. 5603, p. 1, 2005.View at: Google Scholar
R. H. Nelson, “CCD minima for selected eclipsing binaries in 2006,” Information Bulletin on Variable Stars, 5760, #1, 2007.View at: Google Scholar
S. W. Dvorak, “Times of minima for eclipsing binaries 2010,” Information Bulletin on Variable Stars, no. 5974, p. 1, 2011.View at: Google Scholar
R. H. Nelson, “V456 CYG—a detached eclipsing binary,” Information Bulletin on Variable Stars, vol. 5994, pp. 1–4, 2011.View at: Google Scholar
L. Brát, J. Trnka, L. Smelcer et al., “B.R.N.O. contributions #37—times of minima,” Open European Journal on Variable Stars, vol. 137, p. 1, 2011.View at: Google Scholar
A. N. Cox, Allen's Astrophysical Quantities, Springer, New York, NY, USA, 4th edition, 2000.
S. M. Ruciński, “The proximity effects in close binary systems. II. The bolometric reflection effect for stars with deep convective envelopes,” Acta Astronomica, vol. 19, p. 245, 1969.View at: Google Scholar
L. B. Lucy, “Gravity-darkening for stars with convective envelopes,” Zeitschrift für Astrophysik, vol. 65, p. 89, 1967.View at: Google Scholar