## Astrokit-an Efficient Program for High-Precision Differential CCD Photometry and Search for Variable Stars

Artem Y. Burdanov ? , Vadim V. Krushinsky, and Alexander A. Popov

Ural Federal University, Ekaterinburg, 620002, Russia

Received December 24, 2013; accepted May 28, 2014

## ABSTRACT

Having a need to perform di GLYPH&lt;11&gt; erential photometry for tens of thousands stars in a several square degrees field, we developed Astrokit program. The software corrects the star brightness variations caused by variations of atmospheric transparency: to this end, the program selects for each star an individual ensemble of reference stars having similar magnitudes and positions in the frame. With ten or more reference stars in the ensemble, the di GLYPH&lt;11&gt; erences between their spectral types and the spectral type of the object studied become unimportant. Astrokit searches for variable stars using Robust Median Statistics criterion, which allows candidate variables to be selected more e GLYPH&lt;14&gt; ciently than by analyzing the standard deviation of star magnitudes. The software allows very precise automatic analysis of long inhomogeneous sets of photometric observations of a large number of objects to be performed, making it possible to find 'hot Jupiter' type exoplanet transits and low-amplitude variables. We describe the algorithm of the program and the results of its application to reduce the data of the photometric sky survey in Cygnus as well as observations of the open cluster NGC 188 and the transit of the exoplanet WASP-11 b / HAT-P-10b, performed with the MASTER-II-URAL telescope of the Kourovka Astronomical Observatory of the Ural Federal University.

Key words. methods: observational-methods: data analysis-techniques: photometric-stars: variables: general

## 1. Introduction

We developed Astrokit C ++ console application for postprocessing of the results of CCD photometry within the framework of a program for the search for new variable stars and exoplanet transits in the Kourovka Astronomical Observatory of the Ural Federal University. The application is based on an upgraded algorithm of di GLYPH&lt;11&gt; erential photometry using ensembles of comparison stars as described in Everett &amp; Howell (2001).

In this paper we do not consider the sources of errors that influence the precision of CCD photometry. For a detailed discussion and analysis of this problem and the use of other methods of CCD photometry see Howell et al. (1988); Gilliland &amp; Brown (1988, 1992); Newberry (1991); Young et al. (1991); Honeycutt (1992); Gilliland et al. (1993); Merline &amp; Howell (1995); Howell &amp; Everett (2001); Everett et al. (2001). We only recall that the idea of di GLYPH&lt;11&gt; erential photometry consists in determining the di GLYPH&lt;11&gt; erence between the magnitude of the source studied and that or those of one or several reference stars, thereby reducing the influence of time-variable atmospheric e GLYPH&lt;11&gt; ects. In the ideal case the reference star should have similar brightness, color index, and must be located close to the star studied.

The close location of the reference stars to the star studied is especially important in the case of the reduction of widefield images. This is because otherwise local variations of atmospheric transparency would have di GLYPH&lt;11&gt; erent e GLYPH&lt;11&gt; ect on the reference sources and on the object studied, thereby inevitably degrading the resulting photometric precision.

After performing di GLYPH&lt;11&gt; erential photometry Astrokit searches for variable stars. Below we describe the algorithm

? e-mail: burdanov.art@gmail.com

of the program, its implementation, and some of the results obtained.

## 2. Method

Before starting post-processing of the photometric data, photometry proper has to be performed. It makes no di GLYPH&lt;11&gt; erence what software was used to extract fluxes from the CCD frame, the only important thing is to obtain the corresponding magnitudes and their theoretical errors in accordance with the main CCD equation (Howell (1993)).

For the studies carried out at the Astronomical Observatory of Ural Federal University, a dedicated technique of photometric reduction was developed within the framework of IRAF package (Tody (1986)). Before applying IRAF the console version of Astrometry.net application (Lang et al. (2010)) is used to set the correct World Coordinate System parameters in the FITS header of each frame. IRAF package is then used to perform photometric reduction of each frame: subtraction of the dark frame and dividing by the flat-field frame. The PHOT/APPHOT task is then used to perform aperture photometry in each frame with individual aperture and sky background values for each frame. A catalog of objects containing the equatorial coordinates and running numbers of stars is used to this end. The aperture radius used in each particular frame is set equal to 0 8 FWHM, where FWHM is the mean full width : at half maximum value of the stellar PSF in the frame. The resulting data are then transferred to Astrokit , whose general structure is shown in Fig. 1.

Input data for the program are contained in the file created by PDUMP command of IRAF package and include:

Fig. 1: Structure of Astrokit program.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-e3fd22b6-7805-410a-9c76-5b3b7096f060.jpg)

- -identification number of the star ( id ),
- -number of counts (analog-to-digital units) from the star inside the aperture together with the sky background counts ( sum ),
- -aperture area in square pixels ( area ),
- -average sky background in each pixel ( msky ),
- -number of pixels attributed to the sky background ( nsky ),
- -exposure ( itime ).

For the program to operate correctly, it also needs a file containing the equatorial coordinates of the stars. This file may also optionally contain the catalogued color indices which will be later taken into account when selecting ensembles of comparison stars.

Below we give a stage-by-stage description of the algorithm. Let the available set of photometric data consist of j CCDframes each containing i stars from the input catalog.

- (1) The program computes the magnitudes ( m ) and magnitude errors ( merr ) for each star of the input catalog and each frame. These quantities are determined as follows:

$$\ f l u x \ = \ \sum - \ a r e a \times \ m s k y,$$

$$m \quad = \ z \max - 2. 5 \log ( \text{flux} ) + 2. 5 \log ( \text{itime} ),$$

$$\text{merr} \ = \ \frac { 1. 0 8 5 7 } { \text{flux} \times \text{gain} } \left ( \text{flux} \times \text{gain} \quad \ \text{errors of}$$

$$\begin{smallmatrix} & & \text{--} \cdots \cdots \\ & + \text{area} ( 1 + \text{area} / n s k y ) \left ( \text{msky} \times \text{gain} + \text{ron} ^ { 2 } \right ) \right ) ^ { \frac { 1 } { 2 } } \\ & + \text{area} ( 1 + \text{area} / n s k y ) \left ( \text{msky} \times \text{gain} + \text{ron} ^ { 2 } \right ) \right ) ^ { \frac { 1 } { 2 } } \end{smallmatrix}$$

where zmag = 20 is the zero point of the magnitude scale; gain is the CCD gain in e GLYPH&lt;0&gt; / ADU, which is set by the user when beginning to work with the program; ron is the CCD readout noise in e GLYPH&lt;0&gt; / pixel, which is also set by the user.

(2) For each star from the input catalog, an ensemble of reference stars located within a certain radius (the default value is 5 ) is 0

selected whose magnitudes di GLYPH&lt;11&gt; er from that of the star considered by no more than 2 m (this parameter can be changed by the user). The smaller the magnitude di GLYPH&lt;11&gt; erence, the smaller the number of stars in the ensemble, and the farther they are. Another criterion for selecting ensemble stars is based on the di GLYPH&lt;11&gt; erence between the color index of the star considered (i.e., the star for which the ensemble is composed) and candidate comparison stars.

- (3) The weighted average instrumental magnitude h mj i of ensemble stars is computed for each frame of the series

$$\langle m _ { j } \rangle = \frac { \sum _ { k } ^ { K } m _ { k j } \omega _ { k } } { \sum _ { k } \omega _ { k } }, \quad \text{where} \quad \omega _ { k } = \frac { 1 } { \langle \text{merr} _ { k } ^ { 2 } \rangle }, \\ \text{hara $1$ is the minimum number of the $\mathfrak{ta}$ $\mathfrak{ta}$} \, \mathcal { V }$$

where k is the running number of the star in the ensemble, K is the number of stars in the ensemble, j is the frame number, mkj is the magnitude of k -th ensemble star in j -th frame, h merr k i is the mean theoretical error of measured magnitude mk of the ensemble star computed by formula (1) (instead of the error in each frame as in Everett &amp; Howell (2001)).

- (4) The mean magnitude M of all ensemble stars averaged over all frames:

$$M = \frac { \sum _ { j } ^ { N } \langle m _ { j } \rangle } { N }, \\ \sim \epsilon \, \text{$\epsilon$} \, \text{$\epsilon$}$$

where N is the number of frames.

(5) The di GLYPH&lt;11&gt; erence between the weighted average magnitude h mj i of ensemble stars and the mean magnitude M of all ensemble stars averaged over all frames is then subtracted from the observed magnitude of the star for which the e GLYPH&lt;11&gt; ect of the terrestrial atmosphere is determined and for the ensemble stars:

$$m _ { \cot i j } = m _ { i j } - ( \langle m _ { j } \rangle - M ),$$

where m cor i j and mij are respectively the corrected and initial magnitudes of star i in frame j .

(6) After composing the initial ensemble, the standard deviation from the mean magnitude is computed for all stars, and the star with the greatest standard deviation is flagged. If the standard deviation from the mean magnitude is more than twice greater than the mean theoretical photometric error averaged over all frames (we call this the cuto GLYPH&lt;11&gt; ratio of the sigma criterion), which can also be varied, the star is removed from the ensemble, and the procedure is repeated from step 2.

If after all stars with large standard deviations are removed the ensemble contains less than 10 stars, the size of the ensemble domain is increased by 1 0 and all the above steps are repeated. The correction of instrumental magnitudes is thus an iterative process, which is repeated until the ensemble contains more than nine stars or the search radius increases to 30 . 0

The error introduced by the correction of the initial magnitudes using the ensemble stars is determined by the magnitude errors of the ensemble stars:

$$\text{sumsums}. \quad \text{mr} _ { \text{ens} } = \frac { 1 } { \sqrt { \sum _ { k = 1 } ^ { K } \frac { 1 } { \text{mr} _ { k } ^ { 2 } } } }, \\...................$$

where merr k is the measurement error for k -th star of the ensemble. The resulting error for star i is composed of its measurement error and the error introduced by the comparison ensemble:

$$\text{mr} = \sqrt { \text{mr} _ { i } ^ { 2 } + \text{mr} _ { e n s } ^ { 2 } }.$$

The procedure of the formation of ensembles and the correction of instrumental magnitudes is performed for each star of the list. Thus, for each star its individual ensemble of closely located comparison stars is created.

Figure 2 shows schematically the process of the correction of instrumental magnitudes.

Fig. 2: Flowchart of the process of the correction of instrumental magnitudes.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-147380d0-783c-48b5-8d6e-87bbb414111f.jpg)

(7) Variable objects are identified by the algorithm described in Rose &amp; Hintz (2007). For each star the coe GLYPH&lt;14&gt; cient RoMS (Robust Median Statistics) is computed:

$$\eta _ { n } = \frac { \sum _ { i } \frac { | \, m _ { i } - \langle m _ { \text{med} } \rangle \, | } { \sigma _ { \text{rms} } } } { N - 1 },$$

where n is the number of a star; mi is the i -th measurement of the magnitude; h m med i is the median of magnitude measurements of star n ; N is the total number of magnitude measurements for star n ; GLYPH&lt;27&gt; rms is the estimated standard deviation of star n , determined by the least-squares method from the dependence of the standard deviations of the magnitude on the stellar magnitude for stars in the frame.

The RoMS criterion allows estimating the variations of the object brightness. If it exceeds 1, the star is considered to be a suspected variable and is then analyzed more thoroughly.

## 3. Analysis of the technique

To analyze the technique and select the optimum input parameters for Astrokit to ensure the best photometric precision, we used the data of the photometric survey of a sky area in the Milky Way. A total of 3000 50-second R -band frames were taken with the MASTER-II-URAL telescope. In the central 30 0 GLYPH&lt;2&gt; 30 0 region of the frames, we selected 800 stars in the R = 9 m -17 m magnitude interval.

The MASTER-II-URAL is located at the Kourovka Observatory of the Ural Federal University. It is a Hamilton system twin telescope ( D = 40 cm, F = 1000 cm) installed on equatorial mounting and equipped with two Apogee Alta U16M CCD cameras (Lipunov et al. (2010)). The image scale is 1 00 : 8 pixel. / Photometric calibrations are performed using dark frames taken before each observing night and dawn flat-field frames. All observations are performed in automatic mode.

After performing aperture photometry with IRAF, we carried out a series of reduction cycles with Astrokit . We varied such input parameters as the initial radius r of the ensemble, the magnitude di GLYPH&lt;11&gt; erence GLYPH&lt;1&gt; m and color-index di GLYPH&lt;11&gt; erence GLYPH&lt;1&gt; CI of ensemble stars, and the cuto GLYPH&lt;11&gt; ratio ( GLYPH&lt;27&gt; ) of the standard deviation of instrumental magnitude to the theoretical error (the sigma criterion).

We consider the main criterion characterizing the quality of post-processing to be the number of 'good stars, ' i.e., stars with the standard deviation of magnitudes s of less than 0 : m 01 and 0 m : 02 over the entire observing set. We also took into account the minimum computed standard deviation of magnitude for an individual star (hereafter the 'best star').

We first varied the initial radius of the ensemble with the magnitude di GLYPH&lt;11&gt; erence fixed at 1 m (shown by the squares in plots) and the sigma criterion equal to two. The initial ensemble radii were equal to 1 , 2 , 3 , 4 , 5 , 7 , 10 , and 15 . We then counted 0 0 0 0 0 0 0 0 for each case the number of 'good stars' and the minimum standard deviation for each initial radius.

As is evident from Fig. 3, the minimum standard deviation from the mean magnitude is equal to 0 : m 00453 for the initial radii from 1 0 to 5 , 0 increases with further increase of the ensemble radius, and reaches the maximum 5% di GLYPH&lt;11&gt; erence (0 : m 00477) for the ensemble radius of 15 . 0

The number of stars whose standard deviation from the mean magnitude over the entire observing set is less than 0 : m 01 reaches maximum at the ensemble radius of 7 : it is equal to 102 which is 0 7% greater than the corresponding number of stars for the initial ensemble radii from 1 0 to 3 0 (Fig. 4). The size of the subsample of stars with the standard deviation from the mean magnitude of less than 0 : m 02 reaches maximum at the initial ensemble radius of 10 0 (254 stars compared to 245). The increment is equal to 4% (Fig. 5).

Wesimilarly varied the initial ensemble radius with GLYPH&lt;1&gt; m = 2 m (the corresponding curves in the figures are shown by the triangles). In this case the minimum standard deviation from the mean magnitude is 0 : m 00446 and also increases with initial radius (Fig. 3). The number of 'good stars' with s &lt; 0 m : 01 reaches maximum at the initial ensemble radii of 7 0 and 10 0 (it increases by 4%) (Fig. 4). The number of 'good stars' with s &lt; 0 m : 02 reaches maximum at the initial ensemble radius of 5 0 (it increases by 4% compared to the minimum value) (Fig. 5).

In view of the above, we can conclude that the optimum initial radius of the ensemble of comparison stars is 5 -7 0 0 for the magnitude di GLYPH&lt;11&gt; erence of 2 m . In this case the ensemble remains su GLYPH&lt;14&gt; ciently compact while containing a large number of stars. The ensemble permits atmospheric e GLYPH&lt;11&gt; ects to be reduced, and

Fig. 3: Variation of the standard deviation from the mean magnitude for the 'best star' as a function of the initial ensemble radius for GLYPH&lt;27&gt; = 2 and di GLYPH&lt;11&gt; erent GLYPH&lt;1&gt; m values (the dependences for GLYPH&lt;1&gt; m = 1 and GLYPH&lt;1&gt; m = 2 are shown by the squares and triangles respectively).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-aed5621a-5d74-4014-8e71-adafeb5fe25d.jpg)

Fig. 4: Variation of the number of stars with the standard deviation from the mean magnitude of less than 0 : m 01 as a function of the initial ensemble radius for GLYPH&lt;27&gt; = 2 and di GLYPH&lt;11&gt; erent GLYPH&lt;1&gt; m values (the dependences for GLYPH&lt;1&gt; m = 1 and GLYPH&lt;1&gt; m = 2 are shown by the squares and triangles respectively).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-d1c59fe9-69f3-4e87-8bc1-b2b6694437f4.jpg)

this reduction can be expressed in terms of the minimum standard deviation from the mean magnitude for the 'best star' and the maximum number of 'good stars. ' Note that the e GLYPH&lt;11&gt; ect of the varied parameters on the final result is relatively small.

The next stage in the choice of the optimum parameters consisted in varying the sigma criterion with the initial radius and magnitude di GLYPH&lt;11&gt; erence fixed at 5 0 and 2 m respectively.

As is evident from Fig. 6, the optimum value of the sigma criterion is equal to 2. This can be explained by the fact that too 'strict' value decreases the number of stars in the ensemble. A sigma criterion value greater than 2 increases the number of stars in the ensemble by including stars with the greatest standard deviation from the mean magnitude over the entire observing set with inevitable e GLYPH&lt;11&gt; ect on the final precision.

According to postulates of classical di GLYPH&lt;11&gt; erential photometry with a single comparison star and control stars, the best precision can be achieved with a comparison star that is most similar to the object studied both in terms of brightness and spectral type. We

Fig. 5: Variation of the number of stars with the standard deviation from the mean magnitude of less than 0 : m 02 as a function of the initial ensemble radius for GLYPH&lt;27&gt; = 2 and di GLYPH&lt;11&gt; erent GLYPH&lt;1&gt; m values (the dependences for GLYPH&lt;1&gt; m = 1 and GLYPH&lt;1&gt; m = 2 are shown by the squares and triangles respectively).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-651de6a3-34b5-4cb2-ae4a-2d163d5765d1.jpg)

Fig. 6: Variation of the standard deviation from the mean magnitude for the 'best star' (the top plot) and the number of stars with s &lt; 0 m : 01 (the middle plot), s &lt; 0 m : 02 (the bottom plot) as a function of the sigma criterion GLYPH&lt;27&gt; for r = 5 0 and GLYPH&lt;1&gt; m = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-d2b7f474-7a38-4f0e-aca0-5d9007970b1e.jpg)

studied the e GLYPH&lt;11&gt; ect of color index on the precision of photometry. To this end, we varied the di GLYPH&lt;11&gt; erence of the 2MASS J GLYPH&lt;0&gt; H color indices (Skrutskie et al. (2006)) when creating the ensemble of stars with the initial radius of 5 0 and the magnitude di GLYPH&lt;11&gt; erence of 2 m . We set the color-index di GLYPH&lt;11&gt; erence equal to 0 : m 1, 0 : m 2, 0 : m 3, 0 m : 4, 0 : m 5, 0 : m 6, and 0 : m 7.

As is evident from Fig. 7, the closeness of stellar spectral types is not a necessary condition for achieving high precision. However, the classical approach of di GLYPH&lt;11&gt; erential photometry remains the only solution in the case of small fields and insu GLYPH&lt;14&gt; cient number of stars for selecting an ensemble.

What will happen in the case if the ensemble contains the maximum possible number of reference stars? To this end, we varied the cuto GLYPH&lt;11&gt; magnitude for the initial ensemble radius and sigma criterion equal to 5 0 and 2 respectively (Fig. 8).

Note that the 2 m magnitude di GLYPH&lt;11&gt; erence is optimal, because not all possible stars will make it into the ensemble if a smaller value is adopted. Adopting a magnitude di GLYPH&lt;11&gt; erence greater than 2 m

Fig. 7: Variation of the standard deviation from the mean magnitude for the 'best star' (the top plot) and the number of stars with s &lt; 0 m : 01 (the middle plot), s &lt; 0 m : 02 (the bottom plot) as a function of the color-index di GLYPH&lt;11&gt; erence GLYPH&lt;1&gt; CI for r = 5 0 and GLYPH&lt;1&gt; m = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-030c1694-e329-455e-b1fc-6dda643d2a25.jpg)

Fig. 8: Variation of the standard deviation from the mean magnitude for the 'best star' (the top plot) and the number of stars with s &lt; 0 m : 01 (the middle plot), s &lt; 0 m : 02 (the bottom plot) as a function of the magnitude di GLYPH&lt;11&gt; erence GLYPH&lt;1&gt; m for r = 5 0 and GLYPH&lt;27&gt; = 2.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-6ed71705-6ed7-4a06-a0d4-585056804a43.jpg)

would result in the inclusion into the ensemble the stars that are significantly brighter or fainter than the object studied, thereby degrading the final precision.

Our analysis leads us to conclude that in the case of the presence of a su GLYPH&lt;14&gt; cient number of stars in the field, the best photometric precision is achieved using close ensembles with many stars. We thus have the following optimum parameter set for the formation of ensembles of comparison stars: the initial ensemble radius r = 5 -7 , magnitude di 0 0 GLYPH&lt;11&gt; erence GLYPH&lt;1&gt; m = 2 m , and the cuto GLYPH&lt;11&gt; ratio of the standard deviation from the mean magnitude to the theoretical error for the star to be included into the ensemble GLYPH&lt;27&gt; = 2. Note that the closeness of the spectral types of stars is of no importance.

At the first sight, small variations in the number of 'good stars' and the quality of the 'best star' may appear insignificant and not to be worthy of e GLYPH&lt;11&gt; orts to find the optimum parameters. However, first, there is no harm in increasing the precision. Second, the number of 'good stars' may increase significantly in the case of reduction of fields containing thousands of stars, and this may influence the number of candidates for searching for transits of 'hot Jupiter' type exoplanets (e.g., if only stars with s &lt; 0 m : 02 are selected).

## 4. Testing

Astrokit program is used for post reduction of photometric data obtained with wide-field telescopes of the MASTER robotic net. Below we report the results of its use for reducing a photometric survey of a Milky Way area in Cygnus, and the results of observations of the open cluster NGC 188 and the transit of the exoplanet WASP-11 b / HAT-P-10b.

We reduced all data in accordance with the same procedure. We photometrically calibrated the CCD frames using darkcurrent and flat-field frames. We observed five dark frames for the particular set of exposures at dusk hours before observations. Each master dark frame for the particular exposure is made up of five frames via median stacking. We then subtract the master dark frame from the CCD frames with the object studied.

The flat-field frames in the filter considered are acquired automatically at dawn twilight with the telescope tracking turned o GLYPH&lt;11&gt; . The linear range of the CCDs employed is limited by 40 000 ADU, and therefore the exposure in each filter is chosen so that the number of counts in pixels would not exceed this limit. The master dark frame with the appropriate exposure is subtracted from each of the five flat-field frames. The flat-field frames are normalized (the count of each pixel is divided by the median count of all pixels), and then the master flat-field frame is computed as the median of the normalized initial flat-field frames. The CCD frames with the object studied are then divided by the resulting master flat-field frame. The use of flat-field frames allows pixel-to-pixel sensitivity variations and vignetting of the optical system to be taken into account.

After photometric calibration Astrometry.net application is used to write the correct astrometric-calibration (WCS) parameters to the FIST file header. Aperture photometry is then performed using IRAF PHOT/APPHOT task, and the results of this photometry serve as input data for Astrokit .

## 4.1. Photometric Survey in Cygnus

High-precision CCD observations were performed at the Kourovka Astronomical Observatory of the Ural Federal University from May through August 2012 within the framework of a pilot project aimed at the search for exoplanet transits and variable stars. A total of 3600 frames of of a 2 GLYPH&lt;14&gt; GLYPH&lt;2&gt; 2 GLYPH&lt;14&gt; area centered at GLYPH&lt;11&gt; = 20 h 30 m 00 s , GLYPH&lt;14&gt; = 50 30 00 GLYPH&lt;14&gt; 0 00 were acquired with the MASTER-II-URAL telescope, and the photometry of 21 500 stars was performed with a precision of 0 : m 006 to 0 : m 5 in the magnitude interval from 10 m to 18 m . Figure 9 shows the standard deviation from the mean magnitude as a function of magnitude for the entire observing set. The right and left panels show the data before and after reduction with Astrokit .

As is evident from the figure, our program decreases the standard deviation from the mean magnitude by a factor of 10 for some stars.

With the RoMS = 1 threshold, Astrokit selected about 20% of all stars as candidate variables. A visual inspection of the light curves of the candidate variables allowed us to find 360 hitherto unknown variable stars including:

- -139 stars with periods longer than 20 ; d
- -100 stars with periods ranging from 20 d to 0 d : 1;
- -19 stars with periods shorter than 0 d 1; :

Fig. 9: Standard deviation from the mean magnitude as a function of magnitude for the entire observing set: the raw data (right) and the data after reduction with Astrokit (left).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-166a1ed4-6e59-4ecc-a6ab-ef3a3077bee0.jpg)

- -96 eclipsing variables;
- -5 flare stars;
- -1 dwarf nova.

Among the variables found was the outburst of the dwarf nova USNO-B1.01413-0363790 (Burdanov et al. (2012)). Several dozen variable stars were selected as candidate GLYPH&lt;14&gt; Sct stars. The brightness of the star was found to vary with an amplitude as small as 0 : m 005 (Fig. 10).

Fig. 10: Phased light curve of the star 2MASS 20295743 + 5017071 with the period and amplitude of 0 035 and : d 0 m : 005 respectively. The solid line shows the fitted data.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-977107c7-5c1d-4371-9de7-2f33f46be316.jpg)

Other stars were found to exhibit periodic brightness decreases by 0 : m 015. These light variations may be due to transits of 'hot Jupiter' type exoplanets with the radii of 1 10 of the ra-/ dius of the host stars and orbiting with less than one day-long periods (Burdanov et al. (2013)).

## 4.2. Open Cluster NGC 188

The open cluster NGC188 was observed in the R and I filters with the MASTER-II-URAL telescope over five nights in March 2011. A total of 400 frames were acquired, which were post-processed using Astrokit program. We performed aperture R -band photometry of 5513 stars with a precision ranging from 0 m : 006 to 0 : m 05 in a 1 5 GLYPH&lt;14&gt; : GLYPH&lt;2&gt; 1 5 area. Although NGC 188 is a GLYPH&lt;14&gt; :

Fig. 11: Standard deviations from the mean magnitude (left) and coe GLYPH&lt;14&gt; cients RoMS (right) as a function of magnitude. The filled circles show variable stars.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-b8dc069d-761e-452e-9f46-57c4f1be059e.jpg)

well-studied cluster (more than 500 papers in the past 50 years), our algorithm of the search for variable stars found 18 new variables (Popov et al. (2013)). Figure 11 shows the standard deviation from the mean magnitude as a function of magnitude after reduction by our program (left) and the RoMS coe GLYPH&lt;14&gt; cients for all stars (right). Variable stars are shown by the filled circles in both plots. As is evident from Fig. 11, the standard deviation from the mean magnitude for variables often does not di GLYPH&lt;11&gt; er from the corresponding standard deviations for constant stars, whereas the RoMS coe GLYPH&lt;14&gt; cient allows more confident selection of candidate variable objects.

## 4.3. Observations of Exoplanet Transits

During the transits of even the largest 'hot Jupiters' moving in close orbits across the disk of the host star, its brightness dips by about 0 : m 01. To confidently record the very fact of transit and determine its duration and mid-time, in the case of such brightness dips the resulting photometry should have the precision of about 0 m : 001. Achieving the given precision is a di GLYPH&lt;14&gt; cult task for ground-based telescopes and requires special attention during di GLYPH&lt;11&gt; erential photometry.

Below we compare the results obtained with Astrokit and the classical technique of di GLYPH&lt;11&gt; erential photometry with a single comparison star, used to reduce the photometric observations of the transit of a 'hot Jupiter' type exoplanet.

Our R -band observations of the transit of the exoplanet WASP-11b / HAT-P-10b ( R = 11 m : 01, J GLYPH&lt;0&gt; H = 0 46, : H GLYPH&lt;0&gt; K = 0 14) were made with the MASTER-II-Ural telescope : on December 10, 2012.

The light curve of the transit was obtained using the classical method of di GLYPH&lt;11&gt; erential photometry, with a single comparison star 3UC-242-019494 ( R = 12 m , J GLYPH&lt;0&gt; H = 0 55, : H GLYPH&lt;0&gt; K = 0 125), : and a check star 3UC-242-019559 ( R = 11 m : 4, J GLYPH&lt;0&gt; H = 0 175, : H GLYPH&lt;0&gt; K = 0 113). : These stars are located within less than 10 0 from WASP-11 HAT-P-10 ( / R = 11 m : 7) and have close magnitudes and color indices. The light curve obtained from magnitude di GLYPH&lt;11&gt; erence between WASP-11 HAT-P-10 and the compari-/ son star is shown by the triangles in Fig. 12. The standard deviation of the magnitude di GLYPH&lt;11&gt; erence between the comparison and control star was 0 : m 006.

We used Astrokit program with an ensemble of 11 stars to correct photometry for atmospheric transparency variation. The stars are located within 6 0 and their magnitudes di GLYPH&lt;11&gt; er by no

Fig. 12: R -band light curve of the exoplanet WASP-11b / HAT-P-10b transit. The light curve obtained with a single comparison star is shown by the triangles, and that obtained using an ensemble of comparison stars-by the squares (it is shifted for better visualization).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0664v1-figure-a59c3935-e30c-4306-b419-33142e463782.jpg)

more than 2 m from that of the star studied. We estimate the precision of corrected photometry (0 : m 0039) by the standard deviation from the mean magnitude for the control star 3UC-242-019559 mentioned above. The resulting precision is about a factor of 1.5 smaller that the standard deviation of the control-star magnitude in the case the classical method of di GLYPH&lt;11&gt; erential photometry is used, thereby demonstrating the advantage of using ensembles of comparison stars. The resulting transit light curve is shown by the squares in Fig. 12.

This transit is remarkable by the fact that it was observed during total lunar eclipse. We clearly see that the scatter of data points increases toward the end of the transit. The standard deviation of the magnitude of the host star before the transit (the first 11 data points in the light curve) is equal to 0 : m 002 and begins to increase with increasing sky background due to the egress of the Moon out of the Earth shadow. The standard deviation of hoststar magnitude after the transit (the last 11 data points) is four times greater and equal to 0 : m 008.

## 5. Conclusions

The use of a close ensemble of comparison stars while performing di GLYPH&lt;11&gt; erential photometry allows the inhomogeneity of the data series, caused by local variations of atmospheric transparency and sky background, to be taken into account, and reduces the contribution of stellar scintillation to the error budget of the resulting magnitudes (Everett &amp; Howell (2001); Kornilov (2012)). If more than 10 reference stars in a close ensemble are used, the di GLYPH&lt;11&gt; erence between their spectral types and that of the object studied becomes unimportant. However, to achieve the best precision, the magnitude di GLYPH&lt;11&gt; erence should be small (it should not exceed 2 m ) and ensemble stars should be chosen within 5 -7 0 0 of the program star.

The resulting photometric precision after the application of Astrokit allows finding low-amplitude variables and study transits of 'hot Jupiters. '

The use of robust median statistics is approved because of its higher stability against outliers. However, it does not guarantee against finding false variables. A considerable part of suspected variables prove to be constant stars within the photometric errors after further analysis. The number of candidate variables is usually equal to about 10% of the total number of stars in the frame. The final criterion of the variability of a star is the visual inspection of its light curve.

The source code of the program is available at http:// astro.ins.urfu.ru/en/node/1330 .

Acknowledgements. This work was supported by the Russian Foundation for Basic Research (projects No. 2-02-31095 and 14-02-31338). We are also grateful to Ekaterina Avvakumova for her assistance with the preparation of the paper and support. We also thank Kirill Ivanov for testing the program. This research has made use of the VizieR catalogue, SIMBAD database, and Aladin sky atlas, operated at CDS, Strasbourg, France.

## References

Burdanov, A. Y., Krushinsky, V. V., Denisenko, D., et al. 2012, Peremennye Zvezdy Prilozhenie, 12, 24

Burdanov, A. Y., Popov, A. A., Krushinsky, V. V., &amp; Ivanov, K. 2013, Peremennye Zvezdy, 33, 2

Everett, M. E. &amp; Howell, S. B. 2001, PASP, 113, 1428

Everett, M. E., Howell, S. B., &amp; Ousley, D. 2001, in Third Workshop on Photometry, ed. W. J. Borucki &amp; L. E. Lasher, Proceedings of a Workshop held at the SETI Institute in Palo Alto, California, September 24-25, 1998 (Mo GLYPH&lt;11&gt; ett Field, CA 94035-1000: Ames Research Center), 79-83

Gilliland, R. L. &amp; Brown, T. M. 1988, PASP, 100, 754

Gilliland, R. L. &amp; Brown, T. M. 1992, PASP, 104, 582

Gilliland, R. L., Brown, T. M., Kjeldsen, H., et al. 1993, AJ, 106, 2441

Honeycutt, R. K. 1992, PASP, 104, 435

Howell, S. B. 1993, in Stellar Photometry-Current Techniques and Future Developments, ed. C. J. Butler &amp; I. Elliott, Proceedings of the IAU Colloquium No. 136 held in Dublin, Ireland, August 4-7, 1992 (Cambridge: Cambridge University Press), 318

Howell, S. B. &amp; Everett, M. E. 2001, in Third Workshop on Photometry, ed. W. J. Borucki &amp; L. E. Lasher, Proceedings of a Workshop held at the SETI Institute in Palo Alto, California, September 24-25, 1998 (Mo GLYPH&lt;11&gt; ett Field, CA 94035-1000: Ames Research Center), 1-7

Howell, S. B., Warnock, III, A., &amp; Mitchell, K. J. 1988, AJ, 95, 247 Kornilov, V. 2012, MNRAS, 425, 1549

Lang, D., Hogg, D. W., Mierle, K., Blanton, M., &amp; Roweis, S. 2010, AJ, 139, 1782

Lipunov, V., Kornilov, V., Gorbovskoy, E., et al. 2010, Advances in Astronomy, 2010

Merline, W. J. &amp; Howell, S. B. 1995, Experemental Astromomy, 6, 163 Newberry, M. V. 1991, PASP, 103, 122

- Popov, A. A., Krushinsky, V. V., Avvakumova, E. A., et al. 2013, Open European Journal on Variable Stars, 157, 1

Rose, M. B. &amp; Hintz, E. G. 2007, AJ, 134, 2067

Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163

Tody, D. 1986, in Proceedings of the Meeting, Tucson, AZ, March 4-8, 1986, Vol. 627, Instrumentation in astronomy VI, ed. D. L. Crawford (Bellingham, WA: Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series), 733

Young, A. T., Genet, R. M., &amp; Boyd, et al., L. J. 1991, PASP, 103, 221

Translated by A. Dambis
Artem Y. Burdanov Vadim V. Krushinsky Alexander A. Popov
Astrokit -- an Efficient Program for High-Precision Differential CCD Photometry and Search for Variable Stars
Having a need to perform differential photometry for tens of thousands stars
in a several square degrees field, we developed Astrokit program. The software
corrects the star brightness variations caused by variations of atmospheric
transparency: to this end, the program selects for each star an individual
ensemble of reference stars having similar magnitudes and positions in the
frame. With ten or more reference stars in the ensemble, the differences
between their spectral types and the spectral type of the object studied become
unimportant. Astrokit searches for variable stars using Robust Median
Statistics criterion, which allows candidate variables to be selected more
efficiently than by analyzing the standard deviation of star magnitudes. The
software allows very precise automatic analysis of long inhomogeneous sets of
photometric observations of a large number of objects to be performed, making
it possible to find "hot Jupiter" type exoplanet transits and low-amplitude
variables. We describe the algorithm of the program and the results of its
application to reduce the data of the photometric sky survey in Cygnus as well
as observations of the open cluster NGC188 and the transit of the exoplanet
WASP-11 b / HAT-P-10 b, performed with the MASTER-II-URAL telescope of the
Kourovka Astronomical Observatory of the Ural Federal University.