A&A 493, 471-479 (2009)
DOI: 10.1051/0004-6361:200810550
A. Sandage1 - G. A. Tammann2 - B. Reindl2
1 - The Observatories of the Carnegie Institution of Washington,
813 Santa Barbara Street, Pasadena, CA 91101, USA
2 -
Department of Physics and Astronomy, Univ. of Basel
Klingelbergstrasse 82, 4056 Basel, Switzerland
Received 9 July 2008 / Accepted 13 October 2008
Abstract
The photometric data for 460 classical, fundamental-mode Cepheids in
the SMC with
measured by Udalski et al. have been
analyzed for their period-color (P-C) and period-luminosity (P-L)
relations, and for the variation in amplitude across the instability
strip in a similar way to what was done in Papers I and II of this
series.
The SMC Cepheids are bluer in
(B - V)0 at a given period
than for both the Galaxy and the LMC. Their P-C relation in
(B - V)0 is best fit by two lines intersecting at
days. Their break must necessarily exist also in the P-L
relations in B and/or V, but it remains hidden in the magnitude
scatter. An additional pronounced break of the P-L relations in B,
V, and I occurs at
days.
The observed slope of the lines of constant period in the HR diagram
agrees with the theoretical expectation from the pulsation equation.
The largest amplitude Cepheids for periods less than 13 days occur
near the blue edge of the instability strip. The sense is reversed in
the period interval from 13 to 20 days, as in the Galaxy and the LMC.
The SMC P-L relation is significantly flatter than that for the
Galaxy, NGC 3351, NGC 4321, M31, all of which have nearly the
same steep slope. The SMC P-L slope is intermediate between that
of these steep slope cases and the very shallow slope of Cepheids
in the lower metallicity galaxies of NGC 3109 and Sextans A/B,
consistent with the premise that the Cepheid P-L relation varies
from galaxy-to-galaxy as a function of metallicity.
Failure to take the slope differences in the P-L relation into account
as a function of metallicity using Cepheids as distance indicators
results in incorrect Cepheid distances. Part of the 15% difference
between our long distance scale - now independently supported by tip
of the red-giant branch (TRGB) distances - and that of the HST Key
Project short scale is due to the effect of using an inappropriate P-L
relation.
Key words: stars: variables: Cepheids - galaxies: Magellanic Clouds - cosmology: distance scale
The Galaxy, the LMC, and the SMC are the only three galaxies that have reddenings for their Cepheids that are determined by a method that is independent of using a fiducial Cepheid period-color relation that is incorrectly assumed to be universal. Once reliable reddenings are available by such independent methods, and by using data for galaxies where the reddening can be assumed to be nearly zero, differences in the slope and zero points of P-L relations in different galaxies as a function of metallicity follow from the observations (Tammann et al. 2008, hereafter TSR 08a, see their Fig. 4.; Sandage & Tammann 2008, hereafter ST 08, Fig. 4).
In Paper I of this series (Tammann et al. 2003, hereafter TSR 03) we analyzed the CCD photometric data by Berdnikov et al. (2000) for 321 Cepheids in the Galaxy for this purpose. A preliminary comparison of the P-L and P-C relations for the Cepheids in the LMC (593 stars) and SMC (459 stars) was also given there using the CCD photometry by Udalski et al. (1999b,a).
In Paper II (Sandage et al. 2004,
hereafter STR 04) the lines of constant period that
thread the instability strip were derived for the LMC by correlating
the residuals in magnitude at fixed period read from the P-L relation
with those from the P-C relations for individual Cepheids.
We also derived the difference in the slopes and zero points of the
ridge-lines of the instability strip (i.e. the
,
HR diagram) for the Galaxy and the LMC, showing that Cepheids in LMC
are hotter at a given luminosity by about 300 K at
(Fig. 20 of STR 04), but the size of the difference
varies with period (Fig. 3 of STR 04) because of the
slope difference in the two instability strips.
The purpose of the present paper is to continue with a
similar study for the SMC using the data of
Udalski et al. (1999b).
The plan of the paper is this.
The period-color relations in (B-V)0 and (V-I)0 for
SMC in Sect. 2 are similar to Fig. 6 of TSR 03
but with more detail, showing the difference in the P-C relations
between the Galaxy and SMC and emphasizing the break at 10 days for
SMC.
The slope of the SMC P-L relation in Sect. 3 is
significantly flatter than those for the Galaxy and LMC, similar to
Fig. 14b of TSR 03.
Lines of constant period in the HR diagram for SMC are derived in
Sect. 4. The slopes are compared there with that expected
from the pulsation equation.
Section 5 shows the correlation of amplitude with position in
the instability strip for various period and absolute magnitude
ranges.
The SMC instability strip in MV, (B-V)0 and
(V-I)0, with and without the break at 10 days, and the
resulting
instability strip is in
Sect. 6.
Comparison of the SMC P-L relation with those in the high metallicity
galaxies of NGC 3351, NGC 4321 (with slopes that are nearly
identical to the P-L slope for the Galaxy) and the low metallicity
galaxy NGC 3109 is in Sect. 7 showing the significant
difference in the P-L slope as function of metallicity and therefore
that the P-L relation is not universal from galaxy to galaxy.
In the framework of the OGLE program, Udalski et al. (1999b) have
observed UBVI magnitudes of over 2000 variables in 11 strips
covering much of the central parts of SMC (see their Fig. 1).
The magnitudes were corrected for absorption using reddening values of
adjacent red-clump stars. Their adopted reddening-to-absorption ratios
are the same for all practical purposes as those used by
STR 04.
The authors derived intensity-averaged mean magnitudes as well as periods
where possible. The light curves in I are exceptionally well defined
by several hundred epochs; the B and V light curves rest on 15-40 epochs which is sufficient in almost all cases. The authors have
derived also light curve amplitudes and Fourier coefficients R21and
.
These parameters allowed a convincing separation of
classical Cepheids from population II Cepheids and other variables
and to subdivide the classical Cepheids in 1216 fundamental and 833 overtone pulsators.
About half of the fundamental-mode Cepheids have periods shorter
than
(P < 2.5 days). This great abundance of
short-period Cepheids is specific for SMC with its very low
metallicity of [O/H]
(Sakai et al. 2004).
Some such Cepheids are known in other metal-poor galaxies, but they
are rare in young populations with higher metallicity as for instance
in LMC. The light curve parameters like magnitude, color, amplitude,
and Fourier coefficients of the short-period Cepheids in SMC are
continuous extensions of those with longer periods. Therefore they are
most likely classical populations I Cepheids
(see also Sharpee et al. 2002; Bauer et al. 1999), although
taken as a whole they lie below the extrapolated P-L relation defined
by the longer-period Cepheids and they are bluer than the extrapolated
P-C relation defined by the latter. The proliferance of short-period
Cepheids in SMC was explained by Becker et al. (1977) as the effect
of the evolutionary loops of low-mass red giants extending deeper to
the blue side of the CMD, and hence still crossing or at least
penetrating the instability strip, if the metallicity decreases. The
dependence of the loop size on metallicity was first found by
Hofmeister (1967).
In the following we are not concerned with the SMC Cepheids with
because our emphasis lies on the difference of Cepheids
from galaxy to galaxy and on the ensuing effect on distance
determinations. For this purpose the short-period Cepheids are not
helpful because statistically meaningful samples of Cepheids with
are not available in other galaxies. Therefore only SMC
Cepheids with
are considered in the following.
The sample of fundamental-mode SMC Cepheids with
was
cleaned by Udalski et al. (1999b) by removing all objects deviating
by
or more from a common P-L relation. This left 469 Cepheids with V magnitudes. We have removed some additional outliers
resulting in 460 Cepheids with V and I magnitudes, of which 439 objects have also B magnitudes. Their data are the basis of the
following analysis.
![]() |
Figure 1:
The SMC period-color relations in (B-V)0 and
(V-I)0. Small dots are for
|
| Open with DEXTER | |
Least-squares solutions of the (B-V)0 and
(V-I)0 colors of all Cepheids with
selected from
the SMC sample of Udalski et al. yield
For comparison the Galactic P-C relations for
are
repeated here from TSR 03 (the update here in the
equations for the Galaxy compared with the Galactic equations in
TSR 03 and STR 04, obtained by subtracting
Eqs. (16)-(18) there, is insignificant for this comparison):
The individual Cepheids defining the SMC P-C relations of
Eqs. (1) and (2) are plotted in
Fig. 1. Here the (B-V)0 data strongly suggest that
a single-slope P-C relation does not provide an optimum fit, but that
the slope changes around
.
Two separate fits for the
period intervals
lead for
to
For the sake of completeness we also give the separate SMC P-C
relations in (V-I)0 for the period interval
:
A comparison of the Galactic P-C relation in (B-V)0 in
Eq. (3) with the corresponding SMC P-C relations in
Eqs. (5) and (6) reveals a slope
difference, both above and below the 10 day break. It will be shown in
Sect. 5 that these slope differences translate into slope
differences in temperature in the
instability
strip diagram (Figs. 8 and 9 later).
The validity of the color differences between the Galaxy and
SMC in Fig. 1 depends, of course, on the accuracy of the
reddening determinations for both the Galaxy and the SMC. One
might suppose that the observed color difference in Fig. 1
is not real because of errors in the reddenings, although the
methods used by Fernie (1994,1990) and Fernie et al. (1995)
are very powerful. However, this possibility can be disproved by
appeal to the (B-V)0, (V-I)0 color-color diagram for
SMC compared with that for the Galaxy shown in Fig. 7a of
TSR 03. There is an offset in that diagram of
in (B-V)0 at a given (V-I)0,
with SMC Cepheids being bluer than those in the Galaxy. Because
the reddening line is nearly parallel to the observed correlation
line (Panel a of Fig. 7a of TSR 03), errors in the
reddenings will not produce the observed color offsets seen in
Panel c. Furthermore, the observed offset of SMC relative to the
Galaxy is close to what is expected from model atmosphere calculations
by Bell & Tripicco from Table 6 of Sandage et~al. (1999, hereafter
)SBT:99">#, shown in Panel d of Fig. 7a in
TSR 03. Hence, the color offset in Fig. 1 here
of SMC relative to the Galaxy is not due to reddening errors, but is
real.
The SMC P-L relations in B, V, and I from the data by Udalski et al. (1999b) are shown in Fig. 2, using a zero point for the SMC distance modulus of (m-M)0 =18.93 taken from a summary given elsewhere for modulus values determined after 2004 (TSR 08a, Table 7). The symbols are the same as in Fig. 1.
![]() |
Figure 2:
The P-L relations in B, V, and I for the SMC from
the absorption corrected data of Udalski et al. (1999b). The
different symbols for
|
| Open with DEXTER | |
Least squares solutions for the ridge-line P-L equations were
made with and without a break at 10 days using Cepheids with
.
This period restriction is to make the comparison here
with LMC more secure where the same period cut-off was used. The break
equations are shown within the borders of Fig. 2.
The non-break equations are
The P-L relations for the Galaxy from Eqs. (16)-(18) of
STR 04 are shown in Fig. 2 by the dashed
line. The deviation of the SMC P-L relation from the Galaxy is
clear. The SMC Cepheids at
are
brighter
in V than those in the Galaxy at
and are
fainter at
.
The difference is larger in B and smaller in I.
If one distrusts the Galaxy P-L slopes and zero points, either because of questions concerning main sequence fittings and problems with absorption (van Leeuwen et al. 2007) for star clusters and associations, or because of doubts about the reliability of the Baade/Becker/Wesselink moving atmosphere results due to an uncertainty on the velocity projection factor (Fouqué et al. 2007), comparison of the over-all slopes of the SMC P-L relations with those of the LMC again shows a difference, with the SMC slope being flatter in all three colors. The significance of the slope difference decreases from B to I. To avoid confusion, the LMC P-L relation is not overlayed in Fig. 2 to show this, but the LMC P-L Eqs. (7)-(9) of STR 04 differ significantly from the equations set out in Fig. 2 for SMC. The over-all SMC P-L relation compared with LMC is shown later in Fig. 10 in Sect. 7.
The comparison of the P-L relations of SMC and LMC becomes more
complicated if one allows for the break at
.
In the
interval
the SMC slopes are considerably
flatter than those of LMC. The difference is again more
pronounced in B than in I and intermediate in V. On the other
hand the SMC slopes are steeper for
than those in
LMC, the difference being about the same in all three colors.
A very tight P-L relation of SMC is obtained when the individual
Cepheids are slid along the constant-period lines onto the ridge line
of the P-L relation. The procedure is explained in TSR 08a
(Eq. (3)) and the result is illustrated in Fig. 3b there. Because of
the small dispersion of
the P-L relation is very well
defined. The break at
is less pronounced than in LMC,
but still has a significance of
and is of opposite sign!
Any claim that the P-L relations of LMC and SMC in Figs. 3a and 3b of
TSR 08a were the same would be a denial of the evidence.
Lines of constant period in the HR diagram thread from the luminous,
high temperature upper left to the low luminosity, low temperature
lower right. Therefore, Cepheids of a given period near the blue edge
of the instability strip will be brighter than Cepheids near the red
edge. Hence, there should be a correlation between luminosity
residuals in the P-L diagram and the color residuals in the P-C
diagram at fixed periods.
The slope
of the constant-period lines in the instability
strip in the
M, (B-V)0 plane is given by
.
The constant-period lines were well determined for the LMC, with
average slopes of
and
and no significant variation with
period (STR 04, Fig. 9 and Table 4).
The same determination was made there for the SMC and is repeated here
with the results shown in Fig. 3, binned in small intervals
of period. The slopes at given period intervals are similar to those
for LMC. The mean slopes over all periods are
and
for SMC.
![]() |
Figure 3: The slope of the lines of constant period using the ratios of the residuals in (B-V)0 and (V-I)0 from the P-C relations of Fig. 1, and the magnitude residuals from Fig. 2 for the Cepheids in SMC in the period intervals indicated. The sense of the correlations is that Cepheids near the blue border of the instability strip are intrinsically brighter than those near the red border. |
| Open with DEXTER | |
The equation of the family of constant period lines in the
,
HR diagram follows from the
pulsation equation when the mass is
eliminated using a mass-luminosity relation for the Cepheids. It is
shown elsewhere (ST 08) that the resulting equation for
the family is
The predicted slope here of 3.3 for the constant period
lines in V-I agrees well enough with the observed slope of
considering the approximations we have made.
A similar agreement obtains in B-V.
These agreements to within the approximations show that the observed
lines of constant period in Fig. 3 are understood as a
consequence of the Ritter
pulsation condition.
By analogy with the more easily visualized case of the RR Lyrae stars in globular clusters where the horizontal branch cuts the instability strip at nearly constant MV, we make the same correlations for the Cepheids here as is the custom for globular cluster variables. The strong correlation of amplitude with strip position in globular clusters at fixed absolute magnitude (called now the Bailey diagram) shows that the largest amplitudes for RR Lyrae stars occur near the blue edge, decreasing monotonically toward the red edge.
However, the situation is more complicated for the Cepheids because there is no restriction in absolute magnitude as with the RR Lyrae stars. The absolute magnitude in V does not change appreciably for globular cluster variables, hence there is no period-luminosity relation in V. Of course, in other bands than Vwhere the horizontal branch is not horizontal, there is a ``pseudo'' P-L relation caused however by the change of luminosity with color across the strip whereas for Cepheids the vertical change of absolute magnitude with period in the P-L strip is the principal effect.
Said differently, for RR Lyrae stars the principal parameters describing the fine structure of the instability strip are period, color, and amplitude, not luminosity, whereas for Cepheids the principal parameters are period and absolute magnitude, with color and amplitude having only a secondary role. It is the cutting of the strip by the horizontal branch in the globular clusters that restricts the way the parameters appear in the data.
The color-amplitude correlation for the RR Lyrae data is the striking feature of the data for the cluster variables. This translates into the more easily observed period-amplitude Bailey diagram which often dominates the RR Lyrae analyses. The correlation is less noticeable for the Cepheids because of the dominance of the overwhelming luminosity variation with period that is absent in the cluster variables.
Nevertheless there is a color-amplitude effect at fixed
period (or fixed absolute magnitude) for the Cepheids, with
the highest amplitude Cepheids occurring at the blue edge of the
strip, just as for the RR Lyrae stars, at least in the period
range from 3 to 10 days in the LMC (STR 04, Fig. 11) with
the trend reversing for periods from 10 to 20 days, and perhaps
returning again to high amplitudes at the blue edge for
days (STR 04, Fig. 11).
The effect was seen initially in the Galaxy where the reversal of the
sense occurs at
,
and returns to the original sense for
.
This result was originally suggested from a small sample
of Galaxy Cepheids (Sandage & Tammann 1971), and was confirmed in
STR 04 (Sect. 6.3) from the much larger sample of Galaxy
Cepheids in the Berdnikov et al. (2000) sample.
The effect was shown to be the same for the LMC Cepheids in
the period range from 3 to 10 days and for
days, and again
the sense is reversed for
days
(STR 04, Figs. 11 and 12).
We have analyzed the large data sample for the SMC in the same way, with the results shown in Figs. 4 and 5. The sense of the trends are identical with those in the Galaxy and LMC. The reversal of the sense for periods between 10 and 20 days is also evident here for the SMC although the period range is different. For SMC the reversal occurs between 13 and 20 days, while in the Galaxy it is between 7 and 20 days and for the LMC it is between 10 and 20 days. The absolute magnitude range for the change of sense is almost identical in each case at between MV of -4.25and -4.75. The reason for this behavior is not yet understood.
![]() |
Figure 4: Correlation of B band amplitude with color for SMC Cepheids in various period ranges, showing that the largest amplitude variables are near the blue edge of the instability strip except in the period range 13<P<20 days. |
| Open with DEXTER | |
![]() |
Figure 5:
Same as Fig. 4 but binned in absolute magnitude
intervals rather than period intervals. The reversal in the sense
of the correlation of amplitude with color for the absolute
magnitude interval near
|
| Open with DEXTER | |
Because the amplitudes are largest at the blue edge of the
instability strip (except for the period range of 10 to 20 days),
and because the blue edge of the strip is brighter than the red edge
at fixed period (Fig. 3 here and Fig. 9 of
STR 04 for LMC), there must be a period, magnitude,
amplitude correlation in the sense that large amplitude Cepheids are
brightest, seen in Eqs. (36)-(41) of STR 04.
The maximum effect is at the level of
,
and is a
function of the ``amplitude defect'', a parameter introduced by
Kraft (1960), calculated for each Cepheid from the deviation of
a given amplitude from the upper envelope line of the
period-amplitude relation. The bias effect on distances and the
severity of the effect for the extragalactic distance scale is small
as discussed in detail in STR 04, Sect. 6.3, and is not
repeated here. Nevertheless, the existence of the effect leads to the
notion of fine structure in the P-L relation depending on amplitude,
and is to be discussed separately elsewhere.
![]() |
Figure 6:
The Cepheid instability strip for the SMC Cepheids in
the MV-color plane allowing for a break at 10 days (see
text). Variables with periods less than
|
| Open with DEXTER | |
The instability strip in MV vs. (B-V)0 and
(V-I)0 for the SMC Cepheids is in Fig. 6.
The ridge line equations with a break at 10 days, shown at the bottom
of Fig. 6, are obtained by substituting the
term
in the P-L relations at the bottom of Fig. 2 by means of
the P-C relations given in Eqs. (5)-(8).
The ridge lines are shown as full lines in Fig. 6.
The break is suggestive in the (B-V)0 panel but is not as
pronounced in (V-I)0. The blue and red edges of the strip,
shown as dashed lines, are put parallel to the ridge lines using
widths of 0.16 and
,
respectively, for
and
0.14 and
for
.
These values are derived from
the calculated rms spreads in the P-L relations in Fig. 2
and the mean slopes of the lines of constant period from
Sect. 4. The lines of constant period are shown for
of 0.4, 0.7, 1.0, 1.3 and 1.6 using the formulation
in Sect. 4.
The insert diagrams in Fig. 6 repeat the ridge-line equations and are shown as the solid lines with the edges of the instability strip also repeated. The individual points in the inserts are the 53 Cepheids in the Galaxy sample taken from Tables 3 and 4 of TSR 03, as updated in STR 04, Sect. 4.2. The difference in the SMC ridge lines and the position of the Cepheids in the Galaxy, is evident. This difference is the same as seen in Fig. 1 but in a different representation.
Figure 7 is the same as Fig. 6 but with no allowance for the break at 10 days, i.e. the P-L relations in Eqs. (9)-(11) were combined with the P-C relations in Eqs. (1) and (2) to yield the no-break ridge line equations at the bottom of Fig. 7. Again, the differences in the slope and the zero point of the instability strip ridge lines for SMC compared with the Galaxy Cepheids (the individual points) are evident in the insert diagrams.
![]() |
Figure 7: Same as Fig. 6 but the ridge-line solutions are made with no break at 10 days. The insert diagrams again emphasize the difference in the instability strips in SMC (solid lines) compared with the Galaxy (individual points) in both slope and zero point. |
| Open with DEXTER | |
The data in Figs. 6 and 7 have been transformed
to the
,
representation of the instability
strip in the same way as we did for the Cepheids in LMC in Fig. 20 of
STR 04. The (B-V)0 and (V-I)0 colors
are changed into
by interpolations in Table 6 of
SBT 99 using appropriate
values as function of
period
(from the P-L relation and the assumption of the Cepheid
mass-luminosity relation), [A/H] metallicity values (0.0 for the
Galaxy, -0.5 for LMC, and -0.7 for SMC), and a turbulent velocity
of
.
The two values of
for each Cepheid from the
colors were averaged. We also assumed that the bolometric correction
was small enough to be neglected in converting MV to
by
(using
).
The results for the SMC break case are shown in Fig. 8. Those for the no-break solution are in Fig. 9. Also shown are the ridge-line solutions for the Galaxy (dot-dashed line) and the LMC (the heavy solid lines).
![]() |
Figure 8:
The instability strip in
|
| Open with DEXTER | |
Both the SMC and the LMC ridge lines are at higher temperatures than the Galaxy line, and each have different slopes (except for LMC with P<10 days where the ridge line is nearly parallel to that of the Galaxy). The consequences for these different slopes and zero points in the instability strips for differences in the slopes and zero points of the respective P-L relations are discussed in detail elsewhere (ST 08) and are not repeated here. However, the consequences for the P-L relation are shown graphically in the next section.
To emphasize that the data for a number of galaxies now support the proposition that the slopes of the Cepheid P-L relations differ from galaxy-to-galaxy as function of metallicity (low metallicity Cepheids have flatter P-L slopes), Fig. 10 shows a summary of the P-L relations of five galaxies (Galaxy, NGC 3351 plus NGC 4321, LMC, and NGC 3109) compared with the SMC here, illustrating the trend. More complete discussions are given elsewhere (TSR 08a, Table 5 and Fig. 4; ST 08, Fig. 4).
This diagram shows the importance of the SMC data because the difference between the Galaxy (or NGC 3351 plus NGC 4321) and the LMC is so small that a case made only using them is not convincing. However, the SMC P-L slope difference with LMC is large, and is supported more strongly when NGC 3109 is added.
The point to again be made is that the instability strip differences in slope of the ridge lines in Figs. 8 and 9 must show as a difference in the slopes of the P-L relations, as seen in Fig. 10. The case is supported because the instability strip data in Figs. 6-9 have not been used to obtain the slopes of the P-L relations in Fig. 10. These have been determined from the observational data alone using only periods and apparent magnitudes of the Cepheids in each galaxy.
![]() |
Figure 9: Same as Fig. 8 but showing the no-break solution for the SMC. |
| Open with DEXTER | |
![]() |
Figure 10:
Summary of the ridge lines of the P-L relations in 6
galaxies including the data analyzed here for the SMC. The slope
differences are evident. They become flatter as the metallicity of
the host galaxy becomes smaller. The P-L relations shown for SMC
and LMC are single-fit solutions with no allowance for a break at
|
| Open with DEXTER | |
The P-L relation of SMC predicts fainter luminosities than that
of LMC by about
.
For distance determinations the effect is
compensated or even overcompensated by the higher absorption
attributed to the blue SMC Cepheids, such that a Cepheid
sample with wide period coverage yields quite similar distances from
the SMC or LMC P-C and P-L relations. The SMC distance, however, is
larger at a period of 10 days by
.
Of course, for
still higher metallicities the difference can be larger for specific
periods (see TSR 08b, Table 3).
Cepheid distances, consistently reduced as to metallicity and period, are compiled in TSR 08b, Table 2. The resulting ``long'' distance scale is in excellent agreement with the independent tip of the red-giant branch (TRGB) distance scale (TSR 08b), demonstrating the necessity of metal-dependent P-C and P-L relations if they are used for distance determinations.
The existence of a universal P-L relation of classical Cepheids is an only historically justified illusion (see also Romaniello et al. 2008).
Acknowledgements
We thank Dres. Alfred Gautschy and Francoise Combes for helpful informations and comments. A.S. thanks the Carnegie Institution for post retirement support facilities.