We present a byproduct of our long term photometric monitoring of cataclysmic variables. 2MASS J18024395 +4003309 = VSX J180243.9 +400331 was discovered in the field of the intermediate polar V1323 Her observed using the Korean 1m telescope located at Mt. Lemmon, USA. An analysis of the twocolor VR CCD observations of this variable covers all the phase intervals for the first time. The light curves show this object can be classified as an Algoltype variable with tidally distorted components, and an asymmetry of the maxima (the O’Connell effect). The periodogram analysis confirms the cycle numbering of
Andronov et al. (2012)
and for the initial approximation, the ephemeris is used as follows: Min I. BJD = 2456074.4904+0.3348837
E
. For phenomenological modeling, we used the trigonometric polynomial approximation of statistically optimal degree, and a recent method “NAV” (“New Algol Variable”) using local specific shapes for the eclipse. Methodological aspects and estimates of the physical parameters based on analysis of phenomenological parameters are presented. As results of our phenomenological model, we obtained for the inclination
i
=90°,
M
_{1}
=0.745
M
_{◉}
,
M
_{2}
=0.854
M
_{◉}
,
M
=
M
+
M
_{2}
=1.599
M
_{◉}
, the orbital separation
a
=1.65·10
^{9}
m=2.37
R
_{◉}
and relative radii
r
_{1}
=
R
_{1}
/
a
=0.314 and
r
_{2}
=
R
_{2}
/
a=0.360. These estimates may be used as preliminary starting va
lues for further modeling using extended physical models based on the
Wilson & Devinney (1971)
code and it's extensions
1. INTRODUCTION
Chungbuk National University Observatory (CBNUO) is monitoring many cataclysmic variables as a part of an “InterLongitude Astronomy (ILA)” campaign (
Andronov et al. 2010
) in order to study how the physical properties of cataclysmic variables depend on time and luminosity state. CBNUO is also involved in developing an automatic observation system and an analysis program for monitoring the cataclysmic variables (
Yoon et al. 2012
,
2013
). We have published some results from our monitoring data (
Andronov et al. 2011
;
Kim et al. 2004
,
2005
). 1RXS J180340.0+401214 is an intermediate polar, subclass of magnetic cataclysmic variables, which has a magnetic white dwarf accreting from its secondary. This object recently got an official name V1323 Her in the “General Catalogue of Variable Stars” (GCVS) (
Samus et al. 2014
). One eclipsing variable, 2MASS J18024395+4003309, was discovered in the vicinity of this intermediate polar, 1RXS J1803, by
Breus (2012)
. This object was registered in the “Variable Stars index” (VSX,
http://aavso.org/vsx
) and received the name VSX J180243.9 + 400331 (hereafter called VSX1802 for brevity) and an object identifier 282837. Unfortunately, no GCVS name has yet been given to this object. The study of the main target V1323 Her was presented by
Andronov et al. (2011)
.
Due to a relatively large angular distance (14’) from V1323 Her, VSX1802 is seen in the field of V1323 Her only in CCD images with a focal reducer. This was the case for one night reported by
Andronov et al. (2012)
, who noted a sharp profile of the minimum. Additional observations from the Catalina survey (
Drake et al. 2009
) allowed them to determine photometric elements
T
_{0}
= 2456074.4904,
P
=0.3348837±0.0000002
^{d}
. This profile is typical characteristic of systems with nearly equal radii and an inclination of
i
≈90° (i.e., both eclipses are nearly total). Due to incomplete coverage of the phases in the discovery paper by
Andronov et al. (2012)
, more detailed study was needed for VSX1802. Such work can serve as one of the good byproducts of the ILA campaign.
The purpose of this study is to analyse the photometry data with a complete coverage of the phase as well as to carry out phenomenological modeling. The phenomenological model is based on the theory of close binary systems e.g. presented in the classical monographs by
Kopal (1959)
and
Tsessevich (1971)
. Because no detailed photometric or spectroscopic studies have been reported yet, our rough phenomenological model will help to study this object with more detail in the future.
2. OBSERVATIONS
The Korean 1m telescope at Mt. Lemmon in Arizona, USA (LOAO), is equipped with a focal reducer for 2×2 K CCD. The fieldofview is 22.2 square arcminutes, which is very effective for studies not only of the main targets (typically intermediate polars), but also of other variable stars in the field. It should also be noted that another eclipsing binary GSC 0437000206 (now called V442 Cam,
Samus et al. 2014
) has been discovered in the field of MU Cam = 1RXS J062518.2+733433 with this telescope (
Kim et al. 2005
).
In total, we obtained 196 observations in V (range 16.51
^{m}
–17.51
^{m}
) and 242 observations in R (range 15.88
^{m}
–16.77
^{m}
) between 2012 and 2014. The total duration of observations was 45.5 hours during 11 nights in R and 8 nights in the alternatively changing filters VR. The time interval of the observations was JD 2455998–2456722. To improve calibration accuracy, we have used the method of “artificial comparison star” (
Andronov & Baklanov 2004
;
Kim et al. 2004
). As the object is in the field of V1323 Her, we have used the star C1 of
Andronov et al. (2011)
, for which
Henden (2005)
published the magnitudes of V and R filters: V=14.807
^{m}
, Rc=14.436
^{m}
. The original observations (HJD, magnitude) are available upon request.
Fig. 1
is the finding chart of V1323 Her with VSX1802.
Finding chart of V1323 Her (marked with v) and reference stars (marked by numbers). Reference star 22 is VSX 1802. The size of the field is 22’x22’.
3. PERIODOGRAM ANALYSIS
For the periodogram analysis, we have used the trigonometric polynomial fit of a degree s:
where coefficients
C_{α}
,
α
=1..
m,m
=1+2
s
are computed using the least squares method (cf.
Anderson 2003
;
Andronov 1994
,
2003
), and
f
is frequency (in cycles/day) . For fitting, we used the statistical test function
S
(
f
) as
where
σ
_{O}
is the r.m.s. deviation of the observations (O) from the sample mean, C corresponds to calculated values and OC to the deviation of the observed values from the calculated ones (see
Andronov 1994
,
2003
for more details). The periodograms are shown in
Fig. 2
for the filter R and
s
=1 and
s
=2. One can see that the numerous peaks for
s
=1 are much lower than those for
s
=2, indicating that the main signal corresponds not to the true frequency, but to its harmonic, which is characteristic for eclipsing binaries with nearly equal minima (as was shown by
Andronov et al. (2012)
from the photometry from the Catalina survey). The most prominent peak at the periodogram corresponds to
f
=2.986107 cycles/day. A close value of the best frequency,
f
= 2.986114 cycles/day is also seen for the observations in the filter V. The frequencies in the filter V and R are close within error estimates to the value
f
= 2.986111±0.000002 cycles/day (
Andronov et al. 2012
), but differ by a value of 0.5/year from the estimate based on 3 minima by
Parimucha et al. (2012)
. Since our periodogram analysis confirms the cycle numbering of
Andronov et al. (2012)
, for the initial approximation, we use their ephemeris:
4. PHENOMENOLOGICAL MODELING. MULTI – HARMONIC APPROXIMATION
The degree of the trigonometric polynomial
s
is often determined by visual comparison of the light curve with its approximation (cf.
Parenago & Kukarkin 1936
). To determine the statistically optimal degree s of approximation (1), one may use different criteria (see
Andronov 1994
,
2003
for detailed discussion). The classical approach is to use Fischer’s criterion. In
Fig. 2
, the value
L_{s}
= –lg FAP is plotted vs
s
. It is clearly seen that
L
is large for even values of
s
=2
k
, and
L
≈0 for
s
=2
k
+1;
k
=1, 2,.. This can be explained by the good symmetry of the light curve. Adopting the critical value
L_{crit}
=3 (i.e. the false alarm probability FAP=10
^{−3}
), one may suggest
s
=14 for the filter V and
s
=8 for R.
Dependence of the parameter L_{s}= – lg FAP on the degree of the trigonometric polynomial s. The horizontal dashed line corresponds to a critical value L=3.
Another criterion proposed by
Andronov (1994)
is based on the r.m.s. estimate of the accuracy of the smoothing curve
σ
[
x_{c}
] at the moments of the observations. The corresponding dependence is shown in
Fig. 3
. For both filters, the minimum of
σ
[
x_{c}
] corresponds to
s
=6. In this approximation, the frequency
f
is a free parameter determined using differential corrections. The best fit estimate of periods is
P
= 0.3348842
^{d}
±0.0000005
^{d}
(V) , 0.3348845
^{d}
± 0.0000004
^{d}
(R), the difference from the value in (1) is not statistically significant. The corresponding moments of the primary minima are 2456238.9186 ± 0.0005
^{d}
(V), 24 56320.9663±0.0007
^{d}
(R). The difference between these values is due to different cycle numbers close to the mean time of observations in V and R (as required for best accuracy estimates, see
Andronov 1994
for details).
Dependence of the r.m.s. value of the accuracy of the smoothing function σ[x_{c}] on the degree of the trigonometric polynomial s.
The phase curves in the V and R filters for statistically optimal degrees of trigonometric polynomial fit
s
= 6 (best accuracy of the smoothing function) and
s
=14 (Fischer’s criterion) are shown in
Fig. 4
. Comparison of the trigonometric polynomial fits with different degrees
s
shows that, for smaller
s
, the depth of the eclipses is underestimated and there are formal waves at the “outofeclipse” part. These waves are apparently present for larger s, where the accuracy of the smoothing function is larger, but the central part of the eclipses is approximated better.
Phase curves in the V and R filters for statistically optimal degrees of the trigonometric polynomial fit s=6 (best accuracy of the smoothing function) and s=14 (Fischer’s criterion). Filled circles are original observations, the smoothing functions x_{c}(φ) are shown with 1σ and 2σ error corridors.
5. PHENOMENOLOGICAL MODELING. THE NAV ALGORITHM
Phenomenological modeling of the light curves of eclipsing binary stars with relatively narrow eclipses was described in detail by
Andronov (2012)
. The method was called “NAV” (”New Algol Variable”) and applied to a few Algoltype variables (e.g.
Kim et al. 2010b
).
In short, the method may be described as follows. The smoothing function is defined as usual in the linear least squares method
where
C_{α}
are called “coefficients” and
f_{α}
(
t
)  “basic functions”. Contrary to trigonometric polynomials with sines and cosines only, the “New Algol Variable” (NAV) algorithm combines a low – order trigonometric polynomial with a special shape for the eclipses.
Andronov (2012)
and
Tkachenko & Andronov
(2014) compared a few approaches and chose a local shape dependent on a single parameter called
β
:
For
β
=0, the shape is narrow and is physically unrealistic; for
β
=1, the shape at the center of eclipse is triangular; for
β
=2, it is parabolic and when β approaches ∞, the shape tends to be a rectangle. The dimensionless variable
z
is related to phase
ϕ
as
where phase is defined typically for a given initial epoch
T
_{0}
and period
P
:
A discussion of the computation of phase for the case of variable period may be found in
Andronov & Chinarova (2013)
.
Another free parameter is the eclipse halfwidth
Δϕ
=
D
/2, where
D
– the full width of the eclipse, in the GCVS (
Samus et al. 2014
) is expressed in percent of the period.
Other shapes used for determining of the parameters are based on a Gaussian function and its modifications (
Mikulasek et al. 2012
), which are formally of infinite width. An opposite approach is a splitting of the phase interval and approximation of the “outofeclipse” parts by a constant and of the eclipses by a parabola (
Papageorgiou et al. 2014
). The disadvantages of this model are: a) the discontinuity of the smoothing function, b) the underestimation of the depth of the minima because its width is set to a large constant and c) bad fitting of the ascending and descending branches. However, the number of parameters is only 5 (as the width of the eclipse and the phase shifts are set to constants).
Our previous approach was to use 5 parameters for the “constant+parabola” fit (brightness out of eclipse, at the primary and secondary minimum; halfwidth
Δϕ
and shift
ϕ
_{0}
), in order to avoid discontinuity of the smoothing function. Other types of functions were discussed by
Andronov (2005)
and
Andronov & Marsakova (2006)
.
In the NAV algorithm, the basic functions are:

f1(1)=1

f2(ϕ)=cos(φ),φ=2πϕ

f3(ϕ)=cos(2φ)

f4(ϕ)=sin(φ)

f5(ϕ)=sin(2φ)

f6(ϕ)=H((ϕϕ0)/C8,C9)

f7(ϕ)=H((ϕϕ00.5)/C8,C10)
The coefficient
C
_{1}
corresponds to mean stellar magnitude after reducing the observations for the eclipse, effects of reflection, ellipticity and asymmetry. Although the effects of reflection, ellipticity and asymmetry are not strictly sinusoidal, their amplitudes are typically much smaller than those of the eclipses, thus, we use firstorder approximation;
C
_{2}
corresponds to semiamplitude of the reflection effect;
C
_{3}
corresponds to the semiamplitude of the ellipticity effect; the asymmetry (O'Connel effect, usually interpreted by spots) is approximated by terms with
C
_{4}
and
C
_{5}
. The depths of the primary and secondary minimum are described by the coefficients
C
_{6}
and
C
_{7}
, respectively.
To continue numeration of the coefficients, we introduce
C
_{8}
=
Δϕ
,
C
_{9}
=
β
_{1}
(describing the profile of the primary minimum) and
C
_{10}
=
β
_{2}
for the secondary minimum. The phase shift
C
_{11}
=
ϕ
_{0}
may also be added to the set of unknown variables, but we have used a fixed value
ϕ
_{0}
= 0 for the analysis of the present star. Typically the period and the initial epoch determined from longterm observations from different authors are more accurate than those from the smaller set of observations used for the analysis.
For each set of trial values of
C
_{8}
,
C
_{9}
, C
_{10}
, the first 7 coefficients were computed using the method of least squares (
Anderson 2003
;
Press et al 2007
). For the rest of the coefficients, one may choose among various methods (
Cherepashchuk 1993
;
Andronov & Tkachenko 2013a
;
Marquardt 1963
). For our analysis, we used the method of minimizing the test function at a grid.
Usually, the test function for one filter may be written in this way:
where
x_{k}
,
t_{k}
are signal values and times, respectively, and
w_{k}
are possible weight coefficients. The complete theory of statistical properties of the smoothing function in a case of arbitrary functions and (waveletlike) time and scaledependent weight functions was presented by
Andronov (1997)
.
We have computed a dependence
Φ
_{V}
(
C
_{8}
)=min
_{C8,C9}
Φ
_{V}
(
C
_{8}
,
C
_{9}
,
C
_{10}
) for the filter V and, similarly, R. The minimum of
Φ
(
C
_{8}
) are slightly different due to statistical errors. To get a single value, as theoretically expected, we made a scaled sum
The adopted value of the filter halfwidth
C
_{8}
=0.1177, corresponds to the minimum of
Φ
(
C
_{8}
). Other parameters determined by the phenomenological modeling are listed in
Table 1
.
Coefficients of the phenomenological model for our observations in the filters V, R and that for Catalina (the latter published byAndronov et al. 2012).
Coefficients of the phenomenological model for our observations in the filters V, R and that for Catalina (the latter published by Andronov et al. 2012).
The corresponding light curves for the filters V and R are shown in
Fig. 5
.
The “NAV” fit for the phase curve and a corresponding “1σ” corridor, for the filters V (up) and R (bottom). At the phases of eclipses, an additional curve corresponding to a continuation of the “out of eclipse” parts is also shown.
Using the smoothing functions for the two filters, the smoothing function of the color index VR was computed. It is shown in
Fig. 6
.
The smoothing curve for the color index using the “NAV” fits for the phase curve in the filters V and R and a corresponding “1σ” corridor.
6. DISCUSSION
From these coefficients of the “NAV” approximation, the classical phenomenological parameters listed in the GCVS (
Samus et al. 2014
) were determined: Max I=16.567
^{m}
±0.006
^{m}
, Max II =16.592
^{m}
±0.006
^{m}
, Min I = 17.493
^{m}
± 0.014
^{m}
, Min II =17.281
^{m}
±0.008
^{m}
, Min IMax I = 0.926
^{m}
, Min IIMax I=0.714
^{m}
(filter V). The corresponding point at the “depth – depth” diagram (
Fig.1
in
Malkov et al. 2007
) lies close to the line “R” of equal radii, but slightly outside the allowed region. This is caused by the “simplified” model of spherical components without effects of ellipticity, reflection and limb darkening, which
Malkov et al. (2007)
have used.
The parameters
C
_{6}
, and
C
_{7}
correspond to the "reduced depth" of the minimum, i.e., not to the maximum (as in the "General Catalogue of the Variable Stars"), with respect to the continuation of the "outofeclipse" curve, as shown in
Fig. 4
. If one uses them as “depths” of the minima, the corresponding point at the diagram by
Malkov et al. (2007)
is located in the physically reliable region, relatively close to the line of the equal radii.
The mean color index, (VR)=16.674
^{m}
16.105
^{m}
=0.569
^{m}
± 0.069
^{m}
was computed from the values of
C
_{1}
for two filters. In the phase curve of the color index (
Fig. 5
), variations from 0.54
^{m}
± 0.01
^{m}
(at Max I at phase 0.25) to 0.63
^{m}
± 0.03
^{m}
(at Min I at phase 0.00) are present. At Min II at phase 0.50, VR = 0.61
^{m}
± 0.01
^{m}
. The asymmetry of the maxima indicates the presence of the O’Connell effect.
These values of the color index are in the instrumental system, so, before taking into acco unt the color transformation coefficients and determining (VR) in the standard system, we can't estimate more precisely the "color temperature" of the "star1+star2 system" (i.e. weighted mean temperature of two components) and corresponding "spectral class" (also intermediate between that of the components).
Following
Andronov (2012)
, we introduce residual relative intensities
The intensities are in units of sum of intensities of two stars (out of eclipse, but "reduced" for the three effects mentioned above) and the “intensity depth” of minimum (how much light is eclipsed):
d
_{1}
=1I
_{1}
,
d
_{2}
=1I
_{2}
.
Andronov (2012)
introduced a phenomenological parameter defined as
Y
=
d
_{1}
+
d
_{2}
. The case
Y
= 0 corresponds to "no eclipse", 1 to "both eclipses are total", and, for the majority of cases, 0 ≤
Y
≤ 1. For this star, we get values of
In the model of "spherical stars with no limb darkening" (
Shulberg 1971
; Malkov et al. 2005) recently discussed by
Andronov & Tkachenko (2013b
) as a "zeroorder" approximation for determination of physical parameters), the eclipsed part of the light is proportional to the surface of projection at the eclipse
S
; the surface of projection of stars
S
_{1}
=
πr
_{1}
^{2}
,
S
_{2}
=
πr
_{2}
^{2}
Here, we use dimensionless parameters,
r
_{1}
=
R
_{1}
/
a
,
r
_{2}
=
R
_{2}
/
a
, where a is the orbital separation.
Let’s introduce
L
_{1V}
=1
L
_{2V}
as the relative contribution of the emission from the first star in V, similar to other filters. Let’s define that the first star eclipses the second star at the primary (more deep) minimum. Under this definition, the second star has larger surface brightness (temperature).
So one may write a system of 4 equations
with 5 unknowns,
L
_{1V}
,
L
_{1R}
,
S
,
S
_{1}
,
S
_{2}
. Four parameters – the relative surface brightness
F
_{2R}
=
L
_{2R}
/
S
_{2}
and similarly for the star 1 and filter V – are related to that mentioned above. So it is not possible to determine all 5 parameters without additional information (e.g. color index luminosity (or absolute magnitude)radius relation, assuming Main Sequence)  thus we need an "extra" equation.
However, we can estimate
From this, we can derive another combination
i.e., equal to unity within error estimates. A small difference (even if smaller error), could be related to limb darkening.
Other combinations are for surface brightness
These phenomenological values may be used to estimate temperature ratios, e.g. using "Main Sequence" (MS) relations". From these MS assumptions, one may make a sequence of suggestions, which will be discussed below.
The second star (eclipsed at a primary minimum) has a larger surface brightness, consequently a larger radius and luminosity than the first star.
Taking into account Y≈1, one may suggest that the secondary eclipse is total (or close to total). Thus,
S
=
S
_{1}
, and the system of equations may be solved completely:
These values slightly differ by a value of 0.061±0.035=1.7
σ
, and the difference is not statistically significant. The mean weighted value of the ratio of crosssections of the two stars is (
S
_{1}
/
S
_{2}
) = 0.761 ± 0.017, and the ratio of radii is (
R
_{1}
/
R
_{2}
) = 0.872 ± 0.010 (radius is in units of the orbital separation). For the lowmass part of the Main Sequence,
R~M
(e.g.
Faulkner 1971
) so this value may be an estimate of the mass ratio,
q
=
M
_{1}
/
M
_{2}
= 0.872 ± 0.010. From the table of
Allen (1973)
for the Main Sequence one may obtain a statistical relation
R
≈
M
^{0.72}
for a wide range of spectral classes: from A0 to M0. In this case, we get a slightly different estimate
q
=
M
_{1}
/
M
_{2}
≈ (
R
_{1}
/
R
_{2}
)
^{1/0.72}
= 0.826 ± 0.013 . From more recent tables of
Cox (2000)
,
R
≈1.01
M
^{0.80}
, thus
However, the (unknown) deviations from the massradius relation may exceed this maximal 5 percent difference in estimates using different coefficients of the statistical dependence.
From the duration of eclipse, one may suggest an inequality
For such large values, the stars are distorted (and we see this also from large values of
C
_{3}
) and thus close to their Roche lobes.
Although a possible W UMa – type classification may not be completely rejected, these stars are not yet in thermal contact (
Lucy 1976
) based on the difference of mean brightness and we suggest that they are not in contact. Thus, we perfer an EA  type classification with elliptic component(s).
Using the oversimplified form of the MS massradius relation
and, combining with Kepler's third law, we obtain the inclination,
i
=90°,
M
_{1}
= 0.745
M
_{◉}
,
M
_{2}
= 0.854
M
_{◉}
,
M
=
M
_{1}
+
M
_{2}
=1.599
M
_{◉}
the orbital separation
a
=1.65∙10
^{9}
m = 2.37
R
_{◉}
and relative radii
r
_{1}
=
R
_{1}
/
a
= 0.314 and
r
_{2}
=
R
_{2}
/
a
=0.360.
Generally, these are minimal estimates of the radiuses and masses, since for an inclination,
i
, differing from 90°, we’ll get larger estimated values. With the assumption of one total eclipse, one may write cos(
i
)<
r
_{2}

r
_{1}
 = 0.046, thus 87.4° ≤
i
≤ 90°, 0.999 ≤ sin(
i
) ≤ 1, and for VSX 1802, the effects of inclination on estimates of radii and masses are negligible.
Even though the spectral classes of this system have not been reported yet, we can suggest the corresponding spectral classes of these stars to be G8 and K2 according to
Cox (2000)
and the possible Roche model of these VSX 1802 system can be estimated as in
Fig. 7
.
The model of the system VSX 1802: The Roche lobes, the line of centers and circles corresponding to estimated radii of the stars in a sphericallysymmetric approximation.
7. CONCLUSIONS
In this study, we estimated the paramaters only from phenomenological modeling of the phase light curve. The main unjustified suggestions were "no limb darkening" and "similar to MS dependencies" (no real quantitative relations, just "larger stars are brighter"). There are no detailed photometric or spectroscopic studies at present, therefore the true physical parameters of this system cannot be estimated for now.
Our phenomenological modeling could allow rough estimates to make the physical parameters. Although the errors of the phenomenological parameters are not very large (up to few percent), there may be systematic errors up to a dozen percent due to the simplicity of the model. These estimates may be used as preliminary values for further modeling using extended physical models based on the
Wilson & Devinney (1971)
code and it’s extensions (
Wilson 1979
,
1994
,
2012
,
2014
;
Zola et al. 1997
,
2010
;
Bradstreet 2005
;
Kallrath et al, 2009
;
Linnel 2012
;
Reed 2012
;
Rucinski 2010
;
Prsa et al 2012
). Such study of eclipsing binaries with the CBNU observations are reported by
Jeong & Kim (2013)
;
Kim & Jeong(2012)
;
Jeong & Kim (2011)
and
Kim et al. (2010a)
. It should be noted that the lack of spectral data doesn’t allow to study various combination of stars in different evolutionary stages. Some discussions with crude assumptions about the age and metallicity of the system are therefore out of the scope of this study. More detailed photometrical and spectroscopic information will enable more thorough studies.
Acknowledgements
This work was supported by the research grant of Chungbuk National University in 2012. The data acquisition and analysis were partially supported by Basic Science Research Program through the National Research Foundation of Korea (NRF) funded by the Ministry of Education, Science and Technology (20110014954). It is also a part of the “InterLongitude Astronomy” campaign (Andronov et al. 2010) and the project “Ukrainian Virtual Observatory (Vavilova et al. 2012).
Allen CW
1973
Astrophysical quantities
3rd Ed
Athlone Press
London
Anderson TW
2003
An Introduction to Multivariate Statistical Analysis
3rd Ed
Wiley
New York
721 
Andronov IL
1994
(Multi) Frequency Variations of Stars. Some Methods and Results
Odessa Astronomical Publications
Odessa
49 
54
Andronov IL
1997
Method of running parabolae: Spectral and statistical properties of the smoothing function
Astron. Astrophys. Suppl.
http://dx.doi.org/10.1051/aas:1997217
125
207 
217
Andronov IL
,
Sterken C
2003
Multiperiodic versus Noise Variations: Mathematical Methods
Astronomical Society of the Pacific
San Francisco
ASP Conference Series
391 
400
Andronov IL
2003
Multiperiodic versus Noise Variations: Mathematical Methods
Astronomical Society of the Pacific
San Francisco
ASP Conference Series
391 
400
Andronov IL
2012
Phenomenological Modeling of the Light Curves of AlgolType Eclipsing Binary Stars
Astrophys
http://dx.doi.org/10.1007/s1051101292590
55
536 
550
Andronov IL
,
Kim YK
,
Yoon JN
,
Breus VV
,
SmeckerHane TA
2011
TwoColor CCD Photometry of the Intermediate Polar 1RXS J180340.0+401214
J. Korean Astron. Soc.
http://dx.doi.org/10.5303/JKAS.2011.44.3.89
44
89 
96
Andronov IL
,
Antoniuk KA
,
Baklanov AV
,
Breus VV
,
Burwitz V
2010
InterLongitude Astronomy (ILA) Project: Current Highlights and Perspectives.I.Magnetic vs.NonMagnetic Interacting Binary Stars
Odessa Astronomical Publications
Odessa
8 
12
Andronov IL
,
Baklanov AV
2004
Algorithm of the artificial comparison star for the CCD Photometry
Astron School Rep.
5
264 
272
Andronov IL
,
Breus VV
,
Zola S
2012
Determination of Characteristics of Newly Discovered Eclipsing Binary 2MASS J18024395 +4003309 VSX J180243.9+400331
Odessa Astronomical Publications
Odessa
145 
147
Andronov IL
,
Chinarova LL
,
Wszolek B
2013
Method of Running Sines: Modeling Variability in LongPeriod Variables, in Czestochowski Kalendarz Astronomiczny 2014
Astronomia Nova Association
Poland
139 
144
Andronov IL
,
Marsakova VI
2006
Variability of longperiod pulsating stars. I. Methods for analyzing observations
Astrophys
http://dx.doi.org/10.1007/s1051100600378
49
370 
385
Andronov IL
,
Kim Y
,
Yoon JN
,
Breus VV
,
SmeckerHane TA
2011
TwoColor CCD Photometry of the Intermediate Ploar 1RXS J180340.0+401214
J. Korean Astron. Soc.
http://dx.doi.org/10.5303/JKAS.2011.44.3.89
44
89 
96
Andronov IL
,
Tkachenko MG
,
Wszolek B
2013
Comparative Analysis of Numerical Methods for Parameter Determination, Czestochowski Kalendarz Astronomiczny 2014
Astronomia Nova Association
Poland
173 
180
Andronov IL
,
Tkachenko MG
2013
Comparative Analysis of Numerical Methods of Determination of Parameters of Binary Stars. Case of Spherical Components
Odessa Astronomical Publications
Odessa
203 
206
Bradstreet DH
2005
Fundamentals of Solving Eclipsing Binary Light Curves Using Binary Maker 3
Society for Astronomical Sciences
California
The Society for Astronomical Sciences 24th Annual Symposium on Telescope Science
23 
37
Breus VV
2MASS J18024395+4003309
[Internet], reported 2012 May 26, available from:
Cherepashchuk AM
1993
Parametric Models in Inverse Problems of Astrophysics
Astronomicheskii Zhurnal
70
1157 
1176
Cox AN
2000
Allen’s Astrophysical Quantities
AIP Press
New York
Drake AJ
,
Djorgovski SG
,
Mahabal A
,
Beshore E
,
Larson S
2009
First Results from the Catalina RealTime Transient Survey
Astrophys. J.
http://dx.doi.org/10.1088/0004637X/696/1/870
696
870 
884
Faulkner J
1971
UltrashortPeriod Binaries, Gravitational Radiation, and Mass Transfer. I. The Standard Model, with Applications to WZ Sagittae and Z Camelopardalis
Astrophys. J.
http://dx.doi.org/10.1086/180848
170
L99 
L104
Henden A
AAVSO public datasets
uploaded 2005 Nov 8, available form:
Jeong JH
,
Kim CH
2011
WZ Cephei: A Dynamically Active W UMaType Binary Star
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2011.28.3.163
28
163 
172
Jeong MJ
,
Kim CH
2013
The First Comprehensive Photometric Study of the Neglected Binary System V345 Cassiopeiae
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2013.30.4.213
30
213 
221
Kallrath J
,
Milone EF
2009
Eclipsing Binary Stars: Modeling and Analysis
Springer
New York
444 
Kim Y
,
Andronov IL
,
Jeon YB
2004
CCD Photometry Using Multiple Comparison Stars
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2004.21.3.191
21
191 
200
Kim Y
,
Andronov IL
,
Park SS
,
Chinarova LL
,
Baklanov AV
2005
TwoColor VR CCD Photometry of the Intermediate Polar 1RXS J062518.2+733433
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2005.22.3.197
22
197 
210
Kim CH
,
Kim HI
,
Yoon TS
,
Han W
,
Lee JW
2010
SW LyncisAdvances and Questions
J. Astron Space Sci.
http://dx.doi.org/10.5140/JASS.2010.27.4.263
27
263 
273
Kim CH
,
Lee JW
,
Duck HK
,
Andronov IL
2010
Four New Binary Stars in the Field of CL Aurigae. II
Odessa Astronomical Publications
Odessa
62 
64
Kim CH
,
Jeong JH
2012
V700 Cygni: A Dynamically Active W UMatype Binary Star II
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2012.29.2.151
29
151 
161
Kopal Z
1959
Close Binary Systems
Chapman & Hall
London
558 
Linnel AP
2012
The BINSYN Program Package
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2012.29.2.123
29
123 
129
Lucy LB
1976
W Ursae Majoris systems with marginal contact
Astrophys. J.
http://dx.doi.org/10.1086/154265
205
208 
216
Malkov Yu
,
Oblak E
,
Avvakumova EA
,
Torra J
2007
A procedure for the classification of eclipsing binaries
Astron. Astrophys.
http://dx.doi.org/10.1051/00046361:20066103
465
549 
556
Marquardt DA
1963
Method for the Solution of Certain Problems in Least Squares
SIAM J. Appl. Math.
http://dx.doi.org/10.1137/0111030
11
431 
441
Mikulasek Z
,
Zejda M
,
Janik J
2012
Period Analyses Without OC Diagrams
Cambridge University Press
Cambridge
Proceedings IAU Symposium
http://dx.doi.org/10.1017/S1743921311027888
391 
394
Parimucha S
,
Dubovsky P
,
Janík J
,
Kudzej I
,
Solovyova I
2012
Photometric analysis of 8 newly discovered shortperiod eclipsing binaries at Astronomical Observatory at Kolonica Saddle
New Astron
http://dx.doi.org/10.1016/j.newast.2011.06.007
17
93 
100
Papageorgiou A
,
Kleftogiannis G
,
Christopoulou PE
2014
An Automated Search of O'Connell Effect from Surveys of Eclipsing Binaries
Contrib. Astron. Obs. Skalnate Pleso
43
470 
472
Parenago PP
,
Kukarkin BV
1936
The Shapes of Light Curves of Long Period Cepheids
Mit 12 Abbildungen, Zeitschrift fur Astrophysik
11
337 
355
Reed PA
2012
Photometric Observations of Eccentric Accretion in Algoltype Binary Stars
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2012.29.2.141
29
141 
143
Press WH
,
Teukolsky SA
,
Vetterling WT
,
Flannery BP
2007
Numerical Recipes: The Art of Scientific Computing
Cambridge University Press
Cambridge
1193 
Prsa A
,
Guinan EF
,
Devinney EJ
,
Degroote P
,
Bloemen S
2012
Advances in Modeling Eclipsing Binary Stars in the Era of Large AllSky Surveys with EBAI and PHOEBE
Cambridge University Press
Cambridge
Proceeding IAU Symposium
271 
278
Rucinski S
2010
Contact Binaries: The Current State
AIP Conf. Proc.
1314
29 
36
Samus NN
,
Durlevich OV
,
Kazarovets EV
,
Kireeva NN
,
Pastukhova EN
General Catalog of Variable Stars
[Internet], available from:
Shulberg AM
1971
Close Binary Systems with Spherical Components
Nauka
Moscow
246 
Tkachenko MG
,
Andronov IL
Determination of Characteristics of Eclipsing Binaries with Spots: Phenomenological vs Physical Models
Conference of RAGtime
Prague, Czech Republic
Oct 1119, 2014
Tsessevich VP
1971
Eclipsing Variable Stars
Nauka
Moscow
350 
Vavilova IB
,
Pakulyak LK
,
Shlyapnikov AA
,
Protsyuk YI
,
Savanevich VE
2012
Astroinformation Resource of the Ukrainian Virtual Observatory: Joint Observational Data Archive, Scientific Tasks, and Software
Kinem. Phys. Celest. Bodies
http://dx.doi.org/10.3103/S0884591312020067
28
85 
102
Wilson RE
1979
Eccentric Orbit Generalization and Simultaneous Solution of Binary Star Light and Velocity Curves
Astrophys. J.
http://dx.doi.org/10.1086/157588
234
1054 
1066
Wilson RE
1994
Binarystar light Curve Models
Publ. Astron. Soc. Pac.
http://dx.doi.org/10.1086/133464
106
921 
941
Wilson RE
2012
UnderDeveloped and UnderUtilized Eclipsing Binary Model Capabilities
J.Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2012.29.2.115
29
115 
121
Wilson RE
2014
Improving CMD Areal Density Analysis:Algorithms and Strategies Capabilities
J. Astron. Space Sci
http://dx.doi.org/10.5140/JASS.2014.31.2.121
31
121 
130
Wilson RE
,
Devinney EJ
1971
Realization of Accurate CloseBinary Light Curves: Application to MR Cygni
Astrophys. J.
http://dx.doi.org/10.1086/150986
166
605 
619
Yoon JN
,
Lee YS
,
Kim CH
,
Kim Y
,
Yim HS
2012
Statistical Analysis on Weather onditions at Chungbuk National University Observatory in Jincheon, Korea
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2012.29.4.397
29
397 
405
Yoon JN
,
Kim Y
,
Kim DH
,
Yim HS
2013
Development of an AutoFlat program for the acquisition of effective flat images in the automated observation system
J. Astron. Space Sci.
http://dx.doi.org/10.5140/JASS.2013.30.4.327
30
327 
334
Zola S
,
Kolonko M
,
Szczech M
1997
Analysis of a Photoelectric Light Curve of the W UMaType Binary ST Ind
Astron. Astrophys.
324
1010 
1012
Zola S
,
Gazeas K
,
Kreiner JM
,
Ogloza W
,
Siwak M
2010
Physical Parameters of Components in Close Binary Systems  VII
Mon. Not. R. Aston. Soc.
http://dx.doi.org/10.1111/j.13652966.2010.17129.x
408
464 
474