abstract: We measure surface brightness fluctuation (SBF) magnitudes in the F814W filter and $(g_{475}{-}I_{814})$`<!-- -->`{=html} colors for nine bright early-type Fornax cluster galaxies imaged with the *Hubble Space Telescope* Advanced Camera for Surveys (ACS). The goal is to achieve the first systematic SBF calibration for the ACS/F814W bandpass. Because of its much higher throughput, F814W is more efficient for SBF studies of distant galaxies than the ACS/F850LP bandpass that has been used to study nearby systems. Over the color range spanned by the sample galaxies, $1.06<{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}<1.32$ (AB mag), the dependence of SBF magnitude $\overline{m}_{814}$ on $(g_{475}{-}I_{814})$`<!-- -->`{=html} is linear to a good approximation, with slope $\sim2$. When the F850LP SBF distance measurements from the ACS Fornax Cluster Survey are used to derive absolute $\overline{M}_{814}$ magnitudes, the dependence on $(g_{475}{-}I_{814})$`<!-- -->`{=html} becomes extremely tight, with a slope of $1.8\pm0.2$ and scatter of 0.03 mag. The small observed scatter indicates both that the estimated random errors are correct, and that the intrinsic deviations from the SBF–color relation are strongly correlated between the F814W and F850LP bandpasses, as expected. The agreement with predictions from stellar population models is good, both in slope and zero point, indicating that our mean Fornax distance of 20 Mpc is accurate. The models predict curvature in the relation beyond the color limits of our sample; thus, the linear calibration should not be extrapolated naively. In the Appendices, we reconsider the Tonry ground-based and Jensen NICMOS SBF distance catalogues; we provide a correction formula to ameliorate the small apparent bias in the former, and the offset needed to make the latter consistent with other SBF studies. We also tabulate two new SBF distances to galaxies observed in the ACS Virgo Cluster Survey.
author: John P. Blakeslee, Michele Cantiello, Simona Mei, Patrick Côté, Regina Barber DeGraaff, Laura Ferrarese, Andrés Jordán, Eric W. Peng, John L. Tonry, Guy Worthey
title: Surface Brightness Fluctuations in the Hubble Space Telescope ACS/WFC F814W Bandpass
  and an Update on Galaxy Distances

# Introduction

The measurement of galaxy distances has been a central problem in astronomy ever since the discovery of "spiral nebulae" by Lord Rosse in the mid-nineteenth century. Because geometric methods were of no avail, the nature of these objects remained a matter of conjecture for many decades. Only with the work of Henrietta Swan Leavitt did accurate photometric distance estimates to resolved stellar systems become a reality. Her painstaking photographic measurements of the light curves for 1800 variable stars (and hundreds of comparison stars) in the Magellanic Clouds led her to note that for variables of a certain class, the brighter ones had longer temporal periods (Leavitt 1908). The Leavitt Law, relating the luminosities of Cepheid variables to their periods, proved to be the key to the mystery of the spiral nebulae, locating them firmly beyond the spatial extent of our own Galaxy; thus it was also the key to comprehending the true scale of the Universe, eventually accomplished with the work of Hubble (1926, 1929).

Accurate galaxy distances remain of fundamental importance to most astrophysical applications. Errors in distance estimates generally translate into equal or larger fractional uncertainties in derived quantities such as masses (of everything from central black holes to dark matter haloes), linear sizes, dynamical timescales, star formation rates, and ages. Thus, accurate distance estimation is essential for inter-comparing the physical properties of galaxies in the local volume where redshift is not a reliable indicator of distance, and for comparing the properties of nearby galaxies to those at high redshift. Of course, they are also essential for mapping local three-dimensional structures and velocity fields.

There are a variety of distance indicators that can be used within and beyond the limit of classical Cepheid distances; Freedman & Madore (2010) provide a recent overview of the methods. One of the most accurate of these is surface brightness fluctuations (SBF; Tonry & Schneider 1988; for a historical review, see Blakeslee et al. 2009, hereafter ). The SBF method measures the intrinsic variance in a galaxy image resulting from the random variations in the numbers and luminosities of the stars falling within individual pixels of the image. This variance is normalized by the local galaxy surface brightness and then converted to the apparent SBF magnitude $\overline{m}$. The distance modulus, $\overline{m}$$\,{-}$$\overline{M}$ follows once the absolute magnitude $\overline{M}$ is known.

The value of $\overline{M}$ in a given bandpass depends on the stellar population properties. For early-type systems, a single broad-baseline color generally suffices for characterizing the stellar population (e.g., Tonry et al. 1997; Ajhar et al. 1997; Ferrarese et al. 2000; Blakeslee et al. 2001; Mei et al. 2005b; ). The zero-point of the $I$-band $\overline{M}$ calibration has been tied directly to the Cepheid distance scale to an accuracy of 0.08 mag, or $\sim\,$`<!-- -->`{=html}4% in distance, from ground-based SBF measurements in spiral bulges (Tonry et al. 2000). For distance estimation purposes, SBF is best measured in the reddest optical bandpasses. The fluctuations arise mainly from red giant branch stars and thus are brighter at redder wavelengths. However, the stellar population dependence becomes more complicated in the near-infrared (Jensen et al. 2003), which adds intrinsic scatter to that application of the method. Thus, the red end of the optical spectrum is something of a "sweet spot" for SBF.

The SBF method has been an important part of two recent *Hubble Space Telescope* (*HST*) surveys of early-type cluster galaxies carried out with the Advanced Camera for Surveys' Wide Field Channel (ACS/WFC). The 100-orbit ACS Virgo Cluster Survey (ACSVCS; Côté et al. 2004) and the 43-orbit ACS Fornax Cluster Survey (ACSFCS; Jordán et al. 2007) both observed one galaxy per orbit, with the time split between the F475W and F850LP filters. The SBF analysis was performed on the F850LP ($z_{850}$-band) images, while the second bandpass enabled the calibration of the SBF $z_{850}$ magnitude in terms of $(g_{475}{-}z_{850})$`<!-- -->`{=html} color (Mei et al. 2005b; ). These single-orbit observations were adequate at the distances of the Virgo and Fornax clusters, but for galaxies significantly farther away, the relatively low throughput of F850LP can make the required exposure times prohibitive. The F814W filter then becomes more attractive for SBF work simply because the higher total throughput of this passband reduces the exposure time by at least a factor of two. However, although there have been some ACS/WFC measurements of SBF in this bandpass (Cantiello et al. 2007a,b; Barber DeGraaff et al. 2007; Biscardi et al. 2008), there has been no systematic empirical calibration of the SBF method in F814W similar to that in F850LP. The main goal of the present work is to provide such a calibration that will be applicable to the bright early-type galaxies targeted at large distances.

The following section presents our observational sample and basic imaging reductions. Section  describes our SBF and color measurements in detail and compares the results to our previous F850LP measurements. Section  provides a comparison with predictions from stellar population models and discusses the uncertainty in the absolute distance to the Fornax cluster. Section  presents the conclusions from this study. In Appendix , we revisit the ground-based SBF distances from Tonry et al. (2001) and provide a simple correction formula for possible bias in those distances. Finally, Appendix  presents two additional SBF distance measurements from the ACSVCS F475W and F850LP imaging data.

# Observations and Image Reductions

The proximity ($\sim\,$`<!-- -->`{=html}20 Mpc) and compactness of the elliptical-rich Fornax cluster makes it the ideal target for calibrating the stellar population dependence of the SBF method in any bandpass (e.g., Tonry 1991; ). As we are mainly interested in the calibration for luminous early-type galaxies, we selected the eight brightest Fornax galaxies as listed by Jordán et al. (2007), all of which have $M_B<-19$. Because of the precision of ACS/WFC SBF measurements, this sample was deemed suitable for deriving a calibration accurate to about $\pm0.02$ mag over the normal color range for giant ellipticals.

As part of *HST* program GO-10911, we observed six fields with the ACS/WFC targeting large galaxies in the Fornax cluster. Each observation consisted of three dithered exposures in the F814W filter totaling 1224 s and two dithered exposures in F475W totaling 680 s. The galaxy NGC 1375 is separated by only $2\farcm4$ from its neighbor NGC 1374, and we were able to fit them both along the 4 diagonal of a single ACS/WFC field of view. NGC 1374 was placed in the upper left quadrant of the field, on the same CCD chip on which the other giant galaxies were centered, while NGC 1375 was simultaneously in the lower right quadrant. Although NGC 1375 is 1.4 mag fainter than the limit of the other eight, it provides extra leverage for constraining the calibration slope and was available gratis. In addition, we downloaded suitable F814W data from the *HST* archive on the Fornax galaxies NGC 1316 (from GO-9409, PI: Goudfrooij) and NGC 1344 (from GO-9399, PI: Carter) and suitable F475W data on the same two galaxies from GO-10217 (PI: Jordán). Table  summarizes the properties and basic observational details for our nine target galaxies, including the z\_850$z_{850}$`<!-- -->`{=html} SBF distances from the .

The exposures were bias and dark current subtracted and flatfielded by the standard STScI pipeline processing. The calibrated "flt" data were then run through the Apsis pipeline (Blakeslee et al. 2003) to produce geometrically rectified, clean, combined images. Apsis uses Drizzle (Fruchter & Hook 2002) for the geometric transformation of the images, and we use the Lanczos3 interpolation kernel in order to preserve the essential properties of the image power spectra (see Mei et al. 2005a). All photometry reported here is on the AB system, calibrated based on Sirianni et al. (2005) and Bohlin (2007) with the revised zero points from the ACS instrument website[^1] for data obtained before July 2006. In particular, we adopted: $m_1(\hbox{F475W}) = 26.0809$, $m_1(\hbox{F814W}) = 25.9593$, where $m = -2.5\,\log(f) + m_1$, and $f$ is in electrons per second. The ACS/WFC F814W and F475W bandpasses closely approximate the Kron-Cousins $I$ and Sloan Digital Sky Survey $g$ bands, so we refer to magnitudes in these bandpasses as g\_475$g_{475}$`<!-- -->`{=html} and I\_814$I_{814}$`<!-- -->`{=html}, respectively. We corrected for Galactic extinction using the dust maps from Schlegel et al. (1998) with the extinction ratios given by Sirianni et al. (2005). The extinction in this direction is small, $E(B{-}V)\lesssim 0.02$ mag.

# SBF Analysis

## Measurements within Annuli

There have been a variety of implementations of the SBF analysis within different reduction packages (e.g., Tonry et al. 1990; Lorenz et al. 1993; Pahre et al. 1999; Nielsen & Tsvetanov 2000; Liu et al. 2002; Mieske et al. 2003; Mei et al. 2005a; Cantiello et al. 2005; Dunn & Jerjen 2006), but few detailed comparisons between them. One notable exception is the comparison made by Pahre et al. (1999) of their SBF results for one galaxy observed with HST/WFPC2 with the results for the same observation using the Tonry et al. (1990, 1997) software. They found good agreement for this one case. Otherwise, the comparisons have mainly been between sets of published numbers, rather than tests of different analysis procedures applied to the same data.

In order to ensure that our calibration is as general as possible, we performed two complete sets of SBF reductions on all the galaxies. We first used the custom SBF analysis software developed by J. Tonry within the Vista environment[^2] for the ground-based SBF survey (Tonry et al. 1997) and applied most recently to *HST*/ACS data by Barber DeGraaff et al. (2007). We refer to this set of reductions as the "Tonry" SBF analysis; details can be found in the referenced works and several others that have used the same software (e.g., Blakeslee et al. 1997; Jensen et al. 1998). These reductions tend to be interactive, with the user selecting radial limits for the analysis regions within each galaxy and power spectrum fitting limits. Second, we applied the identical SBF analysis code as used in , which was only slightly modified from the analysis routines developed and used by Mei et al. (2005a,b). These routines are written in IDL and the SBF analysis is more automated, which facilitates multiple runs to test for systematics effects. As an example, Mei et al. (2005a) reran the analysis numerous times on simulated images processed with Apsis to assess the effects of different Drizzle interpolation kernels.

While there are no fundamental differences between the two analysis codes likely to cause any significant disagreements, the current data set affords a useful test of the level consistency that can be expected. For both sets of SBF measurements, we model and subtract the galaxy, mask all visible sources, and estimate the residual variance from undetected sources as detailed in our previous works (Mei et al. 2005b; Barber DeGraaff et al. 2007; Jordán et al. 2004, 2007; ). These steps are the same for nearly all implementations of the SBF method, and we used the same galaxy models and background variance corrections for both sets of reductions. For this sample of high-resolution, high-throughput F814W ACS imaging data, the estimated background corrections ranged from 0.011 to 0.019 mag. As in our previous works, we assign 25% uncertainty to these values, meaning that the contribution to the total error in the SBF magnitude is less than 0.005 mag. This contrasts strongly with ground-based analyses, for which the source detection limits tend to be much brighter and the background variance correction is often a major source of uncertainty. We masked dusty regions in the galaxy images as described in Ferrarese et al. (2006) by taking the ratio of the images in the two bands and identifying areas that significantly deviated from the smooth color profile. The dust masks were augmented by hand as needed. The two galaxies in our sample with significant dust features are NGC 1316 and NGC 1380. Thus, this comparison of SBF analyses does not address possible problems with galaxy modeling or residuals due to dust or other contaminants. However, the consistency of the results as a function of radius within a galaxy can be used to check for such problems.

Both the Tonry and IDL analysis procedures are performed in a series of annuli with increasing radii. The inner radius of the innermost annulus was 6 (128 pix) for all galaxies except the dusty post-merger galaxy NGC 1316, for which it was 12. The outer radius of each annulus was typically twice its inner radius, although for the IDL code we used a maximum annular width of 256 pix to be more similar to our past analyses, and only annuli where the galaxy I\_814$I_{814}$`<!-- -->`{=html} surface brightness is at least twice the sky background were used. For the Tonry code analysis, we set the maximum outer radius somewhat arbitrarily at 1000 pix, except in the case of the faintest galaxy NGC 1375, for which it was 512 pix. The two different analysis procedures independently calculate the power spectra and fit the amplitude of the component that has been convolved with the point spread function (PSF), which corresponds to the SBF signal after correcting for the background variance. Section  discusses systematic effects from different PSFs.

In fitting the power spectra, the IDL code omits the lowest wavenumbers, those representing spatial scales larger than 20 pix, since these are affected by the galaxy and large-scale background light subtraction, and it also omits the several highest wavenumbers because of the pixel correlation introduced by the geometric correction (see Mei et al. 2005 for details). However, with the adopted Lanczos3 kernel, the pixel correlation is very mild and the power spectrum fit is not very sensitive to the maximum wavenumber, assuming that the PSF model is accurate. The Tonry code also omits the lowest numbers, although this is done interactively by selecting the point where the fitted amplitude $P_0$ of the PSF component becomes stable. There is no option in this code to omit the highest wavenumbers, so this is another difference in the analyses. Figure  displays power spectra fits for multiple annuli within the program galaxies from the IDL analysis code. The signal is very strong and well fitted in all cases; note the logarithmic scale.

The two analysis procedures both measure the $(g_{475}{-}I_{814})$`<!-- -->`{=html} color in each annulus after applying the same mask as used for the SBF measurement. We then calculate the biweight means of the color and SBF measurements from among the annuli to get the average values for each galaxy. We note that radial color gradients are common in early-type galaxies and SBF magnitudes follow the same trends. SBF and color gradients have been studied in ACS imaging data from a stellar population perspective by Cantiello et al. (2005, 2007a). All of the galaxies in our sample have detectable color gradients except the bluest galaxy NGC 1375. However, as found by , the presence of such gradients is not a problem for our averaging procedure as long as the SBF–color relation is approximately linear over the color range within each galaxy. Figure  shows SBF and color measurements for the individual annuli within the galaxies. The top panel plots the observed values of $\overline{m}_{814}$, while the lower panel plots absolute $\overline{M}_{814}$ after subtracting the z\_850$z_{850}$`<!-- -->`{=html} SBF distance moduli from . There is a clear trend between SBF magnitude and $(g_{475}{-}I_{814})$`<!-- -->`{=html} color, well approximated by a linear relation. The scatter is about $\sim\,$`<!-- -->`{=html}30% lower for $\overline{M}_{814}$ as compared to $\overline{m}_{814}$. The improvement in scatter for a quadratic fit, as compared to linear one, is not significant.

The mean difference in galaxy SBF magnitudes from the Tonry and IDL power spectrum analyses was $0.005\pm0.005$ mag, with an rms scatter of 0.017 mag, which is within the range of the measurement errors. Only NGC 1316 showed an absolute difference as large as 0.03 mag because of the difficulty in selecting "clean" regions in this dusty galaxy; for the other eight, the difference was 0.02 mag or less. This represents the first detailed comparison between the two analysis codes that have been used most extensively for SBF measurement. The close agreement gives us confidence that our SBF calibration is applicable to either one. Here we adopt the results from the IDL analysis simply because it is the same code as used in our past ACS/WFC studies, and the steps have been described in detail by those works. Table  presents the final SBF and color measurements, as well as the distance information derived from the calibration presented below.

## Effect of the Point Spread Function

As noted in the previous section, the SBF measurement involves fitting the amplitude of the component of the power spectrum that has been convolved with the PSF of the image. Thus, the PSF template is an important consideration in the analysis. The ACSVCS and ACSFCS analyses (Mei et al. 2005b; ) both used a composite empirical PSF constructed from archival images by Blakeslee et al. (2006) for deconvolving the profiles of high-redshift cluster galaxies. The PSF variation within ACS/WFC images after distortion-corrected is small (Krist 2003) and occurs on small scales (largest wavenumbers) that do not significantly affect the power spectrum fits (Mei et al. 2005a), although a 2% error in the fitting for different annuli was assumed in  (as well as here, since the same code was used) from the spatial variation in the PSF. However, tests by Mei et al. (2005a) and Cantiello et al. (2005, 2007b) using multiple ACS/WFC PSF stars (generally taken from different images) found variations of 0.02–0.05 mag in the final $\overline{m}_{814}$ values. Similar, though slightly larger, variations of 0.04–0.06 mag were found in tests using multiple PSFs in ground-based data with a variable PSF and in *HST*/NICMOS data (Blakeslee 1999; Jensen et al. 2001). The variations in $\overline{m}$ are due to a combination of effects, including errors in determining the correct normalization of each PSF template (systematic for a given data set), and mismatch between the PSF template and galaxy image (potentially random, if the PSF varies from image to image).

For this study, we searched all of our fields for relatively isolated stars that were bright but unsaturated, and could be used to make a high signal-to-noise model of the PSF power spectrum. The density of stars in these fields, which are at Galactic latitude $b\approx-54^\circ$, is low. We found only one suitable star, which was in the IC 2006 field, and used this as our primary PSF template. We also tried PSFs from other data sets and reran the full IDL power spectrum analysis to assess the effect on $\overline{m}_{814}$. The ACS Instrument Development Team made available the composite PSFs constructed from the white dwarf spectrophotometric standard stars analyzed by Sirianni et al. (2005) and processed with the same interpolation kernel. The composite star for the F814W bandpass has been used for some previous SBF studies (Barber DeGraff et al. 2007; Cantiello et al. 2005, 2007a,b). When we used this template, our final $\overline{m}_{814}$ values became systematically fainter by 0.036 mag. We also tried individual stars from a random ACS F814W image that we had available, and the final $\overline{m}_{814}$ values became brighter by a similar amount.

Consistent with previous studies, these tests indicate variations of $\sim\,$`<!-- -->`{=html}0.04 mag from the PSF. At least for the current, fairly homogeneous, data set, any error of this size resulting from the PSF must be mainly systematic, affecting all the results in the same way, since the scatter in the calibration presented below is less than this, and is consistent with known measurement errors. The low density of Galactic foreground stars makes it difficult to characterize the PSFs in the individual images, but no significant image-to-image variations were found in the mean full widths at half maximum. The fact that we are using a high signal-to-noise template from one of our program fields should ensure that the PSF is accurate for that field, and by extension, given the small scatter that we find for the final calibration (and the quality of the power spectrum fits), the PSF template appears accurate for the full sample. However, it may be beneficial for future SBF studies that use our calibration to compare the results that they obtain with the PSF used here (available from the authors).

## SBF Calibration for F814W

Figure  plots the SBF magnitudes versus $(g_{475}{-}I_{814})$`<!-- -->`{=html} color for our sample of nine Fornax galaxies, similar to Figure , but here using the final averaged quantities for each galaxy. Fitting a linear relation for $\overline{m}_{814}$ in terms of $(g_{475}{-}I_{814})$`<!-- -->`{=html}, including errors in both axes and an additional 0.06 mag of intrinsic (or "cosmic") scatter in $\overline{m}_{814}$ at fixed $(g_{475}{-}I_{814})$`<!-- -->`{=html} due to stellar population effects (), we find $$\ensuremath{\overline{m}_{814}}\;=\; (30.384\pm0.024) \,+ \,(2.10\pm0.33)\,[{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}- 1.2] \,,
\label{eq:appcal}$$ where ${\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}=1.2$ mag is the mean color of our sample. The quoted error bars reflect the fit uncertainties. The rms scatter is 0.064 mag, consistent with measurement errors plus the expected intrinsic scatter and small depth within Fornax ().

We could shift this calibration by the mean distance modulus of $31.54\pm0.02$ mag for our nine sample galaxies, but since we have good z\_850$z_{850}$`<!-- -->`{=html} SBF measurements for each one, and want the calibrations of the two bandpasses to be as consistent as possible, we instead derive absolute $\overline{M}_{814}$ magnitudes using the individual distance measurements. The lower panel of Figure  shows the results, which give the following calibration fit: $$\begin{aligned}
\ensuremath{\overline{M}_{814}}&\;=\;& (-1.168\pm0.013\pm0.092)  \nonumber \\
  & ~ & \;+\;  \,(1.83\pm0.20)\,[{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}- 1.2] \,,
\label{eq:abscal}
\end{aligned}$$ which is valid for AB colors $1.06<{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}<1.32$ mag. The rms scatter of the data points with respect to this linear calibration is 0.029 mag, and the fit has been derived using the $(g_{475}{-}I_{814})$`<!-- -->`{=html} and $\overline{m}_{814}$ measurement errors, plus 0.015 mag of additional error added in quadrature to make $\chi^2_n=1.0$ and ensure that the fit errors are reasonable. It might seem more correct to include the errors in the individual distances used in determining the $\overline{M}_{814}$ values. However, when we do this, we find $\chi^2_n = 0.2$, which has less than a 2% probability of occurrence. The likely explanation relates to the behavior of the intrinsic scatter of 0.06 mag that was estimated by  and has been included in the distance errors. This is the largest component of the error for some galaxies, and apparently it is strongly coupled between the F814W and F850LP bandpasses, so that $\overline{M}_{814}$ and $\overline{M}_{850}$ scatter in the same way (see the following section for more discussion). We have therefore omitted this part of the error in deriving the above fit, and add only 0.015 mag to make $\chi^2_n$ unity. This 0.015 mag represents an allowance for differential intrinsic scatter between the two bands, although it may include a small contribution from some additional measurement error.

The first error bar on the constant coefficient in Equation () reflects the uncertainty from the fit, which is smaller than for Equation () because the scatter of the data points is lower when using $\overline{M}_{814}$. The second error bar gives our estimate of the systematic uncertainty including the 0.08 mag for the tie of the SBF method to the Cepheid scale (Tonry et al. 2000, 2001; Ferrarese et al. 2000), 0.02 mag for the mean distance of these galaxies from , and 0.04 mag for the possible systematic error in the PSF normalization (Sec. ). It does not include the additional uncertainty in the Cepheid distance scale, which may be of order $\sim\,0.1$ mag (Freedman & Madore 2010).

Table  lists the SBF distance moduli $\ifmmode(m{-}M)\else$(m-M)$\fi_{814}$ and linear distances $d_{814}$ derived for the current sample using the calibration given by Equation (). The tabulated distance uncertainties reflect the total random errors, including measurement errors and the estimated 0.06 mag of intrinsic scatter. They do not include the systematic error in the calibration zero-point, which is common to all.

## Comparison with z\_850$z_{850}$`<!-- -->`{=html}-band SBF

It is useful to compare directly our I\_814$I_{814}$`<!-- -->`{=html} SBF magnitudes and distance estimates with the z\_850$z_{850}$`<!-- -->`{=html} results from . We find a mean "SBF color" for these galaxies of $\ensuremath{\overline{m}_{814}}{-}\ensuremath{\overline{m}_{850}}= 0.72$, with an rms of 0.027 mag. However, the SBF color decreases slightly, or gets "bluer," as the SBF magnitudes get fainter in redder galaxies, as shown in Figure , where the fitted relation including errors in both coordinates is given by $$\overline{m}_{814} - \overline{m}_{850} \;=\;
 (0.72\pm0.01) \,- \,(0.33\pm0.18)\,[{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}- 1.2] \,,
%
\label{eq:sbfcolor}$$ with an rms scatter of 0.021 mag. The significance of the slope is slightly less than $2\sigma$, but the same trend occurs as a function of $(g_{475}{-}z_{850})$`<!-- -->`{=html}, with the same scatter. It is not too unusual for the SBF color to become bluer while photometric color becomes redder over some range in population parameters, especially when age is increased at fixed metallicity (Blakeslee et al. 2001). However, SBF and color measurements would be needed for more galaxies in more bandpasses to place strong constraints on the underlying populations. Put another way, the slope of $\overline{m}_{814}$ as a function of $\overline{m}_{850}$ is slightly less than one, as shown in the left panel of Figure , where the plotted linear fit has slope $0.91\pm0.06$. The rms scatter about this fit is 0.022 mag, which is consistent with random measurement errors, as the reduced $\chi^2_n = 0.8$ is reasonable for a linear fit to 9 data points. This supports our conclusion above that any intrinsic deviations in $\overline{m}_{814}$ and $\overline{m}_{850}$ due to stellar population effects must track closely between the two bands.

The right panel of Figure  provides a comparison of the distance moduli from the present work using Equation () with those from . The plotted line indicates unity, and the mean difference between the two sets of measurements is zero by design, since the distance moduli from  went into deriving the $\overline{M}_{814}$ magnitudes used for the I\_814$I_{814}$`<!-- -->`{=html} calibration. The plotted error bars include the 0.06 mag of intrinsic scatter for each band, and naively using these errors, $\chi_n \approx0.1$, which has only a $\sim\,$`<!-- -->`{=html}0.2% probability of occurring. Clearly, the intrinsic scatters in the SBF magnitudes in these two bandpasses are strongly correlated, and the two sets of measurements cannot be considered completely independent, even though they come from different data sets in different filters. However, this comparison is useful in that it shows the random measurement errors are small as expected, $\lesssim 0.03$ mag. Further, the tight correlation ensures consistency of the  distances with those measured here, and thus with future F814W SBF studies that use our calibration.

# Model Comparison and the Distance to Fornax

The SBF distance zero point used in this study is an empirical one, ultimately tied to the metallicity-corrected Cepheid distance scale. More specifically, it assumes a mean distance of 16.5 Mpc to the Virgo cluster based on the Tonry et al. (2001) ground-based SBF distances recalibrated with the final set of *HST* Key Project Cepheid distances of Freedman et al. (2001); see Mei et al. (2005b) and Appendix  of the present work for further details. The mean distance of $20$ Mpc for the Fornax cluster then comes from the measured $0.42\pm0.03$ magnitude offset between the Virgo and Fornax ().

There is also a long history of attempts to calibrate the SBF method from stellar population modeling, thereby making it another "primary" distance indicator not tied to the Cepheids (e.g., Worthey 1993; Buzzoni 1993; Blakeslee et al. 2001; Mei et al. 2001; Cantiello et al. 2003; Raimondo et al. 2005; Marin-Franch & Aparicio 2006). Biscardi et al. (2008) used the Teramo "SPoT" models[^3] (Raimondo et al. 2005) in order to derive a theoretical linear calibration of the absolute SBF magnitude in the ACS/WFC F814W bandpass as a function of $(g_{475}{-}I_{814})$`<!-- -->`{=html} color. This is the same combination of bandpasses as we have calibrated here from an empirical standpoint, although we work in the AB system while Biscardi used the VEGAMAG system. The slope is of course unaffected by this difference, and their quoted value of $2.2\pm0.2$ is close to our empirical results in Equations () and ().

Comparisons between empirical calibrations and those derived from simple stellar population models are complicated because one is never sure which models best represents actual galaxies. But it is a useful exercise to check for general consistency. Figure  provides a direct comparison of our SBF and color observations with the predictions from the SPoT models, as in Biscardi et al. (2008) but transformed to the AB system. The top panel shows the loci for models with metallicities $[\hbox{Fe/H}] = -0.7, -0.3, 0.0, +0.3$, and ages ranging from 3 to 14 Gyr for each. We also show our measurements from Figure , along with the fit from Equation (). The lower panel is similar, but symbols are used to show the model values for 6 different ages at each metallicity. The linear empirical calibration overlaps remarkably well with the locus of the solar metallicity models, which must be at least partly fortuitous, but consistent with the mean metallicities in these galaxies being not too far from solar. The curve in the lower panel is a cubic polynomial fit to the plotted models.

Two things are particularly worthy of note in Figure . First, there is remarkable consistency between the zero points of the observations and models. In fact, the fits to data and models cross very near the mean color of ${\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}= 1.2$. The excellent agreement between two very different avenues for calibrating the F814W SBF magnitudes suggests that neither is likely to be far wrong, and thus the mean Fornax distance of $\sim20$ Mpc appears correct. Second, the fit to the model predictions indicates some curvature in the relation, similar to what we found in  for the $\overline{m}_{850}$ magnitudes of a larger sample of galaxies including dwarfs that extended to much bluer colors. Over the color range of the galaxies in our sample, the deviation between the linear and cubic fits is small, so the linear calibration suffices. However, one should not blindly apply this calibration beyond the fitted range, as it is reasonable to expect that the curvature may be significant.

# Summary

We have measured SBF $\overline{m}_{814}$ magnitudes and $(g_{475}{-}I_{814})$`<!-- -->`{=html} colors for nine bright early-type galaxies in the Fornax cluster from new and archival *HST* ACS/WFC imaging data. We employed two different analysis procedures for measuring the SBF and found excellent agreement between them. For this sample of galaxies, which spans the color range $1.06<{\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}<1.32$ (AB mag), there is a clear trend of $\overline{m}_{814}$ with $(g_{475}{-}I_{814})$`<!-- -->`{=html} that is well approximated by a linear relation of slope $\sim2$. We have used the $\ifmmode z_{850}\else$z\_850$\fi$ SBF distances from  to obtain absolute $\overline{M}_{814}$ for all the galaxies and derive an empirical calibration of the I\_814$I_{814}$`<!-- -->`{=html} SBF method that is consistent with the large sample of existing z\_850$z_{850}$`<!-- -->`{=html} SBF data. Because of the greater efficiency of the F814W bandpass, the linear calibration defined by Equation () will allow accurate SBF distances to be measured with the ACS/WFC at greater distances than can be reached with the z\_850$z_{850}$`<!-- -->`{=html} band in the same amount of time. For instance, 4 orbits in I\_814$I_{814}$`<!-- -->`{=html} suffices for an SBF distance to the Coma cluster at 100 Mpc, whereas $\sim\,$`<!-- -->`{=html}10 orbits would be required in z\_850$z_{850}$`<!-- -->`{=html}.

The observed scatter in the calibration of $\overline{M}_{814}$ as a function of $(g_{475}{-}I_{814})$`<!-- -->`{=html} is only 0.03 mag, significantly smaller than the expected $\sim\,$`<!-- -->`{=html}0.06 mag from intrinsic scatter in the method due to stellar population variations at a given color (Tonry et al. 1997; ). This indicates that there is coherence between the intrinsic scatter in the z\_850$z_{850}$`<!-- -->`{=html} and I\_814$I_{814}$`<!-- -->`{=html} SBF relations, as expected from the similarity, and even partial overlap, of these passbands. Thus, stellar population effects cause the $(\ensuremath{\overline{m}_{814}}{-}\ensuremath{\overline{M}_{814}})$ and $(\ensuremath{\overline{m}_{850}}{-}\ensuremath{\overline{M}_{850}})$ distance moduli to scatter in the same direction, so that they are not fully independent. Direct comparison of the two sets of distances supports this conclusion, which helps to ensure a high degree of consistency between SBF distances measured in these two bands.

Comparison of our measurements with stellar population model predictions provides support for the empirically-based mean distance of 20 Mpc to the Fornax cluster. According to these models, and by analogy to the broader-baseline z\_850$z_{850}$`<!-- -->`{=html} SBF calibration, it is likely that the I\_814$I_{814}$`<!-- -->`{=html} SBF calibration will show curvature beyond the color limits of our sample. Therefore, the calibration should not be linearly extrapolated far beyond the range explored here. However, our sample was designed to have colors similar to those of the giant ellipticals targeted in more distant SBF studies; consequently, there should be little need for extrapolation. At least two different SBF studies in the F814W bandpass are ongoing and can benefit from our calibration. We have emphasized the need for good PSF templates in SBF measurements and the importance of investigating systematic variations due to the PSF. Our results indicate that when properly calibrated and controlled for systematics, the SBF method remains one of the most accurate ways of measuring extragalactic distances.

Support for this work was provided by NASA through grant number GO-10911 from the Space Telescope Science Institute, which is operated by AURA, Inc., under NASA contract NAS 5-26555. M.C. acknowledges support from COFIS ASI-INAF I/016/07/0 and PRIN-INAF 2008 (PI. M. Marconi). A.J. acknowledges support from BASAL CATA PFB-06, FONDAP CFA 15010003 and MIDEPLAN ICM Nucleus P07-021-F. This research made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with NASA.

*Facility:*

# A. SBF Survey Distance Corrections

Tonry et al. (2001; hereafter in this appendix, T01) published distances for 300 galaxies from the ground-based $I$-band SBF Survey of Galaxy Distances (Tonry et al. 1997). The median quoted uncertainty of 0.22 mag, just over 10% in distance, makes this the largest catalogue of nearby galaxy distances with this level of precision. It has proven a useful resource for numerous extragalactic studies. However, following the publication of our combined sample of 134 galaxies with high-quality ACS z\_850$z_{850}$`<!-- -->`{=html}-band SBF distances (median error 0.075 mag) and the direct comparisons for 50 galaxies in common with T01, there has been some confusion over the possible need for correcting the T01 distances. As we have received more than a few queries on this issue, we address here in detail the question of the T01 distances, as well as those of Jensen et al. (2003; hereafter J03).

The ground-based SBF survey was conducted over numerous observing runs at multiple observatories with different types of detectors under variable seeing and photometric conditions during the course of a decade which saw great evolution in astronomical instrumentation and our own reduction software. It targeted galaxies ranging from the Local Group out to the fuzzy, condition-dependent limit of the ground-based SBF method. Although great effort was taken to homogenize the disparate data sets and produce the best possible set of final distances, the limitations of the data were already clear at the time; the authors recommended "further study of the degree of bias present in $[$the T01$]$ data set" (see Sections 4 and 7 of T01).

A separate issue was the zero point. The calibration of the T01 distance catalogue was based on comparisons to six galaxies with *HST* Key Project (KP) Cepheid distances as tabulated by Ferrarese et al. (2000), which yielded an $\overline{M}_I$ zero point of $-1.74$ mag (Tonry et al. 2000). Blakeslee et al. (2002) found that when the calibration was rederived using the revised KP Cepheid distances from Freedman et al. (2001), the resulting zero-point shifted fainter to $-1.68$ mag. If there were no distance-dependent biases, then this simple shift in zero point would mean that all the T01 distance moduli just needed to be revised downward by 0.06 mag. For instance, our ACSVCS papers applied this shift to the mean Virgo distance presented by T01, thus adopting 31.09 mag, or 16.5 Mpc, for Virgo. This was a reasonably "safe" assumption, given that the six galaxies used to tie together the SBF and Cepheid distances extended out to the Virgo cluster, and were generally targeted in above-average observing conditions. Note that some researchers have misinterpreted this 0.06 mag shift as resulting from a change in the adopted Hubble constant, but there is no explicit dependence on the Hubble constant, only a dependence on the zero point of the Cepheid distances.

In addition to rederiving the distance zero point based on the revised KP Cepheid distances, Blakeslee et al. (2002) found some evidence for distance bias in the lowest quality T01 data, in the sense that the relative distances tended to be underestimated for the lower quality data, as suspected by T01. This conclusion was based on comparisons with distances estimated from the fundamental plane and density field-corrected velocity measurements. Note that the absolute data quality generally decreases with distance, e.g., the errors increase with distance in T01, so this could potentially translate into a distance-dependent bias. Because of this combination of zero point and bias issues, Blakeslee et al. (2009) chose to compare the distances for the 50 galaxies in common between the T01 and ACS data sets without first applying any shifts to T01. This contrasted with our approach in Mei et al. (2007), where the comparison was made only for Virgo galaxies and the zero-point shift from Blakeslee et al. (2002) was applied; this difference may have been the cause of some subsequent confusion. Blakeslee et al. (2009) found the same mean Fornax distance modulus to within $\pm0.01$ mag as T01, without any shifting in zero point. However, the T01 Virgo distance modulus is then necessarily too high by $0.06$ mag. As compared to the ACS studies, the data tabulated by T01 therefore give a smaller relative distance modulus of Fornax with respect to Virgo, although the uncertainties overlap. In light of these results, it becomes unclear whether or not a simple shift in the zero point of the T01 distance moduli is useful without a second-order correction that may depend on the distances of the galaxies under consideration. Given the 0.22 mag median error of the T01 distances, it could be argued that these effects are in the noise for any given galaxy. For galaxies within Virgo and Fornax, the tabulated ACS distances from Blakeslee et al. (2009) are much more accurate, and negligibly affected by the background variance corrections that were suggested to bias the ground-based data.

However, for some purposes, it may be desirable to have a general correction formula for the published T01 distances, including both the zero point and second-order bias corrections. Given our high-quality measurements for a sizable fraction of the T01 sample, we have attempted to derive one. We investigated possible correlations of the distance modulus difference $\Delta_{\rm T01} \equiv \ifmmode(m{-}M)\else$(m-M)$\fi_{\rm ACS} - \ifmmode(m{-}M)\else$(m-M)$\fi_{\rm T01}$, where $\ifmmode(m{-}M)\else$(m-M)$\fi_{\rm ACS}$ comes from Blakeslee et al. (2009), with various quantities tabulated by T01 so that the correction is available for that full sample. In general, the correlations are weak, except between $\Delta_{\rm T01}$ and $\ifmmode(m{-}M)\else$(m-M)$\fi_{\rm T01}$; there is no significant correlation with $\ifmmode(m{-}M)\else$(m-M)$\fi_{\rm ACS}$. This is easily understood if the ACSFCS-V distances are taken as the "true" values, or at least unbiased and homogeneous in quality; then the distance measurement error $\Delta_{\rm T01}$ will correlate strongly with the measured distance over any limited range in true distance. But this correlation, caused by measurement error, cannot be used to bias-correct the measured distances for the full sample. There is no correlation of $\Delta_{\rm T01}$ with the quoted measurement error, perhaps because of the relatively small range in the latter quantity and the multitude of effects that contribute to it.

It seemed more promising to investigate correlations with parameters that directly reflect the data quality. It is also necessary for the correction parameter to span reasonably well the same range in the current subsample as in the complete T01 catalogue. T01 supplied two measures of the quality of the data for each galaxy. One such parameter, called $PD$, was the normalized product of the width of the seeing disk in arcseconds and the expected distance based on the galaxy recessional velocity in units of 1000 km s$^{-1}$, thus making it independent of any potential distance biases in the SBF distance measurement itself. $PD$ is proportional to the linear size of the resolution element at the distance of the galaxy, and thus is smaller for better quality data. Another measure of data quality, called $Q$ by T01 and here referred to as $Q_{\rm T01}$, scaled logarithmically with the ratio of the signal-to-noise and $PD$; thus it is larger for better quality data. Blakeslee et al. (2002) investigated bias as a function of $PD$; we have found a slightly more significant trend as a function of $Q_{\rm T01}$. Figure  shows $\Delta_{\rm T01}$ plotted against $Q_{\rm T01}$. The fitted relation suggests that the T01 distance moduli can be bias-corrected according to $$\begin{aligned}
\label{eq:ton01cor}
(m{-}M)_{\rm T01,\,cor} &=& (m{-}M)_{\rm T01,\,raw} \,+\, 0.1 \,-\, 0.03\,\ensuremath{Q_{\rm T01}}\\
                        &=& (m{-}M)_{\rm T01,\,raw} \,-\, 0.03(\ensuremath{Q_{\rm T01}}- 3.3)
                        \,. \nonumber
\end{aligned}$$ Thus, as published without any zero-point shifting, the poorest quality T01 distances with $Q_{\rm T01}\approx1$–2 would be corrected by $\sim\,+$`<!-- -->`{=html}0.05 mag, and the best quality ones with $Q_{\rm T01}\approx7$–8 would be corrected by $\sim\,-$`<!-- -->`{=html}0.12 mag. Given the range of $Q_{\rm T01}$ in Figure , we do not recommend applying corrections any larger than this. The six Cepheid calibrator galaxies have an average $\langle\ensuremath{Q_{\rm T01}}\rangle = 6.0$, giving $\langle\Delta_{\rm T01}\rangle = -0.08$ mag, (or, omitting the Local Group galaxy M31: $\langle\ensuremath{Q_{\rm T01}}\rangle = 5.2$, $\langle\Delta_{\rm T01}\rangle = -0.06$ mag), very close to the expected shift from the Cepheid calibration, which is included within the correction formula. The significance of the fitted slope is only $1.75\sigma$ based on these 50 galaxies, and the relative Fornax–Virgo distance modulus from the sample is unchanged, but if the trend persisted for the full T01 sample of 300 galaxies, the significance would reach $4\sigma$; therefore studies relying on many T01 galaxy distances may benefit from applying a correction such as Equation ().

Finally, we address the J03 *HST*/NICMOS F160W SBF distances for 65 galaxies. These were tied to the KP Cepheid distances from Freedman et al. (2001), but without the espoused $-$`<!-- -->`{=html}0.2 mag dex$^{-1}$ metallicity correction. The evidence for this correction was ambiguous at the time, and J03 found that the agreement with SBF predictions from stellar population models was better when the distance zero point was based on the Cepheid results without metallicity correction. Since then, there has been additional evidence to support a Cepheid metallicity dependence (Sakai et al. 2004; Macri et al. 2006; Scowcroft et al. 2009), and better understanding of the uncertainty in model predictions for near-IR SBF (Raimondo 2009; González-Lópezlira 2010; Lee et al. 2010). Given the high signal-to-noise and good resolution of the J03 NICMOS data, there is no reason to suspect a bias similar to that of T01. However, to bring them into agreement with other work requires putting them on the Cepheid scale that includes the metallicity correction. J03 derived their distance zero point by comparison to 47 T01 galaxies after shifting the distance moduli by $-0.16$ mag, which comes from repeating the Tonry et al. (2000) calibration using the revised Cepheid distances without metallicity correction. As noted above, the change in the T01 moduli when calibrated on the revised Cepheid distances with metallicity correction is $-0.06$ mag. This suggests a shift of $+0.10$ mag to the Jensen et al. distance moduli, although there is some question about possible bias in the particular subsample of T01 that was used. J03 also measured SBF in nine KP Cepheid galaxies and tabulate the resulting $\ensuremath{\overline{M}}_{160}$ values for the distances with and without metallicity correction; the mean difference is 0.08 mag, which is consistent with the shift in the T01 calibration caused by the metallicity correction. We therefore conclude that the J03 distance moduli should be increased uniformly by $+$`<!-- -->`{=html}0.1 mag to bring them into consistency with the corrected T01 distances, our recent ACS $\ifmmode z_{850}\else$z\_850$\fi$ SBF measurements, and the present work.

# B. Additional ACSVCS Dwarf Galaxy Distances

Mei et al. (2007) published SBF measurements for 90 galaxies from the ACSVCS g\_475$g_{475}$`<!-- -->`{=html} and z\_850$z_{850}$`<!-- -->`{=html} imaging data, and Blakeslee et al. (2009) tabulated recalibrated distances for these along with the 43 ACSFCS galaxies and NGC 4697 in the Virgo Southern Extension. No SBF measurements were published for ten of the ACVCS galaxies, mainly owing to irregular morphologies. However, the problem for the galaxies VCC 1192 and VCC 1199 was different: they are both compact elliptical companions of the Virgo brightest cluster galaxy M49 (NGC 4472 or VCC 1226) and the strong background gradient from the giant galaxy's halo caused the photometric measurements for these galaxies to be unreliable. Since the physical distances from M49 of these likely stripped (Ferrarese et al. 2006; Côté et al. 2010) companions are of interest, we have undertaken new measurements of these two galaxies. Although not directly related to the main content of this paper, these data are part of our larger *HST* SBF effort, and we present them here.

We constructed initial isophotal models of the light distributions of VCC1192 and VCC1199 in each filter. We then subtracted these models, masked the background sources, and interpolated two-dimensional cubic spline surfaces on a $16{\times}16$ grid (each grid cell being $\sim270$ pix on a side) to make smooth fits of the background light distribution. Since these galaxies are so compact, and they had been removed by the initial models, they did not perceptibly affect the fits. After subtraction of the fitted surfaces, the backgrounds were extremely flat. We then constructed a new galaxy model for each and proceeded with the SBF and color measurements as usual, using the same code as in the present work. Table  presents our measurements for these galaxies in the same format as for the other ACSVCS galaxies in Table 2 of Blakeslee et al. (2009). This brings the total number of galaxies with measured z\_850$z_{850}$`<!-- -->`{=html} SBF distances to 136.

```latex
\begin{deluxetable}{lrcrccccc}
\tabletypesize{\small}
\tablecaption{Galaxy Data}
 \tablewidth{0pt}
\tablehead{
\colhead{Galaxy} &
\colhead{$B_T$} &
\colhead{Morph} &
\colhead{$v_h$} &
\colhead{Other~} &
\colhead{$T_{\rm F475}$} &
\colhead{$T_{\rm F814}$} &
\colhead{$\ifmmode(m{-}M)\else$(m{-}M)$\fi_{850}$} &
\colhead{Program}
 \\
\colhead{(1)} &
\colhead{(2)} &
\colhead{(3)} &
\colhead{(4)} &
\colhead{(5)} &
\colhead{(6)} &
\colhead{(7)} &
\colhead{(8)} &
\colhead{(9)}
}
\startdata
NGC\,1316 & 9.4 & S0\,pec& 1760 & FCC\,21   & 760 & 4680  & $31.607\pm0.065$ & 10217,9409 \\
NGC\,1344 & 11.3 & E5    & 1169 & NGC\,1340 & 760 &  960  & $31.603\pm0.068$ & 10217,9399 \\
NGC\,1374 & 11.9 & E0    & 1294 & FCC\,147  & 680 & 1224  & $31.459\pm0.070$ & 10911 \\
NGC\,1375 & 13.6 & S0    &  740 & FCC\,148  & 680 & 1224  & $31.499\pm0.072$ & 10911 \\
NGC\,1380 & 11.3 & S0/a  & 1877 & FCC\,167  & 680 & 1224  & $31.632\pm0.075$ & 10911 \\
NGC\,1399 & 10.6 & E1    & 1425 & FCC\,213  & 680 & 1224  & $31.596\pm0.091$ & 10911 \\
NGC\,1404 & 10.9 & E2    & 1947 & FCC\,219  & 680 & 1224  & $31.526\pm0.072$ & 10911 \\
NGC\,1427 & 11.8 & E4    & 1388 & FCC\,276  & 680 & 1224  & $31.459\pm0.068$ & 10911 \\
IC\,2006 & 12.2  & E2    & 1382 & ESO\,359-07& 680 & 1224 & $31.525\pm0.086$ & 10911 \\
\enddata
\label{tab:tab_obs}
\tablecomments{Columns list:  (1)~galaxy name;
  (2)~total $B$~magnitude from the FCC (Ferguson 1989) or NED;
  (3)~morphological type from the FCC or NED;
  (4)~heliocentric velocity from NED (km~s$^{-1}$)
  (5)~alternative galaxy designation;
  (6)~{\it HST}/ACS exposure time for the F475W bandpass (sec);
  (7)~{\it HST}/ACS exposure time for the F814W bandpass (sec);
  (8)~distance modulus and uncertainty from \ifmmode z_{850}\else$z_{850}$\fi-band SBF (ACSFCS-V);
  (9)~{\it HST}\ program IDs for the observations used in this study.
}
\end{deluxetable}
```
```latex
\begin{deluxetable}{lcccc}
\tabletypesize{\small}
\tablecaption{Galaxy Color and SBF Measurements}
\tablewidth{0pt}
\tablehead{
\colhead{Galaxy} &
\colhead{{\ifmmode{(g{-}I)}\else$(g{-}I)$\fi}} &
\colhead{$\overline{m}_{814}$} &
\colhead{$\ifmmode(m{-}M)\else$(m{-}M)$\fi_{814}$} &
\colhead{$d_{814}$~(Mpc)}
 \\
\colhead{(1)} &
\colhead{(2)} &
\colhead{(3)} &
\colhead{(4)} &
\colhead{(5)}
}
\startdata
NGC\,1316 & $1.177\pm0.011$ & $30.428\pm0.020$ & $31.638\pm0.066$ & $21.3\pm0.7$ \\
NGC\,1344 & $1.131\pm0.015$ & $30.315\pm0.039$ & $31.609\pm0.076$ & $21.0\pm0.7$ \\
NGC\,1374 & $1.193\pm0.022$ & $30.286\pm0.015$ & $31.467\pm0.074$ & $19.7\pm0.7$ \\
NGC\,1375 & $1.065\pm0.009$ & $30.116\pm0.037$ & $31.530\pm0.072$ & $20.2\pm0.7$ \\
NGC\,1380 & $1.232\pm0.011$ & $30.460\pm0.016$ & $31.570\pm0.065$ & $20.6\pm0.6$ \\
NGC\,1399 & $1.315\pm0.017$ & $30.662\pm0.022$ & $31.620\pm0.071$ & $21.1\pm0.7$ \\
NGC\,1404 & $1.277\pm0.014$ & $30.516\pm0.018$ & $31.544\pm0.068$ & $20.4\pm0.6$ \\
NGC\,1427 & $1.171\pm0.019$ & $30.213\pm0.017$ & $31.434\pm0.072$ & $19.4\pm0.6$ \\
IC\,2006 & $1.234\pm0.016$ & $30.430\pm0.035$ & $31.536\pm0.075$ & $20.3\pm0.7$ \\
\enddata
\label{tab:tab_result}
\tablecomments{Columns list:  (1)~galaxy name;
  (2)~mean galaxy {\ifmmode{(g_{475}{-}I_{814})}\else$(g_{475}{-}I_{814})$\fi}\ color;
  (3)~mean SBF magnitude \ensuremath{\overline{m}_{814}};
  (4)~distance modulus derived from the linear \ifmmode I_{814}\else$I_{814}$\fi\ SBF calibration presented here;
  (5)~the distance in Mpc.
The distance errors include contributions from SBF and color measurement
errors, as well as the estimated intrinsic (``cosmic'') scatter in the method.}
\end{deluxetable}
```
```latex
\begin{deluxetable}{lrrrrrc}
\tabletypesize{\small}
\tablecaption{Distances for Additional ACS Virgo Cluster Survey Galaxies}
\tablewidth{0pt}
\tablehead{
\colhead{Galaxy} &
\colhead{$(g{-}z)$} &
\colhead{$\ensuremath{\overline{m}_z}$} &
\colhead{\ifmmode(m{-}M)\else$(m{-}M)$\fi} &
\colhead{$d$} &
\colhead{$B_T$} &
\colhead{Name} \\
%
\colhead{(1)} &
\colhead{(2)} &
\colhead{(3)} &
\colhead{(4)} &
\colhead{(5)} &
\colhead{(6)} &
\colhead{(7)}
}
\startdata
VCC1192 & $1.443\pm0.009$ & $29.305\pm0.065$ & $31.078\pm0.091$ & $16.4\pm0.7$ & 15.0 & NGC4467 \\
VCC1199 & $1.469\pm0.030$ & $29.345\pm0.169$ & $31.055\pm0.196$ & $16.3\pm1.5$ & 15.5 & \dots \\
\enddata
\tablecomments{Columns list:  (1) VCC designation (Binggeli {et~al.}\ 1985);
  (2)~mean {\ifmmode{(g_{475}{-}z_{850})}\else$(g_{475}{-}z_{850})$\fi}\ color of the analyzed region;
  (3)~mean F850LP SBF magnitude \ensuremath{\overline{m}_z};
  (4)~distance modulus derived from Equation\,$[6]$ of \hbox{ACSFCS-V};
  (5)~distance in Mpc;
  (6)~total $B$~magnitude from the VCC or NED;
  (7)~common galaxy name.
}
\label{tab:vcc119x}
\end{deluxetable}
```

[^1]: http://www.stsci.edu/hst/acs/analysis/zeropoints

[^2]: http://www.ifa.hawaii.edu/$\sim$jt/soft.html

[^3]: http://www.oao-teramo.inaf.it/SPoT


# INTRODUCTION

The magnetic activity typical for the solar-type stars (spectral classes F8-K2, $T_\textsl{eff} = 6200–4500\;$K) and cooler K–M dwarfs with effective temperatures of up to $2200\;$K determines the diversity of variability observed in different wavelength ranges. The time scale of variability is in the range from seconds (fast flares) to tens of years (activity cycles similar to the solar magnetic cycle); there are also secular cycles of hundreds and thousands of years. On the Sun, manifestations of activity in different layers of the atmosphere are regulated by an 11-year cycle. The development of magnetic activity cycles is determined by the rotation of a star and its internal structure.

According to Nizamov et al. (2017), cycles in G-type dwarfs are established from the epoch when their rotation period reaches 1.1 days, for K and M dwarfs these periods are longer – 3.3 and 7.2 days, respectively. The stars of the corresponding spectral classes, rotating faster, demonstrate the saturated activity regime.

A certain number of stars with cycles longer than the solar one are known, but most of the studied objects were investigated from the modern databases on a limited time interval; therefore, cycles of no longer than 10–15 years were found for these stars (Savanov, 2012; Suárez Mascareño et al., 2016; Boro Saikia et al., 2018). The results concerning stellar cycles require further specification and confirmation from the long-term systematic observations.

Oláh et al. (2016) derived the series of data on time scales up to 36 years for 29 G-K stars from the Mount Wilson HK-Project and determined the duration of their dominant cycles (with a maximum amplitude) and short cycles.

The long-term variations in brightness of several dozens of cool K–M dwarfs were studied from photographic archives and broadband photometry (Bondar', 1995, 2013; Alekseev and Bondar', 1998; Alekseev and Kozhevnikova, 2017); for some of them the activity of photospheric starspots similar to the observed on the Sun has been found or suspected.

The task of our work is to form a sample of stars, whose cycles have been found from the data over several decades; to substantially extent the available photometric series for some G–M dwarfs by the results of recent observations; to refine or determine cycles of their photospheric activity.

Besides, we will compare characteristics of dominant cycles of solar-type stars and cooler dwarfs obtained over long-term series, consider how they are related to rotation periods, reveal features of this connection for different type stars, that is important for the dynamo theory.

# G–M DWARFS WITH CYCLIC ACTIVITY ON LONG-TERM INTERVALS

The first results on activity cycles of F–M stars on or near the main sequence were obtained by Wilson (1978) within his monitoring program of 91 stars in the Ca II H and K emission lines. Based on observations in 1966–1977, for some stars there have been suspected cycles (trends) which require certain evidence from the longer data series. The program was continued in the framework of the HK-Project for 111 stars. In the paper by Baliunas et al. (1995), which became the basis for studying magnetic activity of the solar-type stars, the cycles were searched for and classified from the data series which include the 12-year observations by O. Wilson and new data, and in whole cover the 25-year interval.

Oláh et al. (2016) presented new data about chromospherical activity cycles for 29 HK-Project stars from the 36-year observation series. Most of these stars exhibit cycles with a duration of 7–18 years, i.e. which correspond to the sunspot cycle. Some stars, besides the main cycle, demonstrate one or several shorter cycles (multiple cycles), as well as long-term trends of 25 years or longer.

For our purpose, we take into account the results of a search for cycles on cool dwarfs from the series longer than 20 years, i.e. an investigating time spans must be comparable with those for the solar-type stars in Oláh et al. (2016).

We have used the results derived for 40 dK–dM stars from photographic data series (Bondar', 1995, 2013), from photometric researches by Alekseev and Kozhevnikova (2017) and Messina et al. (2002). For 9 stars (Bondar' 2019) the available data series were added by data from photometric catalogs. Data from the Hipparcos catalog (<http://cds.u-strasbg.fr/cgi-bin/Dic-Simbad?HIP>) cover
the span of 1990–1993, data from the ASAS catalog (Pojmanski, 1997) were obtained in 2002–2009, data from SuperWASP (<https://wasp.cerit-sc.cz/form>) cover the interval between 2004–2008, data obtained in 2008– were taken from the KWS database (Kamogata / Kiso / Kyoto Wide-field Survey, <http://kws.cetus-net.org/~maehara/VSdata.py?object>).

```latex
\begin{table*}[!th]\centering
\vspace*{0.5cm}
\caption{
G- and M-type dwarfs with active cycles determined from new photometric data series.
}
\vspace*{0.5cm}

\begin{tabular}{l|c|l|c|r|c|c|l|l}
Star & B-V & Sp type & $P_\textsl{rot}$, d & $P_\textsl{cyc}$, yr & $A_\textsl{cyc}$, mag & Time span & Variable type & Sources\\
\hline
BE Cet & 0.66 & G2.5 & 7.76 & 6.8 $\pm$ 0.1 & $<$ 0.04 & 1986–2018 & BY Dra & 2, 4, 5\\
DX Leo & 0.77 & G9 & 5.42 & 3.3 $\pm$ 0.1 & 0.03 & 1989–2018 & BY Dra & 1, 2, 4, 5\\
V1005 Ori & 1.37 & M0 & 4.40 & 20.4 $\pm$ 0.2 & 0.12 & 1976–2018 & BY Dra & 1, 2, 4, 7\\
DT Vir & 1.48 & M0 & 2.88 & 31.6 $\pm$ 0.4 & 0.15 & 1911–2018 & fl & 1, 2, 4, 6\\
V647 Her & 1.64 & M3.5 & 1.09 & 76.1 $\pm$ 0.6 & 0.30 & 1939–2018 & fl & 1, 2, 3, 4, 6\\
V577 Mon & 1.69 & M4.5 & 8.00 & 40 $\pm$ 0.4 & 0.50 & 1909–2018 & fl & 1, 2, 4, 6\\
DX Cnc & 2.06 & M6.5 & 0.46 & 2.7 $\pm$ 0.1 & 0.10 & 1984–2018 & fl & 2, 4, 6\\
V374 Peg & 1.26 & M3.5 & 0.44 & var & 0.10 & 1989–2017 & fl & 1, 2, 3, 4\\
HU Del & 1.65 & M4.5 & 0.35 & var & 0.32 & 1940–2017 & fl & 2, 4, 6\\
\end{tabular}
\vspace*{0.5cm}

{
Remarks. Spectral type, B-V and type of variability are taken from the SIMBAD astronomical database; numbers in column "Sources" mark databases Hipparcos (1), ASAS (2), SuperWASP (3), KWS (4); Messina and Guinan, 2002 (5), Bondar', 2013 (6), Alekseev and Kozhevnikova, 2017 (7).
The type of stellar variability is marked according to the SIMBAD astronomical database: BY Dra-type stars (BY Dra), flare stars (fl).
}
\label{Table1}
\end{table*}
```

The list of stars and used data sources are given in Table 1. Two stars from the list – BE Cet and DX Leo – belong to G-dwarfs. On the interval between 1986–2000 Messina and Guinan (2002) found only short-term cycles for these stars. We obtained data series covering a time span twice as long, until 2018; this allows us to consider the stability of short cycles and to search for possible long-term cycles. The cooler stars in our sample have spectral classes M0V–M6.5V; for some of them the long-term photographic measurements of brightness were used from Bondar' (1995). Two stars, V374 Peg and HU Del, are of interest as fast-rotating fully convective dwarfs. Light curves of V374 Peg from the photometry in 1998–2013 were studied by Vida et al. (2016), but the authors found no cycles on this star. The parameters of cycles (the cycle length and amplitude) given in Table 1 were obtained as a result of the analysis of composite light curves formed by us, including recent data of 2017 and 2018. We made a search for the cycles from periodograms derived by the Yurkevich, Hartley and Scargle methods using the AVE program (<http://www.gea.cesca.es>) as described in (Bondar' 2019).

Note that the errors in the $P_\textsl{cyc}$ values show the accuracy of their determination, but not in the duration of cycles, which may vary up to several years at different epochs.

Variations in light curves of V374 Peg and HU Del are not cyclic, and the parameter $A_\textsl{cyc}$ shows the magnitude of the maximum changes in their yearly mean brightness over the studied time interval. No long-term cycles were found on the BE Cet and DX Leo stars. We confirm that only the short cycles detected by Messina and Guinan (2002) have been found from the 30-year photometric series. Considering parameters of the cycle, we took information about the duration of cycles for V833 Tau (78 yr) on the interval of 1899–2009 from Bondar' (2015) and for YZ CMi (27.5 yr) from data on the time span of 1926–2009 (Bondar', 2019).

The rotation period of the star V647 Her was determined according to the SuperWASP data-base obtained with high temporal resolution. A search for the period was performed using the AVE program by the Hartley and Scargle methods on the JD2453128 – 2453160 (2004) interval containing 1292 points for periods of 0.1–5 days and 0.1–10 days. The periodogram in Fig. 1 (the right panel) shows the presence of a peak corresponding to a period of 1.09 days determined by the Scargle method (Scargle, 1982).

We also confirmed the rotation periods of the stars DT Vir and V374 Peg. According to our analysis of data from the ASAS catalog since 2003 to 2008, the rotation period of the star DT Vir $P_\textsl{rot} = 2.88\;$d is consistent with $P_\textsl{rot} = 2.85\;$d derived by Donati et al. (2008); for V374 Peg, according to data in 2004 from the SuperWASP catalog, we found $P_\textsl{rot} = 0.448\;$d to agree with $P_\textsl{rot} = 0.446\;$d acquired by Vida et al. (2016).

# ROTATION PERIODS AND PARAMETERS OF ACTIVITY CYCLES

Here we consider 65 G–M dwarfs, for which a search for cycles has been performed on long-term intervals. A sample includes 36 stars studied from the photometric series and 29 stars, including the Sun, from the Mount Wilson HK-Project (Oláh et al., 2016). The number distribution of objects by spectral types is as follows: 14 G, 28 K, and 23 M dwarfs.

The rotational periods are known for 59 stars of the sample; cycles were detected or suspected for 62 stars. Fig. 2 (left panel) shows the distribution of periods in bins, each of them is equal to 5 d, the spectral classes are marked with different symbols, the dashed line indicates the total number of stars in the bin. Most of G-type dwarfs (10 objects) have periods from 5 to 20 days and only the star HD 81809 (G2V) rotates slower than the Sun, with a period of 39.3 days. Fast-rotating dwarfs with periods of less than 5 days include 25 stars from the sample; mainly, these are cool dwarfs. Three M-stars have rotation periods of less than 0.5 day. 13 K–M dwarfs fall into the range of rotation periods of 30–60 days, the slowest ones are two objects HD 95735 (M2V) and HD 16160 (K3V) with rotation periods of 54 and 57 days, respectively.

The distribution on the length of cycles is shown in Fig. 2 (right panel), bins correspond to 7 years. The cycle durations of 50% of stars are from 7 to 21 years that is comparable to the solar spot cycle. The long cycles, up to 80 years, are detected in fast-rotating K–M dwarfs. The length of dominant cycles for 15 stars is short, less than 7 years. The long-term brightness variations of fully convective M-dwarfs HU Del and V374 Peg with $P_\textsl{rot} <0.5$ d showed no certain periodicity. The light curves for some stars and the results of a search for the cycles are given in (Bondar', 2019).

The amplitude of the dominant cycle in the solar-type stars does not exceed 0.2$^m$, and for the fast-rotating cool dwarfs it can increase up to 0.5$^m$ or more; that indicates on the large areas of active regions. As follows from Fig. 3, K-dwarfs with any rotation periods are distinguished from G and M dwarfs by such parameters as cycle duration and amplitude; this leads to the conclusion that their level of activity is the highest among dwarfs on or near the main sequence.

# FEATURES OF THE DEVELOPMENT OF CYCLES IN G–M DWARFS

The development of activity cycles in the solar-type stars mainly depends on their rotational velocity. On the diagram $P_\textsl{cyc}$ and $P_\textsl{rot}$ (Bondar', 2019) stars are divided into two sequences near a value of $P_\textsl{rot} = 5\;$d. For stars with $P_\textsl{rot} > 5$ d, the cycle length increases if a rotation period becomes slower. This sequence is not homogeneous, since it includes objects with periods from several days to tens of days, i.e. having different ages. For the older stars, rotating slower than the Sun, the above tendency is less pronounced than for the young G and K dwarfs with $5 < P_\textsl{rot} < 20$ days. For these stars cycles of more than 22 years have not been detected yet, but among young K dwarfs there are stars with cycles of about 40 years. The Sun is located between groups of old and young stars.

The youngest stars with $P_\textsl{rot} <5$ d show no any definite connection between the rotation and cycle length. Baliunas et al. (1996) derived a relation between the observed parameters $P_\textsl{cyc}$ and $P_\textsl{rot}$ and the theoretical dynamo number $D$ expressed as $P_\textsl{cyc}/P_\textsl{rot} \sim D^i$. Later, such a relation was considered in several papers for different samples of stars (Oláh et al., 2009; Savanov, 2012; Vida et al., 2014; Katsova et al., 2015; Suárez Mascareño et al., 2016). Vida et al. (2014) have found for all the long and short cycles the value of the slope $i \sim 0.77$, which agrees with Baliunas et al. (1996), but M dwarfs have been excluded from the analysis.

Suárez Mascareño et al. (2016) and Savanov (2012) have considered only M stars; the stellar cycles were determined from observations on intervals of less than 10 years. The obtained large values of $i$ indicate that for the cooler dwarfs such a relation is weaker or absent.

Figure 4 shows the essential difference between the slopes of regression lines for the entire sample of stars (solid line, $i = 1.25 \pm 0.068$) and for the dominant cycles of stars from Oláh et al. (2016) (dashed line, $i = 0.76 \pm 0.15$). The differences in slopes are caused not only by possible uncertainties in the determination of cycle parameters, but are related to a large fraction of rapidly rotating dwarfs in our sample. However, there is still no a sufficient number of stars with $P_\textsl{rot} <5$ d to analyze such a relation for each spectral group.

# CONCLUSION

The stars on or near the main sequence are related by structure (they have convective shells) and the observed manifestations of activity provided by the presence of magnetic fields. Different methods of observations on long-term time intervals reveal slow quasi-periodic changes in brightness of active dwarfs on intervals compared to the duration of the solar cycle. Such a pattern of changes is caused by the evolution of surface spots. For a number of stars, the analogy between the 11-year Schwabe cycle and the dominant cycles of stars is confirmed by the long-term observations. Multiple cycles, the butterfly effects were found in some stars studied in detail. The longer cycles, by several decades, are assumed to be caused by a change in spottedness of the stellar surface. Confident evidence for the evolution of spots on such time spans may be obtained from the long-term observations in different photometric bands. Modern methods for measuring magnetic fields make it possible to detect magnetic spots of different polarities, and from the many-year data one can identify cycles similar to the Hale cycle. The possibility of detecting such cycles on the stars YZ CMi, DT Vir, DX Cnc and others is discussed by Shulyak et al. (2015). Differences in the duration of dominant cycles are associated with differences in the inner structure of stars, their physical parameters, and rotational periods. Certain results on the relation between characteristics of stars and their activity level were obtained for the solar-type stars (Oláh et al., 2016; Fabbian et al., 2017) from many-year observations of indicators of activity in the chromosphere and photosphere. Later stars are mainly investigated on time intervals of about 10 years.

On the basis of long-term observation series, features in the development of activity cycles were considered for 65 G–M dwarfs with rotation periods from fractions of a day to several tens of days. No cycles of more than 40 years were found in the solar-type stars; for most of them, activity cycles are in the range from 7 to 21 years, i.e. compared to the 11-year solar cycle, also having du-ration from 7 to 18 years at different epochs. Cycles of the fast-rotating dwarfs can reach 80 years. So far, such long cycles have been suspected in two stars – V833Tau and V647 Her. About 30% of the considered stars show only short dominant cycles of less than 7 years.

For rapidly rotating dwarfs, there is no definite relationship between the cycle length and the rotational period, which is typical for stars with rotation periods longer than 5 days. Their cycle length increases with a deceleration of rotation. But in this group of stars, the longer cycles are observed in K-dwarfs rotating faster than the Sun.

# ACKNOWLEDGEMENTS

The authors are grateful to the referee for useful remarks and suggestions.

We used data from the following databases: ASAS and SuperWASP catalogs, the International Variable Star Index (VSX) database, the Kamogata – Kiso – Kyoto Wide-Field Survey data. The authors are thankful to all the staff for supporting these databases and providing access to data.

# FUNDING

This work was carried out with the partial support of the Russian Foundation for Basic Research (grants 18-52-06002 Az$_a$ and 19-02-00191a).

# CONFLICT OF INTEREST

The authors declare that they have no conflicts of interest.