Analysis of RR Lyrae Stars in the Northern Sky Variability Survey K. Kinemuchi1 and H.A. Smith Department of Physics & Astronomy, Michigan State University, East Lansing, MI 48824 P. R. Wozniak Los Alamos National Laboratory, Los Alamos, NM 87545 T. A. McKay Department of Physics, 2477 Randall Laboratory, University of Michigan, Ann Arbor, MI 48109 and ROTSE Collaboration ABSTRACT We use data from the Northern Sky Variability Survey (NSVS), obtained from the first generation Robotic Optical Transient Search Experiment (ROTSE-I), to identify and study RR Lyrae variable stars in the solar neighborhood. We initially identified 1197 RRab (RR0) candidate stars brighter than the ROTSE median magnitude V = 14. Periods, amplitudes, and mean V magnitudes are determined arXiv:astro-ph/0606092v1 5 Jun 2006 for a subset of 1188 RRab stars with well defined light curves. Metallicities are determined for 589 stars by the Fourier parameter method and by the relationship between period, amplitude, and [Fe/H]. We comment upon the difficulties of clearly classifying RRc (RR1) variables in the NSVS dataset. Distances to the RRab stars are calculated using an adopted luminosity-metallicity relation with corrections for interstellar extinction. The 589 RRab stars in our final sample are used to study the properties of the RRab population within 5 kpc of the Sun. The Bailey diagram of period versus amplitude shows that the largest component of this sample belongs to Oosterhoff type I. Metal-rich ([Fe/H] > −1) RRab stars appear to be associated 1present address: University of Wyoming, Department of Physics & Astronomy, Dept. 3905, Laramie, WY 82071 –2– with the Galactic disk. Our metal-rich RRab sample may include a thin disk as well as a thick disk population, although the uncertainties are too large to establish this. There is some evidence among the metal-rich RRab stars for a decline in scale height with increasing [Fe/H], as was found by Layden (1995). The distribution of RRab stars with −1 < [Fe/H] < −1.25 indicates that within this metallicity range the RRab stars are a mixture of stars belonging to halo and disk populations. Subject headings: stars: RR Lyrae, Galaxy: thick disk, Oosterhoff dichotomy 1. Introduction RR Lyrae variable stars are versatile objects for astronomical research. Not only are they excellent distance indicators, but they are useful as probes for understanding Galactic evolution and structure. RR Lyrae variables have the virtues of being both luminous and ubiquitous in the Galaxy, tracing the old stellar populations of the bulge, disk, and halo components. In this paper, we present results from a new survey of Bailey ab-type RR Lyrae stars (RRab or RR0) in the solar neighborhood based upon data in the Northern Sky Variability Survey (Wo´zniak et al. 2004a), hereafter NSVS. Bailey type c (RRc or RR1) variables in the NSVS will be discussed in a future paper. To get a clearer picture of the Galaxy, we need as complete a database as possible of both field and globular cluster RR Lyrae stars. The General Catalog of Variable Stars (GCVS) (Kholopov et al. 1996) summarizes much previous work on the discovery of RR Lyrae stars in the Galactic field, but is significantly incomplete, especially for RR Lyrae stars of low amplitude (Akerlof et al. 2000). More recent CCD surveys of RR Lyrae stars in the Galactic bulge (OGLE: Wo´zniak et al. (2002), MACHO: Alcock et al. (1996)) and halo (SDSS: Ivezi´c et al. (2000), QUEST: Vivas et al. (2004)) have extended the coverage of the GCVS, but do not emphasize the discovery of brighter RR Lyrae stars in the solar neighborhood. The Northern Sky Variability Survey (Wo´zniak et al. 2004a) is based on photometric observations obtained with the first generation telescope of the Robotic Optical Transient Search Experiment (ROTSE-I) (see section 2). Within the NSVS, RR Lyrae stars as faint as 15th magnitude in V can be detected, extending to a distance of about 7-9 kpc from the Sun (Akerlof et al. 2000). This region includes part of the inner halo and thick disk components of the Galaxy. Hence, this survey is an excellent complement to the bulge and outer halo surveys, bridging the portions of the Galaxy that they cover. The NSVS data complement those of the QUEST survey (Vivas et al. 2004), which includes RR Lyrae stars from 4 to 60 –3– from the sun, but principally samples RR Lyrae stars more distant than those in the NSVS. Another survey that will complement the NSVS is the All Sky Automated Survey (ASAS) (Pojmanski 1997). The ASAS database surveys the southern sky up to δ = +28 and has a magnitude limit of V = 14. In this paper, we will limit ourselves to a consideration of those NSVS RRab stars that are brighter than 14th magnitude in the ROTSE-I system. These stars have light curves sufficiently well defined that the periods are secure, and the classification of Bailey type is clear. Within certain limits, which we shall describe later, we expect this survey to provide a more complete and unbiased catalog of the RRab stars in the solar neighborhood than is provided by the GCVS. Because of this magnitude restriction, our paper mainly concerns RRab stars to a distance of 5 kpc from the Sun. Section 2 of the paper briefly describes the ROTSE-I telescope and the NSVS database. In section 3, we discuss the selection criteria used to identify RRab stars, the determination of periods, the photometric calibration, and methods for determining photometric metallicities. In section 4, the properties of the RRab sample and its distribution in space are discussed. Results are summarized in section 5. 2. NSVS Database The NSVS database1 is constructed from observations taken from the first generation Robotic Optical Transient Experiment (ROTSE-I). The ROTSE-I telescope, located in Los Alamos, New Mexico, was originally designed to find the optical counterparts to gamma ray bursters. While waiting for a trigger from a gamma ray burster, the telescope was set to a patrol mode, which scanned the sky. The NSVS is based on a year’s worth of observations, beginning in April of 1999 and ending in March 2000. The telescope, now retired, was comprised of 4 smaller telescopes mounted on a rapidly slewing platform. Each telescope of the 2 × 2 array of telescopes was equipped with a 2048 × 2048 CCD. The main optical element of each telescope was a Canon FD 200mm f/1.8 telephoto lens (Kehoe et al. 2001). The plate scale was 14.4′′ per pixel, which corresponded to a 16.4◦ × 16.4◦ field of view. Each individual telescope had a field of view of 8.2◦ × 8.2◦. Thus, the sky was divided into 206 fields (see Figure 1 of Akerlof et al. 2000b). Each night the visible fields were observed in pairs of 80 second exposures, while in patrol mode. In some cases, multiple pairs of observations were taken for some fields in a single night. 1See http://skydot.lanl.gov –4– The ROTSE-I telescope did not have a standard filter set installed. The instrumental photometric bandpass therefore does not exactly correspond to those of any of the Johnson- Cousins filters, although the peak sensitivity of the system is similar to that of the Cousins R band. The calibration of instrumental magnitudes to Johnson V will be discussed in section 3.4. The basic data reductions, i.e. bias/dark subtraction and flat fielding, source detection, centroiding and preliminary calibrations of positions and magnitudes were performed using the legacy ROTSE-I software (Akerlof et al. 2000). The ROTSE-I pipeline performs source extraction and flux measurements using the SExtractor package (Bertin & Arnouts 1996). The astrometric and photometric calibration of each image is based on approximately 1500 field stars from the Tycho catalog. NSVS magnitudes are unfiltered ROTSE-I magnitudes converted to a V-band scale of a mean Tycho star. The details of the light curve formation, the photometric corrections, and data quality indicators are given in Wo´zniak et al. (2004a). The NSVS data processing flags capture a number of potential problems with the measure- ments and helped to identify observations suitable for our analysis. In our work, we have converted the modified Julian dates listed in the NSVS to heliocentric Julian dates. 3. RR Lyrae Data The data used for this work were obtained from a version of the NSVS database at the University of Michigan. At the University of Michigan, a cross correlation of the variable star candidates with the 2MASS database had been completed. This gave us additional information for the selection of RR Lyrae candidates. The results presented in this paper is based on the data obtained from this version of the NSVS database. 3.1. Selection Criteria To identify whether a star in the NSVS database was a variable star, a variability index was used. The index is based on previous work outlined in Akerlof et al. (2000). The variability index itself is based on the work of Welch & Stetson (1993) and Stetson (1987, 1996). The main variability index used in this work is Stetson’s L index, which is constructed from Stetson’s J and K indices, which are defined below. This particular method of finding variable stars was chosen because the L variability index is optimized for pairs of observations. This method also allows us to give a lower weight to spuriously bright or faint measurements and to avoid misidentification of non-variable stars. This index helps exclude transient events such as flare stars, which, though interesting in their own right, are not the focus of this paper. –5– Three separate calculations must be completed to arrive at the final L index.