next previous
Up: V 700 Cygni

4. Photometric solution

The light curve analysis of both systems is quite difficult for the following reasons: (a) no spectroscopic mass-ratio is known, (b) the light curves show significant anomalies with O'Connell effect, and (c) the systems exhibit partial eclipses.

The light curves of V 700 Cyg reveal brightness fluctuations mainly around the maxima. The O'Connell effect is observed in both light curves and the magnitude difference between the two maxima is MaxII - MaxI = tex2html_wrap_inline1716 on the average for both light curves. A small deficiency of light is noted in the phase interval 0.13 - 0.24, mainly in the B band. The difference between the two minima is tex2html_wrap_inline1722 in B and tex2html_wrap_inline1726 in V.

The light curves of AW Vir show brightness disturbances not only around the maxima, but also at other phases. A decrease in brightness is present in both B and V light curves in the phase interval 0.60 - 0.78 and a small excess of light is seen in both curves in the phase interval 0.39 - 0.46. Other minor light variations are present in other phase intervals. A small O'Connell effect is observed in both light curves. The magnitude difference between the two maxima is MaxII - MaxI = tex2html_wrap_inline1740 on the average for both light curves. The difference between the two minima is about tex2html_wrap_inline1742 in B and tex2html_wrap_inline1746 in V.

With such anomalies in the light curves, the analysis of both systems is quite difficult. We tried to model the light curves by invoking spots on one or both components. The need to place cool (sunspot-like) and/or hot spots on the components of contact binaries to explain the light curve anomalies and the O'Connell effect, has been suggested by several investigators (e.g. Binnendijk 1960; Hilditch 1981; Linnell 1982; Van Hamme & Wilson 1985; Milone et al. 1987; van't Veer & Maceroni 1988, 1989; Maceroni et al. 1990).

4.1. Unspotted solution

The light curve analysis was carried out by using the most recent (1993) version of the Wilson-Devinney (Wilson 1990) synthetic light curve code, which has the capacity of automatically adjusting the spots. In order to reduce computational time and to smooth the scatter around the maxima, normal points were formed from the individual observations and assigned weights equal to the number of observations per normal. The normal points of V 700 Cyg are given in Table 3 and those of AW Vir in Table 4 (both tables are available from cdsarc.u-strasbg.fr). The mean standard deviation for normal points in B and V are tex2html_wrap_inline1766 and tex2html_wrap_inline1768 for V 700 Cyg, and tex2html_wrap_inline1770 and tex2html_wrap_inline1772 for AW Vir, respectively. Care has been taken to ensure a faithful resemblance of the normal points to the actual shape of the minima.

We proceeded to the unspotted solution by assuming that there are no spots on the components of the systems. Therefore the unperturbed parts of the light curves were used. In accordance with the light curve anomalies mentioned before, the observations in the phase intervals 0.13 - 0.24 and 0.61 - 0.82 were not included in the case of V 700 Cyg, while those in the interval 0.60 - 0.78 were omitted in the case of AW Vir.

In the following, the subscripts h and c refer to the hotter and cooler component, respectively. The shape of the light curve minima shows partial eclipses, so that we do not know a priori whether the systems belong to A or W type. In both cases the deeper primary minimum indicates that the hotter star is eclipsed at primary minimum. Therefore, the subscripts 1 and 2 in the DC program (with phase of conjunction tex2html_wrap_inline1780) are identical with h and c, respectively. A preliminary set of input parameters for the DC program was obtained by the Binary Maker 2.0 program (Bradstreet 1993). The DC program was used in mode 3. In the subsequent analysis the following assumptions were made:

Case V 700 Cyg: A mean surface temperature tex2html_wrap_inline1786 according to the spectral type G5V; bolometric albedos tex2html_wrap_inline1788 and gravity darkening coefficients tex2html_wrap_inline1790 were assigned values typical for stars with convective envelopes; limb darkening coefficients tex2html_wrap_inline1792 in B and tex2html_wrap_inline1796 in V were taken from Al-Naimiy's (1978) tables; bolometric linear limb darkening coefficients tex2html_wrap_inline1800 were taken from Van Hamme (1993); third light was assumed to be tex2html_wrap_inline1802.

Case AW Vir: A mean surface temperature tex2html_wrap_inline1804 according to the spectral type F8V; bolometric albedos and gravity darkening coefficients were assigned the same values as for V 700 Cyg; limb darkening coefficients tex2html_wrap_inline1806 in B and tex2html_wrap_inline1810 in V were taken from Al-Naimiy's tables; bolometric linear limb darkening coefficients tex2html_wrap_inline1814 were taken again from Van Hamme (1993); third light was assumed to be tex2html_wrap_inline1816.

The adjustable parameters were in both cases: the phase of conjunction tex2html_wrap_inline1818, the inclination i, the temperature tex2html_wrap_inline1822, the non-dimensional potential tex2html_wrap_inline1824, the monochromatic luminosity tex2html_wrap_inline1826 and the mass-ratio tex2html_wrap_inline1828. The quantity tex2html_wrap_inline1830 was adjusted only in the first few iterations, since it showed no tendency to vary significantly. The lack of a spectroscopic mass-ratio for both systems led us to search for the solution with several fixed values for the mass-ratio q in the range 0.2 - 4. The values of q<1 correspond to a transit at the primary (A-type system) and those of q>1 to an occultation (W-type system). The lowest values of the sum tex2html_wrap_inline1840 of the weighted squares of the residuals occured at q=0.8 and q=1.6 for the case of V 700 Cyg and at q=0.8 and q=1.2 for AW Vir, with almost equal values in both cases for each system. Figures 3 (click here) and 4 (click here) show the fit parameter tex2html_wrap_inline1850 as a function of the mass-ratio q. The range in q for V 700 Cyg is 0.4 - 2.8, since for values of q outside this range no convergent acceptable solution could be obtained.

  figure305
Figure 3: V 700 Cyg: The fit parameter tex2html_wrap_inline1860 as a function of the mass-ratio q

  figure310
Figure 4: AW Vir: The fit parameter tex2html_wrap_inline1864 as a function of the mass-ratio q

In order to find the final unspotted solution we continued the analysis by applying the DC program for the above values of q, treating q as a free parameter. The solutions for V 700 Cyg converged to q=0.8895 and q=1.5699 with values of tex2html_wrap_inline1876 0.0330 and 0.0319, respectively. In the case of AW Vir the solutions converged to q=0.8006 and q=1.2222 and the corresponding values of tex2html_wrap_inline1882 were found to be 0.0210 and 0.0202. Of these two solutions (A and W-type for each system), we finally adopted the W-type solutions by taking into account (a) the better fit of the W-type solution and (b) the large mass-ratio of the A-type solution, which would be unusual for an A-type system. In all DC solutions, the method of subsets (Wilson & Biermann 1976) was used because of the strong correlation between the adjustable parameters. The parameters i and q were always taken in separate subsets, since in many cases q and i can be combined in such a way that the quality of the solution becomes good for a large range of mass-ratios. The results of the unspotted solution for both systems are given in Table 5 (click here). In Figs. 5 (click here) and 6 (click here), the corresponding theoretical light curves are shown as dashed lines.

  figure319
Figure 5: Normal points and theoretical B and V light curves of V 700 Cyg. Dashed lines: unspotted solution; Solid lines: spotted solution

  figure324
Figure 6: Normal points and theoretical B and V light curves of AW Vir. Dashed lines: unspotted solution; Solid lines: spotted solution

4.2. Spotted solution

 

V 700 Cyg V 700 Cyg AW Vir AW Vir AW Vir
Parameter unspotted spotted unspotted spotted unspotted solution
solution solution solution solution (Lapasset et al. 1996)
tex2html_wrap_inline1916 tex2html_wrap_inline1918 tex2html_wrap_inline1920 tex2html_wrap_inline1922 tex2html_wrap_inline1924
i (degrees) tex2html_wrap_inline1928 tex2html_wrap_inline1930 tex2html_wrap_inline1932 tex2html_wrap_inline1934 tex2html_wrap_inline1936
tex2html_wrap_inline1938 0.tex2html_wrap_inline1940 0.tex2html_wrap_inline1942 0.tex2html_wrap_inline1944 0.tex2html_wrap_inline1946 0.tex2html_wrap_inline1948
tex2html_wrap_inline1950 (K) tex2html_wrap_inline1952 tex2html_wrap_inline1954 tex2html_wrap_inline1956 tex2html_wrap_inline1958 tex2html_wrap_inline1960
tex2html_wrap_inline1962 (K) tex2html_wrap_inline1964 tex2html_wrap_inline1966 tex2html_wrap_inline1968 tex2html_wrap_inline1970 tex2html_wrap_inline1972
tex2html_wrap_inline1974 tex2html_wrap_inline1976 0.tex2html_wrap_inline1978 0.tex2html_wrap_inline1980 0.tex2html_wrap_inline1982 0.tex2html_wrap_inline1984 0.tex2html_wrap_inline1986
tex2html_wrap_inline1988 tex2html_wrap_inline1990 tex2html_wrap_inline1992 tex2html_wrap_inline1994 tex2html_wrap_inline1996 tex2html_wrap_inline1998 tex2html_wrap_inline2000
tex2html_wrap_inline2002 tex2html_wrap_inline2004tex2html_wrap_inline2006 tex2html_wrap_inline2008 tex2html_wrap_inline2010 1.48
tex2html_wrap_inline2012 (B) tex2html_wrap_inline2016 tex2html_wrap_inline2018 tex2html_wrap_inline2020 tex2html_wrap_inline2022 0.585
tex2html_wrap_inline2024 (V) tex2html_wrap_inline2028 tex2html_wrap_inline2030 tex2html_wrap_inline2032 tex2html_wrap_inline2034 0.584
tex2html_wrap_inline2036 tex2html_wrap_inline2038 (B) 0.tex2html_wrap_inline2042 0.tex2html_wrap_inline2044 0.tex2html_wrap_inline2046 0.tex2html_wrap_inline2048 0.tex2html_wrap_inline2050
tex2html_wrap_inline2052 tex2html_wrap_inline2054 (V) 0.tex2html_wrap_inline2058 0.tex2html_wrap_inline2060 0.tex2html_wrap_inline2062 0.tex2html_wrap_inline2064 0.tex2html_wrap_inline2066
tex2html_wrap_inline2068 tex2html_wrap_inline2070 (bolo) 0.tex2html_wrap_inline2072 0.tex2html_wrap_inline2074 0.tex2html_wrap_inline2076 0.tex2html_wrap_inline2078
% overcontact 26% 27% 9% 8% 2.1%
tex2html_wrap_inline2080 (pole) tex2html_wrap_inline2082 tex2html_wrap_inline2084 tex2html_wrap_inline2086 tex2html_wrap_inline2088 0.332
tex2html_wrap_inline2090 (side) tex2html_wrap_inline2092 tex2html_wrap_inline2094 tex2html_wrap_inline2096tex2html_wrap_inline2098 0.348
tex2html_wrap_inline2100 (back) tex2html_wrap_inline2102 tex2html_wrap_inline2104 tex2html_wrap_inline2106 tex2html_wrap_inline2108 0.383
tex2html_wrap_inline2110 (pole) tex2html_wrap_inline2112 tex2html_wrap_inline2114 tex2html_wrap_inline2116 tex2html_wrap_inline2118 0.391
tex2html_wrap_inline2120 (side) tex2html_wrap_inline2122 tex2html_wrap_inline2124 tex2html_wrap_inline2126tex2html_wrap_inline2128 0.414
tex2html_wrap_inline2130 (back) tex2html_wrap_inline2132 tex2html_wrap_inline2134 tex2html_wrap_inline2136 tex2html_wrap_inline2138 0.446
tex2html_wrap_inline2140 0.0319 0.0338 0.0202 0.0157
tex2html_wrap_inline2142 0.60 0.84
tex2html_wrap_inline2144 0.92 1.11
tex2html_wrap_inline2146 0.86 0.95
tex2html_wrap_inline2148 1.04 1.08
tex2html_wrap_inline2150 -0.13 0.08
tex2html_wrap_inline2154 -0.08 0.10
tex2html_wrap_inline2158assumed.
Table 5: Light curve solutions of V 700 Cyg and AW Vir

 

Inspection of Figs. 5 (click here) and 6 (click here) reveals that the theoretical ``unspotted'' light curves do not fit satisfactorily the observations (normal points). The disagreement for V 700 Cyg is relatively small in the phase regions before MaxI and around MaxII, while the disagreement for AW Vir is severe around MaxII (O'Connell effect). In order to take care of the decrease of brightness in the above phase intervals in both systems, a spotted solution was carried out by adopting the simplest spot model with a physical meaning, i.e. by assuming that the systems have cool spots on the primary component of the same nature as solar magnetic spots (Mullan 1975). Two cool spots were placed on the primary (larger, more massive and cooler) component of V 700 Cyg and one cool spot on the primary (larger, more massive and cooler) component of AW Vir. The Binary Maker 2.0 program was used to obtain the best fit by adjusting the spot parameters: the latitude b, the longitude l, the angular radius R and the temperature factor tex2html_wrap_inline2166

Once the best fit was obtained, the DC program was used to derive the final solution. At this stage of the solution the DC program (1993 version) was used to adjust some of the spot parameters (longitude, angular radius and temperature factor) together with the system parameters i, tex2html_wrap_inline2170, tex2html_wrap_inline2172, q and tex2html_wrap_inline2176. Because of convergence problems with a simultaneous adjustment of all parameters, we employed the method of subsets in the sense: system parameters, spot parameters, system parameters etc. (Van Hamme, private communication). The differential corrections were computed until the corrections became smaller than their probable errors. A few more iterations were always performed after obtaining a convergent solution to ensure its stability. The final results of the spotted solution of V 700 Cyg and AW Vir are given in Table 5 (click here), and the corresponding theoretical light curves are shown as solid lines in Figs. 5 (click here) and 6 (click here), respectively. The final spot parameters for the two systems are given in Table 6 (click here). The final spotted solution for V 700 Cyg (where the spots were adjusted) places the two spots on the primary in positions where both are visible at phase 0.0 and hence should introduce perturbations in that phase region. These perturbations are quite small, of the order of the statistical error lines (see Sect. 5 for a discussion), because the temperature factor of the spots are very close to unity. The tex2html_wrap_inline2178 differences between the observed and calculated points for the unspotted and spotted solutions for the two systems are shown in Figs. 7 (click here) and 8 (click here), respectively.

  figure746
Figure 7: The light curve tex2html_wrap_inline2180 residuals for V 700 Cyg in B and V band. Crosses refer to unspotted solution; asterisks refer to spotted solution

  figure751
Figure 8: The light curve tex2html_wrap_inline2186 residuals for AW Vir in B and V band. Crosses refer to unspotted solution; asterisks refer to spotted solution

 

V 700 Cyg V 700 Cyg AW Vir
parameter spot 1 spot 2 spot 1
b (degrees) 90 90 80
l (degrees) 241.77tex2html_wrap_inline2192 6.41 123.20tex2html_wrap_inline2194 4.58 309.51tex2html_wrap_inline21964.98
R (degrees) 20.16tex2html_wrap_inline2192 0.84 19.60tex2html_wrap_inline2194 1.07 14.56tex2html_wrap_inline21960.32
tex2html_wrap_inline2210 0.945tex2html_wrap_inline2192 0.005 0.952tex2html_wrap_inline2194 0.006 0.720tex2html_wrap_inline2196 0.028
Table 6: Spot parameters

 


next previous
Up: V 700 Cygni

Copyright by the European Southern Observatory (ESO)
web@ed-phys.fr