Hostname: page-component-cd9895bd7-dk4vv Total loading time: 0 Render date: 2024-12-22T21:32:49.532Z Has data issue: false hasContentIssue false

Solar Space Density of the Red Clump Stars and the Scale-Length of the Thin Disc

Published online by Cambridge University Press:  08 April 2015

E. Yaz Gökçe*
Affiliation:
Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, 34119, Beyazıt, Istanbul, Turkey
S. Karaali
Affiliation:
Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, 34119, Beyazıt, Istanbul, Turkey
Ş. Duran
Affiliation:
Department of Astronomy and Space Sciences, Graduate School of Science and Engineering, Istanbul University, 34116, Beyazıt, Istanbul, Turkey
S. Bilir
Affiliation:
Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, 34119, Beyazıt, Istanbul, Turkey
A. Yalçınkaya
Affiliation:
Mithatpaşa cd. 16/4, 35320, Narlıdere, İzmir, Turkey
S. Ak
Affiliation:
Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, 34119, Beyazıt, Istanbul, Turkey
T. Ak
Affiliation:
Faculty of Science, Department of Astronomy and Space Sciences, Istanbul University, 34119, Beyazıt, Istanbul, Turkey
M. López-Corredoira
Affiliation:
Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain Departamento Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain
A. Cabrera-Lavers
Affiliation:
Instituto de Astrofísica de Canarias, E-38205 La Laguna, Tenerife, Spain Departamento Astrofísica, Universidad de La Laguna, E-38206 La Laguna, Tenerife, Spain
Rights & Permissions [Opens in a new window]

Abstract

We estimated the scale-length of the thin disc with the J and W1 magnitudes of the most probable red clump stars in the Galactic plane, − 0°.5 ⩽ b ⩽ +0°.5, in 19 equal sized fields with consecutive Galactic longitudes which cover the interval 90° ⩽ l ⩽ 270°. Our results are constrained with respect to the solar space density (D* = 5.95), which indicates that the radial variation of the density is lower for higher Galactocentric distances. The scale-length of the thin disc is 2 kpc for the fields in the Galactic anticentre direction or close to this direction, while it decreases continuously in the second and third quadrants reaching to a lower limit of h = 1.6 kpc at the Galactic longitudes l = 90° and l = 270°. The distribution of the scale-length in 19 fields is consistent with the predictions from the Galaxia model and its variation with longitude is probably due to the inhomogeneity structure of the disc caused by the accreted material or other features such as warp and flare.

Type
Research Article
Copyright
Copyright © Astronomical Society of Australia 2015 

1 INTRODUCTION

The photometrical and spectroscopical data obtained from the systematical sky survey programs have an important impact on the investigation of the structure, formation and evolution of the Galaxy. By the beginning of 2000s, the data from the near-infrared sky surveys have contributed to the studies of the Galactic disc; for instance, the relation between the star formation regions and spiral arms, distribution of open clusters through the Galactic disc, formation and evolution of the open clusters, three-dimensional distribution of the interstellar matter in the Galactic disc, spatial distribution of the stars from different spectral types.

Estimation of the Galactic model parameters, such as the scale-length, the scale-height, and the relative space densities of the Galactic populations, is a part of the investigation of the Galactic structure. The researchers used different methods for estimation of the Galactic model parameters for the thin disc, thick disc, and halo of our Galaxy. There has been an improvement for numerical values of the model parameters over 30 years, starting from the work of Gilmore & Reid (Reference Gilmore and Reid1983). We present those for only the thin disc in Table 1. There are many striking points for the numerical values of the Galactic model parameters for a given Galactic component, thin and thick discs and halo, as stated in the following: (a) the range of values is large. For example, Chen et al. (Reference Chen2001) and Siegel et al. (Reference Siegel, Majewski, Reid and Thompson2002) give the local space density of the thick disc as 6.5 − 13.5% and 6–10 %, respectively. (b) the Galactic model parameters are absolute magnitude dependent (Karaali, Bilir, & Hamzaoǧlu Reference Karaali, Bilir and Hamzaoǧlu2004; Bilir, Karaali, & Gilmore Reference Bilir, Karaali and Gilmore2006a), (c) the Galactic model parameters are Galactic longitude/latitude dependent (Buser, Rong, & Karaali Reference Buser, Rong and Karaali1998, Reference Buser, Rong and Karaali1999; Bilir et al. Reference Bilir, Karaali, Ak, Yaz and Hamzaoğlu2006b, Reference Bilir, Karaali, Güver, Karataş and Ak2006c, Reference Bilir, Cabrera-Lavers, Karaali, Ak, Yaz and López-Corredoira2008; Cabrera-Lavers et al. Reference Cabrera-Lavers, Bilir, Ak, Yaz and López-Corredoira2007a; Ak et al. Reference Ak, Bilir, Karaali and Buser2007; Yaz & Karaali Reference Yaz and Karaali2010), (d) the Galactic model parameters change with limiting distance of completeness (Karaali et al. Reference Karaali, Bilir, Yaz, Hamzaoǧlu and Buser2007). We should add that the scale-length also depends on the wavelength (cf. Beckman et al. Reference Beckman, Peletier, Knapen, Corradi and Gentet1996; Tuffs et al. Reference Tuffs, Popescu, Vlk, Kylafis and Dopita2004).

Table 1. Scale-height and scale-length values of thin disc in the literature.

The observed features cited above can be explained by disc flares and warps. The disc of our Galaxy is not radially smooth and uniform. While the warp bends the Galactic plane upwards within the longitudes 0° < l ⩽ 180° and downwards within the longitudes 180° < l ⩽ 360°, the flare changes the scale-height as a function of radial distance. Warp has been detected in the dust, gas and star components of the Galaxy (Drimmel & Spergel Reference Drimmel and Spergel2001; López-Corredoira et al. Reference López-Corredoira, Cabrera-Lavers, Garzón and Hammersley2002; Momany et al. Reference Momany, Zaggia, Gilmore, Piotto, Carraro, Bedin and de Angeli2006; Marshall et al. Reference Marshall, Robin, Reylé, Schultheis and Picaud2006). The stellar and gaseous flares are also consistent with each other (Momany et al. Reference Momany, Zaggia, Gilmore, Piotto, Carraro, Bedin and de Angeli2006; Kalberla et al. Reference Kalberla, Kerp, Dedes and Haud2014).

The overdensity regions also do affect the numerical values of the Galactic model parameters for a given Galactic component, thin and thick discs, and halo. This is the case for very large fields such as in the work of Jurić et al. (Reference Jurić2008), 6500 deg2. Such fields are contaminated by remnants of mergers which cause different numerical values for a specific model parameter, relative to the ones estimated in a relatively homogeneous region of the Galaxy. However, one should pay attention in explaining the origin of the overdensity regions in the Galaxy, i.e., such a region may originate either from a merger or from the affect of a warp of the Galaxy. The Canis Major overdensity region can be given as a typical example. Martin et al. (Reference Martin, Ibata, Conn, Lewis, Bellazzini, Irwin and McConnachie2004) interpreted the overdensity region with Galactic coordinate (l, b) = (240°, −7°) as the core of a satellite galaxy. However, Momany et al. (Reference Momany, Zaggia, Gilmore, Piotto, Carraro, Bedin and de Angeli2006), López-Corredoira (Reference López-Corredoira2006) and López-Corredoira et al. (Reference López-Corredoira, Momany, Zaggia and Cabrera-Lavers2007) confirmed its coincidence with the Southern stellar maximum warp occurring near l ~ 240° and R ~ 7 kpc.

The red clump (RC) stars can be used for the investigation of the galactic structure. RC stars are core helium-burning giants. Their position in the colour-magnitude diagram (CMD) of open clusters as well as globular clusters is discrete relative to the main sequence stars and red giants of these clusters. RC stars are also abundant in the solar neighbourhood which provide accurate absolute magnitudes due to their parallaxes. The small scattering of this population in the CMD has been used as an advantage to adopt a mean absolute magnitude for them in different optical and infrared (IR) bands. The final aim of the researchers was to use the mean absolute magnitude in distance estimation; could the RC stars be used as distance indicators?

Keenan & Barnbaum (Reference Keenan and Barnbaum1999) claimed that the absolute magnitudes in the optical range of the RC stars within the spectral range of G8 III–K2 III lie within 0.7 ⩽ MV ⩽ 1 mag. Also, we should note that the Alves (Reference Alves2000) sample covered a metallicity range of − 0.8 < [Fe/H] < +0.2 dex. The absolute magnitude computed by Grocholski & Sarajedini (Reference Grocholski and Sarajedini2002) in the same band based on 14 open clusters with metallicity and age restrictions is compatible with the one of Alves (Reference Alves2000).

An absolute magnitude in any optical or IR band free of metallicity and age would be the ideal one for the researchers. This was the case for the I band in Udalski (Reference Udalski1998). However, a modest variation in MI with colour and metallicity has been claimed in the literature (Paczyński & Stanek Reference Paczyński and Stanek1998; Stanek & Garnavich Reference Stanek and Garnavich1998; Sarajedini Reference Sarajedini1999; Zhao et al. Reference Zhao, Qiu and Zhang2001; Kubiak et al. Reference Kubiak, McWilliam, Udalski and Gorski2002). The dependence of the MI absolute magnitude on metallicity has been confirmed also theoretically (Salaris & Girardi Reference Salaris and Girardi2002).

The numerical values claimed that the absolute magnitude M Ks in recent works are compatible with the former ones. van Helshoecht & Groenewegen (Reference van Helshoecht and Groenewegen2007) showed that M Ks = −1.57 ± 0.05 mag can be used for the determination of the distances of the open clusters with metallicities − 0.5 < [Fe/H] ⩽ +0.4 dex and age 0.31 ⩽ t ⩽ 7.94 Gyr. A value slightly different than the one just cited is the one in Groenewegen (Reference Groenewegen2008), i.e., M Ks = −1.54 ± 0.03 mag, which is based on the re-reduced Hipparcos catalogue (van Leeuwen Reference van Leeuwen2007). Groenewegen (Reference Groenewegen2008) determined also the IR absolute magnitude for this sample as MI = −0.22 ± 0.03 mag.

There are even more recent works related to the absolute magnitude of RC stars. Laney , Joner, & Pietrzyski (Reference Laney, Joner and Pietrzyski2012) determined the absolute magnitude M Ks = −1.613 ± 0.015 mag for RC stars in the solar neighbourhood and applied it to estimate the distance of the Large Magellanic Cloud. Karaali, Bilir & Yaz Gökçe (Reference Karaali, Bilir and Yaz Gökçe2013) fitted the M Ks absolute magnitudes of the RC stars to a linear polynomial in terms of VKs colour, and Bilir et al. (Reference Bilir, Önal, Karaali, Cabrera-Lavers and Çakmak2013a) calibrated the M Ks, MV, MJ and Mg absolute magnitudes in terms of colours corresponding to the photometric systems where the Ks, V, J and g bands are defined. Finally, we quote the paper of Bilir et al. (Reference Bilir, Ak, Ak, Yontan and Bostancı2013b) where the MV absolute magnitude had been calibrated in terms of metallicity and BV colour. The data are based on a large sample of open and globular clusters with a wide range of metallicity. A mean value for the M Ks absolute magnitude with weak metallicity dependence can be used for distance estimation which is rather useful for the Galactic structure and interstellar extinction. Examples can be found in López-Corredoira et al. (Reference López-Corredoira, Cabrera-Lavers, Garzón and Hammersley2002, Reference López-Corredoira, Cabrera-Lavers, Gerhard and Garzón2004), Cabrera-Lavers, Garzón, & Hammersley (Reference Cabrera-Lavers, Garzón and Hammersley2005); Cabrera-Lavers et al. (Reference Cabrera-Lavers, Bilir, Ak, Yaz and López-Corredoira2007a, Reference Cabrera-Lavers, Hammersley, González-Fernández, López-Corredoira, Garzón and Mahoney2007b, Reference Cabrera-Lavers, González-Fernández, Garzón, Hammersley and López-Corredoira2008) and Bilir et al. (Reference Bilir, Karaali, Dağtekin, Önal, Ak, Ak and Cabrera-Lavers2012) and references therein.

The absolute magnitudes of the RC stars has been investigated also with the Wide-field Infrared Survey Explorer (WISE; Wright et al. Reference Wright2010). Yaz Gökçe et al. (Reference Gökçe, Bilir, Öztürkmen, Duran, Ak, Ak and Karaali2013) used the RC stars in the solar neighbourhood and determined a mean absolute magnitude on each band of W1 (3.4 μm) and W3 (12 μm) as well as a mean colour, i.e., M W 1 = −1.635 ± 0.026, M W 2 = −1.606 ± 0.024, and (W1 − W3)0 = −0.028 ± 0.001 mag. WISE can go at least one magnitude deeper than the Two Micron All Sky Survey (2MASS; Skrutskie et al. Reference Skrutskie2006). Also, WISE can be affected from the Galactic extinction less than the other photometric systems, including 2MASS. Hence, one expects more accurate results for a larger sample of RC stars investigated with WISE.

In this paper, we estimate the scale-length of the thin disc in the Galactic plane. We used the RC stars with Galactic latitude − 0°.5 ⩽ b ⩽ +0°.5 and Galactic longitude 90° ⩽ l ⩽ 270° and evaluated the space density functions for equally spaced 19 fields. The Galactocentric radial distance range in our fields is 8 ⩽ R ⩽ 16 kpc. Additionally, the Galactic latitude range, − 0°.5 ⩽ b ⩽ 0°.5, is rather small. Hence, we considered the flare as negligible and fitted the observational space density profiles to a Galaxy model which involves only the exponential factor with scale-length. The evaluation of the solar space density of the RC stars by using the Hipparcos catalogue (van Leeuwen Reference van Leeuwen2007) is given in Section 2. Section 3 is devoted to the estimation of scale-length, and finally a discussion is presented in Section 4.

2 SOLAR SPACE DENSITY OF THE RC STARS

The RC star sample (the main sample) in the WISE All-Sky Data Release catalogue of Cutri et al. (Reference Cutri2013) used for the estimation of the scale-length in the Galactic plane is not complete within distances 0 < r ⩽ 1.5 kpc. Here, r is the distance from the Sun. Hence, we used a second sample taken from the re-reduced Hipparcos catalogue (van Leeuwen Reference van Leeuwen2007) just to fix the position of the RC stars in the solar neighbourhood. Thus, the solar space density of the second sample would be used to avoid any probable degeneracy in the space density functions of the RC stars in the main sample.

RC stars consist of a sub-sample of Hipparcos catalogue which contains about 118 000 objects (van Leeuwen, Reference van Leeuwen2007). We omitted those with relative parallax errors σπ/π ⩽ 0, so we used 113 596 of them with σπ/π > 0 to obtain more accurate data. The estimation of the interstellar extinction is based on the following density law for dust in Sharma et al. (Reference Sharma, Bland-Hawthorn, Johnston and Binney2011) which can be reduced to an equation with respect to the Galactic coordinates and distance of a star:

(1) \begin{eqnarray} \rho _{\text{dust}}(R,z)=\frac{\rho _0}{k_{fl}} \times \exp \left[\frac{R_0-R}{h_R}-\frac{z-z_w}{k_{fl}h_z}\right], \end{eqnarray}

where R 0 = 8 kpc is the distance of the Sun to the Galactic centre, hR (4.2 kpc) and hz (88 pc) are the scale-length and the scale-height of the dust, ρ0 = 0.54 mag kpc− 1, and kfl and zw describe the flaring and warping of the gas, respectively (Sharma et al. Reference Sharma, Bland-Hawthorn, Johnston and Binney2011):

(2) \begin{eqnarray} k_{fl}(R) & = & 1+\gamma _{fl} \times \min (R_{fl} ,R-R_{fl}), \nonumber\\ z_w(R,\Phi ) & =&\gamma _{w} \times \min (R_w ,R-R_w) \times \sin \Phi . \end{eqnarray}

Here, γfl, γw, Rfl, and Rw values are 0.0054, 0.18, 1.12R 0 and 8.4 kpc, respectively. Φ is the Galactocentric azimuth that increases in the direction of the Galactic rotation and place the Sun at Φ = 0.

Following the comments of Arce & Goodman (Reference Arce and Goodman1999) and Schlafly & Finkbeiner (Reference Schlafly and Finkbeiner2011), we reduced the colour excess at infinite distance in Schlegel, Finkbeiner, & Davis (Reference Schlegel, Finkbeiner and Davis1998), which will be labelled as $E_{\infty } (B-V)_{\text{initial}}$ hereafter, to the final one by the following equation (Binney et al. Reference Binney2014):

(3) \begin{eqnarray} E_{\infty } (B-V)_{\text{final}}&=&E_{\infty } (B-V)_{\text{initial}} \nonumber\\ &&\times C(E(B-V)), \end{eqnarray}

where,

(4) \begin{eqnarray} C(E(B-V)) &=& 0.6 + 0.2 \nonumber\\ &&\times \left[1-\tanh \left(\frac{E_{\infty } (B-V)_{\text{initial}}-0.15}{0.3}\right)\right]. \end{eqnarray}

Then, the colour excess of a star at distance r can be evaluated by the following equation:

(5) \begin{eqnarray} E_r(B-V)=E_{\infty } (B-V)_{\text{final}} \times \frac{\int _0^r \! \rho (s)ds}{\int _0^{\infty } \! \rho (s)ds}, \end{eqnarray}

where ρ(s) is the reduced density law in Equation (1) in terms of distance relative to the Sun, r, and Galactic coordinates of a given star.

In our calculations, the distance r of a sample star is evaluated by means of its parallax taken from the Hipparcos catalogue, and the integrals in Equation (5) are replaced by two sums, i.e., $\displaystyle \sum _{s=0}^{r} \rho (s)$ and $\displaystyle \sum _{s=0}^\infty \rho (s)$, where the density ρ(s) has been evaluated in steps of 5 pc of the distance. We used the classical equation AV/E(BV) = 3.1 for evaluation of the total absorption, AV, and combined it with the V apparent magnitude of the star in question, taken from the ESA (1997), to estimate its absolute magnitude:

(6) \begin{eqnarray} M_V=V-A_V+5-5\log r. \end{eqnarray}

The absolute magnitude-colour diagram of the star sample is plotted in Figure 1. The position of the RC stars is prominent in this diagram. However, we do not exclude the possibility of contamination of the red giants (RGs). As stated in Nataf, Cassisi, & Athanassoula (Reference Nataf, Cassisi and Athanassoula2014), there is only a small colour difference between the RC stars and the red giant branch (RGB), i.e., $(V - I)_{\text{RGB}}-(V - I)_{\text{RC}}\approx $0.15 mag, for star sample with a given metallicity and luminosity. However, we limited the size of the region occupied by the RC stars to decrease the possible contamination of the RGB stars.

Figure 1. Absolute magnitude-colour diagram for the stars in the re-reduced Hipparcos catalogue. The location of the most probable RC stars is indicated with a rectangular.

We used 9 789 of them in the condense region with colour 0.8 ⩽ (BV)0 ⩽ 1.2 mag and absolute magnitude 0 ⩽ MV ⩽ 2 in our calculations. This constraint is rather useful for our purpose, because our aim in this chapter is to fix a space density for the RC stars in the solar neighbourhood but not to obtain a density distribution for a distance interval. RC stars satisfying this constraint lie within the distances 100 < r ⩽ 300 pc. The space density of these stars has been evaluated by the formula D(r) = N/V(r), where V(r) is the volume between two spheres with radii r 1 = 100 and r 2 = 300 pc, i.e., $V(r) = \frac{4}{3} \pi (r^3_2-r^3_1)$, and N is the number of RC stars in this volume. The logarithmic value of the space density of the RC stars thus evaluated is D* = logD(r) + 10 = 5.954 ± 0.004.

3 SCALE-LENGTH OF THE THIN DISC IN THE GALACTIC PLANE

3.1. Identification of the RC stars in the (J, JW1) CMD

The RC stars used for the estimation of the scale-length of the thin disc in the Galactic plane are taken from the WISE catalogue (Cutri et al. Reference Cutri2013). We plotted the (J, JW1) CMD of the stars with Galactic latitudes − 0°.5 ⩽ b ⩽ +0°.5 in 19 fields with consecutive Galactic longitude intervals, i.e., 85° < l ⩽ 95°, 95° < l ⩽ 105°, . . ., 265° < l ⩽ 275° and identified the RC stars in each field using the Galaxia model (Sharma et al. Reference Sharma, Bland-Hawthorn, Johnston and Binney2011) as explained in the following. We applied the Galaxia model to each field defined in our study and obtained the (J, JW1) CMD for a set of synthetic stars. Then, we adopted the following constrains to separate the RC ones in each field: 1.8 < logg(cgs) ⩽ 2.6 and $4\,500 < T_{\text{eff}} (K) \le 4\,900$. As stated in Section 2, the colour difference between the RC stars and RGs is small. Then, the same case holds also for the temperatures for two sub-samples. Additionally, Nataf et al. (Reference Nataf, Cassisi and Athanassoula2014) claimed the surface gravity of the red giant branch bump (RGBB) as log g ≈ 2.6 (cgs) which is close to that of RC stars we adopted. Hence, we expect a possible contamination of the RG stars in our study.

The CMD for the field with 175° < l ⩽ 185° is plotted in Figure 2 as an example. It turned out that the right mode in the right panel of this figure corresponds to the RC stars. Then, the same case should hold for the observed RC stars, i.e., the right mode in the left panel consists of RC stars we are interested in, while the left one is occupied by the dwarfs. Although the limiting apparent J magnitude of the RC stars in our study is J = 17 mag, we could consider only those with J ⩽ 14.5 mag which corresponds to r ⩽ 8 kpc.

Figure 2. Colour-magnitude diagrams for the observed (left panel) and synthetic (right panel) stars in the field with Galactic longitude 175° < l ⩽ 185° and Galactic latitude − 0°.5 ⩽ b ⩽ +0°.5. Blue solid lines show intrinsic colour, (JW1)0 = 0.665 mag of RC stars (Yaz Gökçe et al. Reference Gökçe, Bilir, Öztürkmen, Duran, Ak, Ak and Karaali2013), light blue-dashed lines display the borders of the most probable RC stars.

Figure 2 shows that there is a gap between the dwarfs and RC stars with apparent magnitude J ⩽ 13.5 mag. We could extend the borders separating the two star categories down to the apparent magnitude J = 14.5 mag by checking the consistency of the space densities of the observed RC stars with the ones of the synthetic RC stars, for a given field.

We applied the constraint ‘data quality flags’ to the WISE catalogue (Cutri et al. Reference Cutri2013) to select the best possible data. Hence, we expect the best possible data and best possible result in scale-length estimation of the thin disc in our study.

3.2. Estimation of the scale-length of the thin disc

The scale-lengths of the thin disc in the Galactic plane have been estimated by fitting the observed space density functions of RC stars in 19 fields to the Galaxy models. The distances of the RC stars relative to the Sun used in the density functions have been evaluated as follows. The E(JW1) colour excess of a RC star has been adopted as the difference between its observed JW1 colour index and the mean true colour index (JW1)0 = 0.665 ± 0.032 mag of the RC sample (Yaz Gökçe et al. Reference Gökçe, Bilir, Öztürkmen, Duran, Ak, Ak and Karaali2013), while the equation for the total absorption, AJ = 1.217 × E(JW1), is adopted from Bilir et al. (Reference Bilir, Karaali, Ak, Daǧtekin, Önal, Yaz, Coşkunoǧlu and Cabrera-Lavers2011). Then, the distance r could be calculated from the following well-known formula:

(7) \begin{eqnarray} J-M_J-A_J=5\log r-5, \end{eqnarray}

here, MJ = −0.970 ± 0.016 mag (Yaz Gökçe et al. Reference Gökçe, Bilir, Öztürkmen, Duran, Ak, Ak and Karaali2013).

We assumed that the RC stars in the Galactic latitude range − 0°.5 ⩽ b ⩽ +0°.5 are all thin disc stars. Hence, the logarithmic space densities D* = logD + 10 evaluated for 19 star fields belong only to the thin disc population of our Galaxy. Here, D = NV 1, 2; $\Delta V_{1,2} = (\frac{\pi }{180})^2 \frac{A}{3}(r^3_2-r^3_1)$; A denotes the size of the field (10 deg2); r 1 and r 2 are the lower and upper limiting distances of the volume ΔV 1, 2; N is the number of RC stars in this volume; r* = [(r 31 + r 32)/2]1/3 is the centroid distance of the volume ΔV 1, 2.

We adopted the density law of Basle group (Buser et al. Reference Buser, Rong and Karaali1998, Reference Buser, Rong and Karaali1999) for the thin disc. Disc structures are usually parametrised in cylindrical coordinates using radial and vertical exponentials:

(8) \begin{eqnarray} D(x,z)=n \times \exp \left(\frac{-|z|}{H}\right) \times \exp \left(\frac{-(x-R_0)}{h}\right), \end{eqnarray}

where z = z + rsinb, r is the distance from the Sun, b is the Galactic latitude, z is the vertical distance of the Sun from the Galactic plane (25 pc; Jurić et al. Reference Jurić2008), x is the projection of the Galactocentric distance on the Galactic plane, R 0 is the solar distance from the Galactic centre (8 kpc), H and h are the scale-height and scale-length, respectively, and n is the normalised local density at the solar radius. In our case, we assume that the space density does not vary in the vertical direction due to the small range of the fields, i.e., − 0°.5 ⩽ b ⩽ +0°.5. Then, assuming that exp(− |z|/H) ≈ 1, Equation (8) can be reduced to the following form:

(9) \begin{eqnarray} D(x)=n \times \exp \left(\frac{-(x-R_0)}{h}\right). \end{eqnarray}

We fitted the space density functions derived from the observations to the density law adopted for the thin disc and estimated the model parameters n and h for each field. We used the classical $\chi _{\text{min}}^2$ statistic for the estimation of the model parameters (Phleps et al. Reference Phleps, Meisenheimer, Fuchs and Wolf2000, Reference Phleps, Drepper, Meisenheimer and Fuchs2005; Chen et al. Reference Chen2001; Siegel et al. Reference Siegel, Majewski, Reid and Thompson2002; Du et al. Reference Du2003, Reference Du, Ma, Wu and Zhou2006; Jurić et al. Reference Jurić2008; Ak et al. Reference Ak, Bilir, Karaali and Buser2007; Bilir et al. Reference Bilir, Cabrera-Lavers, Karaali, Ak, Yaz and López-Corredoira2008; Yaz & Karaali Reference Yaz and Karaali2010). The comparison of the logarithmic density functions derived from the observations and the analytical density law for the thin disc is given in Figure 3 or four fields with Galactic longitude l = 90°, 140°, 180° and 270°, as example, while the results are tabulated in Table 2 for all fields. The errors have been evaluated by changing a Galactic model parameter until an increase or decrease in Δχ2 by 1 was achieved (Phleps et al. Reference Phleps, Meisenheimer, Fuchs and Wolf2000).

Figure 3. Logarithmic space density functions for the most probable RC stars in the fields with Galactic longitudes (a) l = 90°, (b) 140°, (c) 180°, (d) 270°. Filled circles indicate the observed data, while the open ones correspond to the synthetic data in Galaxia.

Table 2. Galactic model parameters for the thin disc estimated by using the observed RC stars in WISE and the synthetic ones in Galaxia for 19 star fields. N indicates the number of star for a given field.

We used the synthetic RC stars identified by the Galaxia model (Sharma et al. Reference Sharma, Bland-Hawthorn, Johnston and Binney2011) and estimated a second set of Galactic model parameters which are also given in Table 2. Also, the logarithmic density functions for the same fields, l = 90°, 140°, 180° and 270°, are presented in Figure 3 as example, as well as for comparison with the corresponding ones of observed RC stars. Table 2 shows that the solar normalisations estimated for two sets of data, i.e., for the observed RC stars and for the synthetic ones, for a given field are compatible. Also, they confirm the solar space density, D*(0) = 5.954, estimated in Section 2 by means of Hipparcos RC stars. We plotted the scale-lengths estimated for the same data in Figure 4, versus Galactic longitudes. Table 2 and Figure 4 show also the agreement of the scale-lengths corresponding two sets of data for a given longitude.

Figure 4. Scale-lengths for the observed (●) and synthetic (⋄) data versus Galactic longitude.

The distribution of the scale-length for the thin disc in Figure 4 is not flat. It makes a peak, h = 2 kpc, at the Galactic longitude l = 180° and it decreases symmetrically in the Galactic longitude intervals 90° ⩽ l < 180° and 180° < l ⩽ 270°, reaching to a minimum value of h = 1.6 kpc at the terminal Galactic longitudes, l = 90° and l = 270°.

4 DISCUSSION

We estimated the scale-length of the thin disc by means of the space density functions of the RC stars in 19 star fields with consecutive Galactic longitudes, l = 90°, 100°, . . ., 270° and with Galactic latitudes − 0°.5 ⩽ b ⩽ +0°.5. We used the solar space density of the RC stars in the Hipparcos catalogue, logD(0) + 10 = 5.954 ± 0.004, evaluated in this study in fitting the observed space density functions to the Galaxy models. We used the data in the WISE catalogue (Cutri et al. Reference Cutri2013) and identified the RC stars by comparing their positions in the (J, JW1) CMDs with the synthetic ones in Galaxia (private communication and Sharma et al. Reference Sharma, Bland-Hawthorn, Johnston and Binney2011). The RC stars are adopted as the objects which satisfy the following conditions: (i) 1.8 < logg(cgs) ⩽ 2.6 and (ii) $4\,500 < T_{\text{eff}} (K) \le 4900$. Then, the right mode in the (J, JW1) CMD consists of RC stars, while the left one belongs to the dwarfs (Figure 2). We used the gap between the dwarfs and RC stars with apparent magnitudes J ⩽ 13.5 mag as a clue and separated these star categories down to the apparent magnitude J = 14.5 mag by checking the consistency of the space densities of the observed RC stars with the ones of the synthetic RC stars, for a given field.

The scale-lengths estimated for the thin disc with the RC stars in different fields are not equal to each other. The largest scale-length, h = 2137+ 114− 58 pc, corresponds to the Galactic longitude l = 180°. The scale-lengths estimated for the fields with Galactic longitudes close to l = 180°, i.e., l = 150°, 160°, 170°, 190° are rather close to h = 2137+ 114− 58 pc. However, they deviate from this value systematically in the fields with Galactic longitudes l ⩽ 150° and l > 190°, and they make two flat distributions at the two terminal longitudes, l = 90° and l = 270°. The smallest scale-length, h = 1574+ 94− 141 pc, in the second quadrant belongs to the field with Galactic longitude l = 110°, while in the third quadrant it corresponds to l = 250°, h = 1549+ 119− 153 pc. The distribution of the observed scale-length is confirmed by the synthetic RC stars, hence it is real. We estimated the scale-length of the thin disc for three fields, l = 90°, 180°, and 270° by adopting a scale-height, H = 350 pc (Siegel et al. Reference Siegel, Majewski, Reid and Thompson2002) and compared them with the ones in Table 2. The differences for two sets of data are less than 0.1 kpc for three fields, which are within the errors given in Table 2.

The difference between the scale-lengths estimated for different longitudes can be explained as follows. The first reason is that in the Galactic anticentre direction, we reach in average higher values of Galactocentric distance (R). A second explanation can be done by the overdensities in the direction perpendicular to the Galactic anticentre direction or close to this direction, such as the fields with Galactic longitudes l = 90°, 100°, 110°, 120°, 240°, 250°, 260°, 270°. The scale-lengths in the fields with such Galactic longitudes are relatively small, while small number of stars in the fields in the anticentre direction exhibit larger scale-lengths (Table 2). Additionally, Figure 5 shows that the distance histograms for three Galactic longitudes, l = 90°, 180°, and 270° are different. The one for the longitude l = 180° is flat, while the number of RC stars for the Galactic longitudes l = 90° and 270° are monotonously increasing with distance. As the volume increases with distance, the flat distribution results a decreasing space density for the Galactic longitude l = 180° (Figure 3(c)), while the increasing number of stars for the Galactic longitudes l = 90° and 270° supply relatively flat space density functions (Figure 3(a)).

Figure 5. Distance histograms for the most probable RC stars in three different star fields: (a) l = 90°, (b) l = 180°, and (c) l = 270°.

The scale-length for the thin disc appeared in the literature lies in the range 2 ⩽ h ⩽ 4 kpc, the least ones being those of Robin et al. (Reference Robin, Haywood, Crézé, Ojha and Bienayme1996); Ojha et al. (Reference Ojha, Bienaymé, Mohan and Robin1999); Robin, Reylé, & Crézé (Reference Robin, Reylé and Crézé2000); Chen et al. (Reference Chen2001); Siegel et al. (Reference Siegel, Majewski, Reid and Thompson2002), i.e., 2 ⩽ h ⩽ 2.5 kpc. The most recent scale-length is claimed by López-Corredoira & Molgó (Reference López-Corredoira and Molgó2014), h = 2 kpc. Their observations were carried out in an off-plane low Galactic regions with size 1400 deg2 and they were fitted to a model of flared thin and thick disc. If we normalise the scale-length, estimated for the RC stars in the fields in the opposite direction of the Galactic centre or close to this direction, to one decimal we obtain exactly the same value, h = 2 kpc. This is a good result which can be used as an argument to claim that the value h = 2 kpc can be adopted as a standard scale-length for the thin disc for the regions of the Galaxy which are not contaminated by the accreted material or other features such as warp and flare.

In this study, we did not consider the effect of the contamination of the RGs. However, this contamination shows its effect more or less in the studies related to the RC stars due to several reasons. First, the RC stars and the RGB fall in close colour intervals, i.e., the difference between their intrinsic colours in (VI)0 is only ≈ 0.15 mag (Nataf et al. Reference Nataf, Cassisi and Athanassoula2014) which complicates their separation due to the errors in colour estimation and interstellar reddening. Second, the number of RC stars is larger than the one of RGs in a given sample, and their ratio is an increasing function of metallicity. Nataf et al. (Reference Nataf, Cassisi and Athanassoula2014) give this ratio as 1.75 and 2.32 for the metallicities [M/H] = −1.27 and [M/H] = 0.40 dex, respectively. Finally, as cited in Nataf et al. (Reference Nataf, Cassisi and Athanassoula2014) the fainter RC stars and the RGBB overlap at the distance modulus 14.80 mag (their Figure 2).

In our study, the probable effect of the first reason is unvoidable, while we expect minimum effect due to the second reason, because our RC sample covers only the disc stars with solar metallicity in the Galactic plane. The distance modulus 14.80 mag just mentioned corresponds to the apparent magnitude J 0 ≈ 13.80 mag in our study which indicates that only the faintest stars, 13.50 < J 0 < 14.50 mag, can be affected from the overlapping in question. However, the faint segments of our space density functions show a smooth extension relative to the bright ones. Hence, we do not expect a large contamination due to this reason. All these statements are valid for the Galaxia data, because they are based on a set of atmospheric parameters which can be valid for RC stars as well as RGs.

5 CONCLUSIONS

The scale-length of the thin disc estimated with the most probable RC stars in the Galactic plane is 2 kpc for the fields in the Galactic anticentre direction or close to this direction, while it decreases continuously in the second and third quadrants reaching to a lower limit of h = 1.6 kpc at the Galactic longitudes l = 90° and l = 270°.

ACKNOWLEDGEMENTS

We would like to thank anonymous referee for the useful comments and corrections. Also, we would like to thank Dr. Sanjib Sharma for providing an updated version of Galaxia model. This study has been supported in part by the Scientific and Technological Research Council (TÜBİTAK) 112T120. Part of this work was supported by the Research Fund of the University of İstanbul, Project Number Numbers: 48483. MLC was supported by the Spanish Ministry of Economy and Competitiveness (MINECO) under the grant AYA2012-33211. This research has made use of NASA’s Astrophysics Data System Bibliographic Services, the VizieR catalogue access tool, CDS, Strasbourg, France and NASA/IPAC Infrared Science Archive and Extragalactic Database (NED), which are operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This publication makes use of data products from the Wide-field Infrared Survey Explorer and the Two Micron All Sky Survey, which are a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.

References

REFERENCES

Ak, S., Bilir, S., Karaali, S., & Buser, R. 2007, AN, 328, 169Google Scholar
Alves, D. R. 2000, ApJ, 539, 732CrossRefGoogle Scholar
Arce, H. G., & Goodman, A. A. 1999, ApJ, 512, 135CrossRefGoogle Scholar
Beckman, J. E., Peletier, R. F., Knapen, J. H., Corradi, R. L. M., & Gentet, L. J. 1996, ApJ, 467, 175CrossRefGoogle Scholar
Bilir, S., Cabrera-Lavers, A., Karaali, S., Ak, S., Yaz, E., & López-Corredoira, M. 2008, PASA, 25, 69CrossRefGoogle Scholar
Bilir, S., Karaali, S., Ak, S., Daǧtekin, N. D., Önal, Ö., Yaz, E., Coşkunoǧlu, B., & Cabrera-Lavers, A. 2011, MNRAS, 417, 2230Google Scholar
Bilir, S., Ak, T., Ak, S., Yontan, T., & Bostancı, Z. F. 2013b, NewA, 23, 88CrossRefGoogle Scholar
Bilir, S., Karaali, S., Ak, S., Yaz, E., & Hamzaoğlu, E. 2006b, NewA, 12, 234Google Scholar
Bilir, S., Karaali, S., Dağtekin, N. D., Önal, Ö., Ak, S., Ak, T., & Cabrera-Lavers, A. 2012, PASA, 29, 121CrossRefGoogle Scholar
Bilir, S., Karaali, S., & Gilmore, G. 2006a, MNRAS, 366, 1295CrossRefGoogle Scholar
Bilir, S., Karaali, S., Güver, T., Karataş, Y., & Ak, S. G. 2006c, AN, 327, 72Google Scholar
Bilir, S., Önal, Ö., Karaali, S., Cabrera-Lavers, A., & Çakmak, H. 2013a, Ap&SS, 344, 417Google Scholar
Binney, J., et al. 2014, MNRAS, 437, 351CrossRefGoogle Scholar
Bovy, J., Rix, H.-W., Liu, C., Hogg, D. W., Beers, T. C., & Lee, Y. S. 2012, ApJ, 753, 148CrossRefGoogle Scholar
Buser, R., Rong, J., & Karaali, S. 1998, A&A, 331, 934Google Scholar
Buser, R., Rong, J., & Karaali, S. 1999, A&A, 348, 98Google Scholar
Cabrera-Lavers, A., Bilir, S., Ak, S., Yaz, E., & López-Corredoira, M. 2007a, A&A, 464, 565Google Scholar
Cabrera-Lavers, A., Garzón, F., & Hammersley, P. L. 2005, A&A, 433, 173Google Scholar
Cabrera-Lavers, A., González-Fernández, C., Garzón, F., Hammersley, P. L., & López-Corredoira, M. 2008, A&A, 491, 781Google Scholar
Cabrera-Lavers, A., Hammersley, P. L., González-Fernández, C., López-Corredoira, M., Garzón, F., & Mahoney, T. J. 2007b, A&A, 465, 825Google Scholar
Chang, C.-K., Ko, C.-M., & Peng, T.-H. 2011, ApJ, 740, 34CrossRefGoogle Scholar
Chen, B., et al. 2001, ApJ, 553, 184CrossRefGoogle Scholar
Cheng, J. Y., et al. 2012, ApJ, 752, 51CrossRefGoogle Scholar
Cutri, R. M., et al. 2013. WISE All-Sky Data Release, CDS/ADC Electronic Catalogues, 2328Google Scholar
Drimmel, R., & Spergel, D. N. 2001, ApJ, 556, 181CrossRefGoogle Scholar
Du, C., et al. 2003, A&A, 407, 541Google Scholar
Du, C., Ma, J., Wu, Z., & Zhou, X. 2006, MNRAS, 372, 1304CrossRefGoogle Scholar
ESA 1997, The Hipparcos and Tycho Catalogues, CDS/ADC Electronic Catalogues, 1239Google Scholar
Gilmore, G. 1984, MNRAS, 207, 223CrossRefGoogle Scholar
Gilmore, G., & Reid, N. 1983, MNRAS, 202, 1025CrossRefGoogle Scholar
Grocholski, A. J., & Sarajedini, A. 2002, AJ, 123, 1603CrossRefGoogle Scholar
Groenewegen, M. A. T. 2008, A&A, 488, 935Google Scholar
Jurić, M., et al. 2008, ApJ, 673, 864CrossRefGoogle Scholar
Kalberla, P. M. W., Kerp, J., Dedes, L., & Haud, U. 2014, ApJ, 794, 90CrossRefGoogle Scholar
Karaali, S., Bilir, S., & Hamzaoǧlu, E. 2004, MNRAS, 355, 307Google Scholar
Karaali, S., Bilir, S., Yaz, E., Hamzaoǧlu, E., & Buser, R. 2007, PASA, 24, 208CrossRefGoogle Scholar
Karaali, S., Bilir, S., & Yaz Gökçe, E. 2013, Ap&SS, 346, 89Google Scholar
Keenan, P. C., & Barnbaum, C. 1999, ApJ, 518, 859CrossRefGoogle Scholar
Kubiak, M., McWilliam, A., Udalski, A., & Gorski, K. 2002, AcA, 52, 159Google Scholar
Laney, C. D., Joner, M. D., & Pietrzyski, G. 2012, MNRAS, 419, 1637CrossRefGoogle Scholar
Larsen, J. A. 1996, PhD thesis, University of MinnesotaGoogle Scholar
Larsen, J. A., & Humphreys, R. M. 2003, AJ, 125, 1958CrossRefGoogle Scholar
López-Corredoira, M., Cabrera-Lavers, A., Garzón, F., & Hammersley, P. L. 2002, A&A, 394, 883Google Scholar
López-Corredoira, M. 2006, MNRAS, 369, 1911CrossRefGoogle Scholar
López-Corredoira, M., Cabrera-Lavers, A., Gerhard, O. E., & Garzón, F. 2004, A&A, 421, 953Google Scholar
López-Corredoira, M., & Molgó, J. 2014, A&A, 567A, 106Google Scholar
López-Corredoira, M., Momany, Y., Zaggia, S., & Cabrera-Lavers, A. 2007, A&A, 472, L47Google Scholar
Marshall, D. J., Robin, A. C., Reylé, C., Schultheis, M., & Picaud, S. 2006, A&A, 453, 635Google Scholar
Martin, N. F., Ibata, R. A., Conn, B. C., Lewis, G. F., Bellazzini, M., Irwin, M. J., & McConnachie, A. W. 2004, MNRAS, 355, 33CrossRefGoogle Scholar
Momany, Y., Zaggia, S., Gilmore, G., Piotto, G., Carraro, G., Bedin, L. R., & de Angeli, F. 2006, A&A, 451, 515Google Scholar
Nataf, D. M., Cassisi, S., & Athanassoula, E. 2014, MNRAS, 442, 2075CrossRefGoogle Scholar
Ojha, D. K. 2001, MNRAS, 322, 426CrossRefGoogle Scholar
Ojha, D. K., Bienaymé, O., Mohan, V., & Robin, A. C. 1999, A&A, 351, 945Google Scholar
Ojha, D. K., Bienayme, O., Robin, A. C., Creze, M., & Mohan, V. 1996, A&A, 311, 456Google Scholar
Paczyński, B., & Stanek, K. Z. 1998, ApJ, 494, 219CrossRefGoogle Scholar
Phleps, S., Drepper, S., Meisenheimer, K., & Fuchs, B. 2005, A&A, 443, 929Google Scholar
Phleps, S., Meisenheimer, K., Fuchs, B., & Wolf, C. 2000, A&A, 356, 108Google Scholar
Polido, P., Jablonski, F., & Lépine, J. R. D. 2013, ApJ, 778, 32Google Scholar
Robin, A. C., Haywood, M., Crézé, M., Ojha, D. K., & Bienayme, O. 1996, A&A, 305, 125Google Scholar
Robin, A. C., Reylé, C., & Crézé, M. 2000, A&A, 359, 103Google Scholar
Salaris, M., & Girardi, L. 2002, MNRAS, 337, 332CrossRefGoogle Scholar
Sarajedini, A. 1999, AJ, 118, 2321CrossRefGoogle Scholar
Schlafly, E. F., & Finkbeiner, D. P. 2011, ApJ, 737, 103CrossRefGoogle Scholar
Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525CrossRefGoogle Scholar
Sharma, S., Bland-Hawthorn, J., Johnston, K. V., & Binney, J. 2011, ApJ, 730, 3CrossRefGoogle Scholar
Siegel, M. H., Majewski, S. R., Reid, I. N., & Thompson, I. B. 2002, ApJ, 578, 151Google Scholar
Skrutskie, M. F.et al. 2006 AJ, 131, 1163CrossRefGoogle Scholar
Stanek, K. Z., & Garnavich, P. M. 1998, ApJ, 494, 219Google Scholar
Tuffs, R. J., Popescu, C. C., Vlk, H. J., Kylafis, N. D., & Dopita, M. A. 2004, A&A, 419, 821Google Scholar
Udalski, A. 1998, AcA, 48, 383Google Scholar
van Helshoecht, V., & Groenewegen, M. A. T. 2007, A&A, 463, 559Google Scholar
van Leeuwen, F. 2007, A&A, 474, 653Google Scholar
Wright, E. L., et al. 2010, AJ, 140, 1868CrossRefGoogle Scholar
Xu, Y., Voronkov, M. A., Pandian, J. D., Li, J. J., Sobolev, A. M., Brunthaler, A., Ritter, B., Menten, K. M. 2009, A&A, 507, 1117Google Scholar
Yamagata, T., Yoshii, Y. 1992, AJ, 103, 117Google Scholar
Yaz, E., Karaali, S. 2010, NewA, 15, 234CrossRefGoogle Scholar
Yaz Gökçe, E., Bilir, S., Öztürkmen, N. D., Duran, Ş., Ak, T., Ak, S., Karaali, S. 2013, NewA, 25, 19CrossRefGoogle Scholar
Zhao, G., Qiu, H., Zhang, H. 2001, ChA&A, 25, 59Google Scholar
Figure 0

Table 1. Scale-height and scale-length values of thin disc in the literature.

Figure 1

Figure 1. Absolute magnitude-colour diagram for the stars in the re-reduced Hipparcos catalogue. The location of the most probable RC stars is indicated with a rectangular.

Figure 2

Figure 2. Colour-magnitude diagrams for the observed (left panel) and synthetic (right panel) stars in the field with Galactic longitude 175° < l ⩽ 185° and Galactic latitude − 0°.5 ⩽ b ⩽ +0°.5. Blue solid lines show intrinsic colour, (JW1)0 = 0.665 mag of RC stars (Yaz Gökçe et al. 2013), light blue-dashed lines display the borders of the most probable RC stars.

Figure 3

Figure 3. Logarithmic space density functions for the most probable RC stars in the fields with Galactic longitudes (a) l = 90°, (b) 140°, (c) 180°, (d) 270°. Filled circles indicate the observed data, while the open ones correspond to the synthetic data in Galaxia.

Figure 4

Table 2. Galactic model parameters for the thin disc estimated by using the observed RC stars in WISE and the synthetic ones in Galaxia for 19 star fields. N indicates the number of star for a given field.

Figure 5

Figure 4. Scale-lengths for the observed (●) and synthetic (⋄) data versus Galactic longitude.

Figure 6

Figure 5. Distance histograms for the most probable RC stars in three different star fields: (a) l = 90°, (b) l = 180°, and (c) l = 270°.