Planetary Nebula Photometry and Distances

Planetary Nebulae in Face-On Spiral Galaxies.
I. Planetary Nebula Photometry and Distances

Kimberly A. Herrmann11affiliation: Visiting Astronomer, Cerro Tololo Inter-American Observatory. CTIO is operated by AURA, Inc. under contract to the National Science Foundation. 22affiliation: Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the National Optical Astronomy Observatories. , Robin Ciardullo11affiliation: Visiting Astronomer, Cerro Tololo Inter-American Observatory. CTIO is operated by AURA, Inc. under contract to the National Science Foundation. 22affiliation: Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the National Optical Astronomy Observatories. Department of Astronomy & Astrophysics, The Pennsylvania State University
525 Davey Lab, University Park, PA 16802,
John J. Feldmeier22affiliation: Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the National Optical Astronomy Observatories. Department of Physics & Astronomy, Youngstown State University, Youngstown, OH 44555-2001 Matt Vinciguerra22affiliation: Visiting Astronomer, Kitt Peak National Observatory, National Optical Astronomy Observatories, which is operated by the Association of Universities for Research in Astronomy, Inc. (AURA) under cooperative agreement with the National Science Foundation. The WIYN Observatory is a joint facility of the University of Wisconsin-Madison, Indiana University, Yale University, and the National Optical Astronomy Observatories. Department of Astronomy & Astrophysics, The Pennsylvania State University
525 Davey Lab, University Park, PA 16802

As the first step to determine disk mass-to-light ratios for normal spiral galaxies, we present the results of an imaging survey for planetary nebulae (PNe) in six nearby, face-on systems: IC 342, M74 (NGC 628), M83 (NGC 5236), M94 (NGC 4736), NGC 5068, and NGC 6946. Using Blanco/Mosaic II and WIYN/OPTIC, we identify 165, 153, 241, 150, 19, and 71 PN candidates, respectively, and use the Planetary Nebula Luminosity Function (PNLF) to obtain distances. For M74 and NGC 5068, our distances of  Mpc and  Mpc are the first reliable estimates to these objects; for IC 342 ( Mpc), M83 ( Mpc), M94 ( Mpc), and NGC 6946 ( Mpc) our values agree well with those in the literature. In the larger systems, we find no evidence for any systematic change in the PNLF with galactic position, though we do see minor field-to-field variations in the luminosity function. In most cases, these changes do not affect the measurement of distance, but in one case the fluctuations result in a 0.2 mag shift in the location of the PNLF cutoff. We discuss the possible causes of these small-scale changes, including internal extinction in the host galaxies and age/metallicity changes in the underlying stellar population.

distance scale — galaxies: distances and redshifts — planetary nebulae: general

1 Introduction

Planetary Nebulae (PNe) are dying low-mass () stars composed of an expanding rarefied shell of gas surrounding a soon-to-be white dwarf (e.g., Kwok et al., 2003). The central stars of PNe can be among the brightest stellar objects in a galaxy with maximum luminosities above 6000  (Vassiliadis & Wood, 1994). Due to the physics of photoionization, up to 10% of this flux can be reprocessed into a single emission line of doubly ionized oxygen at 5007 Å, transforming the objects into very bright monochromatic stars (Jacoby, 1989; Jacoby & Ciardullo, 1999; Marigo et al., 2004). Moreover, because all stars with main-sequence masses between and eventually evolve into PNe, these objects are plentiful in all stellar populations with ages between and  yr.

The above properties make PNe extremely useful probes of (1) the internal kinematics of spiral (Merrett et al., 2006; Ciardullo et al., 2004; Douglas et al., 2000), elliptical (Douglas et al., 2007; De Lorenzi et al., 2008), and interacting (Durrell et al., 2003; Peng et al., 2004) galaxies, (2) the chemical evolution of the Local Group (e.g., Dopita et al., 1997; Richer et al., 1998; Magrini et al., 2007), (3) the stellar populations of early-type systems (Buzzoni et al., 2006; Ciardullo, 2006), (4) the extragalactic distance scale (e.g., Ciardullo et al., 2002; Ciardullo, 2005), and (5) the dynamical evolution of galaxy clusters (Feldmeier et al., 2004; Arnaboldi et al., 2004; Gerhard et al., 2007). We seek to use PNe as kinematic test particles in face-on spiral galaxies in order to use their motions to measure the distribution of mass in spiral disks. (This analysis will be the subject of the second paper in this series, hereafter referred to as Paper II.) The first step in doing this is to identify large samples of PNe for analysis. With these data, we can also measure the PNe’s photometric properties, use their luminosity function to estimate distance, and explore how the PNLF behaves in a range of star-forming environments.

In this paper, we describe an [O iii] and H survey for planetary nebulae in six nearby, face-on spiral galaxies. In §2 we present the details of our observing runs and image reduction. In §3 we explain how we identified our PN candidates and measured their magnitudes. We discuss the discrimination of possible contaminants in §4, including how we used photometry and follow-up spectroscopy to eliminate H ii regions, supernova remnants, possible background galaxies, and four asteroids from the sample. In §5, we show how we determined the completeness limits of our samples and also compare our narrow-band selected PNe sample in M94 to a sample of planetaries found via slitless spectroscopy. In general, we find good agreement between the two data sets, although the additional information provided by our H photometry does allow us to eliminate a few compact H ii regions from consideration. In §6, we derive the PNLF distance to each of our six galaxies and compare these new measurements to estimates in the literature. In all cases, our values are consistent with previous distance determinations, and in three cases, our measurements are significantly more precise than those of previous studies. In §7 we study the shape of the PNLF, and test for differences between various PN subsamples. We show that small variations in the shape of the PNLF are not uncommon, although these changes rarely affect our ability to measure distance. We conclude by discussing the possible implications of these fluctuations.

2 Observations and Image Reduction

We chose for our study six nearby spiral galaxies, all with low recessional velocities ( km s), all relatively face-on (), and, with one exception, all with large optical angular diameters (). These galaxies are listed in Table 1, along with their photometric properties and metallicities. The two southernmost objects in the sample, M83 and NGC 5068, were imaged with the Mosaic II camera on the CTIO 4 m Blanco telescope. This instrument is composed of 8 SITe CCDs, giving it a pixel scale of pixel and a field-of-view large enough to encompass the entire galaxy (i.e.,  or 10 disk scale lengths for M83 and more for NGC 5068). The remaining galaxies were imaged with OPTIC, an orthogonal transfer CCD camera on the WIYN 3.5 m telescope at Kitt Peak (Howell et al., 2003; Tonry et al., 2002). OPTIC, which consists of two Lincoln Lab CCID28 chips separated by , produces a pixel scale of pixel and a field-of-view, , that is slightly smaller than the typical size of our galaxies. As a result, our surveys with this instrument extended over 1.2, 1.8, 3.5, and 7 disk scale lengths for IC 342, NGC 6946, M74, and M94, respectively. Although OPTIC is designed to track rapid image motion by real-time clocking of charge in both the vertical and horizontal direction, the lack of suitably bright field stars prevented us from using this feature. Table 2 presents a log of our observations.

For our survey, we followed the standard observing procedure for imaging extragalactic planetary nebulae (e.g., Jacoby et al., 1992; Feldmeier et al., 1997). First we chose a narrow-band filter to isolate the [O iii] emission line at the redshift of the target galaxy. For M83 and NGC 5068, this filter was custom made, with a central wavelength of and a full-width-half-maximum (FWHM) of 46 Å in the converging f/2.87 beam of the telescope. For the remaining galaxies, standard Kitt Peak 50 Å FWHM filters were used. These filters are listed in Table 2; models for their transmission curves at dome temperature are displayed in Figure 1 along with the wavelength of [O iii] shifted to the systemic velocity of each target galaxy.

Figure 1: Transmission curves for the [O iii] filters used in this survey. These curves represent the filter at the outside temperature and in the converging beam of the telescope. The wavelengths of [O iii] at the systemic velocity of each galaxy are marked.

Each galaxy was imaged via a series of 30 minute exposures, with dithering sequences sometimes used to fill in the gaps between chips. For M83 and NGC 5068, corresponding 1 minute exposures through a standard filter served as off-band images; for the WIYN observations, continuum observations were obtained via a set of 7.5 min exposures through the WIYN 32 filter (, FWHM = 230 Å). In general, these off-band exposures were scaled to go at least 0.2 mag deeper than their on-band counterparts to avoid false detections at the frame limit. To help discriminate PNe from (lower-excitation) H ii regions, each galaxy was also imaged in H. For M83 and NGC 5068, the H filter had a central wavelength of and a FWHM of 80 Å, while for the other four galaxies, and FWHM = 46 Å. Finally, for IC 342, M74, M83, and NGC 5068, images were also taken through a standard filter. These served as off-band images for the H exposures.

All images were reduced using packages in IRAF111IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. (Valdes, 1998). The Mosaic II exposures of M83 and NGC 5068 were reduced entirely within the mscred package: the data were trimmed, bias subtracted, and flatfielded via ccdproc, placed on an astrometric system by applying mscmatch to 2800 stars from the USNO-B1.0 catalog (Monet et al., 2003), projected onto the tangent plane with mscimage, sky subtracted with mscskysub, aligned with mscimatch, and co-added (with cosmic-ray removal) with mscstack. During this process, special care was taken to define the images’ world coordinate system, as we found it necessary to run msccmatch separately on each individual readout section of each chip (i.e., 16 runs per exposure) to obtain a proper stacking and an astrometric precision of . We found the NOAO Deep Wide-Field Survey MOSAIC Data Reductions website222 to be very useful.

The OPTIC images of IC 342, M74, M94, and NGC 6946 were reduced using a combination of IRAF’s imred and mscred packages. First, all individual exposures (including biases and flats) were converted to Multi-Extension FITS (MEF) format using mkmsc; this enabled us to use the mscred package to correct for the four different overscan regions and the gap between chips. The OPTIC reduction notes333 available off the WIYN website were very useful for this first step. Next ccdproc was used to trim and flatfield as above. Again, special care was needed to define the frames’ world coordinate system, as each readout section (4 per frame) required a manual initialization with a centrally located star (via ccsetwcs) and then a refinement (with ccfind and ccmap). For M94, only 20 USNO-B1.0 stars (Monet et al., 2003) per quadrant were available for the astrometric solution; for M74, this number was closer to 40. (Since IC 342 and NGC 6946 are both close to the Galactic plane, more than 100 stars per quadrant were available for their solutions.) In all cases, our final astrometry attained a precision of . Finally, mscimage was run to project the frames onto the tangent plane, linmatch was used to determine the scaling between exposures, and imcombine was used to remove cosmic rays and stack the images.

3 Planetary Nebula Candidate Identification and Photometry

PN candidates were identified by examining the stacked [O iii] images, their H and off-band counterparts, and difference images created by scaling, then subtracting the off-band images from their on-band frames. PN candidates had to (1) have an image profile consistent with that of a point source, (2) be present in [O iii] but invisible in the off-band image (or extremely weak in ), and (3) be invisible or weak in H (and when available). This latter condition served to eliminate H ii regions and supernova remnants from the sample (see below).

Photometric measurements of all the candidates in [O iii] and H were performed with the IRAF version of DAOPHOT (Stetson, 1987; Stetson et al., 1990). First, the point spread function (PSF) fitting routines in allstar were used to define the PN magnitudes on the difference images relative to those of field stars on the on-band frames. These results were then checked by using DAOPHOT’s substar option to subtract off a scaled-PSF representation of each PN from its position on the on-band and difference frames. Although the majority of these subtractions produced little or no residual, a few objects, usually in or near the galaxy’s spiral arms, were found to be significantly under- or over-subtracted, sometimes by 0.5 mag. The cause of this error was an incorrect estimate of the background: in regions of active star formation, bright emission from H ii regions, supernova remnants, and the diffuse ISM can play havoc with the underlying “sky,” particularly in H. When this occurred, the PN magnitudes were manually adjusted until the residuals on the star-subtracted images appeared reasonable.

Next our relative [O iii] and H magnitudes were placed on the standard AB system by comparing large-aperture measurements of our field stars to similar measurements of spectrophotometric standards. For our southern hemisphere measurements, these standards came from Stone & Baldwin (1983); the stars in the north were drawn principally from the list of Stone (1977). Finally, the narrow-band AB magnitudes were converted to monochromatic fluxes using models for the filter transmission curves in the converging beams of the telescopes (Figure 1), knowledge of the galaxies’ systemic velocities (listed in Table 1), and the photometric procedures for emission-line objects described by Jacoby et al. (1987). In the case of our [O iii] measurements, these fluxes were converted to [O iii] magnitudes using


where is given in ergs cm s (Jacoby, 1989). At this time, we also used the photometric uncertainties derived by DAOPHOT to determine the mean photometric error of our measurements as a function of [O iii] brightness. This photometric error function, which is given in Table 3, was used to help determine the PNLF distances of §4.

4 Discriminating Pne From Contaminants

Planetary nebulae are not the only sources that are detectable via our on-band/off-band survey technique. Deep images through an [O iii] filter will also detect H ii regions, supernova remnants, background [O ii] emitting galaxies (at ), high-redshift Ly sources (at ), and even asteroids. We discuss each of these below.

4.1 H ii Regions

At the typical distance of our galaxies (5 Mpc), corresponds to 25 pc. Consequently, at groundbased resolution, compact H ii regions will be unresolved, and can be misclassified as planetaries. Fortunately, a quantitative way exists to eliminate most of these contaminants using two of a PN’s brightest emission lines. Ciardullo et al. (2002) have shown that bright PNe populate a distinctive region in [O iii]-H emission-line space. Specifically, when [N ii]) is plotted against absolute [O iii] magnitude, true PNe occupy a cone defined by


This equation, which is valid over the top 4 mag of the PN luminosity function, is equivalent to the statement that the ratio of [O iii] to depends on the luminosity of the planetary, via


where  ergs s is the absolute [O iii] luminosity of a PN at the bright-end of the PNLF. Practically speaking, this means that PNe in the top 1.5 mag of the [O iii] luminosity function always have [O iii] brighter than H. This contrasts with the line ratios of the vast majority of H ii regions, which have H much brighter than [O iii] (e.g., Shaver et al., 1983; Kniazev et al., 2005; Peña et al., 2007).

For three of our galaxies, our narrow-band photometry was sufficient to locate each PN candidate in the [O iii] -H emission line diagram. Candidates which failed to satisfy the criterion given by equation (2) within the uncertainties were excluded from further consideration. Unfortunately, due to clouds and poor seeing, the H images of the other three galaxies (IC 342, M74, and NGC 6946) did not go deeply enough to detect the H emission from all the possible PNe. Specifically, in IC 342 only 55% of our [O iii] point sources were detected in H; in M74, the H recovery fraction was just 20%, and in NGC 6946, none of our bright PN candidates were detected in H. In the case of IC 342, its nearness ( Mpc; see §4) and extreme extinction ( mag) helped to enhance H with respect to [O iii] , somewhat lessening the shallow H imaging problem.

For IC 342 and M74, we compensated for the lack of H photometry using spectroscopy. Approximately 68% of the PN candidates in IC 342 have spectra obtained from the Hydra spectrograph on the WIYN telescope. (A full description of these data is presented in Paper II.) For M74, the fraction of PNe with spectra is 65%. Although these data were obtained for the purpose of obtaining radial velocities and only covered the wavelength range 4500 Å to 5500 Å, it is possible to use them to estimate the H flux of each candidate.

To do this, we began by subtracting the continuum from each spectrum, and measuring the relative line fluxes of [O iii] , [O iii] , and H. Appropriately, the ratio of the two oxygen lines was (on average), three-to-one, indicating that the response function within this small wavelength range was approximately flat. We then examined the throughput of the spectrograph using observations of a standard star. Again, the data indicated only a slight (8%) decrease in efficiency between the wavelengths of [O iii] and H. (This is close to the 5% drop expected from the blaze function of the grating.) Using this knowledge, we derived the true [O iii] -H ratio for each object, then scaled these values to H, using an estimate of foreground Galactic extinction (Schlegel et al., 1998), a Cardelli et al. (1989) reddening law with , and an expected intrinsic H to H ratio of 2.86 (Brocklehurst, 1971). For cases where the H flux was not detected photometrically, this spectroscopic measurement was used to estimate and eliminate compact H ii regions from the sample. We note that this technique may overestimate the [O iii]-to-H flux ratio for PNe that are heavily reddened by circumstellar dust. However, these objects are likely to fall well below the PNLF’s bright-end cutoff, and not be an important part of the sample.

In the extreme case of NGC 6946, where no spectroscopy was available, the best we could do is place an upper limit on the H brightness of our PN candidates. To do this, we began by using the brightest point source magnitude returned by allstar during our analysis of the H frame. We then scaled the frame’s PSF to this magnitude, and used substar to subtract off this “brightest star” from the location of each PN. We then visually inspected the resulting image. As expected, due to differences in the underlying background, this procedure occasionally resulted in an oversubtraction. When this occurred, we manually adjusted the limiting allstar magnitude until the subtraction residual showed no sign of a flux deficit. This magnitude was then used to define the upper limit to the object’s H flux. Unfortunately, even with these limits, it was impossible to discriminate PNe from H ii regions with sizes under 25 pc. Thus, all the PN identifications in this galaxy must be considered tentative.

4.2 Supernova Remnants

Although compact supernova remnants (SNRs) are much rarer than PNe, they may present a problem at the very bright end of the PNLF. This is particularly true for the more distant galaxies, where groundbased 1″ seeing limits our spatial resolution to 25 pc. The difficulty is well-illustrated by Roth et al. (2004), who used spatially resolved spectrophotometry to show that object 276 in the list of M31 bulge PNe presented by Ciardullo et al. (1989) is actually a supernova remnant. Fortunately, most SNRs have line ratios that clearly distinguish themselves from planetary nebulae. For example, of the 60 M31 supernova remnants surveyed by Galarza et al. (1999), all but one has [O iii] fainter than H.

More specifically, we can compare the photometric properties of PNe directly with those of supernova remnants using the optical sample of M83 SNRs created by Blair & Long (2004). Their catalog contains 71 supernova remnants, selected on the basis of strong [S ii] relative to H. Of these, only three fall within of one of our PN candidates, and in each case, there is a better, more distinctive SNR candidate nearby. None of the other 68 SNRs are within of a PN. Blair & Long (2004) also presented 18 additional [O iii]-bright objects, which they believe might be extremely young SNRs. Interestingly, all of these sources are significantly fainter than the prototype of the class, the Milky Way remnant Cas A. They are, however, comparable in brightness to the well-known SMC remnant, 1E 01027219, and in theory, such objects can be confused with PNe. Seventeen of the Blair & Long (2004) [O iii]-bright sources are more than from any of our candidates and thus not on our list. The remaining source (O8) is coincident with our M83 PN candidate 10, which satisfies our selection criteria in every respect. Though we targeted this object with the Hydra spectrograph, the resulting spectrum has very strong H emission. While there is virtually no H emission at the location of PN 10, there is a strong H emitter within . We believe our spectrum is contaminated by this nearby H emitter due to a slight misplacement of the fiber by Hydra.

In addition, there are 17 historical supernovae in our target galaxies444The Central Bureau for Astronomical Telegrams (CBAT) List of Supernovae website ( was used to obtain the locations of the 17 historical SNe described here.. An examination of our images demonstrates that none of the remnants from these sources are in our PN sample. We did find an emission line object within of the location of both SNe in M74 (SN 2003gd and 2002ap) but found nothing of interest at the locations of SNe 1983N, 1968L, and 1957D in M83 and only well-resolved emission-line sources within 5″ of the locations of SNe 1950B, 1945B, and 1923A (also in M83). Of the nine historical SNe in NGC 6946, we found emission line objects near three (1980K, 1948B, and 1939C), some H ii regions and blue stars near three others (2002hh, 1969P, and 1917A), and only star-like objects around the remaining three (2008S, 2004et, and 1968D).

4.3 Background Galaxies

Star-forming galaxies at have their [O ii] doublet redshifted into the bandpass of our [O iii] filter. Consequently, it is theoretically possible for small ( kpc) star-bursting objects to be mistaken for planetaries. In practice, however, such a misidentification is extremely unlikely. At , our survey volume is 8500 Mpc for the CTIO 4 m Mosaic II data, and only 700 Mpc for the galaxies observed with OPTIC. According to the luminosity functions of Hogg et al. (1998), Gallego et al. (2002), and Teplitz et al. (2003), this implies that our M83 and NGC 5068 fields contain less than 200 [O ii] galaxies above our completeness limit, and that our OPTIC images contain no more than 20 [O ii] emitters. Less than 2% of these objects will have the extreme rest frame equivalent widths (60 Å) needed to appear on our on-band frame, but be completely invisible in the continuum (Hogg et al., 1998). Moreover, many of these objects will be (at least) marginally resolved. Thus, the number of background [O ii] interlopers in our sample must be negligible.

In principal, Ly emitting galaxies (LAEs) are another source of concern. These objects can have extremely large (100 Å) observer-frame equivalent widths, and, in the absence of deep continuum images, can easily be mistaken for PNe. Large numbers of such objects have been discovered in deep [O iii] surveys of Virgo’s intracluster space (Feldmeier et al., 2004), and one has even been found masquerading as a PN in the M51/NGC 5195 system (Durrell et al., 2003). Fortunately, the surveys described here do not reach the flux levels needed to sample this population. None of the 162 LAEs found by Gronwall et al. (2007) in their deep 0.28 deg survey of the Extended Chandra Deep Field South would have been detected on our images of IC 342, M94, NGC 5068, and NGC 6946, and only their two brightest sources (one an AGN) would have been seen in M83. Our survey of M74 does go deep enough to detect high-redshift Ly emitters, but this galaxy was surveyed with the small field-of-view of OPTIC. Consequently, the luminosity function of Gronwall et al. (2007) predicts that only 1 of our PN candidates should be a LAE. The spectroscopy of Paper II confirms this result: Ly emitters are not an important source of contamination in our survey.

4.4 Asteroids

Since both M74 and M83 are located at low ecliptic latitudes (5 and 18, respectively), one might expect to find asteroids in the images of these two galaxies. Asteroids show up as trails on the long on-band exposures since they move quickly (0 min) across the field of view and as nearly point sources on the short off-band exposures. This makes them very easily distinguishable from PNe. We found no asteroids on our 10 M74 images but found four on our 37 images of M83. The asteroid trails in the stacked on-band image correspond to asteroid sky position angles and speeds of 323 at 0 min, 260 at 0 min, 312 at 0 min, and 299 at 0 min, from brightest to faintest. A cursory check of the IAU Minor Planet Center website suggests the brightest asteroid on our images might be minor planet (16522) Tell, which moves 0 min with a PA of and is expected to have been 45″ from the observed location on each of our four frames and have virtually the same motion. However, a more detailed analysis would be necessary to verify this possibility.

5 Completeness of the Final Samples

After the application of the procedures of §4, we were left with the following numbers of strong PN candidates in each galaxy: 165 in IC 342, 153 in M74, 241 in M83, 150 in M94, 19 in NGC 5068, and 71 in NGC 6946. The positions of these objects are displayed in Figure 2; Figure 3 shows the location of the objects in the [O iii] -H+[N ii] emission-line space. For the latter diagram, we have plotted the photometric line ratio whenever possible; in cases of unreliable H photometry, either the spectroscopic estimate or a lower limit has been displayed. A table of PN positions, fluxes, line-ratios, and radial velocities will appear in Paper II. Note that it is very difficult to find PNe near the centers of galaxies and within the spiral arms, due to the bright stellar background and diffuse emission associated with regions of active star formation. Since the ease of PN discovery increases as one moves to the lower-surface brightness outer regions of the galaxies, many of our PN candidates are at large galactocentric radii: for example, the percentages of PNe beyond 2 -band scale lengths are 43% in M74, 93% in M83, and 73% in M94.

Figure 2: Our [O iii] images of the six galaxies studied in this survey. North is up and east is to the left; the positions of the PN candidates are colored such that red, orange, yellow, green, and blue represent objects in the top 0.5, 1.0, 1.5, 2.0, and 2.5 mag of the PNLF, respectively. The major and minor axes are shown as white lines. For IC 342, M74, M83, and M94, the inner ellipse marks the high surface-brightness region of the galaxy within which PN identifications were severely incomplete, while the outer ellipse shows the dividing line between our inner and outer samples. For M94, the green rectangles indicate the approximate field of view of the Douglas et al. (2000) slitless spectroscopy survey, the black crosses indicate candidates in both samples, and the yellow diamonds indicate their candidates that are not present in our sample.
Figure 3: The flux ratio between [O iii] and H+[N ii] for the planetary nebula candidates of our six galaxies. These ratios have been corrected for foreground Galactic reddening using the extinction estimates of Schlegel et al. (1998), but the effects of internal and circumstellar reddening are still present in the figure. The dashed lines indicate values determined from [O iii] -H spectroscopy; arrows indicate lower limits. The black cone-shaped outline illustrates the locus of Local Group PNe (Ciardullo et al., 2002). Since the vast majority of H ii regions have H brighter than [O iii] , a diagram such as this is extremely useful for discriminating PNe from low-excitation interlopers.

The next step in our analysis was to determine our level of photometric completeness. To do this, we performed a series of artificial star experiments using the addstar task within DAOPHOT. First we examined each image and found the locations where [O iii] emission from H ii regions, supernova remnants, and shocks made PN detections problematic. These regions were masked out and excluded from further analysis. At this time, we also excluded bad CCD columns, detector edges, and the regions between the two CCDs of OPTIC. This removed 25% of the disk area of M74, M83, and M94 within the galaxies’ inner three disk scale lengths, and 16% of the disk area of IC 342 and NGC 6946 within one scale radius. (For the small galaxy NGC 5068, no exclusion regions were necessary.) Next, we scattered a small number of artificial stars of varying [O iii] brightness throughout our on-band image. By adding no more than 10% of the original number of detected planetaries to each frame, we ensured that image crowding was not a factor in our analysis (Fleming et al., 1995). Then we searched the frame for these artificial stars using the daofind task within DAOPHOT and a threshold value of three times the background sky. Each recovery (and non-recovery) was noted, and the process was repeated several thousand times until the completeness fraction, as a function of [O iii] magnitude and position in the galaxy, was well determined. Finally, we fit this function to an equation of the form


(Fleming et al., 1995). For purposes of this paper, we define our limiting magnitude for PN detections to be the 90% completeness limit. This number is listed in Table 2.

5.1 Comparison with Previous PN Samples: M94

One of the galaxies in our sample has previously been surveyed for planetary nebulae. M94 was observed with a dual-beam, “counter-dispersed imaging” spectrograph, designed to create H images with one arm, and measure PN velocities via [O iii] slitless spectroscopy with the other (Douglas et al., 2000). Although PN photometry through an instrument such as this is difficult and subject to a range of systematic errors, it is still meaningful to compare the two samples.

Douglas et al. (2000) studied two regions along M94’s major and minor axes, and found 53 and 14 PN candidates, respectively. Their fields are displayed in Figure 2. Of the 67 PN candidates found via slitless spectroscopy, 44 are located within of one of our candidates, 16 are likely H ii regions and not included in our sample, one is outside our field-of-view, one appears to be a duplicate, and five are below the limit of our completeness and not present in our sample. Conversely, 22 of the PNe found in our survey are not cataloged by Douglas et al. (2000). Of these, eight are near the edge of the spectrograph’s field, and thus may have been missed for that reason. Of the remaining 14, two are within 0.5 mag of the bright end of the PNLF, nine are between 1 and 1.5 mag down the luminosity function, and three are extremely faint.

A comparison of the objects’ [O iii] magnitudes shows generally good agreement. As described above, 44 PNe were detected in both surveys. Of these, four have slitless spectroscopy magnitudes that are more than 0.8 mag fainter than our measurements, and two have no measured magnitude at all. However, the remaining 38 have magnitudes that are consistent with our own. The scatter between the two data sets is  mag, and the systematic zero-point offset between the two samples is only  mag, with our values being slightly fainter.

Finally, when we examine the objects that we tagged as likely H ii regions, we again find good agreement. Of the 16 sources, 10 were classified by Douglas et al. (2000) as H ii regions based on their extremely bright [O iii] magnitudes. Conversely, one object that Douglas et al. (2000) considered too bright to be a PN is classified as a planetary in our sample. This disagreement arises from the fact that Douglas et al. (2000) assumed a distance of 6 Mpc in their analysis. As we describe below, our distance to the galaxy is substantially smaller than this value.

Of the 44 candidates in common between the two samples, we have spectra of 37. A comparison of the velocities from these two samples will be included in Paper II.

6 Pnlf Distances and Comparison to Existing Estimates

The Planetary Nebula Luminosity Function is an excellent standard candle: observations in 50 elliptical, spiral, and irregular galaxies demonstrate that the bright-end of this function is well represented by an equation of the form


where is the absolute magnitude of the most luminous planetary (Ciardullo et al., 1989; Jacoby et al., 1992). An external calibration using galaxies with well-determined Cepheid distances defines the nominal value of to be , although this cutoff appears to fade in low-metallicity objects (Ciardullo et al., 2002).

To derive the PNLF distances, we began by convolving the empirical PNLF of equation (5) with the photometric error functions defined in Table 3. We then fit those PNe brighter than our 90% completeness limit to the convolved function via the method of maximum likelihood (Ciardullo et al., 1989). (Since the PNLF of M83 appears to deviate from this law at faint magnitudes, we truncated its fit at .) Foreground Galactic extinction was taken into account by adopting the extinction estimates of Schlegel et al. (1998) and a Cardelli et al. (1989) reddening law, i.e., . This analysis yielded the galaxies’ distances, along with their formal fitting errors. To these errors, we then added (in quadrature) the possible systematic uncertainties associated with our photometric standard star measurements (0.02 mag), the filter response function (0.02 mag), and the foreground extinction estimate (; see Table 4). Finally, a Kolmogorov-Smirnov test was used to determine whether the best-fit empirical function was indeed a valid representation of the observed data. In all cases, this is confirmed to be true: the worst fit is for M74, which has a strange dearth of PNe between 0.2 and 0.4 mag below the PNLF cutoff (see §7). However, even for this galaxy, the null hypothesis can only be excluded at the 90% confidence level. Figure 4 displays the observed luminosity functions for each galaxy, along with the best fit models.

Figure 4: The [O iii] PNLFs of the six galaxies analyzed in this study. The curves represent the best-fitting empirical PNLFs convolved with the photometric error function and shifted to the most likely distance of the galaxy. No corrections for reddening have been applied here. The open circles indicate points fainter than the completeness limit.

We note that our analysis of the PNLF ignores the effects of internal extinction. This omission is justified for several reasons. First, Giovanelli et al. (1994) have found that Sc galaxies have very little dust extinction beyond two or three scale lengths of the galactic center. As we mentioned above, many of our PN candidates in M74, M83, and M94 have galactocentric radii larger than this limit. Moreover, Feldmeier et al. (1997) modeled the effects of internal extinction on the PNLF by assuming that the vertical scale heights of dust and PNe in other spiral galaxies are similar to those observed in the Milky Way. They found that the larger scale height of the PNe make the PNLF method relatively insensitive to internal extinction. Specifically, derived PNLF distances are never more than 0.1 mag less than their nominal value, even when the total amount of extinction is large. In theory, this question can be addressed directly, via spectroscopy of the PNe’s H and H emission lines, though the faintness of H and the effects of atmospheric dispersion make the observations challenging.

In the sections below, we compare our PNLF distances to values determined from the analyses of Cepheids, the Tip of the Red Giant Branch (TRGB), Surface Brightness Fluctuations (SBF), the Brightest Blue Super Giants (BBSG), the Brightest (Red and Blue) Supergiants (BSM), and the Expanding Photospheres of Supernovae (EPM). Though some of the galaxies also have Tully-Fisher measurements, the fact that we selected our targets to have low inclination makes these determinations unreliable. Figure 5 and Table 5 compare our distance moduli to measurements in the literature. Note that the zero point for all our PNLF distances is tied to that of Cepheids via Ciardullo et al. (2002). Thus, following Freedman et al. (2001), our distance scale assumes an LMC distance modulus of .

Figure 5: Final PNLF distances compared to existing distances in the literature. In most cases, only the distances with the smallest uncertainties are shown, but a few measurements with large uncertainties and one without quoted errors are provided for comparison. See Table 5 for more information on the distances, including the reference for each value.

6.1 Ic 342

IC 342, properly dubbed the Hidden Galaxy, would be one of the brightest galaxies in the sky except for its location near the Galactic plane () and the large amount of foreground extinction (). The galaxy is part of the Maffei group, a system that garnered some attention due to the suggestion of a past gravitational encounter with the Local Group (Zheng et al., 1991; Valtonen et al., 1993; Peebles, 1994). More recently, multiple studies have presented evidence against this possibility (Krismer et al., 1995; Karachentsev et al., 2003a; Fingerhut et al., 2007). Although the group of 16 galaxies was originally thought to have contained only two large galaxies (IC 342 and the giant elliptical Maffei 1), a recent analysis by Fingerhut et al. (2007) has resulted in the addition of another member, the heavily extinguished (5.6 mag) spiral Maffei 2.

A number of different distance indicators have been used for IC 342. Saha et al. (2002) used two-passband observations of 20 Cepheids to estimate a distance modulus of and a total reddening of . With an updated period-luminosity relation and metallicity measurement, this number increases slightly to (Fingerhut et al., 2007). For comparison, the magnitudes of the three brightest blue and red supergiants gives a distance modulus of (Karachentsev & Tikhonov, 1993, revised by Karachentsev et al. (1997)).

Out of the 165 PNe found in IC 342, the 65 in the top 0.9 mag of the PNLF qualify as our complete sample. Our application of the maximum likelihood method then gives a true distance modulus of ( Mpc). Note that the error in our distance is completely dominated by the uncertainty in the foreground extinction. (Without this term, our error would only be 0.05 mag!) Since the extinction derived by Fingerhut et al. (2007) is completely consistent with that determined from the COBE satellite (Schlegel et al., 1998), we adopt their 0.16 mag uncertainty in our calculation as well.

Figure 5 compares our distance estimate to that of the Cepheids and supergiant stars. Our value is statistically identical to the former value, but substantially larger than that associated with the galaxy’s brightest stars.

6.2 M74

The Perfect Spiral, M74 (NGC 628), is perhaps the most face-on galaxy in the nearby universe, with (Kamphuis & Briggs, 1992). This grand-design system dominates a small group which includes three dwarf irregular objects (Sharina et al., 1996), and plays host to both our nearest Type Ic “hypernova” (Mazzali et al., 2002) and an extremely variable ultraluminous x-ray source (Krauss et al., 2005).

Though M74 is well-studied, it lacks a reliable distance estimate. Values in the literature range from 7 to 20 Mpc (Sohn & Davidge, 1996) and this uncertainty propagates into all aspects of its analysis, including luminosity estimates for its supernovae and x-ray sources. M74’s three least uncertain distances come from the photometry of its three brightest blue supergiants (; Sharina et al., 1996), the photometry of its brightest supergiants of all colors (; Hendry et al., 2005), and the expansion velocity-absolute magnitude relation of the Type II-P supernova SN 2003gd (; Hendry et al., 2005).

Our complete sample of planetaries consists of 82 objects out of a total population of 153. By fitting the top 1.1 mag of the luminosity function, we obtain a distance modulus of , or Mpc. Though this distance differs by only 10% from the values obtained from the Type II-P supernova and the galaxy’s supergiants, it has considerably less uncertainty.

6.3 M83 and NGC 5068

M83 (NGC 5236), the barred Southern Pinwheel, is one of the two central galaxies of the nearby Cen A (NGC 5128) group, an 87 galaxy complex which also includes NGC 5068. Though often viewed as a single system, Karachentsev et al. (2002, 2007) have recently shown that the group actually contains two spatially separated associations, one centered on the interacting elliptical NGC 5128, and the other clustered about M83. M83 itself is one of the great supernova factories in the sky, with six recorded events in the last century.

To date, the two most accurate distances to M83 have come from Cepheids and stars at the tip of the red giant branch. The former method, based on and -band observations of 12 Cepheids with the VLT, gives a distance modulus of (Thim et al., 2003; Saha et al., 2006); the latter, which comes from HST F606W and F814W photometry of halo stars, yields (Karachentsev et al., 2007). A group distance, determined via TRGB measurements to 10 of the galaxy’s surrounding dwarfs, is .

Our PN survey with the CTIO 4-m telescope detected 241 PNe in M83; 164 of these fall in the top 1.3 mag of the PNLF and comprise the complete sample before the luminosity function turns over. Our maximum-likelihood fit to the data yields a distance modulus of , or Mpc. Interestingly, this value lies between the Cepheid and TRGB distances to the galaxy, and is closer to the mean group distance than either of the two direct measurements.

For NGC 5068, 11 of our 19 PNe are above our completeness limit and in the statistical sample. Formally, our fit to these data produces a distance modulus of , a value that places the galaxy in back of the main Cen A/M83 complex. However, Ciardullo et al. (2002) have noted that for systems more metal-poor than the LMC, the PNLF cutoff fades, leading to systematically larger distances. To correct for this effect, Ciardullo et al. (2002) suggest using the theoretical relation derived by Dopita et al. (1992) for the behavior of with metallicity, i.e.,


where the solar abundance of oxygen is assumed to be (Grevesse et al., 1996). With this relation and the observed low oxygen abundance of the galaxy, (Pilyugin et al., 2004), NGC 5068 moves closer to M83, to  Mpc. This is the first direct distance measurement to this object.

6.4 M94

One of the many galaxies in the Canes Venatici I cloud, M94 (NGC 4736) is surprisingly isolated with no visible companions down to a -band surface brightness limit of 25 mag arcsec (Karachentseva & Karachentsev, 1998) within 3 degrees or 230 kpc. This Sab spiral has a small, but bright bulge, a Low Ionization Nuclear Emission Region (Heckman, 1980), and an extremely bright H ii ring encircling its nucleus at a distance of 1.1 kpc (Muñoz-Tuñón et al., 2004).

The distance to M94 is not well established. Although Karachentsev et al. (2003b) obtained a value of from HST images of the Tip of the Red Giant Branch, and Tonry et al. (2001) derived from Surface Brightness Fluctuations of the galaxy’s bulge, a significantly larger distance (), based simply on the object’s Hubble Flow and  km s Mpc is often used in the literature.

We identified 150 PNe in M94, with 71 in a complete sample which extends 1.3 mag down the PNLF. Our distance modulus from fitting these data is , or  Mpc. This easily fits within the uncertainty of the TRGB measurement but is inconsistent with the more distant SBF value. However, Ciardullo et al. (2002) have shown that the zero point of the SBF distance scale as described by Tonry et al. (2001) is systematically larger than that of the PNLF by  mag. (Half of this discrepancy has since been resolved by Jensen et al. (2003) via a scale adjustment to the SBF calibrators.) If this correction is taken into account, then the SBF distance becomes statistically indistinguishable from our value.

6.5 Ngc 6946

NGC 6946, the Fireworks Galaxy, has a nuclear starburst, significant star formation activity throughout its spiral arms, and an impressive number of recent supernovae (9 in all). It is relatively isolated, having only 10 low surface brightness objects in its group (Karachentsev et al., 2000), and it is located close to the Galactic plane, (). It therefore falls behind a relatively large amount of Galactic extinction (; Schlegel et al., 1998).

There are few reliable distances to NGC 6946. Sharina et al. (1997) used ground-based and images to identify three of the galaxy’s brightest supergiants, and estimated a distance modulus of . Similarly, Karachentsev et al. (2000) analyzed the brightnesses of blue supergiants in the galaxy’s surrounding dwarfs and inferred . Both these values imply a distance that is smaller than that found by Schmidt et al. (1992) by applying the Expanding Photosphere Method to the Type II-L supernova 1980K ().

Because the spiral arms of NGC 6946 are filled with bright H ii regions, identifying PNe in this galaxy is particularly difficult. As a result, our PN sample consists of only 71 candidates, and only 11 of these are brighter than our 90% completeness limit. Moreover, as detailed in §4.1, our H observations did not extend deep enough to confirm the high-excitation nature of these objects, and no follow-up spectroscopy was available in lieu of this limitation. Still, the observed shape of our PNLF is consistent with that expected from the empirical model, and our fit of the top 0.3 mag of this function yields ( Mpc). We note that despite the limited number of PNe in the fit, our distance uncertainty is still dominated by the error associated with the galaxy’s foreground extinction, 0.2 mag. We also note that, although some of our PN candidates may in fact be H ii regions, it is likely that our PNLF distance will remain largely unaffected. Specifically, a careful examination of our images suggests that four of our final 71 PN candidates might be slightly resolved. If we eliminate these objects from our sample, our derived distance to NGC 6946 increases by less than 2%.

7 Variations in the Pnlf

The PNLF distance method has been tested as well as any general purpose extragalactic standard candle. Despite these tests, a few lingering doubts remain. In the HST Key Project’s calibration of the Population II distance scale, Ferrarese et al. (2000) chose not to include the PNLF results in their final distance determination, since its distances were systematically shorter than those found from Surface Brightness Fluctuations. More recently, Sambhus et al. (2006) claimed the detection of two kinematically distinct PN populations in the E6 elliptical galaxy NGC 4697, with each following its own luminosity function. Moreover, surveys of Coma and Virgo have presented evidence for the possible existence of overluminous planetary nebulae in intergalactic space (Arnaboldi et al., 2007, 2008).

The best way to test the robustness of the PNLF method is to examine subsets of PNe within a single galaxy. Population gradients in spiral and elliptical galaxies are common (Zaritsky, et al., 1994; Peletier et al., 1990), so if the PNLF is observed to be constant across the face of a galaxy, that presents strong evidence that neither age nor metallicity is affecting the luminosity function. Unfortunately, to date, only a few PN populations have been studied in this way. Analyses of the PNe across M31 (Merrett et al., 2006), M33 (Magrini et al., 2000; Ciardullo et al., 2004), M81 (Magrini et al., 2001), M104 (Ford et al., 1996), and NGC 5128 (Hui et al., 1993) have found no radial change in the PNLF, while the gradients found in M87 (Ciardullo et al., 1998) and NGC 4526 (Feldmeier et al., 2007) have been attributed to the presence of foreground intracluster stars.

The planetary nebula samples of IC 342, M74, M83, and M94 are large enough for us to examine the systematics of subsamples. We therefore divided the PNe in these galaxies in three different ways: radially (i.e., by defining an inner and an outer sample), by the direction of rotation (assuming that all PNe are rotating along with the galactic disk), and by the orientation of the galaxy (i.e., near side versus far side). The first of these divisions should be sensitive to population changes within the program galaxy, the second to wavelength-dependent variations in the transmission of the narrow-band filter, and the third to host galaxy internal extinction. Figure 2 illustrates these divisions.

Figure 6 compares the cumulative luminosity functions of our various subsamples via a Kolmogorov-Smirnov test. In general, there is excellent consistency within the galaxies. However, as the statistic shows, a few anomalies do exist. In IC 342, the PNLF derived from the sample of planetaries on the approaching side of the galaxy differs from that derived from PNe on the receding side at the 90% confidence level. Interestingly, as the differential plot of Figure 7 illustrates, both sets of objects yield identical values for the PNLF distance. However, while the planetaries on the receding side obey the luminosity function given in equation 5, the approaching PNe deviate from this law with 90% confidence. The PNe of M74 show a similar effect, except this time the gradient lies in the radial direction. Again, the PNLF distances inferred from both samples are nearly identical, but the detailed shapes of the two functions differ with 95% confidence. In this case, the problem can be traced to a deficit of objects 0.3 mag down from . Finally, in M83 and M94, the PNe on the near side of the galaxy follow a slightly different luminosity function than those on the far side (at 90% confidence). In M94, this difference does not affect the determination of , since it is primarily caused by a deficit of objects 0.5 mag below the cutoff in the sample of planetaries north of the galaxy’s major axis. However, in M83, is affected: although both samples of planetaries are well fit by equation (5), the PNLF for objects on the northwest side of the galaxy yields an apparent distance modulus that is 0.18 mag farther than that derived from the southeast sample. This discrepancy is 3 times the formal error of the fits.

Figure 6: The cumulative [O iii] luminosity functions for subsamples of PNe within four galaxies, with the Kolmogorov-Smirnov test statistic marked. In four out of the 12 comparisons, this statistic excludes the null hypothesis with more than 90% confidence. The approaching PNe of IC 342, the inner PNe of M74, and the southwest PNe in M94 all deviate from the empirical law with more than 90% confidence. Although the near-side and far-side samples of PNe in M83 both follow the empirical law, the cutoff magnitude changes by 0.18 mag.
Figure 7: The differential [O iii] luminosity function for the four systems where the Kolmogorov-Smirnov statistic indicates a change in the PNLF across the galaxy. In the case of IC 342, M74, and M94, the most-likely distance modulus is not affected by these fluctuations. In M83 however, the best-fit value for changes by 0.18 mag.

There are two possible explanations for the inconsistencies seen in the figure. The first is internal extinction in the host galaxies. In the Milky Way, the scale height of PNe is only slightly larger than that of the dust (Pandey & Mahra, 1987; Zijlstra & Pottasch, 1991; Phillips, 2001). Such an intermingling can depress the magnitudes of some objects, but not others, and create a luminosity function very similar to that seen in the inner regions of M74 or the receding side of IC 342. Of course, to explain all the variations seen in Figure 6, one needs the extinction to be patchy, or at least have some azimuthal dependence. Moreover, it is well known that the amount of dust in a galaxy declines with radius (Giovanelli et al., 1994; Kylafis et al., 2001), so that if internal extinction is truly responsible for these effects, one would expect to see the PNLF change more with radius than with position angle.

An alternative explanation involves assuming that the PNLF’s shape is population dependent, and that the differences seen in Figure 6 are due to subtle changes in the age or metallicity distribution function of the progenitor stars. To some extent, such changes should not be surprising. We know that the faint end of the PNLF is population dependent: while old systems such as M31’s bulge have exponentially increasing PNLFs, star-forming systems such as the Small Magellanic Cloud exhibit a “dip” in the luminosity function starting 2 mag below the cutoff (Ciardullo et al., 2002; Jacoby & De Marco, 2002). In fact, Figure 4 suggests that the same dip seen in M33 and the SMC may be present in M83 as well. However, there has never before been any evidence for small-scale features that distort at the bright-end of the PNLF.

It is becoming increasingly clear that the term “planetary nebula” refers to a heterogeneous collection of objects. PNe from single stars will have luminosities which depend on the mass of their central star (Vassiliadis & Wood, 1994). Since central star mass depends on progenitor mass through the initial mass-final-mass relation (Weidemann, 2000; Williams, 2007), populations of different ages will produce different PNLFs. Binary stars also produce planetaries via common envelope interactions, and their contribution to the PNLF may be very different from that of single stars (Soker, 2003; Moe & De Marco, 2006). Finally, Ciardullo et al. (2005) have proposed that the descendants of blue stragglers may contribute greatly to the bright-end of the planetary nebula luminosity function. The observed shape of a population’s PNLF may involve a complex mixture of contributions from all these sources. The resiliency of the PNLF technique may simply be due to the processes which place a hard upper limit on the emitted flux at [O iii] .

8 Discussion

The results of our analysis demonstrate that the PNLF is still an excellent standard candle. The distances presented to the six galaxies discussed in this paper are in good agreement with previous measurements, and the subsample comparisons performed in §7 indicate that there are no major population gradients in the PNLF cutoff. However, minor unexplained variations in the luminosity function do exist, and in some cases, this can introduce an additional 0.2 mag error into the distances. At this time, it is impossible to determine whether these variations are intrinsic to the stellar population, or caused by internal extinction in the host galaxies. The fact that these are star-forming galaxies with large amounts of dust suggests that extinction is the cause, but the data from the Sambhus et al. (2006) survey of the elliptical galaxy NGC 4697 and the observation of Virgo’s intracluster stars (Arnaboldi et al., 2008) suggest otherwise. In theory, one can discriminate between these (non-exclusive) hypotheses via spectrophotometry of the objects’ H/H ratio. In practice, however, such a study is extremely difficult. In a typical bright planetary, H is 10 times weaker than [O iii] , so obtaining accurate measurements is a challenge. Moreover, PNe are subject to circumstellar extinction, as well as Galactic extinction. Disentangling the two components will require the analysis of large samples of objects, in many different galactic environments.

We would like to thank the KPNO and CTIO personnel for friendly travel, telescope, and instrumental support and specifically acknowledge Heidi Schweiker for measuring the transmission curves of our filters. We would also like to thank the anonymous referee for useful comments. This research has made extensive use of the USNOFS Image and Catalogue Archive operated by the United States Naval Observatory, Flagstaff Station (, NASA’s Astrophysics Data System, and the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This work was supported by NSF grant AST 06-07416 and a Pennsylvania Space Grant Fellowship. Facilities: Blanco:Mosaic II, WIYN:OPTIC.


  • Arnaboldi et al. (2008) Arnaboldi, M., Doherty, M., Gerhard, O., Ciardullo, R., Aguerri, J.A.L., Feldmeier, J.J., Freeman, K.C., & Jacoby, G.H. 2008, ApJ, 674, L17
  • Arnaboldi et al. (2004) Arnaboldi, M., Gerhard, O., Aguerri, J.A.L., Freeman, K.C., Napolitano, N.R., Okamura, S., & Yasuda, N. 2004, ApJ, 614, L33
  • Arnaboldi et al. (2007) Arnaboldi, M., Gerhard, O., Okamura, S., Kashikawa, N., Yasuda, N., & Freeman, K.C. 2007, PASJ, 59, 419
  • Blair & Long (2004) Blair, W.P., & Long, K.S. 2004, ApJS, 155, 101
  • Brocklehurst (1971) Brocklehurst, M. 1971, MNRAS, 153, 471
  • Buzzoni et al. (2006) Buzzoni, A., Arnaboldi, M., & Corradi, R.L.M. 2006, MNRAS, 368, 877
  • Cardelli et al. (1989) Cardelli, J.A., Clayton, G.C., & Mathis, J.S. 1989, ApJ, 345, 245
  • Ciardullo (2005) Ciardullo, R. 2005, in AIP Conf. 804, Planetary Nebulae as Astronomical Tools, ed. R. Szczerba, G. Stasińska, & S.K. Górny (Melville: American Institute of Physics), 277
  • Ciardullo (2006) Ciardullo, R. 2006, in IAU Symp. 234, Planetary Nebulae in Our Galaxy and Beyond, ed. M.J. Barlow & R.H. Méndez (Cambridge: Cambridge University Press), 325
  • Ciardullo et al. (2004) Ciardullo, R., Durrell, P.R., Laychak, M.B., Herrmann, K.A., Moody, K., Jacoby, G.H., & Feldmeier, J.J. 2004, ApJ, 614, 167
  • Ciardullo et al. (2002) Ciardullo, R., Feldmeier, J.J., Jacoby, G.H., Kuzio de Naray, R., Laychak, M.B., & Durrell, P.R. 2002, ApJ, 577, 31
  • Ciardullo et al. (1998) Ciardullo, R., Jacoby, G.H., Feldmeier, J.J., & Bartlett, R.E. 1998, ApJ, 492, 62
  • Ciardullo et al. (1989) Ciardullo, R., Jacoby, G.H., Ford, H.C., & Neill, J.D. 1989, ApJ, 339, 53
  • Ciardullo et al. (2005) Ciardullo, R., Sigurdsson, S., Feldmeier, J.J., & Jacoby, G.H. 2005, ApJ, 629, 499
  • De Lorenzi et al. (2008) De Lorenzi, F., Gerhard, O., Saglia, R.P., Sambhus, N., Debattista, V.P., Pannella, M., & Mendez, R.H. 2008, ArXiv e-prints, 802, arXiv:0802.1726
  • Dopita et al. (1992) Dopita, M.A., Jacoby, G.H., & Vassiliadis, E. 1992, ApJ, 389, 27
  • Dopita et al. (1997) Dopita, M.A., et al. 1997, ApJ, 474, 188
  • Douglas et al. (2000) Douglas, N.G., Gerssen, J., Kuijken, K., & Merrifield, M.R. 2000, MNRAS, 316, 795
  • Douglas et al. (2007) Douglas, N.G., et al. 2007, ApJ, 664, 257
  • Durrell et al. (2003) Durrell, P.R., Mihos, J.C., Feldmeier, J.J., Jacoby, G.H., & Ciardullo, R. 2003, ApJ, 582, 170
  • Feldmeier et al. (1997) Feldmeier, J.J., Ciardullo, R., & Jacoby, G.H. 1997, ApJ, 479, 231
  • Feldmeier et al. (2004) Feldmeier, J.J., Ciardullo, R., Jacoby, G.H., & Durrell, P.R. 2004, ApJ, 615, 196
  • Feldmeier et al. (2007) Feldmeier, J.J., Jacoby, G.H., & Phillips, M.M. 2007, ApJ, 657, 76
  • Ferrarese et al. (2000) Ferrarese, L., et al. 2000, ApJ, 529, 745
  • Fingerhut et al. (2007) Fingerhut, R.L., Lee, H., McCall, M.L., & Richer, M.G. 2007, ApJ, 655, 814
  • Fleming et al. (1995) Fleming, D.E.B., Harris, W.E., Pritchet, C.J., & Hanes, D.A. 1995, AJ, 109, 1044
  • Ford et al. (1996) Ford, H.C., Hui, X., Ciardullo, R., Jacoby, G.H., & Freeman, K.C. 1996, ApJ, 458, 455
  • Freedman et al. (2001) Freedman, W.L., et al. 2001, ApJ, 553, 47
  • Gallego et al. (2002) Gallego, J., García-Dabó, C.E., Zamorano, J., Aragón-Salamanca, A., & Rego, M. 2002, ApJ, 570, L1
  • Galarza et al. (1999) Galarza, V.C., Walterbos, R.A.M., & Braun, R. 1999, AJ, 118, 2775
  • Gerhard et al. (2007) Gerhard, O., Arnaboldi, M., Freeman, K.C., Okamura, S., Kashikawa, N., & Yasuda, N. 2007, A&A, 468, 815
  • Giovanelli et al. (1994) Giovanelli, R., Haynes, M.P., Salzer, J.J., Wegner, G., da Costa, L.N., & Freudling, W. 1994, AJ, 107, 2036
  • Grevesse et al. (1996) Grevesse, N., Noels, A., & Sauval, A.J. 1996, in ASP Conf. Ser. 99, Cosmic Abundances, ed. S.S. Holt & G. Sonneborn (San Francisco: ASP), 117
  • Gronwall et al. (2007) Gronwall, C., et al. 2007, ApJ, 667, 79
  • Heckman (1980) Heckman, T.M. 1980, A&A, 87, 152
  • Hendry et al. (2005) Hendry, M.A., et al. 2005, MNRAS, 359, 906
  • Hogg et al. (1998) Hogg, D.W., Cohen, J.G., Blandford, R., & Pahre, M.A. 1998, ApJ, 504, 622
  • Howell et al. (2003) Howell, S.B., Everett, M.E., Tonry, J.L., Pickles, A., & Dain, C. 2003, PASP, 115, 1340
  • Hui et al. (1993) Hui, X., Ford, H.C., Ciardullo, R., & Jacoby, G.H. 1993, ApJ, 414, 463
  • Jacoby (1989) Jacoby, G.H. 1989, ApJ, 339, 39
  • Jacoby & Ciardullo (1999) Jacoby, G.H., & Ciardullo, R. 1999, ApJ, 515, 169
  • Jacoby & De Marco (2002) Jacoby, G.H., & De Marco, O. 2002, AJ, 123, 269
  • Jacoby et al. (1987) Jacoby, G.H., Quigley, R.J., & Africano, J.L. 1987, PASP, 99, 672
  • Jacoby et al. (1992) Jacoby, G.H., et al. 1992, PASP, 104, 599
  • Jensen et al. (2003) Jensen, J.B., Tonry, J.L., Barris, B.J., Thompson, R.I., Liu, M.C., Rieke, M.J., Ajhar, E.A., & Blakeslee, J.P. 2003, ApJ, 583, 712
  • Kamphuis & Briggs (1992) Kamphuis, J., & Briggs, F. 1992, A&A, 253, 335 (KB92)
  • Karachentsev et al. (1997) Karachentsev, I., Drozdovsky, I., Kajsin, S., Takalo, L.O., Heinämäki, P., & Valtonen, M. 1997, A&AS, 124, 559
  • Karachentsev et al. (2003a) Karachentsev, I.D., Sharina, M.E., Dolphin, A.E., & Grebel, E.K. 2003, A&A, 408, 111
  • Karachentsev et al. (2000) Karachentsev, I.D., Sharina, M.E., & Huchtmeier, W.K. 2000, A&A, 362, 544
  • Karachentsev & Tikhonov (1993) Karachentsev, I.D., & Tikhonov, N.A. 1993, A&AS, 100, 227
  • Karachentsev et al. (2002) Karachentsev, I.D., et al. 2002, A&A, 385, 21
  • Karachentsev et al. (2003b) Karachentsev, I.D., et al. 2003, A&A, 398, 467
  • Karachentsev et al. (2007) Karachentsev, I.D., et al. 2007, AJ, 133, 504
  • Karachentseva & Karachentsev (1998) Karachentseva, V.E., & Karachentsev, I.D. 1998, A&AS, 127, 409
  • Kniazev et al. (2005) Kniazev, A.Y., Grebel, E.K., Pustilnik, S.A., Pramskij, A.G., & Zucker, D.B. 2005, AJ, 130, 1558
  • Krauss et al. (2005) Krauss, M.I., Kilgard, R.E., Garcia, M.R., Roberts, T.P., & Prestwich, A.H. 2005, ApJ, 630, 228
  • Krismer et al. (1995) Krismer, M., Tully, R.B., & Gioia, I.M. 1995, AJ, 110, 1584
  • Kwok et al. (2003) Kwok, S., Dopita, M., & Sutherland, R. 2003, IAU Symp. 209, Planetary Nebulae: Their Evolution and Role in the Universe (San Francisco: ASP)
  • Kylafis et al. (2001) Kylafis, N.D., Misiriotis, A., Papamastorakis, J., & Xilouris, E.M. 2001, Ap&SS, 276, 531
  • Magrini et al. (2007) Magrini, L., Corbelli, E., & Galli, D. 2007, A&A, 470, 843
  • Magrini et al. (2000) Magrini, L., Corradi, R.L.M., Mampaso, A., & Perinotto, M. 2000, A&A, 355, 713
  • Magrini et al. (2001) Magrini, L., Perinotto, M., Corradi, R.L.M., & Mampaso, A. 2001, A&A, 379, 90
  • Marigo et al. (2004) Marigo, P., Girardi, L., Weiss, A., Groenewegen, M.A.T., & Chiosi, C. 2004, A&A, 423, 995
  • Mazzali et al. (2002) Mazzali, P.A., et al. 2002, ApJ, 572, L61
  • Merrett et al. (2006) Merrett, H.R., et al. 2006, MNRAS, 369, 120
  • Moe & De Marco (2006) Moe, M., & De Marco, O. 2006, ApJ, 650, 916
  • Monet et al. (2003) Monet, D.G., et al. 2003, AJ, 125, 984
  • Muñoz-Tuñón et al. (2004) Muñoz-Tuñón, C., Caon, N., & Aguerri, J.A.L. 2004, AJ, 127, 58
  • Pandey & Mahra (1987) Pandey, A.K., & Mahra, H.S. 1987, MNRAS, 226, 635
  • Peebles (1994) Peebles, P.J.E. 1994, ApJ, 429, 43
  • Peletier et al. (1990) Peletier, R.F., Davies, R.L., Illingworth, G.D., Davis, L.E., & Cawson, M. 1990, AJ, 100, 1091
  • Peña et al. (2007) Peña, M., Stasińska, G., & Richer, M.G. 2007, A&A, 476, 745
  • Peng et al. (2004) Peng, E.W., Ford, H.C., & Freeman, K.C. 2004, ApJ, 602, 685
  • Phillips (2001) Phillips, J.P. 2001, PASP, 113, 839
  • Pilyugin et al. (2004) Pilyugin, L.S., Vílchez, J.M., & Contini, T. 2004, A&A, 425, 849
  • Richer et al. (1998) Richer, M., McCall, M.L., & Stasińska, G. 1998, A&A, 340, 67
  • Roth et al. (2004) Roth, M.M., Becker, T., Kelz, A., & Schmoll, J. 2004, ApJ, 603, 531
  • Saha et al. (2002) Saha, A., Claver, J., & Hoessel, J.G. 2002, AJ, 124, 839
  • Saha et al. (2006) Saha, A., Thim, F., Tammann, G.A., Reindl, B., & Sandage, A. 2006, ApJS, 165, 108
  • Sambhus et al. (2006) Sambhus, N., Gerhard, O., & Méndez, R.H. 2006, AJ, 131, 837
  • Schlegel et al. (1998) Schlegel, D.J., Finkbeiner, D.P., & Davis, M. 1998, ApJ, 500, 525 (SFD98)
  • Schmidt et al. (1992) Schmidt, B.P., Kirshner, R.P., & Eastman, R.G. 1992, ApJ, 395, 366
  • Sharina et al. (1996) Sharina, M.E., Karachentsev, I.D., & Tikhonov, N.A. 1996, A&AS, 119, 499
  • Sharina et al. (1997) Sharina, M.E., Karachentsev, I.D., & Tikhonov, N.A. 1997, AstL, 23, 373a
  • Shaver et al. (1983) Shaver, P.A., McGee, R.X., Newton, L.M., Danks, A.C., & Pottasch, S.R. 1983, MNRAS, 204, 53
  • Sohn & Davidge (1996) Sohn, Y.-J., & Davidge, T.J. 1996, AJ, 111, 2280
  • Soker (2003) Soker, N. 2003, in IAU Symp. 209, Planetary Nebulae: Their Evolution and Role in the Universe, ed. S. Kwok, M. Dopita, & R. Sutherland (San Francisco: ASP), 223
  • Stetson (1987) Stetson, P.B. 1987, PASP, 99, 191
  • Stetson et al. (1990) Stetson, P.B., Davis, L.E., & Crabtree, D.R. 1990, in ASP Conf. Ser. 8, CCDs in Astronomy, ed. G.H. Jacoby (San Francisco: ASP), 289
  • Stone (1977) Stone, R.P.S. 1977, ApJ, 218, 767
  • Stone & Baldwin (1983) Stone, R.P.S., & Baldwin, J.A. 1983, MNRAS, 204, 347
  • Teplitz et al. (2003) Teplitz, H.I., Collins, N.R., Gardner, J.P., Hill, R.S., & Rhodes, J. 2003, ApJ, 589, 704
  • Thim et al. (2003) Thim, F., Tammann, G.A., Saha, A., Dolphin, A., Sandage, A., Tolstoy, E., & Labhardt, L. 2003, ApJ, 590, 256
  • Tonry et al. (2001) Tonry, J.L., Dressler, A., Blakeslee, J.P., Ajhar, E.A., Fletcher, A.B., Luppino, G.A., Metzger, M.R., & Moore, C.B. 2001, ApJ, 546, 681
  • Tonry et al. (2002) Tonry, J.L., Luppino, G.A., Kaiser, N., Burke, B.E., & Jacoby, G.H. 2002, Proc. SPIE, 4836, 206
  • Valdes (1998) Valdes, F.G. 1998, in ASP Conf. Ser. 145, Astronomical Data Analysis Software and Systems VII, ed. R. Albrecht, R.N. Hook, & H.A. Bushhouse (San Francisco: ASP), 53
  • Valtonen et al. (1993) Valtonen, M.J., Byrd, G.G., McCall, M.L., & Innanen, K.A. 1993, AJ, 105, 886
  • Vassiliadis & Wood (1994) Vassiliadis, E., & Wood, P.R. 1994, ApJS, 92, 125
  • Weidemann (2000) Weidemann, V. 2000, A&A, 363, 647
  • Williams (2007) Williams, K.A. 2007, in ASP Conf. Ser. 372, 15th European Workshop on White Dwarfs, ed. R. Napiwotzki & M.R. Burleigh (San Francisco: ASP), 85
  • Zaritsky, et al. (1994) Zaritsky, D., Kennicutt, R.C., Jr., & Huchra, J.P. 1994, ApJ, 420, 87
  • Zheng et al. (1991) Zheng, J.-Q., Valtonen, M.J., & Byrd, G.G. 1991, A&A, 247, 20
  • Zijlstra & Pottasch (1991) Zijlstra, A.A., & Pottasch, S.R. 1991, A&A, 243, 478
Galaxy Type V Size PA
(km s) (deg) (deg)
IC 342 Scd 31 5.58 8.85 0.558
M74 Sc 657 9.76 8.68 0.070
M83 SBc 513 8.79 8.79 0.066
M94 Sab 308 8.75 8.60 0.018
NGC 5068 SBd 668 10.09 8.32 0.102 60
NGC 6946 Scd 48 7.78 8.70 0.342

References. – a: NASA/IPAC Extragalactic Database (NED); b: Pilyugin et al. 2004; c: various sources (see Paper II)

Table 1: Basic Galaxy Properties
[O iii] Filter Exp Time (min) Limiting
Galaxy Name Bandpass Date [O iii] H Seeing
M83 5001/46 2004 May 28-30 180 180 25.9
NGC 5068 5001/46 2004 May 28-30 120 150 25.4
M94 KP1467 4996/45 2005 May 14-16 150 75 25.1
NGC 6946 KP1467 4996/45 2005 May 14-16 150 30 25.9
IC 342 KP1467 4996/45 2005 Nov 5-6 210 30 26.1
M74 KP1590 5020/51 2005 Nov 5-6 360 90 26.5
Table 2: Observing Log
Photometric Uncertainty
Magnitude IC 342 M74 M83 M94 NGC 5068 NGC 6946
24.0 0.040
24.2 0.027 0.043
24.4 0.031 0.051
24.6 0.037 0.064
24.8 0.046 0.076 0.039
25.0 0.055 0.092 0.045
25.2 0.049 0.064 0.108 0.050
25.4 0.056 0.072 0.127 0.057
25.6 0.066 0.045 0.083 0.149 0.063 0.078
25.8 0.080 0.050 0.095 0.175 0.074 0.093
26.0 0.097 0.058 0.114 0.210 0.087 0.114
26.2 0.115 0.066 0.138 0.254 0.104 0.134
26.4 0.135 0.076 0.163 0.310 0.155
26.6 0.155 0.093 0.184 0.180
26.8 0.175 0.110 0.205 0.207
27.0 0.198 0.135 0.229 0.232
27.2 0.225 0.163 0.255 0.269
27.4 0.260 0.192 0.325
27.6 0.306 0.225
Table 3: Photometric Error versus Magnitude
Photometric Filter
Galaxy Zero Point Response Extinction
IC 342 0.015 0.010 0.160
M74 0.028 0.010 0.038
M83 0.015 0.010 0.038
M94 0.040 0.030 0.010
NGC 5068 0.010 0.010 0.056
NGC 6946 0.037 0.010 0.191

References. – a: Schlegel et al. 1998; b: Cardelli et al. 1989

Table 4: Possible Systematic Uncertainties (in mags)
Galaxy (in Mpc) Method Reference
IC 342 Revision of WIYN Cepheids Fingerhut et al. (2007)
IC 342 PNLF This study
IC 342 27.29 2.87 Blue Supergiants corrected for H ii regions Karachentsev et al. (1997)
M94 PNLF This study
M94 SBF Tonry et al. (2001); Jensen et al. (2003)
M94 HST TRGB Karachentsev et al. (2003b)
M83 HST TRGB of group Karachentsev et al. (2007)
M83 PNLF This study
M83 revised VLT Cepheids Saha et al. (2006)
M83 HST TRGB Karachentsev et al. (2007)
NGC 5068 Metallicity corrected PNLF This study
NGC 6946 Brightest Blue Supergiants of group Karachentsev et al. (2000)
NGC 6946 Brightest Supergiants Sharina et al. (1997)
NGC 6946 PNLF This study
NGC 6946 Expanding Photospheres Method Schmidt et al. (1992)
M74 PNLF This study
M74 Brightest Blue Supergiants Sharina et al. (1996)
M74 Brightest Supergiants Hendry et al. (2005)
M74 SN II-P expansion rate Hendry et al. (2005)

Note. – Galaxies are listed by increasing distance and distance moduli are listed by increasing uncertainty.

Table 5: Comparison of Distance Moduli
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description