THE ALFALFA "ALMOST DARKS" CAMPAIGN: PILOT VLA HI OBSERVATIONS OF FIVE HIGH MASS-TO-LIGHT RATIO SYSTEMS

, , , , , , , , , , , , , , , and

Published 2015 January 22 © 2015. The American Astronomical Society. All rights reserved.
, , Citation John M. Cannon et al 2015 AJ 149 72 DOI 10.1088/0004-6256/149/2/72

1538-3881/149/2/72

ABSTRACT

We present new Very Large Array (VLA) H i spectral line imaging of five sources discovered by the ALFALFA extragalactic survey. These targets are drawn from a larger sample of systems that were not uniquely identified with optical counterparts during ALFALFA processing, and as such have unusually high H i mass to light ratios. The candidate "Almost Dark" objects fall into four broad categories: (1) objects with nearby H i neighbors that are likely of tidal origin; (2) objects that appear to be part of a system of multiple H i sources, but which may not be tidal in origin; (3) objects isolated from nearby ALFALFA H i detections, but located near a gas-poor early type galaxy; (4) apparently isolated sources, with no object of coincident redshift within ∼400 kpc. Roughly 75% of the 200 objects without identified counterparts in the α.40 database (Haynes et al. 2011) fall into category 1 (likely tidal), and were not considered for synthesis follow-up observations. The pilot sample presented here (AGC193953, AGC208602, AGC208399, AGC226178, and AGC233638) contains the first five sources observed as part of a larger effort to characterize H i sources with no readily identifiable optical counterpart at single dish resolution (3farcm5). These objects span a range of H i mass [7.41 < log(MHi ) < 9.51] and H i mass to B-band luminosity ratios (3 < MHi /LB < 9). We compare the H i total intensity and velocity fields to optical imaging drawn from the Sloan Digital Sky Survey and to ultraviolet imaging drawn from archival GALEX observations. Four of the sources with uncertain or no optical counterpart in the ALFALFA data are identified with low surface brightness optical counterparts in Sloan Digital Sky Survey imaging when compared with VLA H i intensity maps, and appear to be galaxies with clear signs of ordered rotation in the H i velocity fields. Three of these are detected in far-ultraviolet GALEX images, a likely indication of star formation within the last few hundred Myrs. One source (AGC208602) is likely tidal in nature, associated with the NGC 3370 group. Consistent with previous efforts, we find no "dark galaxies" in this limited sample. However, the present observations do reveal complex sources with suppressed star formation, highlighting both the observational difficulties and the necessity of synthesis follow-up observations to understand these extreme objects.

Export citation and abstract BibTeX RIS

1. INTRODUCTION

ALFALFA, the Arecibo Legacy Fast ALFA extragalactic H i survey, has obtained a census of H i-bearing galaxies over 7000 square degrees of high Galactic latitude sky visible from Arecibo between $0{}^\circ \;\lt $ δ $\lt +36{}^\circ $ (Giovanelli et al. 2005; Haynes et al. 2011). Its combination of sensitivity, bandwidth, and sky coverage make ALFALFA the first wide area H i survey to sample a cosmologically significant volume, delivering robust determinations of the H i mass, width, and diameter functions, the H i correlation function, and the cosmic density of H i at $z\sim 0$ (Dowell 2010; Martin et al. 2010; Papastergis et al. 2011; Martin et al. 2012). When completed, the final ALFALFA database will contain more than 30,000 high signal-to-noise ratio (S/N) H i detections.

As discussed in Haynes et al. (2011), the 40% ALFALFA catalog (hereafter "α.40") contains 15,855 extragalactic H i sources. Less than 1.5% of these sources were not clearly identified with an associated stellar population (hereafter referred to as an "optical counterpart") in an optical survey (SDSS or DSS). Of these, only ∼50 were not readily associated with a known structure or nearby H i source (that is, roughly 75% of these "dark" candidates are likely tidal in origin). The nature of the remaining objects remains unclear; their optical associations are often ambiguous given the ALFALFA pointing accuracy (∼26'' for reasonable S/N sources) and unresolved given the ALFALFA beam (3farcm5). As discussed in detail below, there are only a few viable candidate "dark" galaxies in the current literature; the rich population of sources measured by ALFALFA is a promising place to look for more.

To probe the nature of these candidate "dark" objects, we have initiated the ALFALFA "Almost Darks" program: a comprehensive, multi-wavelength follow-up observing campaign to characterize the nature of ALFALFA H i detections that were not identified with an unambiguous optical counterpart in initial ALFALFA processing. We anticipate that detailed studies of these objects will result in the confirmation of associated optical counterparts for the vast majority of systems. Figure 1 shows these objects in comparison to the full α.40 catalog; the 11,300 ALFALFA sources with good SDSS photometry (see Haynes et al. 2011) are plotted as either individual black points, or as colored contours at 10% levels (between 10% and 90%) of the number density of sources. Optical luminosities for these sources are calculated using SDSS catalog magnitudes from the ALFALFA-SDSS cross match presented in Haynes et al. (2011). The 52 candidate "dark" sources are shown in red, with the five objects selected for the pilot sample studied in detail in the present manuscript shown in green. These 52 sources are those that have no clear optical counterpart assigned in the initial ALFALFA processing, and that were not deemed as likely to be tidal in nature. Red points represent upper limits in luminosity estimated at the limiting magnitude for a source 10'' in diameter in downloaded SDSS images. Green points that are not upper limits are calculated using the SDSS catalog magnitudes of the optical counterparts identified after Very Large Array (VLA) observations; prior to source localization with the VLA these objects were upper limits similar to the points plotted in red.

Figure 1.

Figure 1. Plot of the logarithm of the neutral hydrogen to optical luminosity ratio, M $_{{\rm HI}}$/L $_{B}$, versus the logarithm of the H i mass, for all members of the α.40 catalog (Haynes et al. 2011). All black points and areas within contours represent galaxies with reliable SDSS optical counterparts. The contours show number densities in intervals of 10%, from 10% to 90% of the number density maximum (which occurs at log(M $_{{\rm HI}}$) = 9.8; see Haynes et al. 2011). The candidate "Almost Dark" sources are shown in red or green (studied in this work). These sources have high M $_{H}$/L but lack obvious or uncertain optical counterparts. The green point at log(M $_{{\rm HI}}$) = 9.51 is AGC233638, a massive H i source with a discrepant SDSS redshift whose association was uncertain based on ALFALFA alone; see detailed discussion in Section 3.3.2. For consistency, all M $_{{\rm HI}}$/L values calculated here use SDSS catalog magnitudes.

Standard image High-resolution image

While Figure 1 was not used for the selection of candidate "Almost Dark" sources, it does clearly demonstrate that these systems harbor high H i mass to stellar light ratios (hereafter, MHi /L) compared to the vast majority of α.40 detections. Note that some sources with elevated MHi /L ratios are not classified as candidate "Almost Dark" systems: these sources have clearly associated optical counterparts. Further, while no strict minimum value of M $_{{\rm HI}}$/L was enforced to create the sample of "Almost Dark" systems, an empirical threshold (which does have significant uncertainty at low luminosity and low surface brightness) appears around M $_{{\rm HI}}$/L $\simeq $ 10.

In this manuscript we present pilot VLA H i spectral line observations of five candidate "Almost Dark" sources. This represents the first step in a major program to understand the nature and origin of the most extreme M $_{{\rm HI}}$/L systems within the context of various scenarios for their existence: as tidal features, as galaxies with extremely low surface brightness stellar populations, or as genuine "dark" galaxies. The ALFALFA "Almost Dark" objects fall into four broad categories: (1) objects with nearby H i neighbors that are likely of tidal origin; (2) objects that appear to be part of a system of multiple H i sources, but which may not be tidal in origin; (3) objects isolated from nearby ALFALFA H i detections, but located near a gas-poor early type galaxy; (4) apparently isolated sources, with no object of coincident redshift within ∼400 kpc. The five sources presented here (AGC193953, AGC208602, AGC208399, AGC226178, and AGC233638) populate each of the last three categories (sources with a clear tidal origin from ALFALFA data products were not targeted for interferometric follow-up). Table 1 summarizes the basic properties of the five sources, including position, distance, and H i properties (derived from ALFALFA measurements). This limited sample is not meant to be representative of the "Almost Dark" galaxy sample as a whole; these specific sources were selected in part to fit into pre-approved LST observing slots for VLA scheduling (see Section 2.1).

Table 1.  Basic and Global H i Properties of the Candidate "Almost Dark" Systems

AGC R.A.a Decl.a Offsetb S $_{{\rm HI}}$c V $_{{\rm sys}}$ W $_{50}$ Dd log(M $_{{\rm HI}}$)
Number (J2000) (J2000) (arcminutes) (Jy km s−1) (km s−1) (km s−1) (Mpc) (M $_{\odot }$)
193953e 143.6444 12.394 1.32 0.43 2591 32 39.9 8.21
208399e 160.0417 4.910 0.36 1.06 747 31 10.1 7.41
208602f 163.6558 17.635 0.92 1.45 1093 50 18.5 8.07
226178e,g 191.6782 10.369 1.33 0.62 1581 28 16.5 7.60
233638e 202.8267 16.167 0.05 1.10 7587 62 112 9.51

aR.A. and Decl. of the H i centroid measured in the VLA images. bAngular offsets between ALFALFA centroid coordinates and VLA centroid coordinates. cTotal integrated H i line flux density measured by ALFALFA in Jy km s−1. dDistances based on the local flow model by Masters (2005) and Martin et al. (2010). eVery low surface brightness optical counterpart. fPositional offset from early type galaxy NGC 3457. gMultiple possible optical counterparts.

Download table as:  ASCIITypeset image

There has been widespread interest in identifying truly "dark" galaxies, if they exist. With the advent of sensitive, wide-area H i surveys such as ALFALFA, and with the increased availability of mapping observations of nearby groups with single dish telescopes, a number of candidate "dark galaxies" have appeared in the literature. These are systems whose H i properties are representative of galaxies, but that lack a clear optical counterpart. The VIRGO HI 21 cloud was originally postulated to be a true dark galaxy (Minchin et al. 2005, 2007). However, subsequent observations (Haynes et al. 2007) and modeling (Duc & Bournaud 2008) have shown that this source is likely tidal in nature. More recently, single-dish telescopes have been used in mapping mode to produce deep targeted H i images of nearby groups. In such a survey of the M101 group, Mihos et al. (2012) discovered a gas cloud labeled as GBT 1355 + 5439. In that same work, very deep ground-based optical imaging revealed no obvious stellar counterpart to a limiting surface brightness of ${{\mu }_{V}}$ = 29 mag arcsec−2. Subsequent interferometric H i observations by Oosterloo et al. (2013) confirm the extreme characteristics of GBT 1355 + 5439; both a tidal interpretation and an interpretation as a truly "dark" object remain viable. Numerous other sources are interesting and relevant in this context (see discussion and references in Walter et al. 2005; Kent et al. 2009; English et al. 2010; Kent 2010; Matsuoka et al. 2012; Serra et al. 2013; Adams et al. 2015; Lee-Waddell et al. 2014).

We organize this manuscript as follows. In Section 2 we discuss the details of the H i data reductions and image production, as well as the techniques applied to derive optical and ultraviolet (UV) properties from SDSS and GALEX images. Section 3 presents the H i , optical, and UV images of each of the five sources; using these data, we can localize the H i gas with confidence, assign possible optical counterparts (by comparison with optical imaging), and search for ordered motions that are indicative of galaxy-like rotation. In Section 4 we further discuss the five "Almost Dark" sources in the context of the ALFALFA database, the low-mass galaxies within it, and also among the (growing) list of "Almost Dark" candidates. We draw our conclusions in Section 5; primary among these is that none of the five candidate objects appears to qualify as a bona fide "dark galaxy."

2. OBSERVATIONS AND DATA HANDLING

2.1. H i Observations

H i spectral line imaging of the five sources presented in this work was obtained with the National Radio Astronomy Observatoryʼs Karl G. Jansky VLA12 in 2013 April, for program 13A-028 (legacy ID AC1116; P.I. Cannon). The five sources selected for observing were the best fits based on the LST oversubscription of the VLA during the 13A semester. The observations were acquired during four separate 2 hr observation blocks and during one 2.5 hr block (for source AGC193953) with the observatory in the "D" (most compact) configuration. The WIDAR correlator was used to divide an 8 MHz bandwidth into 1024 channels. The resulting channel separation of 7.812 kHz ch−1 delivers a native spectral resolution of 1.65 km s−1 ch−1 at the rest frequency of the H i spectral line. Reductions followed standard prescriptions and used both the Astronomical Image Processing System (AIPS) and the Common Astronomy Software Applications packages.

Radio frequency interference (RFI) and bad data were removed manually from each data set. The RFI contamination of the AGC193953 data set was significant in a few channels close to those containing the source; the loss of sensitivity in these channels has affected our recovery of flux from this source (see further discussion below). Bandpass and flux calibrations used the standard calibrator 3C286 for all five observing blocks. For each source, a nearby quasar was observed every ∼20 minutes and used to calibrate the phases. We estimate our absolute flux calibration to be accurate at the ∼10% level.

The imaging of the calibrated uv databases followed standard prescriptions. To increase S/N, each database was spectrally averaged by a factor of two, to produce a velocity resolution of 15.625 kHz ch−1 (3.3 km s−1 ch−1) prior to imaging. Each cube was inverted with minimal cleaning, and source(s) were identified. Using these full cubes, regions were identified that excluded H i line sources and bandpass edge effects. Using these regions (typically 50 channels or more, depending on the complexity of emission within each cube) as fitting points, the continuum was then subtracted in the uv plane using a first-order polynomial function.

The final spectral line cubes used in this analysis were produced using the IMAGR task in AIPS, with a weighting ROBUST factor of 0.5. Using clean boxes that tightly encompass the high S/N line emission, each cube was deeply cleaned to one half of the rms level found in line-free channels. Residual flux rescaling was enforced during the imaging process (Jorsater & van Moorsel 1995). The moment zero maps shown below were produced following the THINGS procedures as described in Walter et al. (2008). The original-resolution cubes were convolved to the nearest integer circular sizes in arcseconds (so, for example, the 48farcs 97 × 43farcs 00 cube for AGC208602 was convolved to a circular beam size of 49'') and were then used for all subsequent analyses. A second convolution was enacted to produce cubes with circular beam sizes ∼10% larger; these smoothed cubes were blanked at the 2.5σ level, closely examined by eye to retain only emission that is coherent in velocity space, and then used as masks against the original (circular beam size) cubes. Moment zero maps (representing total H i emission) were produced by collapsing these masked cubes. The resulting images of each source are discussed in detail below. The total H i flux integral for each of the five sources is less than the flux integral derived from the ALFALFA observations; we discuss the recovered flux from each source in the sections that follow.

The velocity structure of the H i gas in each source plays an important role in the interpretation of these data; the presence of ordered gas motion is very strong evidence for classifying a candidate "Almost Dark" source as a bona fide galaxy. As discussed in detail for each individual system below, whenever possible we have used Gaussian fitting to the full (unblanked) datacubes in order to extract the velocity fields of the "Almost Dark" systems presented below. Specifically, we use the GIPSY task XGAUFIT to fit Gaussian functions to emission above the 3σ level in regions of each datacube immediately surrounding four of the five "Almost Dark" sources; these Gaussian-fitted velocity fields are shown in the four-panel images presented in Section 3. For AGC208602, the incoherent velocity structure of the source precluded meaningful Gaussian fitting. For those "Almost Dark" sources that have neighboring H i-rich galaxies that are detected in the datacubes (AGC208602, AGC226178, and AGC233638), we also show full-field images including all detected galaxies using the more traditional moment one approach (i.e., collapsing of the blanked datacube) for creating the velocity fields.

2.2. Optical Observations

To quantify the optical properties of these five "Almost Dark" sources, we extract photometry from the Sloan Digital Sky Survey 10th data release products (Ahn et al. 2014). Mosaic frames in the SDSS-g and the SDSS-r bands, using primary covering fields, were created using the SDSS online interface. Integrated magnitudes were derived for each of the "Almost Dark" sources using the LAMBDAR source extraction code (A. Wright et al. 2015, in preparation13 ). Elliptical apertures, centered on the H i centroid positions (see Table 1) and with semimajor and semiminor axis lengths as shown in Table 2, were used to extract the SDSS-g magnitudes and (g−r) colors of each source. The aperture sizes and orientations were set by hand in order to enclose all of the optical flux from each source when the SDSS images were displayed at the highest contrast levels.

Table 2.  Observed Properties of the Candidate "Almost Dark" Systems

AGC Aperturea m $_{{\rm FUV}}$b m $_{{\rm NUV}}$c m $_{g,{\rm SDSS}}$d m $_{g}$e m $_{r}$f
Number Size (mag) (mag) (mag) (mag) (mag)
193953 16'' × 11'' N/A 20.0 ± 0.3 20.8 ± 0.1 20.2 ± 0.3 20.1 ± 0.3
208399 26'' × 18'' 21.6 ± 0.3 19.9 ± 0.3 20.00 ± 0.03 18.6 ± 0.2 18.3 ± 0.2
208602 15'' × 15'' >21.4 20.3 ± 0.3 >20.5 >21.0
226178 42'' × 8'' 19.6 ± 0.2 20.2 ± 0.3 22.4 ± 0.1 19.4 ± 0.3 19.2 ± 0.3
233638 25'' × 18'' 19.9 ± 0.2 19.3 ± 0.2 18.72 ± 0.01 18.0 ± 0.2 17.8 ± 0.2

aSemi-major and semi-minor axis lengths of photometric aperture. bApparent integrated magnitude of the source in the GALEX FUV band. No Galactic absorption correction has been applied. cApparent integrated magnitude of the source in the GALEX NUV band. No Galactic absorption correction has been applied. dSDSS catalog model magnitude. eApparent integrated magnitude of the source in the SDSS g band. No Galactic absorption correction has been applied. fApparent integrated magnitude of the source in the SDSS r band. No Galactic absorption correction has been applied.

Download table as:  ASCIITypeset image

The photometry in Table 2 is used to derive the physical parameters of the five "Almost Dark" sources shown in Table 3. After correcting for foreground Galactic absorption (Schlafly & Finkbeiner 2011), the SDSS g magnitude is converted to a Johnson-Cousins B magnitude via the prescriptions in Table 2 of Blanton & Roweis (2007). These (absorption corrected) apparent magnitudes are converted to absolute magnitudes using our assumed distances from Table 1, and then converted to solar luminosities using the values from Table 1 of (Blanton & Roweis 2007). The resulting g and B-band luminosities of each "Almost Dark" source are tabulated and compared with the total H i masses (using the ALFALFA flux integrals) in Table 3. This table reveals that the M $_{{\rm HI}}$/L ratios span more than a factor of 10 across this small sample. It is important to stress that all of our measurements of the optical properties of these sources are brighter than the pipeline SDSS values; this is a result of the more careful aperture placement and background subtraction done here than by the standard SDSS pipeline measuring tools.

Table 3.  Derived Properties of the Candidate "Almost Dark" Systems

AGC L $_{g}$a L $_{B}$b M $_{\star }$c M $_{{\rm HI}}$/L $_{g}$d M $_{{\rm HI}}$/L $_{B}$e M $_{{\rm HI}}$/M $_{\star }$f SFR$_{{\rm FUV}}$g
Number (L $_{\odot }$) (L $_{\odot }$) (M $_{\odot }$) (M $_{\odot }$/L $_{\odot }$) (M $_{\odot }$/L $_{\odot }$) (M $_{\odot }$/M $_{\odot }$) (M $_{\odot }$ yr−1)
193953 1.6 × 107 1.9 × 107 2.2 × 106 10 9 75 ...
208399 4.4 ×106 5.0 × 106 1.7 × 106 6 5 15 5.1 × 10−5
208602 <2.6 × 106 <2.8 × 106 <2.5 × 106 >48 >45 >49 <2.6 × 10−4
226178 5.7 × 106 6.7 × 106 1.3 × 106 7 6 30 0.0017
233638 9.6 × 108 1.1 × 109 2.2 × 108 3 3 14 0.047

ag-band luminosity, corrected for Galactic extinction via Schlafly & Finkbeiner (2011), derived from g band magnitudes extracted from downloaded SDSS images using LAMBDAR. bB-band luminosity, corrected for Galactic extinction via Schlafly & Finkbeiner (2011), derived from g and r band magnitudes extracted from downloaded SDSS images using LAMBDAR. cStellar mass, derived using the LAMBDAR r-band magnitudes, g-r colors, and the formalism in Zibetti et al. (2009). dRatio of H i mass in M $_{\odot }$ to SDSS g-band luminosity in L $_{\odot }$. eRatio of H i mass in M $_{\odot }$ to B-band luminosity in L $_{\odot }$. fRatio of H i mass in M $_{\odot }$ to stellar mass in M $_{\odot }$. gFar-UV star formation rate.

Download table as:  ASCIITypeset image

We use the measured r-band magnitude, (g−r) color, and the relations from Zibetti et al. (2009) to derive the stellar mass and M $_{{\rm HI}}$/M $_{\star }$ ratio of each candidate "Almost Dark" source shown in Table 3. We stress that these representative stellar masses are uncertain at the ∼50% level due to systematic uncertainties alone. Nonetheless, these values can be compared to those for the full α.40 catalog using the same technique; we explore this comparison later in this manuscript.

A detailed investigation of the optical properties of the ALFALFA "Almost Darks" sample is currently underway with the One Degree Imager (ODI14 ) on the WIYN 3.5 m telescope.15 The results of this endeavor will be presented in a forthcoming manuscript (Janowiecki et al. 2015). As of this writing, preliminary images of two members of the present sample (AGC193953 and AGC233638) have been obtained as part of this effort. In the interests of uniformity, we only derive quantitative results for the members of the present pilot sample from the SDSS images; the preliminary derivations of SDSS g-band magnitudes for AGC193953 and AGC233638 from the partially filled focal plane ODI images agree within errors with the values derived from the original SDSS images.

2.3. Ultraviolet Observations

The UV properties of the five "Almost Dark" sources were derived using archival GALEX 16 imaging. The depth of the GALEX images is heterogeneous across the sample; two sources (AGC193953 and AGC233638) have only All-Sky Imaging Survey depth imaging (∼100 seconds per integration), and one of these (AGC193953) is a non-detection. The other three sources have deeper imaging (∼1600 seconds per integration); one source (AGC208602) is a non-detection at this sensitivity level. Table 2 shows the magnitudes in the GALEX near and far-UV channels, measured in identical apertures to those used to extract the optical photometry shown in Table 2. We use the far-UV fluxes to derive the recent (∼100 Myr) star formation rates (SFR$_{{\rm FUV}}$) of the sources using the prescriptions in Salim et al. (2007); these quantities are discussed in detail below.

3. THE GASEOUS AND STELLAR COMPONENTS OF THE "ALMOST DARK" SOURCES

In this section we present the H i optical and UV images of each of the "Almost Dark" sources. We group these five systems into three categories: isolated systems (no known sources with coincident redshift within a projected distance of 400 kpc), systems with an apparent tidal origin, or systems near other objects.

3.1. Isolated Sources

We classify two of the five "Almost Dark" sources as isolated; no known sources with coincident redshift are located within a projected distance of 400 kpc of either source.

3.1.1. AGC193953

AGC193953 has the weakest ALFALFA H i line flux in this small sample (S $_{{\rm HI}}$ = 0.43 Jy km s−1; see Table 1). Located at ∼40 Mpc, the H i line width (W $_{50}$ = 32 km s−1, where W $_{50}$ is the width of the H i line at 50% of maximum) and H i reservoir of 1.6 × 108 M $_{\odot }$ are indicative of an average dwarf galaxy. This source was initially cataloged as a "dark" candidate due to the extremely low surface brightness nature of the optical counterpart, which is significantly displaced from the center of the ALFALFA beam. As discussed above, RFI compromised a number of channels in the datacube close to those that contain emission from the source. As such, the total H i flux in the VLA data (0.13 Jy km s−1) only amounts to ∼30% of the ALFALFA flux; the H i properties (e.g., the more narrow range of velocities detected in the VLA data compared to the ALFALFA data) should only be considered to be representative based on these data.

Figure 2 shows the H i, SDSS-r band, and GALEX near-UV images of AGC193953. At D configuration angular resolution the source is unresolved. The H i gas has a single surface density concentration that reaches a maximum column density of only 6.5 × 1019 cm−2. This surface density maximum is slightly offset spatially (though within one beam size) from the very low surface brightness optical counterpart, which is clearly visible in the SDSS-r band image. The Gaussian-fitted velocity field shows coherent and ordered rotation within the beam. As panels (b) and (d) of Figure 2 show, the isovelocity contours are parallel over ∼7 km s−1(∼2 velocity resolution elements).

Figure 2.

Figure 2. H i, SDSS R-band, and GALEX images of AGC193953. Panel (a) shows H i column density contours at the (0.4, 0.5, 0.6) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (b) shows H i isovelocity contours from 2596 to 2603 km s−1 in intervals of 1 km s−1, overlaid on a color representation of the velocity field image. Panel (c) shows the GALEX near-UV image, while panel (d) shows the SDSS r-band image (displayed in a logarithmic stretch to highlight low surface brightness emission); the contours in (c) are the same as in panel (a), while the contours in (d) are the same as in panel (b). The velocity field was derived by fitting Gaussian functions to the full data cube. The H i beam size is 50''; the source is unresolved at this angular resolution.

Standard image High-resolution image

The optical counterpart of AGC193953 is of very low surface brightness. As Tables 2 and 3 show, the SDSS integrated magnitudes imply a total B-band luminosity of ∼1.9 × 107 L $_{\odot }$ and a total stellar mass of 2.2 × 106 M $_{\odot }$. This faint stellar component is detected in the shallow GALEX near-UV image; while we have no direct constraints on the nature of recent star formation (which requires far-UV imaging), the faint nature of the source in the near-UV argues for quiescence over the most recent ∼100 Myr. Using the ALFALFA total H i flux and the SDSS-g and derived B-band luminosities, we derive M $_{{\rm HI}}$/L $_{g}$ = 10, M $_{{\rm HI}}$/L $_{B}$ = 9, and M $_{{\rm HI}}$/M $_{\star }$ = 75. While our analysis of the SDSS images have recovered more flux than the SDSS pipeline products, the source is close to the detection limits of the SDSS imaging; it is likely that deeper optical imaging would recover more low surface brightness emission, thus decreasing the M $_{{\rm HI}}$/L and ${{M}_{{\rm HI}}}$/M $_{\star }$ ratios.

3.1.2. AGC208399

AGC208399 is a very low surface brightness source that was included in the Huang et al. (2012) study of gas, stars, and star formation in ALFALFA dwarfs. AGC208399 has the lowest recessional velocity of the five "Almost Dark" sources studied here; its flow model distance (see Table 1) is 10.1 Mpc. The flux integral (S $_{{\rm HI}}$ = 1.06 Jy km s−1) implies an H i mass of 2.6 × 107 M $_{\odot }$. By comparison with the properties in the Huang et al. (2012) sample, AGC208399 is an extreme low surface brightness dwarf galaxy.

The H i, optical, and UV images of AGC208399 shown in Figure 3 reveal a centrally concentrated H i distribution that is spatially coincident with a diffuse, very low surface brightness optical counterpart in the SDSS r-band image. The VLA images recover ∼85% of the ALFALFA flux integral. The H i column density peaks at 2.1 × 1020 cm−2 at 52'' angular resolution, which marginally resolves the source. There is bulk rotation of the H i gas, extending across nearly 20 km s−1. Most of this projected rotation is coherent, as evidenced by the nearly parallel isovelocity contours over most of the disk. However, there are some departures from circular motion in the outer regions, especially on the approaching side of the disk.

Figure 3.

Figure 3. H i, SDSS R-band, and GALEX images of AGC208399. Panel (a) shows H i column density contours at the (0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (b) shows H i isovelocity contours from 760 to 776 km s−1 in intervals of 2 km s−1, overlaid on a color representation of the velocity field image. Panel (c) shows a color representation of the GALEX UV images, while panel (d) shows the SDSS r-band image (displayed in a logarithmic stretch to highlight low surface brightness emission); the contours in (c) are the same as in panel (a), while the contours in (d) are the same as in panel (b). The velocity field was derived by fitting Gaussian functions to the full data cube. The H i beam size is 52''; the source is marginally resolved at this angular resolution. The two galaxies in panel (d) located close to AGC208399 (to the W and SW) are background galaxies with SDSS redshifts of 0.065 and 0.071, respectively.

Standard image High-resolution image

The GALEX imaging shown in Figure 3(c) shows a very low surface brightness blue component in close angular proximity to the H i surface density maximum. It is difficult to discern if this faint emission arises from the stellar population associated with AGC208399 (and thus likely indicates recent star formation within the past ∼100 Myr) or with unrelated background sources; indeed, the SDSS-r image shown in panel (d) of Figure 3 reveals some small point sources that could be background star-forming galaxies. With this caution in mind, the UV fluxes are extracted using the same apertures as for the SDSS photometry (see Table 2). If the far-UV flux is associated with AGC208399, then the implied SFR$_{{\rm FUV}}$ = 5 × 10−5 M $_{\odot }$ yr−1.

AGC208399 shares many physical characteristics with the most massive dwarf galaxies in SHIELD (Cannon et al. 2011), a complementary ALFALFA follow-up observing campaign to study the physical properties of extremely low-mass systems cataloged by ALFALFA. Specifically, the SDSS colors, H i line widths, total H i masses, and distances of the most massive SHIELD galaxies are similar to those of AGC208399. Interestingly, however, the majority of the SHIELD galaxies show evidence for recent star formation. The gas-rich (M $_{{\rm HI}}$/L $_{g}$ = 6; M $_{{\rm HI}}$/L $_{B}$ = 5; M $_{{\rm HI}}$/M $_{\star }$ = 15) but quiescent (SFR$_{{\rm FUV}}$ = 5 × 10−5 M $_{\odot }$ yr−1) nature of AGC208399 warrants further detailed exploration.

3.2. Systems With A Tidal Origin

3.2.1. AGC208602

AGC208602 was classified as a source with H i in the vicinity of, but offset from, an early type galaxy (NGC 3457). A comparison of the ALFALFA detections near AGC208602 with SDSS imaging revealed that this candidate "Almost Dark" source is located within 30' of at least four known H i sources, including NGC 3443, NGC 3454, NGC 3455, and UGC 6022. This loose collection is part of the NGC 3370 group of galaxies, which is in turn part of the larger Leo II structure, which contains multiple gas-rich spirals and irregular galaxies. It also contains at least one early type galaxy; NGC 3457 is in close angular proximity to AGC208602, and its optical velocity (V $_{{\rm sys}}$= 1148 km s−1; Cappellari et al. 2011) is in fair agreement with those of other members of the group (see below).

Figure 4 shows VLA H i and SDSS optical images of the region surrounding the "Almost Dark" candidate, AGC208602. The size of the field of view (FOV) (∼36') is roughly equal to that of the primary beam of the VLA at H i. Within this field, five other H i line sources are detected in the velocity range immediately surrounding AGC208602 (specifically, between ∼985 and ∼1210 km s−1); these five sources are noted in each panel of Figure 4. In addition to NGC 3443, NGC 3454, NGC 3455, and UGC 6022, we also detect the neutral hydrogen in a small system labeled as SDSS J105506.66 + 172745.4. Each of these five sources shows coherent signatures of rotation.

Figure 4.

Figure 4. H i and SDSS R-band images of the AGC208602 field. Panel (a) shows H i column density contours overlaid at the 3.5 × 1019 cm−2 level (red, highlighting the low surface brightness nature of AGC208602) and at the (2.5, 5, 10) × 1020 cm−2 levels (black, highlighting nearby sources in the field), overlaid on a grayscale representation of the image. Panel (c) shows H i isovelocity contours from 985 to 1210 km s−1 in intervals of 35 km s−1, overlaid on a color representation of the velocity field image; the color bar is set to encompass the full velocity range of H i gas in the six detected sources. Panels (b) and (d) show the SDSS r-band image of the AGC208602 field, displayed in a logarithmic stretch to highlight low surface brightness emission; the contours in (b) are the same as in panel (a), while the contours in (d) are the same as in panel (c). The six gas-bearing systems in the field are labeled, as is the early type galaxy NGC 3457. The H i beam size is 49'', as shown in the bottom left of panel (a).

Standard image High-resolution image

Figure 5 shows a closer view of the detected H i from AGC208602, and its relationship to NGC 3457. The neutral hydrogen in this source is of very low surface brightness, characterized by flocculent patches of emission that are at times disjointed in velocity space. While the VLA images recover ∼65% of the ALFALFA flux, the integrated column densities at 49'' resolution only reach the ∼5.5 × 1019 cm−2 level. The velocity structure of this source is sufficiently incoherent that attempts to fit Gaussians to the emission in the cube failed. Thus, the moment one maps shown in Figures 4 and 5 were derived by standard routines. Emission is detected over ∼50 km s−1; there is no convincing evidence of coherent rotation of this source in the velocity field or in the three-dimensional cube.

Figure 5.

Figure 5. H i, SDSS R-band, and GALEX images of AGC208602. Panel (a) shows H i column density contours at the (0.35, 0.45, 0.55) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (b) shows H i isovelocity contours from 1085 to 1130 km s−1 in intervals of 5 km s−1, overlaid on a color representation of the velocity field image. Panel (c) shows a color representation of the GALEX UV images, while panel (d) shows the SDSS r-band image (displayed in a logarithmic stretch to highlight low surface brightness emission); the contours in (c) are the same as in panel (a), while the contours in (d) are the same as in panel (b). The early type galaxy NGC 3457 is labeled in panel (c); the systemic velocity of this galaxy (V = 1148 km s−1; Cappellari et al. 2011) suggests a physical association with the other galaxies in the AGC208602 field. The velocity field was derived via a standard first moment analysis, using the same pixels that contribute to the column density (moment zero) image; the incoherent velocity structure of AGC208602 precluded meaningful Gaussian fitting to the full data cube. The H i beam size is 49'', as shown in panel (a).

Standard image High-resolution image

The SDSS and GALEX images of AGC208602 presented in Figure 5 reveal no detectable optical counterpart. The magnitudes derived in Table 2 are thus listed as limits; the corresponding luminosities, M $_{{\rm HI}}$/L ratios (M $_{{\rm HI}}$/L $_{g}$ > 48; M $_{{\rm HI}}$/L $_{B}$ > 45), and M $_{{\rm HI}}$/M $_{\star }$ ratio (>49) should likewise be considered limiting values. Deeper optical imaging of this source might be fruitful, but the interpretation may be complicated by the outer halo of NGC 3457.

AGC208602 stands out markedly in Table 3; it has the largest M $_{{\rm HI}}$/L ratios of any of the five sources in the present work. Further, the irregular morphology and incoherent velocity structure are unlike those found in any of the other systems studied here or in other low-mass galaxies. Based on these factors, and on the sourceʼs membership in the NGC 3370 Group, we conclude that this feature may have a tidal origin. It shares some physical characteristics with other ALFALFA-detected groups with extended tidal structures (e.g., Lee-Waddell et al. 2014 and references therein).

It is not clear which sources in the NGC 3370 Group might have interacted to give rise to this object. At the adopted distance of 18.5 Mpc, the 1.2 × 108 M $_{\odot }$ of H i associated with AGC208602 spans more than 20 kpc. This remarkable cloud is located roughly 65 kpc from UGC 6022 and SDSS J105506.66 + 172745.4, the two gas-rich systems with the closest angular projection on the sky. Further, AGC208602 is located between NGC 3443 (∼126 kpc projected separation) and the early type galaxy NGC 3457 (∼15 kpc projected separation); it may be the case that the feature was associated with NGC 3443, and that a recent close encounter with NGC 3457 has stripped the H i away from NGC 3443. A very deep single-dish map of the entire NGC 3370 Group might be fruitful; this would provide sensitivity to possible low surface brightness tidal structure.

3.3. Systems Near Other Galaxies

The remaining two "Almost Dark" sources are located in sufficiently close angular proximity to other known sources that the assignment of an unambiguous optical counterpart was not possible using the ALFALFA data alone.

3.3.1. AGC226178

AGC226178 is an H i source with no optical counterpart at the ALFALFA centroid. The system is a possible member of the Virgo Cluster; Kent et al. (2008) catalog its H i properties, but find no optical counterpart. They note the possible association with AGC229166, an optically faint, low surface brightness object (not detected in H i by ALFALFA, but cataloged as an AGC galaxy based on the extended optical emission which was apparent when examining the ALFALFA data) located at (α,δ) = (12h46m41fs7, 10°23'09''), roughly 2.3' from the ALFALFA centroid for AGC226178. The VLA images presented here show the H i centroid to be 1' south of AGC229166, and suggest that the H i is associated with a second possible optical/UV counterpart with an irregular morphology.

AGC226178 is in close angular and velocity proximity to other H i-bearing sources in our spectral line data. As shown in Figure 6, two other Virgo dwarf galaxies are detected at high significance: VCC 2034 (AGC221001; V $_{{\rm sys}}$ = 1504 km s−1; angular separation $\simeq $910'' $\simeq $72.8 kpc) and VCC 2037 (AGC221004; V $_{{\rm sys}}$ = 1142 km s−1; angular separation $\simeq $720'' $\simeq $57.5 kpc). We do not find clear evidence for an ongoing interaction between any of these sources at our current sensitivity level.

Figure 6.

Figure 6. H i and SDSS R-band images of the AGC226178 field. Panel (a) shows H i column density contours at the (0.5, 1, 2) × 10 20 cm−2 levels, overlaid on a grayscale representation of the image. Panel (c) shows H i isovelocity contours at the (1145, 1148, 1495, 1500, 1505, 1510, 1515, 1586, 1588, 1590) km s−1 levels overlaid on a color representation of the image; the color bar is set to encompass the full velocity range of H i gas in the three detected galaxies. Panels (b) and (d) show the SDSS r-band image of the AGC226178 field; panel (b) is displayed in a logarithmic stretch to highlight low surface brightness emission, while panel (d) has been smoothed by a Gaussian kernel (see discussion in Section 3.2). The contours in (b) are the same as in panel (a), while the contours in (d) are the same as in panel (c). The very low optical surface brightness source AGC229166 is also labeled in panels (b) and (d). The H i beam size is 49'', as shown in the bottom left of panel (a).

Standard image High-resolution image

At SDSS imaging depths, the optical counterpart of the AGC226178 H i source is ambiguous. To highlight the low surface brightness nature of both AGC226178 and AGC229166, panel (d) of Figure 6 shows the SDSS-r band image from panel (b), smoothed by a Gaussian kernel with σ = 3 pixels. The H i is coincident with a very faint irregular optical counterpart candidate, but offset from AGC229166. Figure 7 shows the H i, optical, and UV properties in a smaller region containing only these two sources; note that the SDSS-r band image in this figure is the same as shown in panel (d) of Figure 6 (i.e., it has been smoothed by a Gaussian kernel). The optical morphology of the possible stellar counterpart is highly irregular, appearing as a number of low surface brightness knots of emission oriented from southwest to northeast. Remarkably, this same morphology is apparent in the color GALEX image presented in panel (c); the blue color is indicative of a significant far-UV flux. Note that AGC229166 is not detected with significance in the GALEX images.

Figure 7.

Figure 7. H i and SDSS R-band images of AGC226178. Panel (a) shows H i column density contours at the (0.5, 0.7, 0.9) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (b) shows H i isovelocity contours from 1584 to 1592 km s−1 in intervals of 1 km s−1, overlaid on a color representation of the velocity field image. Panel (c) shows a color representation of the GALEX UV images, while panel (d) shows the SDSS r-band image (convolved by a Gaussian kernel to highlight low surface brightness emission); the contours in (c) are the same as in panel (a), while the contours in (d) are the same as in panel (b). The velocity field was derived by fitting Gaussian functions to the full data cube. The H i beam size is 49''; the source is unresolved at this angular resolution.

Standard image High-resolution image

The H i associated with AGC226178 shows properties that are akin to those of low-mass dwarf galaxies. The neutral hydrogen is centrally concentrated, reaching a maximum column density of 1 × 1020 cm−2. While the source is unresolved by these VLA observations, there is a clearly defined velocity gradient spanning ∼10 km s−1. The total ALFALFA flux integral implies an H i mass of ∼4 × 107 M $_{\odot }$. Our VLA images recover ∼40% of this flux; there is minor contamination from RFI, but it is likely that some of the H i detected by Arecibo is distributed over larger angular scales than those probed by the present VLA observations.

The optical and UV properties of AGC226178 are extracted at the position of the H i centroid and using the parameters given in Table 2. Note that the elliptical aperture size closely matches the UV and optical morphology of the source. Using the SDSS images, we derive M $_{{\rm HI}}$/L $_{g}$ = 7, M $_{{\rm HI}}$/L $_{B}$ = 6, and ${{M}_{{\rm HI}}}$/M $_{\star }$ = 30. The far-UV flux implies a recent SFR$_{{\rm FUV}}$ = 0.0017 M $_{\odot }$ yr−1.

We conclude that the H i source AGC226178 is a dwarf galaxy. We base this conclusion on two primary lines of evidence. First, there is excellent positional agreement between the H i centroid and the UV-bright stellar component. Second, the UV morphology of the putative stellar component suggests an inclined disk morphology. Assuming circular rotation, the orientation of this disk is perpendicular to the isovelocity contours shown in Figure 7(d); this is exactly what would be expected in a rotating disk seen nearly edge-on.

If the H i source AGC226178 were instead associated with the very low surface brightness object labeled as AGC229166 in Figures 6 and 7, then significant interpretive challenges would arise. First, the excellent positional and kinematic agreements noted above would need to occur by chance with an unrelated background source. Second, the ∼1' angular offset between the H i source and the centroid of AGC229166 corresponds to a minimum ∼5 kpc physical separation (assuming AGC226178 and AGC229166 are at the same distance); moving ∼107 M $_{\odot }$ of H i over such a distance would require energies that could likely only arise during a gravitational interaction or ram pressure stripping event.

An alternative but speculative interpretation is one in which AGC226178 and AGC229166 are in fact the same physical system. There are very weak suggestions of spiral structure in the smoothed SDSS-r band image of AGC229166 shown in Figure 7(d). It is conceivable that this is a nearly face-on, gas-poor system. In such a scenario, AGC226178 could be an outer spiral arm that has experienced some sort of infall or interaction episode that has ignited recent star formation. While such a scenario seems physically unlikely, it cannot be completely ruled out with the present data.

3.3.2. AGC233638

AGC233638 was classified as a pair of galaxies with very low optical surface brightness. It is the most distant source in our pilot sample. The total ALFALFA H i flux integral (of which ∼65% is recovered in our VLA observations) implies an H i mass of 3.3 × 109 M $_{\odot }$. Figure 8 shows that the VLA observations resolve these two sources (AGC233638 and AGC743043), which are separated by ∼3farcm3 (∼108 kpc at the adopted distance of 112 Mpc) and by less than 100 km s−1. The SDSS-r band imaging allows a clear identification of an associated optical counterpart with each H i peak. AGC233638 was included in the "Almost Dark" sample due to the ambiguous nature of multiple possible optical counterparts; the higher spatial resolutions presented here allow us to characterize the source with confidence.

Figure 8.

Figure 8. H i and SDSS R-band images of the AGC233638 field. Panel (a) shows H i column density contours at the (0.5, 1, 2) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (c) shows H i isovelocity contours at the (7330, 7335, 7340, 7345, 7400, 7405, 7410, 7415) km s−1 levels overlaid on a color representation of the image; the color bar is set to encompass the full velocity range of H i gas in the two detected galaxies. Panels (b) and (d) show the SDSS r-band image of the AGC226178 field, displayed in a logarithmic stretch to highlight low surface brightness emission; the contours in (b) are the same as in panel (a), while the contours in (d) are the same as in panel (c). The two systems in the field are labeled. The H i beam size is 60''; at this angular resolution, AGC743043 is unresolved, while AGC233638 is marginally resolved.

Standard image High-resolution image

As Figure 9 shows, the candidate "Almost Dark" source is a star-forming dwarf irregular galaxy. There is coherent and well-ordered rotation spanning a projected velocity width of ∼40 km s−1. The H i surface density maximum is exactly co-spatial with the GALEX centroid. While the GALEX images are shallow, they nonetheless allow a robust measurement of the UV fluxes in both bands (see SDSS and GALEX photometry in Table 2). Based on the far-UV flux, this source has undergone significant star formation during the most recent ∼100 Myr (SFR$_{{\rm FUV}}$ = 0.05 M $_{\odot }$ yr−1). The M $_{{\rm HI}}$/L and M $_{{\rm HI}}$/M $_{\star }$ ratios for AGC233638 are typical of those found in star-forming galaxies: M $_{{\rm HI}}$/L $_{g}$ = 3, M $_{{\rm HI}}$/L $_{B}$ = 3, M $_{{\rm HI}}$/M $_{\star }$ = 14. We conclude that AGC233638 is a very intriguing star-forming dwarf irregular galaxy; it has an absolute magnitude (M $_{B}$ $\simeq $ −17) roughly equal to that of the Large Magellanic Cloud, a total H i mass (M $_{{\rm HI}}$ = 3.2 × 109 M $_{\odot }$) roughly equal to that of the Milky Way, and a stellar mass (M $_{\star }$ = 2.2 × 108 M $_{\odot }$) typical of those found in dwarf galaxies.

Figure 9.

Figure 9. H i, SDSS R-band, and GALEX images of AGC233638. Panel (a) shows H i column density contours at the (0.5, 0.75, 1.0, 1.25, 1.5, 1.75, 2.0, 2.25) × 1020 cm−2 levels, overlaid on a grayscale representation of the image. Panel (b) shows H i isovelocity contours from 7385 to 7420 km s−1 in intervals of 5 km s−1, overlaid on a color representation of the velocity field image. Panel (c) shows a color representation of the GALEX UV images, while panel (d) shows the SDSS r-band image (displayed in a logarithmic stretch to highlight low surface brightness emission); the contours in (c) are the same as in panel (a), while the contours in (d) are the same as in panel (b). The velocity field was derived by fitting Gaussian functions to the full data cube. The H i beam size is 60''; AGC233638 is slightly resolved at this angular resolution.

Standard image High-resolution image

A bright nuclear point source is coincident with the optical counterpart of AGC233638; an SDSS spectrum of this bright central knot is available. The SDSS pipeline software proposes a redshift of z = 0.940. Inspection of the SDSS spectrum suggests that this high-redshift classification may be erroneous; no obvious emission or absorption lines agree with this redshift. Further, this would be highly inconsistent with the present observations of the velocity of the neutral interstellar medium. If this compact source is not associated with AGC233638 itself, then the optical luminosity and resulting stellar mass of AGC233638 will both decrease, raising the ${{M}_{{\rm HI}}}$/L and M $_{{\rm HI}}$/M $_{\star }$ ratios. Follow-up optical spectroscopy of this compact source would be very valuable.

4. DISCUSSION

4.1. The "Almost Dark" Sources in Context

The five "Almost Dark" sources studied in this work are representative of the sample of ALFALFA-discovered systems that are not uniquely identified with optical counterparts using ALFALFA data alone. As a result, most harbor unusually high H i mass to light ratios. In the present work, we use multi-wavelength follow-up observations to determine the nature of these objects: four sources are likely galaxies, while one source is likely a tidal remnant.

Keeping in mind the uncertainties that arise because different works have measured optical properties in different ways, the systems studied here can be compared to both the sample of ALFALFA dwarf galaxies in Huang et al. (2012) and to the full α.40 database presented in Haynes et al. (2011). In Figure 10 we thus show the logarithm of the M $_{{\rm HI}}$/M $_{\star }$ ratio versus the logarithm of the stellar mass for all of these systems. There is a clear trend for higher M $_{{\rm HI}}$/M $_{\star }$ ratios in systems of fainter absolute magnitude (lower stellar mass), although the scatter is significant (becoming more pronounced as stellar mass decreases). With the exception of AGC208399, the five sources studied here preferentially populate the extreme upper left region of this plot; these are gas-rich systems with low stellar masses.

Figure 10.

Figure 10. Plot of the logarithm of the neutral hydrogen to stellar mass ratio, ${{M}_{{\rm HI}}}/{{M}_{\star }}$, versus the logarithm of the stellar mass, for all members of the α.40 catalog (Haynes et al. 2011). All small black points and areas within contours represent galaxies with reliable SDSS optical counterparts in the α.40 catalog. The contours show number densities in intervals of 10%, from 10% to 90% of the number density maximum (which occurs at log(M $_{{\rm HI}}$) = 9.8; see Haynes et al. 2011). The dwarf galaxies from the Huang et al. (2012) sample are shown as thick black points; the five "Almost Dark" sources studied in this work are shown in green. With the exception of AGC233638 (which may have contaminated photometry; see detailed discussion in Sections 1 and 3.3.2), these objects preferentially populate the extreme upper left region of this plot.

Standard image High-resolution image

A second comparison can be made against targeted studies of gas-rich, low surface brightness dwarf galaxies. For example, van Zee et al. (1997) computes the M $_{{\rm HI}}$/L ratios in a sample of seven low surface brightness dwarfs and a comparison sample of five dwarfs with more typical hydrogen mass to light ratios. All twelve of these systems have M $_{{\rm HI}}$/L $_{B}$ values near unity (the lowest ratio being 0.6 in UGC 891 and the highest ratio being 6.0 in UGC 5716). The M $_{{\rm HI}}$/L $_{B}$ values of the four "Almost Dark" systems in the present study that appear to be bona fide galaxies (see Table 3) range from 3 (AGC233638) to 9 (AGC193953). The galaxies in the present work thus harbor M $_{{\rm HI}}$/L $_{B}$ that are comparable to the most extreme systems in van Zee et al. (1997).

Further comparisons can made with studies of low surface brightness galaxies that display elevated M $_{{\rm HI}}$/L ratios (e.g., Longmore et al. 1982; van Zee et al. 1995; de Blok et al. 1996; Matthews & Gallagher 1997; O'Neil et al. 2000; Doyle et al. 2005). While some physical characteristics of the present "Almost Dark" objects are similar to those found for objects in these previous studies, the ALFALFA "Almost Dark" sources constitute a different class of object. For example, all of the low surface brightness sources in the de Blok et al. (1996) sample that are within the ALFALFA footprint are detected in H i, and all of them are assigned an optical counterpart by comparison with SDSS imaging. In contrast, the ALFALFA "Almost Dark" sources represent the most extreme objects (<0.5% of the full ALFALFA catalog) for which an optical counterpart remains elusive, even when comparing the sensitive ALFALFA and SDSS/DSS2B survey products.

4.2. The Physical Characteristics of the "Almost Dark" Sources

It is tempting to study the nature of the recent (within ∼100 Myr) star formation in those systems in which star formation is apparent based on a measurable GALEX far-UV flux (AGC208399, AGC226178, AGC233638). However, with the angular resolution afforded by these D configuration observations, all of these sources are only marginally resolved (if at all). The resulting H i column densities and corresponding mass surface densities will necessarily be lower limits that will increase upon observation at higher angular resolution. Taken at face value, the observed column density maxima (N $_{{\rm HI}}$ = (1–2) × 1020 cm−2; see Figures 3, 7, and 9) are significantly lower than the peak column densities observed in the more nearby galaxies of the van Zee et al. (1997) sample (N $_{{\rm HI}}$ = (7–30) × 1020 cm−2). Higher spatial resolution observations of these "Almost Dark" sources would be valuable for studying the nature of the star formation law in the extreme M $_{{\rm HI}}$/L regime by allowing us to resolve the H i surface densities in the disks and correlating these values with well-understood tracers of ongoing star formation.

The elevated M $_{{\rm HI}}$/L and M $_{{\rm HI}}$/M $_{\star }$ ratios in the four non-tidal "Almost Dark" candidates raise interesting questions about galaxy evolution. The H i masses are not unusually low for any of these sources; as with all known non-spheroidal galaxies, an H i reservoir larger than 107 M $_{\odot }$ appears to be at least partially converted into stars. However, in these sources, this conversion appears to be extremely inefficient. Using the UV-based SFRs tabulated in Table 3, the gas consumption timescales exceed the Hubble time for all sources. These "Almost Dark" galaxies are an extreme case of inefficient star formation, similar to, but more dramatic than, what is found in the outer regions of spiral disks and in gas-dominated dwarfs (see Bigiel et al. 2008).

The angular resolution of the H i observations presented here results in physical resolution elements ranging from 2.5 (AGC208399) to 32.6 kpc (AGC233638). This precludes rotation curve work and a detailed mass decomposition. It is interesting to note in this context that the objects identified as galaxies do show irregular kinematics and morphology compared to dynamically relaxed systems. For example, AGC208399 and AGC226178 show lopsided kinematics; AGC193953 shows a misalignment between the total H i intensity and the kinematical axes; AGC223638 may be warped or show symmetric (bar-like) streaming motions. Higher resolution observations of the nearest of these sources would facilitate an examination of the dark and baryonic components in these systems.

5. CONCLUSIONS AND OUTLOOK

We have presented new H i spectral line observations of five sources extracted from the ALFALFA catalog as being representative members of the "Almost Dark" population: systems for which an optical counterpart cannot be assigned using ALFALFA data alone. These sources harbor extreme hydrogen mass to optical luminosity ratios. With the optical, UV, and interferometric H i imaging presented in this work, we now have a much more clear understanding of the nature of these objects; none of them appears to be a bona fide "dark galaxy."

The H i images of these sources localize the H i gas and provide a coarse representation of the rotation and dynamics of the neutral interstellar medium. We compare these H i results with SDSS and GALEX imaging. In four of the five systems, a stellar counterpart is detected that is co-spatial with the H i centroid. Each of these objects shows signatures of coherent rotation; we classify these systems as galaxies.

One system in this study (AGC208602) lacks a detectable stellar population in SDSS and GALEX imaging. This source has the largest M $_{{\rm HI}}$/L ratio in our sample, has an irregular H i morphology, and lacks clear signs of coherent rotation. Given its location in the NGC 3370 group, we conclude that this feature has a tidal origin.

The gaseous and stellar components of AGC226178 make it a very intriguing system. It is located within a projected separation of ∼75 kpc from VCC 2034 and VCC 2037, both of which are detected in our H i imaging. It is also in very close angular projection (∼1', or ∼5 kpc at the adopted distance) to AGC229166, a source with very low optical surface brightness that is not detected in H i by either ALFALFA or the present observations or in the UV in GALEX imaging. Surprisingly, these GALEX images reveal a significant UV flux from an irregular-morphology source that is coincident with the H i centroid of AGC226178. We interpret this source as the stellar counterpart of AGC226178 and thus classify it as a dwarf galaxy. The nature of AGC229166, and its possible association with AGC226178, remains unclear with the present data.

The remaining three sources in this study (AGC193953, AGC208399, and AGC233638) appear to be galaxies with higher than average M $_{{\rm HI}}$/L ratios in comparison with those found in the samples of van Zee et al. (1997), Haynes et al. (2011), and Huang et al. (2012). Of these sources, AGC193953 and AGC208399 are the most promising for more detailed investigation. With the highest M $_{{\rm HI}}$/L ratios of the galaxies in this study, and very low UV-based star formation rates, these sources can facilitate an exploration of the star formation threshold in a sparsely explored region of parameter space. AGC233638 is a comparatively distant source with a lower mass to light ratio; higher resolution observations of this source will be less valuable.

As of this writing, there are roughly 50 viable "Almost Dark" candidate sources in the ALFALFA catalog. The present pilot study of five sources from this larger sample demonstrates the complexity of the members of this heterogenous class of objects: the higher angular resolution VLA images localize the H i gas (thereby facilitating an association with a stellar component for most of these objects) and reveal ordered motions indicative of galaxy-like objects. The discovery potential of a dedicated follow-up observing campaign to study the gaseous components of all candidate "Almost Dark" systems is thus exceptionally high.

Our team is pursuing similar H i observations of the full "Almost Dark" sample, as well as deep optical images of the sources with large FOV instruments. For example, preliminary results from these initiatives demonstrate significant similarities between AGC226178 (studied here) and AGC229385 and its neighbors (Janowiecki et al. 2015). These future works hold significant promise for furthering our understanding of the evolution of extremely gas-rich systems. They will allow us to explore the nature of star formation in the high M $_{{\rm HI}}$/L regime; the preliminary results from the present work suggest a very inefficient or suppressed star formation process. Most importantly, these endeavors will allow us to quantify the frequency of truly "dark" galaxies (if they exist), using the statistically robust ALFALFA database with well-understood sensitivity and completeness limits.

J.M.C. is supported by NSF grant 1211683. J.M.C., R.A.K., C.P.M., N.N., and P.T. acknowledge NSF grant AST-1211005. The ALFALFA work at Cornell is supported by NSF grant AST-1107390 and by the Brinson Foundation. K.L.R. acknowledges support from NSF grant AST-0847109. J.M.C. would like to thank the Instituto Nazionale di Astrofisica and the Osservatorio Astronomico di Padova for their hospitality during a productive sabbatical leave. L.L. would like to thank Angus Wright for technical support using his code, LAMBDAR. G.J. started to participate in this project as an employee of ASTRON, Dwingeloo, The Netherlands (Nederlandse Organisatie voor Wetenschappelijk Onderzoek, NWO). This investigation has made use of 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, and NASAʼs Astrophysics Data System. This investigation has made use of data from the Sloan Digital Sky Survey (SDSS). Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.

Footnotes

  • 12 

    The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

  • 13 
  • 14 
  • 15 
  • 16 

    Based on observations made with the NASA Galaxy Evolution Explorer. GALEX is operated for NASA by the California Institute of Technology under NASA contract NAS5-98034.

Please wait… references are loading.
10.1088/0004-6256/149/2/72