A publishing partnership

Articles

NEW EVIDENCE SUPPORTING CLUSTER MEMBERSHIP FOR THE KEYSTONE CALIBRATOR DELTA CEPHEI

, , and

Published 2012 February 24 © 2012. The American Astronomical Society. All rights reserved.
, , Citation D. Majaess et al 2012 ApJ 747 145 DOI 10.1088/0004-637X/747/2/145

0004-637X/747/2/145

ABSTRACT

New and existing UBVJHKs, spectroscopic, NOMAD, Hubble Space Telescope, and revised Hipparcos observations are employed to determine properties for δ Cep and its host star cluster. The multi-faceted approach ensured that uncertainties were mitigated (σ/d ∼ 2%). The following fundamental parameters were inferred for δ Cep: E(BV) = 0.073 ± 0.018(σ), log τ = 7.9 ± 0.1, and $d=272\,\pm\, 3(\sigma _{\bar{x}})\, \pm\, 5 (\sigma)$ pc. The cluster exhibits a turnoff near B6 (M*/M ∼ 5), and the brightest host cluster members are the supergiants ζ Cep (K1.5Ib) and δ Cep. To within the uncertainties, the two stars share common astrometric parameters (π, μα, μδ, RV ∼ −17 km s−1) and are tied to bluer members via the evolutionary track implied by the cluster's UBVJHKs color–color and color–magnitude diagrams. The cluster's existence is bolstered by the absence of an early-type sequence in color–magnitude diagrams for comparison fields. NOMAD data provided a means to identify potential cluster members (n ∼ 30) and double the existing sample. That number could increase with forthcoming precise proper motions (DASCH) for fainter main-sequence stars associated with classical Cepheids (e.g., δ Cep), which may invariably foster efforts to strengthen the Galactic Cepheid calibration and reduce uncertainties tied to H0.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

Cepheid variables are utilized to establish extragalactic distances and constrain cosmological models (Macri & Riess 2009; Shappee & Stanek 2011; Riess et al. 2011). However, the reliability of the derived parameters is invariably tied to the Cepheid calibration. Freedman et al. (2001) noted that ambiguities related to the zero point of the calibration account for a sizable fraction of the total uncertainty associated with H0. The uncertainty hinders efforts to constrain dark energy, since the parameter is acutely dependent on an accurate Hubble constant ($\sigma _{\rm w}\sim 2 \sigma _{\rm H_0}$; Macri & Riess 2009). The next generation follow-up to the Hubble Space Telescope (HST) key project to measure H0 (the Hubble Carnegie project; Freedman & Madore 2010) aims to mitigate that problem by relying on LMC and Galactic calibrators (Benedict et al. 2002, 2007; Turner 2010; Storm et al. 2011). Consequently, bolstering the Galactic calibration should support efforts by the Hubble Carnegie, Araucaria, and SH0ES projects to determine H0 to within 2%–4% (Gieren et al. 2005; Riess et al. 2011).

In this study, UBVJHKs, spectroscopic, NOMAD, HST, and revised Hipparcos (HIP) observations for stars physically associated with δ Cep are employed to constrain its fundamental parameters: age (log τ), color excess (EBV), distance, progenitor mass, and absolute Wesenheit magnitude ($W_{VI_c,0}$).

2. ANALYSIS

2.1. Revised Hipparcos Observations for Cep OB6

In their comprehensive study, de Zeeuw et al. (1999) discovered that δ Cep was a member of a group subsequently denoted as Cep OB6 (see also Hoogerwerf et al. 1997). Twenty stars identified by de Zeeuw et al. (1999) as Cep OB6 members are highlighted in Table 1. Stars exhibiting spectral types inconsistent with cluster membership based on their UBVJHKs color–color and color–magnitude positions were flagged as non-members. For example, HIP 110459 was previously assigned a membership probability of 100% (Table 1), yet the star exhibits JHKs photometry and a late-type temperature class (K5; Skiff 2010) indicative of a field red clump giant (Figure 2). UBVJHKs photometry was obtained from Mermilliod (1991) and Two Micron All Sky Survey (2MASS; Cutri et al. 2003). Spectral types were assigned to stars featured in the Catalogue of Stellar Spectral Classifications (Skiff 2010).

Table 1. Cep OB6 Member List (Z99)

HIP ID HIP πa V07 πa μα, μδb m.p. (Z99)c m.p.d E(BV)e
  (mas) (mas) (mas yr−1)      
109426 3.8 ± 0.7 3.6 ± 0.5 16.6 ± 0.6, 4.1 ± 0.5 94 m  ⋅⋅⋅
109492 4.5 ± 0.5 3.9 ± 0.1 13.3 ± 0.4, 4.4 ± 0.3 97 m  ⋅⋅⋅
110266 3.9 ± 0.6 4.4 ± 0.3 19.0 ± 0.5, 5.1 ± 0.5 96 m 0.070
110275 4.0 ± 1.0 4.0 ± 0.9 14.7 ± 1.0, 5.7 ± 0.9 89 m  ⋅⋅⋅
110356 3.4 ± 0.7 2.9 ± 0.5 11.7 ± 0.6, 3.2 ± 0.6 100 m 0.085
110497 3.8 ± 0.6 3.2 ± 0.4 17.4 ± 0.5, 4.8 ± 0.5 98 m 0.060
110648 3.9 ± 1.0 3.3 ± 0.9 16.2 ± 1.0, 6.2 ± 0.8 84 m  ⋅⋅⋅
110807 4.0 ± 0.6 3.5 ± 0.4 16.1 ± 0.5, 5.4 ± 0.4 92 m 0.060
110925 4.3 ± 0.9 5.1 ± 0.8 21.3 ± 1.1, 4.8 ± 0.8 86 m 0.060f
110988 3.4 ± 0.6 3.7 ± 0.5 16.4 ± 0.7, 4.7 ± 0.7 100 m 0.085f
111060 5.0 ± 0.8 5.3 ± 0.7 17.2 ± 0.7, 4.5 ± 0.7 100 m  ⋅⋅⋅
112141 3.2 ± 0.8 3.4 ± 0.6 14.6 ± 0.7, 2.1 ± 0.6 89 m  ⋅⋅⋅
113255 4.3 ± 0.7 4.5 ± 0.5 19.8 ± 0.6, 3.0 ± 0.6 99 m 0.095
113316 3.2 ± 0.7 3.6 ± 0.4 14.1 ± 0.7, 3.2 ± 0.6 99 m 0.100
δ Cep 3.3 ± 0.6 3.8 ± 0.2 16.4 ± 0.6, 3.5 ± 0.6 89 m  ⋅⋅⋅
110459 4.1 ± 0.9 4.5 ± 0.7 16.2 ± 0.9, 5.2 ± 0.7 100 nm  ⋅⋅⋅
111069 3.2 ± 0.9 3.1 ± 0.8 15.1 ± 0.8, 6.9 ± 0.7 79 nm  ⋅⋅⋅
112473 3.6 ± 0.8 5.0 ± 0.7 13.9 ± 0.8, 2.8 ± 0.8 97 nm  ⋅⋅⋅
112998 2.5 ± 0.6 2.6 ± 0.3 12.7 ± 0.6, 2.2 ± 0.5 98 nm  ⋅⋅⋅
113993 3.8 ± 0.7 3.7 ± 0.4 14.3 ± 0.6, 3.9 ± 0.5 79 nm  ⋅⋅⋅

Notes. aParallaxes from Perryman et al. (1997, HIP) and van Leeuwen (2007, V07). bProper-motion data from NOMAD (Zacharias et al. 2004). cMembership probability assigned by de Zeeuw et al. (1999, Z99). dMembership inferred from UBVJHKs and spectroscopic observations (Figures 1 and 2). eReddenings derived from the UBV color–color analysis (Figure 1). fStars in close proximity to δ Cep.

Download table as:  ASCIITypeset image

Revised HIP parallaxes (van Leeuwen 2007) were tabulated for 15 stars in Table 1 which were identified as probable cluster members. The revised HIP parallaxes exhibit a ∼30% reduction in uncertainties relative to existing data (Perryman & ESA 1997), and the parallaxes for δ Cep and HD 213307 (r ∼ 0farcm7) were increased from π = 3.32  ±  0.58:3.43  ±  0.64 to π = 3.77  ±  0.16:3.69  ±  0.46 mas. A mean distance for the revised cluster sample outlined in Table 1 is $d=271\,{\pm}\, 11(\sigma _{\bar{x}})\,{\pm}\, 42(\sigma)$ pc (see also de Zeeuw et al. 1999, and their Appendix B).

2.2. Reddening

UBVJHKs color–color analyses permit an assessment of the sample's extinction properties. UBV data are particularly efficient at identifying early-type stars owing to the U-band's sensitivity to the Balmer decrement (e.g., Turner 1989; Carraro et al. 2006).

Cluster members featuring UBV photometry (Mermilliod 1991) are offset from the intrinsic UBV color–color relation by E(BV) = 0.077 ± 0.016(σ) (Table 1 and Figure 1). The findings support the field reddening determined for δ Cep by Benedict et al. (2002) and are consistent with those established from spectroscopic and JHKs observations. However, uncertainties associated with the latter hamper a precise assessment. The minimal spread, E(BV) ∼ 0.06–0.10 (Table 1), may be indicative of marginal differential reddening, rotation, binarity, or photometric uncertainties.

Figure 1.

Figure 1. Color–color diagram for stars in Table 1 which are associated with δ Cep and possess UBV photometry. The sample is offset from the intrinsic relation (dotted line) by E(BV) = 0.077 ± 0.016(σ) (solid line). The result confirms the reddening established for δ Cep by Benedict et al. (2002). The intrinsic relation and reddening law for the region were adopted from Turner (1976, 1989).

Standard image High-resolution image

2.3. Age

UBVJHKs color–color analyses imply a cluster turnoff near B5–B7 (M*/M ∼ 5), according to intrinsic relations (Padova models; Straižys & Lazauskaitė 2009; Turner 1976, 1989, 2011). A cluster age of log τ = 7.9 ± 0.1 was inferred from the turnoff determined (Figure 1 and Table 1), and since the corresponding Padova evolutionary track (Figure 2) matches bluer and redder evolved members (δ and ζ Cep). The result agrees with the Cepheid's predicted age (Turner 1996; Bono et al. 2005). The temporal match is pertinent evidence in support of cluster membership for δ Cep.

Figure 2.

Figure 2. Left to right: JH and BV color–magnitude diagrams for the Table 1 and NOMAD samples (panels 1 and 2), and comparison fields (panels 3 and 4). Small dots denote calibration stars from Majaess et al. (2011a), which were employed to tie the cluster distance to a geometrically anchored scale (van Leeuwen 2009; Majaess et al. 2011a). Large dots characterize stars exhibiting μα = 11–19 and μδ = 2–7 mas yr−1. Open circles are likely field stars (Table 1). HIP 110459 (circled dot) was previously identified as a cluster member, yet BVJHKs photometry implies the object is a field star. Panel 2: a Padova log τ = 7.9 isochrone was applied. The brightest cluster members are likely the supergiants ζ Cep (K1.5Ib) and δ Cep (amplitude variation indicated). An early-type cluster sequence is absent from the comparison fields (HIP data for the cluster μαδ).

Standard image High-resolution image

2.4. Cluster Distance

A precise distance may be established since two of four principal parameters associated with isochrone fitting were constrained by the UBVJHKs color–color and spectroscopic analyses, namely, the reddening and age (spectral type at the turnoff). δ Cep exhibits solar abundances, and hence the remaining parameter is the shift required in magnitude space to overlay the intrinsic relation upon the data. The resulting distance is d = 277 ± 15 pc (Figure 2). The zero point is tied to seven benchmark open clusters (d < 250 pc) which exhibit matching JHKs and revised Hipparcos distances (the Hyades, α Per, Praesepe, Coma Ber, IC 2391, IC 2609, and NGC 2451; van Leeuwen 2009; Majaess et al. 2011a). A redetermination of the HST parallax for the Hyades supports that scale (McArthur et al. 2011). Isochrones, models, and the distance scale should be anchored (and evaluated) using clusters where consensus exists, rather than the discrepant case (i.e., the Pleiades). A ratio of total to selective extinction RJ was adopted from Majaess et al. (2011b; see also Bonatto et al. 2004), whereas a value for RV was adopted from Turner (1976). An advantage of employing JHKs observations is that the cluster reddening is negligible in that part of the spectrum (EJH ∼ 0.3 × EBV, Majaess et al. 2008; Bonatto et al. 2004, and references therein), which consequently mitigates the impact of uncertainties in Rλ (J0 = JEJH × RJ).

The distance derived from the cluster color–magnitude diagrams (Figure 2) is tied to additional potential cluster members identified using NOMAD (Zacharias et al. 2004). That database was queried for stars exhibiting μα = 11–19 and μδ = 2–7 mas yr−1 (Table 1). Stars fainter than J ∼ 9.2 were culled from the resulting sample to reduce field contamination. Proper motions may be less reliable for such stars, and spectroscopic and UBV observations are typically unavailable. Stars redder than JH ∼ 0.4 were likewise removed to mitigate contamination from field red clump giants. In addition, stars featuring anomalous positions in the multiband color–color and color–magnitude diagrams were removed. The remaining stars double the number of existing potential cluster members and are highlighted in Table 2. The stars in Tables 1 and 2 are potential members pending further evidence. The NOMAD proper motions are consistent with estimates from the PPMXL catalog (Table 2). UBVJHKs photometry was sought from Mermilliod (1991), 2MASS (Cutri et al. 2003), and pertinent resources. For example, new observations acquired from the Bright Star Monitor, which is part of the AAVSO's robotic telescope network, provided BV photometry for HD 239949: V = 10.013 ± 0.031 and BV = 0.392 ± 0.048. The cluster reddening was redetermined (EBV = 0.073 ± 0.018(σ)) after including four earlier type stars highlighted in Table 2 which possess UBV photometry.

Table 2. Additional Potential Cluster Members

ID V B − V U − B J H Ks μα, μδa μα, μδb
              (mas yr−1) (mas yr−1)
HD 210071c 6.39 −0.10 −0.45 6.49 6.57 6.58 16.3 ± 0.4, 2.5 ± 0.4 16.3 ± 0.4, 2.4 ± 0.4
HD 209636c 7.01 −0.05 −0.23 7.10 7.17 7.16 15.0 ± 0.5, 2.7 ± 0.5 15.0 ± 0.5, 2.6 ± 0.5
HD 214259 8.72 0.15 0.10 8.27 8.24 8.24 18.1 ± 1.6, 3.7 ± 1.6 18.4 ± 1.3, 4.3 ± 1.3
HD 214512 8.74 0.15  ⋅⋅⋅ 8.30 8.29 8.29 17.6 ± 1.8, 6.0 ± 1.8 −21.5 ± 1.6, 53.1 ± 1.6  
HD 212093 8.25 −0.02 −0.29 8.31 8.36 8.35 15.0 ± 1.6, 3.3 ± 1.6 14.1 ± 1.2, 2.5 ± 1.2
HD 210480 8.71 0.15 0.08 8.41 8.35 8.31 14.5 ± 1.6, 2.2 ± 1.6 14.9 ± 1.3, 2.3 ± 1.3
HD 211226 8.65 0.07 0.03 8.45 8.49 8.47 17.3 ± 1.3, 6.9 ± 1.3 15.4 ± 1.2, 6.7 ± 1.2
HD 215879 8.98 0.12  ⋅⋅⋅ 8.65 8.62 8.59 13.8 ± 1.6, 2.6 ± 1.6 13.4 ± 1.3, 4.4 ± 1.2
HD 212711 9.25 0.23  ⋅⋅⋅ 8.78 8.72 8.66 17.3 ± 2.2, 3.6 ± 2.1 16.7 ± 1.9, 5.2 ± 1.9
HD 240052 9.44 0.30  ⋅⋅⋅ 8.79 8.74 8.71 18.1 ± 1.6, 2.8 ± 1.6 17.3 ± 1.3, 4.1 ± 1.2
HD 212137 9.19 0.09  ⋅⋅⋅ 8.95 8.95 8.94 12.4 ± 1.1, 2.8 ± 1.1 12.9 ± 1.2, 1.8 ± 1.2
BD+54°2675 9.48 0.28  ⋅⋅⋅ 9.04 9.01 8.94 11.9 ± 1.6, 4.5 ± 1.6 12.2 ± 1.8, 7.5 ± 1.7
HD 239949 10.01 0.39  ⋅⋅⋅ 9.11 9.10 9.06 16.8 ± 2.5, 4.5 ± 2.3    5.4 ± 1.6, −0.8 ± 1.5
BD+59°2523 9.75 0.27  ⋅⋅⋅ 9.16 9.07 9.11 14.8 ± 2.3, 3.3 ± 2.2 15.3 ± 1.6, 4.2 ± 1.6

Notes. aNOMAD proper motions (Zacharias et al. 2004). bThe PPMXL catalog (Roeser et al. 2010). cπ = 5.06 ± 0.33: 5.54 ± 0.39 mas (van Leeuwen 2007).

Download table as:  ASCIITypeset image

Cluster members appear to aggregate near J2000 coordinates of $22^{\rm h} 22\mbox{$.\!\!^{\mathrm m}$}5 +56\mbox{$^\circ $}34\mbox{$^\prime $}$ (Figure 3). δ Cep lies at the periphery of the density enhancement, and within the confines of the cluster since the corona extends further (Kholopov 1969; Turner 1985). The Cepheid is r ∼ 2° from the aforementioned coordinates, which is equivalent to a linear projected separation of ∼9 pc. The revised HIP and HST parallaxes for δ Cep and HD 213307 (r ∼ 0farcm7 from δ Cep), in tandem with their apparent positions, suggest that the distance to the cluster center and Cepheid is analogous to within the uncertainties.

Figure 3.

Figure 3. Top: J2000 R.A./decl. positions for probable (squares) and non-members (open circles) highlighted in Table 1, and new potential members (filled circles) outlined in Table 2. Bottom: dotted open circles represent all HIP stars near the (approximate) cluster center.

Standard image High-resolution image

2.5. Mean Distance to δ Cep

The mean HIP parallax for the cluster (Section 2.1) agrees with the HIP parallax for δ Cep (π = 3.77 ± 0.16 mas), the HST parallax for δ Cep (π = 3.66 ± 0.15 mas; Benedict et al. 2002), and the distance inferred for the host cluster from UBVJHKs and spectroscopic observations (Figures 1 and 2). Assigning equal weight to each method, the mean of the four distances for δ Cep is: $d=272\pm 3(\sigma _{\bar{x}}) \pm 5 (\sigma)$ pc. That result agrees with the Storm et al. (2011) determination from the infrared surface brightness technique (IRSB; Fouque & Gieren 1997). The associated standard error and deviation provide a realistic estimate for the systematic uncertainty, which is often difficult to characterize.

The resulting VIc Wesenheit magnitude for δ Cep is $W_{VI_c,0}=-5.12$ ($R_{VI_c}=2.55$). That is consistent with results established for classical Cepheids displaying similar pulsation periods: CV Mon, V Cen, Y Sgr, and CS Vel (Benedict et al. 2007; Turner 2010; Majaess et al. 2011b).

3. CONCLUSION AND FUTURE RESEARCH

The evidence presented bolsters the assertion by de Zeeuw et al. (1999) that δ Cep is a constituent of an intermediate-age cluster. The brightest cluster member is the K1.5Ib supergiant ζ Cep. δ and ζ Cep share similar HIP parallaxes (π = 3.77 ± 0.16:3.90 ± 0.10 mas), proper motions, radial velocities (RV ∼ −17 km s−1), and evolutionary histories (Figure 2). In tandem with the aforementioned evidence, the cluster's existence is supported by the absence of early-type stars from the comparison fields (Figure 2). NOMAD data were employed to identify additional potential cluster members (Tables 1 and 2).

The Cepheid exhibits parameters of: E(BV) = 0.073 ± 0.018(σ), log τ = 7.9 ± 0.1, and $d=272\pm 3(\sigma _{\bar{x}}) \pm 5 (\sigma)$ pc (Table 1, and Figures 1 and 2). The results are tied in part to spectroscopic and UBVJHKs observations, and may be adopted to refine classical Cepheid period-color, period-age, period-mass, period-luminosity, and period-Wesenheit relations (e.g., Turner 2010).

Potential future research entails establishing precise proper motions for fainter stars near Cepheids using photographic plates stored at the CfA (Grindlay 2007; DASCH),3 thereby extending the astrometric coverage provided by HIP. The plates stored at the CFA offer unmatched multi-epoch observations over a ∼100 year temporal baseline. A concurrent venture pertains to employing new VVV JHKs observations for young clusters in Galactic spiral arms to calibrate adjacent long-period classical Cepheids (Minniti et al. 2010; Moni Bidin et al. 2011; Majaess et al. 2011b). The aforementioned initiatives, in harmony with the analysis presented here, shall complement a suite of diverse efforts aimed at reducing uncertainties associated with H0 in order to constrain cosmological models (e.g., Benedict et al. 2002, 2007; Feast 2008; Gerke et al. 2011; Ngeow 2011).

The agreement between the distances inferred for δ Cep from cluster membership and the IRSB technique suggests that the systematic uncertainties have been marginalized. An analysis of n = 20 cluster Cepheids featuring revised IRSB distances (Storm et al. 2011) yields a mean fractional difference of −3% ± 3%. That result is reassuring, and subsequent research on the discrepant Cepheid calibrators (e.g., S. Vul & S. U. Cas, in preparation) may reduce the remaining offset. Further research is likewise required on the candidate cluster members associated with δ Cep (Tables 1 and 2).

D.M. is grateful to the following individuals and consortia whose efforts lie at the foundation of the research: F. van Leeuwen and M. Perryman (HIP), F. Benedict (HST), P. de Zeeuw, R. Hoogerwerf, J-C. Mermilliod, B. Skiff, NOMAD, 2MASS, AAVSO (T. Krajci, A. Henden, and M. Simonsen), and staff at the CDS, arXiv, and NASA ADS. W.G. gratefully acknowledges financial support for this work from the Chilean Center for Astrophysics FONDAP 15010003, and from the BASAL Centro de Astrofisica y Tecnologias Afines (CATA) PFB-06/2007.

Footnotes

Please wait… references are loading.
10.1088/0004-637X/747/2/145