Long-slit Spectroscopy of Edge-on Low Surface Brightness Galaxies
We present long-slit optical spectra of 12 edge-on low surface brightness galaxies (LSBGs) positioned along their major axes. After performing reddening corrections for the emission-line fluxes measured from the extracted integrated spectra, we measured the gas-phase metallicities of our LSBG sample using both the [N ii]/H and the diagnostics. Both sets of oxygen abundances show good agreement with each other, giving a median value of 12 + log(O/H) = 8.26 dex. In the luminosity-metallicity plot, our LSBG sample is consistent with the behavior of normal galaxies. In the mass-metallicity diagram, our LSBG sample has lower metallicities for lower stellar mass, similar to normal galaxies. The stellar masses estimated from -band luminosities are comparable to those of prominent spirals. In a plot of the gas mass fraction versus metallicity, our LSBG sample generally agrees with other samples in the high gas mass fraction space. Additionally, we have studied stellar populations of 3 LSBGs which have relatively reliable spectral continua and high signal-to-noise ratios, and qualitatively conclude that they have a potential dearth of stars with ages Gyr instead of being dominated by stellar populations with ages Gyr. Regarding the chemical evolution of our sample, the LSBG data appear to allow for up to 30% metal loss, but we cannot completely rule out the closed-box model. Additionally, we find evidence that our galaxies retain up to about 3 times as much of their metals compared with dwarfs, consistent with metal retention being related to galaxy mass. In conclusion, our data support the view that LSBGs are probably just normal disk galaxies continuously extending to the low end of surface brightness.
Low surface brightness galaxies (LSBGs) were theorized to exist by Disney (1976) and then verified by the discovery of Malin 1 (Bothun et al., 1987). They refer to galaxies which have lower surface brightness than the background (Impey Bothun, 1997) and are traditionally defined as galaxies with central surface brightness fainter than some threshold which varies between 22.0 and 23.0 mag arcsec (e.g., Impey et al., 2001; Ceccarelli et al., 2012). Recent studies indicate that LSBGs contribute 20% of the dynamical mass of galaxies in the Universe (Minchin et al., 2004) and 30%–60% of the number density of local galaxies (McGaugh, 1996; Bothun et al., 1997; O’Neil et al., 2000; Trachternach et al., 2006; Haberzettl et al., 2007), suggesting that the contribution of LSBGs to the Universe can not be negligible. However, selection effects cause LSBGs to be underrepresented in the present-day galaxy catalogs.
Compared with disk galaxies having high surface brightness (HSB), LSBGs generally tend to be dust-poor (e.g., Matthews Wood, 2001) and H i-rich, with diffuse, low-density stellar disks and metallicity solar abundance (e.g., McGaugh, 1994). This suggests that LSBGs have been inefficient star formers over their lifetimes and have low-level star-formation activity at present (e.g., Das et al., 2009; Galaz et al., 2011). Indeed, H i surface densities in LSBGs are often observed to be well below the critical threshold for star formation (e.g., van der Hulst et al., 1993). Nonetheless, LSBGs typically contain some signatures of ongoing star formation, including modest amounts of H emission from a very small number of H ii regions Schombert et al. (e.g., 1992, 2013) and blue colors indicative of young stellar populations (e.g., McGaugh Bothun, 1994; Matthews Gallagher, 1997).
Metallicities of stars and gas within a galaxy provide both a fossil record of its star formation history (SFH) and evolution. For galaxies, metallicities are traditionally represented by the oxygen abundance, 12 + log(O/H), which directly reflects the metal content of the gas in galaxies and indirectly reveals the metal content of the stars in galaxies that ionize the gas. For LSBGs, many previous works have demonstrated that LSBGs are metal-poor. The LSBG sample of McGaugh (1994) has been measured to have low oxygen abundances, 12 + log(O/H) dex. Burkholder et al. (2001) found a mean value 12 + log(O/H) = 8.22 dex for 17 LSBGs. Lam et al. (2015) studied KKR17, an LSBG dominated by H i, and measured an average 12 + log(O/H) = 8.0 dex. The oxygen abundances of de Naray (2004) are 12 + log(O/H) –8.2 dex. Apparently, most LSBGs have 12 + log(O/H) lower than the solar value of oxygen abundance of 12 + log(O/H) = 8.87 dex (Grevesse et al., 1996), and even lower than the transition value (12 + log (O/H) = 8.4 dex) between metal-poor and metal-rich branches for the method (see §4.3). To derive gas-phase metallicities of LSBGs, spectra should be required. However, it is not easy to derive spectra of LSBGs due to their low surface brightness. Fortunately, we have carried out spectroscopic observations for a sample of LSBGs. Later in this paper, we will introduce this observation and then study properties which are mainly relevant to metallicities on basis of the observational data.
Given the faint nature of LSBGs, with peak surface brightness mag arcsec, it is extremely difficult to obtain spectra of their stellar continua with relatively high signal-to-noise ratio () because the sky background always dominates the signal. Studies of the stellar populations of LSBGs are currently underrepresented and are based mainly on spectra of their H ii regions (e.g., Schombert et al., 2013) and central portions (e.g., Morelli et al., 2012). But the bright H ii regions contain a large fraction of O- and B-type stars, so stellar populations studied through spectra of the H ii regions would be biased to younger ages. Compared to observing face-on LSBGs (Du et al., 2015, e.g.,) or fiber spectroscopy, the spectral qualities can be highly improved by observing edge-on LSBGs, especially those having thin disks, with a single slit of appropriate length and width, positioned along the major axis. This scheme makes it possible, although still very difficult, to obtain reasonable-quality data. Very fortunately, spectra of several few LSBGs of our sample have relatively reliable stellar continua which can be used for stellar population studies, so we additionally study stellar populations of these few LSBGs and to investigate their star formation histories as byproducts of our observations (§6.1).
The outline of this paper is as follows. Section 2 describes our sample, and detailed accounts of the observations and data reduction are given in §3 and §4. We study the physical properties of our edge-on LSBGs in §5, including oxygen abundances in §5.1, the luminosity-metallicity relation in §5.2, the mass-metallicity relation in §5.3, and the behavior of gas mass fraction versus metallicity in §5.4. In Section 6, we discuss stellar populations in §6.1, stellar mass estimation in §6.2, chemical evolution in §6.3, and the true nature of our LSBG sample in §6.4. Finally, we conclude in §7.
2 The Sample
A catalog of edge-on disk galaxies in SDSS DR7 (EGIS; Bizyaev et al., 2014) consists of 5747 genuine edge-on galaxies in the Sloan Digital Sky Survey (SDSS; Abazajian et al., 2009). By analyzing photometric profiles in each of the -, -, and -band images, the EGIS catalog provides estimates for structural parameters of the stellar disks of the galaxies, such as the disk scale height (), radial scale length (), and central surface brightness (). To further study the EGIS catalog in combination with H i data, we correlated it with the .40 catalog (Haynes et al., 2011), which lists H i line sources from 40% of the area of the final blind 21 cm H i ALFALFA (Arecibo Legacy Fast Arecibo L-band Feed Array) survey for local galaxies. Using the transformation formula of Smith et al. (2002), we derived the central surface brightness in the band, , from those in the and bands. Then, according to a widely used definition of a LSBG, there are 287 galaxies in the EGIS-.40 sample that can be identified as LSBGs with mag arcsec.
Out of the EGIS-.40 LSBG sample, only 21 galaxies (blue filled circles in Figure 1(a),(b)) are in the Fall sky region () which matched our available observing time. Of the 21 galaxies, 14 galaxies have mag (dashed black line in Figure 1), sufficiently blue to be considered as the final targets for our 3 nights in late-September 2015. Unfortunately, only 50% of the total time was useful for observations because of stormy and cloudy conditions; we ultimately observed 9 galaxies. Comparing with the 14 prepared targets (magenta lines), the observed 9 galaxies (red) are generally consistent in distributions of , magnitude, color, and as shown in Figure 1(c)–(f). Furthermore, comparing with the entire EGIS-.40 LSBG sample (black), these 9 galaxies (red) sample almost the full range in from 22.5 to 24.5 mag arcsec, but they sample only the fainter part in ( mag) and only the blue part in ( mag).
On the other hand, in February 2015, we had already spectroscopically observed in a similar way several galaxies (not confined to be LSBGs) in the EGIS-.40 catalog, but with a different telescope using spare time from another program. From these objects, there are 3 galaxies with mag arcsec. We overplot them (blue) in Figure 1 and find that they are comparable to our observed (Fall 2015) 9 LSBGs in distributions of main properties. So, we add these 3 LSBGs into our 9-galaxy LSBG sample in order to enlarge the overall observed sample in this paper.
We display the SDSS images of the 12 edge-on LSBGs in Figure 2 and list their main properties in Table 1. Statistically, values of the 12 LSBGs range from 22.6 to 24.4 mag arcsec. There is one object (AGC 228010; see Figure 2) with ; it is a super-thin disk galaxy, previously defined to have an axial ratio between the disk scale length () and the scale height () in the band larger than 10.0 (van der Kruit, 1981, 1988; van der Kruit et al., 2001; de Grijs, 1998). Additionally, our 12 galaxies have little (or no) optical bulges or central light concentrations, with in the SDSS band, but they are not the only ones with (see Figure 1(f)).
|Name||(J2000)||(J2000)||aaRadial velocity from Haynes et al. (2011).||Distance||bb-band central surface brightness from Bizyaev et al. (2014).|
|(deg)||(deg)||(km s)||(Mpc)||(mag arcsec)||(mag)||(mag)||(arcsec)||(arcsec)|
|PA||ccBulge-to-total luminosity ratio of the galaxy from the SDSS catalog.||log((H i)/M)|
|(deg)||(km s)||(Jy km s)|
Long-slit spectroscopy of three LSBGs in the Spring sky were obtained on 2015 Feb. 18 UT with the Kast double spectrograph (Miller Stone, 1993) at the Cassegrain focus of the 3-m Shane telescope at Lick Observatory. A long slit 2 wide was positioned along the major axis of each target, and the seeing was . For the blue side, a 452 lines mm grating blazed at 3306 Å was used in first order, yielding a wavelength range of 3150–6300 Å with a spectral resolution of 1.41 Å pixel. Use of a dichroic allows simultaneous observations at longer wavelengths with the red side; however, Lick spectroscopy was obtained only with the blue side.
The other nine LSBGs were observed on 2015 Sep. 21 UT (seeing ) and Sep. 23 UT (seeing ) with the 5-m Hale telescope at Palomar Observatory. We used the Double Spectrograph (DBSP; Oke & Gunn, 1982) and a -long slit of width 2, positioned along the major axis of each target. A D68 dichroic was used to divide the wavelength coverage into two regions, 3671–6727 Å for the blue side (600 lines mm grating blazed at 4000 Å) and 6223–9577 Å for the red side (600 lines mm grating blazed at 10,000 Å). The respective spectral resolutions are 1.078 Å pixel and 0.819 Å pixel, so both gratings resolve all possible line features of interest to us.
The total exposure times at Lick and Palomar for a galaxy were 3600–7200 s, subdivided into 2–4 subexposures in order to remove cosmic rays and enhance the signal-to-noise ratio (S/N) of the combined spectra. See Table 2 for a log of observations.
3.2 Light-Loss Problem
The atmosphere is a dispersive medium. At nonzero zenith angles, a stellar image will be spread out into a small spectrum, and this spread can be significant compared to the size of the seeing disk. Under these conditions, if the spectrograph slit is not oriented along the direction of atmospheric dispersion (the parallactic angle), starlight at some wavelengths will not enter the narrow slit, while starlight at other wavelengths will pass into the instrument relatively unhindered.
During our observations, the slit was actually oriented along the major axis of the target object (position angle) instead of the parallactic angle, and the object was centered in the slit at Å. Therefore, at nonzero zenith angles, object light from either wavelength end (blue or red) may be selectively lost. For example, as seen in Table 2, the spectrum of AGC 191770 at airmass 1.08 was obtained with a 128 2 slit on a night when the seeing was 1.5. Then, following the calculations of Filippenko (1982), if AGC 191770 was centered in the slit at a wavelength of 6000 Å, it would be off center by 0.61 at 3600 Å and by 0.05 in the opposite direction at 6500 Å . The parallactic angle was east of north, but the slit was actually aligned along a position angle of east of north. So, if AGC 191770 was a star with a Gaussian point-spread function, the percentage of light entering the slit at 3600 Å , 6000 Å , and 6500 Å would be (respectively) 83.5, 88.3, and 88.3. Thus, the light loss at 3600 Å would be 4.8% greater than that at 6000 Å.
Differential losses are a bigger concern for point sources than for extended objects such as galaxies. For extended sources, especially if the surface brightness is roughly constant over an area larger than the slit width, blue-end light loss from the position being observed can be partly gained back from neighboring positions within the source. The effect on the relative flux calibration will not be very large for our targets because of the following reasons. (1) Most spectra of our targets were observed at small airmasses (Table 2). (2) Our targets are edge-on galaxies with extended areas larger than the slit width (see Figure 2), so the blue and red light losses are partly gained back from the neighboring positions within the objects. For completeness, we do the calculations above for all our targets (at the midpoint of the exposures) and list the results in Table 2.
|Name||UT Date||Exp.||blue range||red range||Tel. & Instr.||airmass||Position Angle|
|YYYY-MM-DD||(s)||(Å)||(Å)||deg (East of North)|
|deg (East of North)||arcsec||arcsec||arcsec|
4 Data Reduction
Data reduction was carried out using standard IRAF procedures. After correction for overscan, bias, and flat field, we reject cosmic rays using the IRAF task L.A.Cosmic (van Dokkum et al., 2001), which identifies and removes cosmic rays robustly from astronomical images via a Laplacian algorithm. Then, we performed wavelength calibration on the two-dimensional (2D) spectral images of objects and standard stars. Sky subtraction will be described more detailedly in next subsection. After sky subtraction, we defined the appropriate aperture, traced the flux peaks along the dispersion direction, and interactively extracted one-dimensional (1D) spectrum by using the APALL task of IRAF. Finally, we averagely combined spectra from multiple exposures for an object into one final spectrum.
In Figure 3, we show the combined blue-side spectra for the three galaxies observed with the Lick Shane telescope (they have no red-side observations). For the nine galaxies observed with the Palomar Hale telescope, both the blue- and red-side spectra are shown (with blue and red colors, respectively), after scaling the red side to match the flux level of the blue side in the small overlap region. Strong emission lines are marked. Generally, the night-sky lines are not clearly subtracted in the faint LSBG spectra, so we simply mask four wavelength regions having the strongest night-sky lines (green).
4.2 Sky Subtraction
Sky subtraction quality is very important for LSBGs as they are naturally fainter than the sky background in surface brightness. So, we describe the three methods of sky subtraction we have tried in data reduction for our observations below.
The first one is called the “median-filter method.” Before performing this method, we make the spatial and dispersion directions perpendicular to each other through rectification, such as by performing 2D wavelength calibration on spectral images of objects. Then, on the rectified images, assuming the direction of night-sky lines is along the column, we perform a median filter represented by a array along each column to generate a sky image. Here, must be larger than at least two times the width (major axis in pixels) of the object projected on the images. This method is good at removing sky lines but may overestimate the sky continuum blended in the object.
The second one is the sky-fitting method adopted in the APALL task of IRAF, which uses a second-order polynomial to fit the background from the substantial sky region along the slit. This method can do well in estimating sky continuum blended in the object but does not do very well in removing sky lines.
The third one is the Kelson (2003) method which makes full use of the raw image data (before the optical distortions and spectral line curvatures have been rectified) to model the sky background. This method is not very good for our faint objects, because it either underestimates the sky continuum or leaves sky-line residuals in our images.
None of the three methods above is sufficiently good at removing both sky continuum and sky lines clearly from spectra of our faint LSBGs. However, we have no other good ideas to implement it ideally. So, for Lick data which have no obvious distortions of raw images, the first “median-filter” method has been used to simultaneously subtract the sky lines and sky continuum. However, for Hale data which have relatively serious distortions of raw images, this “median-filter” method does not do well in sky continuum subtraction, so we finally used the second method to improve the sky continuum subtraction a little better, but we warn the reader that the final spectral continuum shapes for Hale spectra are still probably not accurate.
4.3 Problems in the Spectra
Examining the extracted spectra, it appears that the reductions of the Hale spectra are not good at wavelengths Å. Unexpectedly, many absorption lines appear in the observed spectra of our standard stars in this blue-end region due to unresolved reasons, so we are not able to determine an accurate continuum shape for fitting the flux-calibration curve in this blue-end region. Fortunately, except for AGC 123165, our LSBG galaxies all have redshift , so the bluest useful line ([O ii] 3727) falls at Å.
Although we have extracted the red-side Palomar Hale spectra, their continuum S/N is quite low, so we abandoned using the red-side spectral continuum. However, the emission lines in the red spectra, such as [N ii] 6584, H, and [S ii] 6716, 6731, are reasonably strong, and they could therefore be used for dust-extinction corrections and oxygen-abundance measurements.
Moreover, as mentioned above, the spectral continuum shapes of the nine galaxies observed with the Hale telescope are not reliable, so we will not use them to draw any firm conclusions.
5.1 Oxygen Abundance
The gas-phase metallicity (oxygen abundance, log(O/H)) of a galaxy is a powerful probe of the past star formation and growth. It can also constrain the current stage of galaxy evolution as it is intimately connected to the fueling and regulation of star formation. Usually, flux ratios of strong optical emission lines are used as the abundance-sensitive diagnostics, such as , [N ii]/[O ii], [N ii]/H, [N ii]/[S ii], and so on (Pagel et al., 1979; Kewley Dopita, 2002). Based on extracted spectra of our LSBGs, we measure the observed fluxes of strong emission lines in the blue-side ([O ii] 3727, H, [O iii] 4959, 5007) and red-side ([N ii] 6548, 6584, H, [S ii] 6716, 6731).
An intrinsic H/H flux ratio of 2.86 is assumed for Case B recombination at K (Osterbrock, 1989) because the relative intensities of the hydrogen Balmer lines are nearly independent of density and temperature. Therefore, for observed line fluxes measured from the 9 Hale spectra which have both blue- and red-side, we use the Balmer decrement (H/H) and the CCM extinction law with (Cardelli et al., 1989) to correct dust extinction. Based on these intrinsic fluxes, we compute several useful flux ratios of [N ii]/[O ii], [N ii]/[O ii], [N ii]/H, and [where ([O ii] 3727 + [O iii] 4959, 5007)/H], among which [N ii]/H and are individually used as metallicity diagnostics. For the [N ii]/H diagnostic, both the PP04 (Pettini Pagel, 2004) and D02 (Denicolo et al., 2002) calibrations are used to measure the oxygen abundances. For the diagnostic, both the M91 (McGaugh, 1991) and Z94 (Zaritsky et al., 1994) calibrations are employed to measure the oxygen abundances. Thus, we have four sets of oxygen abundances for each of the 9 Hale-observed galaxies except for AGC 123165 which have no –based metallicities because its spectra has some problem in the blue-end wavelength range covering [O ii] 3727 (see §4.3). It should be noted that when using the M91 calibration, we need to determine which branch is appropriate for each galaxy because the versus O/H relation is double-valued (Kewley Dopita, 2002; de Naray, 2004). Distinguishing a branch (high- or low-) can be accomplished with the [N ii] 6584 line: strong (weak) nitrogen emission indicates that the high- (low-) branch is appropriate. The dividing line is at log([N ii]/[O ii]) (Kewley Dopita, 2002). However, for observed line fluxes measured from the 3 Lick spectra which have no observed red-side spectra (neither H nor [N ii]), no reddening correction is performed and only diagnostic and –based metallicities are measured. When using the M91 calibration, no available log([N ii]/[O ii]) can be used for distinguishing branches, so we empirically adopt the low- branch to measure metallicities as LSBGs are generally acknowledged to be metal-poor.
In Table 3, we report measurement results above and also show the average () of both -based metallicities ( and ) and the average () of both [N ii]/H-based metallicities ( and ) for each LSBG. In statistic, our measurements show that the -based metallicities () are 8.00–8.62 dex, with a mean of 8.30 dex and a median of 8.26 dex. The [N ii]/H-based metallicities () are 7.67–8.54 dex, with a mean of 8.20 dex and a median of 8.26 dex. Generally, the two sets of metallicities are consistent with each other and are both subsolar, with the log(O/H) value being 8.87 dex for the Sun (Grevesse et al., 1996). Our metallicities for the LSBG sample agree with the results for LSBGs from de Naray (2004) (see Figure 4).
|AGC No.||(H/H)||log([N ii]/[O ii])||log([N ii]/H)||bb is the average value of and , which are both based on the [N ii]/H diagnostic.||ccMF refers to mass fraction to the total stellar mass.|
5.2 The Luminosity-Metallicity Relation
The strong correlation between galaxy luminosity () and galaxy metallicity () is one of the most significant phenomenological results in chemical evolution studies of galaxies. Previous work found that the correlations between and for irregular galaxies (Lequeux et al., 1979; Skillman et al., 1989), disk galaxies (predominantly high surface brightness galaxies, HSBGs; Garnett et al. 1987), and elliptical galaxies (Faber, 1973; Brodie Huchra, 1991) are very similar. Thus, it is instructive to plot our LSBG sample onto the familiar – relation defined by previous work on other galaxy types to validate the consistency of our measurements of LSBGs.
In the – plot in Figure 4, we show both -based (red filled circles) and the [N ii]/H-based (black filled circles) metallicities for our LSBGs. Here, the -band absolute magnitude for our LSBGs are computed from the galaxy distances (provided by 40 catalog) and the -band apparent magnitudes (which are transformed from magnitudes of - and -bands measured by our own surface photometry (Du et al., 2015)). Data of LSBGs from de Naray (2004) (purple triangles) and dwarf irregulars (dIrrs; dark-green squares) from Lee et al. (2003) are also plotted in this – figure for comparisons. Besides, we show a variety of recent fits to the – relation determined respectively for the dIrr sample of Lee et al. (2003) (green line), a large sample of SDSS star-forming galaxies from Tremonti et al. (2004) (orange line), a large sample of starbursting emission-line galaxies drawn from the KISS (KPNO International Spectroscopic Survey) sample from Melbourne Salzer (2002) (cyan line), and a large sample of 2dFGRS (Two-Degree Field Galaxy Redshift Survey) star-forming galaxies from Lamareille et al. (2004) (blue line) from the literature. Among these relations, it seems that there is some disparities, which are probably caused by the utilization of different samples and different diagnostics and calibrations for measuring abundances.
It is not easy to constrain a quantitatively reliable – relation for our LSBG data; they do not have sufficient leverage because of their small coverage in luminosity space. As shown in Figure 4, if all 12 galaxies (red filled circles) are considered, our sample covers only mag of absolute magnitude, and if only the 9 galaxies (black filled circles) having both blue- and red-side spectra are considered, it covers only mag. Instead of being quantitative, we thus carry out our studies qualitatively. From Figure 4, it appears that our LSBG sample with both sets of metallicities, especially with the [N ii]/H–based metallicities, seems to generally follow the similar distribution as the LSBGs from de Naray (2004) in their common covering region in luminosity and metallicity space. Both LSBG samples are consistent with the general distribution trend of other samples (normal galaxies, dIrrs) in this – plot, although with some degree of scatter.
Additionally, we checked for quantitative consistency between our sample and the LSBG sample from de Naray (2004). We did the test 50 times. Each time, we randomly picked a subsample of size equal to our small LSBG sample with [N ii]/H-based metallicities from the LSBG sample of de Naray (2004), and then checked to see whether this subsample gives a distribution consistent with our small LSBG sample by using a K-S test. Our criteria for each K-S test to decide whether the subsample is consistent with being drawn from the same distribution as our small LSBG sample is the probability value being larger than 0.90 (P0.90). The result is that in the total of 50 trials, 32% of the subsamples from the de Naray (2004) LSBG sample give distributions that are consistent with our limited LSBG sample in the – plot. Any fraction more than 10% and less than 90% is basically inconclusive. So, such a fraction shows that the samples are too small to provide a conclusive result. If our sample could be improved to a larger size and wider luminosity coverage, such quantitative tests might give a more conclusive statistical fraction for consistency.
5.3 The Mass-Metallicity Relation
5.3.1 Stellar Mass Estimate
For our LSBGs, we derived two sets of stellar masses. One set is the “-based” stellar masses which are derived by using their -band luminosity and color transformed from measured in (Du et al., 2015). The other set is the “-based” stellar masses in the SDSS MPA-JHU catalog which are computed from the SDSS -band luminosities in synergy with the -band stellar derived by the SDSS multi-band SED fitting method (Kauffmann et al., 2003). Resultantly for our LSBGs, the “-based” stellar masses are – M, with a mean of M and a median of M. The “-based” stellar masses are – M, with a mean of M and a median of M. We report the results in Table 4. Apparently, the set of “-based” stellar masses is on average order of magnitude lower than the set of “-based” stellar masses. In this paper, we prefer the “-based” stellar masses for our LSBGs, and the reason will be given in §6.2.
5.3.2 – plot
We show our LSBG data in the mass-metallicity (–) plot in Figure 5. For metallicity, both -based (black) and [N ii]/H-based (red) oxygen abundances are displayed. For stellar mass, both the “-based” (open circles) and “-based” (filled circles) stellar masses are shown. So, in this – plot, our LSBG sample appears as four individual sets of distributions, among which the two sets represented by red and black open circles are systematically biased to lower stellar mass because the “-based” stellar masses may be underestimated (as mentioned in §5.3.1 and discussed in §6.2). So, comparing with distributions represented by the red and black open circles, we prefer distributions based on the “-based” stellar mass estimations (red and black filled circles). For comparison, the LSBG sample (purple triangles) of de Naray (2004), the dIrr sample (dark-green squares) of Lee et al. (2003), and the massive spiral galaxy sample (brown asterisks) of Garnett et al. (1987) in the very local universe (including M101, M31, NGC 1365, M81, NGC 2997, and so on) are also included in this – plot. Additionally, the – relation determined for SDSS star-forming galaxies of Tremonti et al. (2004) is shown. The black open diamonds represent the median in bins of 0.1 dex in stellar mass that include at least 100 data points. The magenta line shows the polynomial fit of Tremonti et al. (2004) to the median data. The cyan blue dashed and solid lines are the contours that enclose 68% and 95% of the data.
Obviously, the distributions of our LSBG sample with “-based” stellar mass are generally consistent with the LSBG sample from de Naray (2004) at their common stellar mass space. Additionally, for normal star-forming galaxies and other types, there is a tendency that galaxies with lower stellar masses generally have lower metallicities. Specifically, for normal star-forming galaxies with stellar masses between and M, Tremonti et al. (2004) find the correlation between stellar mass and metallicity to be roughly linear. Owing to the small coverage in stellar mass space, our LSBG sample is not sufficient for fitting a certain – relation, but galaxies in our sample are consistent with the relatively low abundances for their stellar masses.
5.4 The Gas Mass Fraction vs. Oxygen Abundance
The gas consisting mostly of hydrogen and helium is the raw material out of which stars and metals are formed. In LSBGs, the largest constituent of gas is assumed to be neutral atomic hydrogen. The molecular gas remains unknown in quantity but, fortunately, it is not expected to contribute greatly to the total mass of gas in LSBGs due to their low metallicities and low dust content (Schombert et al., 1990; de Blok, 1998). As for ionized gas, it contributes negligibly to the mass of LSBGs (Schombert et al., 2001).
The gas mass fraction () refers to the fraction of baryons in gaseous form. It can be simply measured by the ratio of the total mass in gas (, including molecular, atomic, ionized, and metals) to the combined total mass in gas and stars ():
So, taking the contribution of helium and metals into account, the total gas mass of LSBGs in units of solar mass can be given by
where is the fraction of the gas mass in the form of hydrogen, adopted to be 0.733 (de Blok, 1998).
Following the formula above, for our edge-on LSBG sample can be easily computed by multiplying (Table 1) which is available from the 40 catalog (Haynes et al., 2011) by 1/0.733. For ,we have previously derived two sets of stellar masses for our LSBG sample in §5.3.1 (Table 4). Consequently, we have two sets of gas mass fraction, , for our sample.
In Figure 6, we show oxygen abundance versus inverse gas mass fraction for our LSBG sample. Both sets of gas mass fractions are shown for our sample. The filled circles represent for being computed from the “-based ” while the open circles are for being computed from the “-based” . Besides, two sets of metallicities are also presented with red circles showing -based metallicities and black circles showing [N ii]/H-based metallicities. For comparisons, the LSBG sample from de Naray (2004) (purple triangles) and the dIrr sample from Lee et al. (2003) (green squares) are also plotted.
For our LSBG sample, data points with “-based” stellar mass (red and black open circles) or with “-based” stellar mass but -based metallicities (red filled circles) have systematically biased to much larger gas mass fraction, compared with other samples (e.g., the LSBG sample from de Naray 2004; the dIrr sample from Lee et al. 2003). However, data points with “-based” stellar mass and [N ii]/H-based metallicities (black filled circles) show general consistency with other samples in the common high- region. The systematically larger gas mass fractions are more likely to be caused by different methods used for estimating stellar mass and metallicity. Again, we prefer the “-based” stellar masses for our sample because LSBGs are acknowledged to be dominated by low-mass stars that would radiate light predominantly at red and near-infrared wavelengths. Using -band luminosity method is very likely to underestimate stellar masses of LSBGs and then overestimate their gas mass fractions. More details on this issue will be discussed in §6.2.
6.1 Stellar Populations
6.1.1 Coarse Stellar Populations Results
Stellar populations are the fossils of the past star formation and evolution of galaxies. However, the sky background makes it very difficult to study the faint stellar continuum and absorption lines in the entire LSBG unless very large telescopes and extraordinarily long exposure times are used. Because of such difficulties, most spectroscopic studies of LSBG stellar populations have been based on spectra of the brighter bulges (e.g., Morelli et al. 2012) or H ii regions. A slit of appropriate width and length positioned along the major axis can generally capture light from almost all regions of an edge-on galaxy, especially of a super-thin galaxy. Thus, the spectral properties may in principle be much better determined for edge-on LSBGs compared with face-on ones.
If we expect to derive reliable stellar populations with full-spectrum fitting stellar population synthesis method, we must use high-quality spectra, especially with well-calibrated spectral continuum shapes. Given the uncertainties associated with continuum shapes of the Palomar Hale spectra, we decided to use only the three Lick Shane spectra (of AGC 191770, AGC 228010, and AGC 253926) for our stellar-population study of LSBGs. All these three spectra have S/N in the continuum wavelength interval 4730–4780 Å and meet our requirements.
For stellar population synthesis, STARLIGHT (Cid Fernandes et al., 2005) is a widely used tool. It can derive the age and metallicity of a galaxy, as well as its extinction and velocity dispersion, by fitting its observed spectrum with a superposition of populations of various ages and metallicities from stellar population models such as STARLIGHT (Bruzual Charlot, 2003, hereafter BC03). We used STARLIGHT to perform stellar population synthesis for the integrated spectra of the three edge-on LSBGs observed at Lick. For the reddening, both the Calzetti et al. (2000, hereafter CAL) law and the Cardelli et al. (1989, hereafter CCM) law were separately used during the fit, although the problem of dust attenuation is minimal in LSBGs because they are acknowledged to be dust poor.
We show the fit plot in Figure 7. As said by Cid Fernandes et al. (2005), the individual stellar population components from STARLIGHT are very uncertain because the existence of multiple solutions is an old known problem in stellar population synthesis. Thus, we should avoid focusing in detail on any individual stellar populations. Instead, we jump straight to results based on more robust descriptions of the resultant stellar populations. As suggested in the STARLIGHT paper (Cid Fernandes et al., 2005), a coarse but robust description of the star-formation history of our galaxies may be obtained by binning the resultant individual stellar populations onto “young” (Y; age Gyr) and “old ” (O; age Gyr) components. In order to show the respective contributions of the two components, we provide the “light fraction” (LF) out of the total light and the “mass fraction” (MF) out of the total mass for Y and O components in Table 5. Results from both reddening laws (CAL and CCM) for each LSBG seem to be consistent. Qualitatively, for both light and mass, the contributions from the O components seem to be much larger than those from the Y components for these three LSBGs. Thus, rather than a quantitative conclusion, we can draw a qualitative conclusion that these three LSBGs may have a potential dearth of stars with ages Gyr. Instead, they are more likely to be dominated by stars with ages Gyr. This supports the general notion that LSBGs have a very low level of current star-formation activity.
Comparing with the previously measured subsolar gas-phase metallicities for the first three galaxies in Table 3, their stellar metallicities (from stellar population synthesis; Table 5) are about solar. On one hand, this inconsistency is likely to be caused mainly by the large uncertainties of the stellar metallicities. As can be deduced from the amount of work devoted to this topic of stellar population synthesis, recovering the stellar content of a galaxy with accurately quantitative ages and metallicities from its observed integrated spectrum is not an easy task. It needs a high-quality integrated spectrum, as well as models which are made and calibrated to mimic nature to the greatest degree. Any stellar population model, such as STARLIGHT (BC03), is based on a stellar library. Usually, the stellar library does not sample a sufficiently wide range in metallicity, making stellar metallicity notoriously more difficult to access than other properties. Also, to get the metallicity correct, one must have correct ages, as they affect one another in these integrated spectral fits. Since our derived stellar ages have large uncertainties, the stellar metallicities are necessarily also poorly constrained. So, stellar metallicities determined here from stellar population synthesis (Table 5) carry large uncertainties, and are not reliable as specific values for comparisons and scientific studies; the gas-phase metallicities (Table 3) from the emission-line flux ratios are much more reliable.
6.1.2 Effect of Blue Light Loss on Conclusions
We use the blue-side spectra for studying stellar populations by full-spectrum fitting, so the accuracy of the relative flux calibration becomes an important factor for the reliability of the stellar population results. If the slit is not positioned along the atmospheric dispersion, the light loss at blue-end wavelengths would be greater than that at red-end wavelengths. As shown in Table 2 in Section 3.2, if a star was being centered at 6000 Å and observed with the same slit positioning schemes and observing configurations as those in our program, it will be off center by 0.61 at 3600 Å and 0.05 in the opposite direction at 6500 Å for AGC 191770, by 0.8 at 3600 Å and 0.06 in the opposite direction at 6500 Å for AGC 228010, and by 0.8 at 3600 Å and 0.06 in the opposite direction at 6500 Å for AGC 253926. The light loss at 3600 Å is (respectively) 4.8%, 15.2%, and 0.6% greater than that at 6500 Å for AGC 191770, AGC 228010, and AGC 253926.
Obviously, the relative blue-light loss of AGC 228010 is the largest among the three galaxies. However, such an amount of light loss is calculated under the assumption that our targets were stars. Actually, as seen in Figure 2, they are galaxies with extended areas larger than the slit width. So, the lost blue light is partly gained back from the neighboring pixels, and the effect on relative flux calibration of blue-end spectra is not that significant, although it still exists. The probable light loss at blue-end wavelengths ( Å) would cause an underestimate of late A-type stars. Even with some blue-light loss, we could still obtain the same qualitative conclusion that these three LSBGs have a potential dearth of stars with ages Gyr because late A-type stars typically have ages older than Gyr.
6.1.3 Limitations Caused by Very Small Sample
Limited by the reliability of a spectral continuum, the number of spectra that we can use for a possibly reliable conclusion from full-spectrum fitting is only three. Undoubtedly, this is from a representative sample of the whole edge-on LSBG sample, and not even representative of our final 12 observing targets. Thus, conclusions from the three galaxies may not apply more generally to LSBGs. To get more universal conclusions on stellar populations of LSBGs, we will seek sufficient observing time in the future, obtaining stellar continuum spectra of good quality for a larger and more representative sample of edge-on LSBGs.
|AGC 191770||AGC 228010||AGC 253926|
|aa is the average value of and , which are both based on the diagnostic.||0.030||0.030||0.050||0.050||0.037||0.040|
|MFccMF refers to mass fraction to the total stellar mass.||0.0||0.0||0.0||0.0||0.04||0.1|
6.2 Stellar Mass
In Figure 6, the gas mass fractions computed from two different sets of stellar masses for our sample have obvious discrepancies. The gas mass fractions from the “-based” stellar masses are apparently higher than those from the “-based” stellar masses.
Which set of stellar mass is more appropriate for LSBGs? Here, we prefer the “-based” stellar masses because of the following reasons: (1) The luminosity at relatively longer wavelengths (the red or near-infrared bands) is more ideal for tracing the stellar mass because it is less affected by extinction and other factors that affect shorter-wavelength bands. (2) LSBGs have few blue, massive stars formed in the past 1 Gyr; the stellar mass is predominantly in low mass stars which radiate most of their light at red and near-infrared wavelengths. We checked the WISE satellite 3.4 m and 4.6 m images of our LSBGs and found that all of them are detected in both bands, especially in the bulge regions. Thus, the stellar masses of LSBGs would be underestimated from their -band luminosity. (3) The “-based” stellar masses of our LSBG sample are – M, which agrees with the conclusion of other studies that LSBGs are mostly comparable in mass with prominent spirals that define the Hubble sequence (e.g., McGaugh, 1994).
6.3 Chemical Evolution of LSBGs
In the closed-box model of the chemical evolution of galaxies, the initially metal-free gas is permanently locked into stars with neither inflow nor outflow and only evolves “in isolation” with time. Other assumptions include an invariant stellar initial mass function and instantaneous recycling. The relation between the metallicity () and the gas mass fraction () for a system that evolves in isolation is
where is the “yield,” the ratio of the mass of newly formed metals to the mass of gas permanently locked in stars. Here, is best expressed as an elemental abundance in logarithmic form. For oxygen, following Lee et al. (2003), this equation can be rewritten as
In a plot of oxygen abundance, 12 + log(O/H), against inverse gas mass fraction, conveyed by log log(), the closed-box model predicts a slope of unity. Deviations from this slope may be an indication for the inflow or outflow of gas. The oxygen yield, , can be derived from the intercept of the plot.
Following Lee et al. (2003), for the field dIrrs, there is an excellent correlation represented by a line with a slope of unity (green solid line in Figure 6) which is consistent with the closed-box model of the chemical evolution of galaxies. The standard deviation of the fit is 0.126 dex in log(O/H). From the intercept of the line, the effective oxygen yield by mass () for the field dIrrs from Lee et al. (2003) is 2.2210. Obviously, the LSBG data (purple open triangles) from de Naray (2004) are generally consistent with the field dIrr data (green open squares) and can also be fitted by the closed-box model fit (green solid line) of Lee et al. (2003), albeit with more scatter.
We cannot tell for sure whether our LSBG sample can be explained with the closed-box model by fitting; our sample is too small and also too narrow in the coverage of both gas mass fraction and gas-phase metallicity to show a true relation, even if some essential relations do indeed exist. However, to compare with the field dIrrs from Lee et al. (2003) and the LSBGs from de Naray (2004), we still forced our LSBG data (all filled circles with “-based” stellar mass) to be fitted by a line with a slope of unity (red solid line). The standard deviation of this fit is 0.27 dex in log(O/H). From this artificial fitting, the effective oxygen yield by mass for our LSBG sample is under the assumption of a closed-box model. The value of true yield is often considered to be the solar oxygen abundance by mass. For the Sun, the oxygen abundance is log(O/H) dex (Grevesse et al., 1996; Lee et al., 2003), which is by mass. Apparently, the effective oxygen yield for our LSBG sample is larger than that for the dIrr sample, and is 69% of the value for the true oxygen yield. The difference between the effective and true yields roughly represents the fraction of metal lost (up to 30%) up to now (Garnett, 2002).
Such a value of metal loss can be reasonable since it is often argued that LSBGs have shallow potential wells, making it easier for mass loss to occur. However, given the substantial uncertainties in both metallicity and gas mass fraction for our LSBG sample, this mass-loss value of 30% is only a rough estimate, and we cannot completely rule out a closed-box scenario. Additionally, judging from the offset between the two closed-box lines in Figure 6, we find evidence that our LSBGs retain up to times as much of their metals compared with dwarf irregulars (green squares), consistent with metal retention being related to galaxy mass. In any case, although the closed-box scenario cannot be completely ruled out for our LSBG sample, there is still substantial room for considering other, more complicated chemical evolution models with outflow or even infall (e.g., Edmunds, 1990; Chang et al., 2010) for our LSBG sample.
6.4 Are LSBGs a Fundamentally Different Type?
In the -based stellar masses (Table 4), our LSBG sample is estimated to be – M, showing no fundamental differences from prominent spiral galaxies. Furthermore, comparing with previously published luminosity-metallicity plots for LSBGs and other galaxy types (Figure 4), our LSBG sample agrees well with the LSBG sample of de Naray (2004) and also generally the sample of normal galaxies. In the mass-metallicity plot (Figure 5), our LSBG sample is also consistent with the LSBG sample of de Naray (2004) and with the relatively low abundance for their stellar mass, which is appropriate for normal galaxies. Therefore, from the behavior in the – and – plots, our LSBG sample does not seem to be fundamentally different from normal disk galaxies.
Regarding the plot of gas mass fraction versus oxygen abundance (Figure 6), we prefer the -based stellar masses (black or red filled circles). Our LSBG sample is not likely to be explained by simple closed-box models, although the galaxies evolve in relative isolation. Systematically biased to larger gas mass fractions, our LSBG sample has higher oxygen abundances for smaller gas mass fractions, which is also shown by other samples. Such qualitative agreement in gas mass fraction versus metallicity distribution makes the LSBGs seem consistent with other galaxy types.
We therefore suggest that LSBGs are probably not fundamentally different from normal galaxies, but are just the continuous extension to the lower surface-brightness end of normal disk galaxies.
We have obtained long-slit spectra of 12 edge-on LSBGs selected from the EGIS-40 sample. Observations were made with the Kast spectrograph on the 3-m Shane telescope at Lick Observatory for 3 galaxies and with the Double Spectrograph on the 5-m Hale telescope at Palomar Observatory for 9 galaxies. We adopted a median-filter method to estimate the sky background of objects observed at Lick Observatory and a second-order polynomial to fit the sky background of objects observed at Palomar Observatory, which can be comparable to or even exceed the galaxy surface brightnesses.
The stellar populations of three of the LSBGs were studied with the program STARLIGHT, as their spectral quality is sufficiently high. We find that our LSBGs are dominated by stellar populations with ages Gyr and have a potential dearth of stars with ages Gyr. This indicates that LSBGs probably have a very low level of star-formation activity over the past 1 Gyr, which may be the most likely explanation for their low surface brightness.
We have measured oxygen abundances (12 + log(O/H)) using both the and [N ii]/H diagnostics after performing reddening corrections for the observed fluxes of emission lines. The -based metallicities are 8.00–8.62, with a mean of 8.30 dex and a median of 8.26 dex. The [N ii]/H-based metallicities are 7.67–8.54, with a mean of 8.20 dex and a median of 8.26 dex. The two sets of oxygen abundances are consistent with each other. Comparing with the luminosity-metallicity relation established in previous work on LSBGs, normal spirals, and other galaxy types, our LSBGs tend to be consistent with normal galaxies.
We derive two sets of stellar masses for our LSBG sample. One set is the -based stellar masses which are – M, and the other set is the -based stellar masses which are – M, 0.5 order of magnitude lower. We argue that the -based stellar masses are better for our LSBG sample because such galaxies primarily consist of low mass stars which radiate their light mostly at red and infrared wavelengths. Stellar masses derived from the band might be underestimated. According to the -based stellar masses, our LSBGs are comparable to normal spiral galaxies in stellar mass, but they have considerably lower abundances.
In the plot of gas mass fraction versus oxygen abundance, we find that like other galaxy types, our LSBGs with relatively high “-based” stellar masses have lower metallicities for higher gas mass fractions. Regarding the chemical evolution of our sample, the LSBG data appear to allow for up to 30% metal loss, but we cannot completely rule out the closed-box model. Additionally, we find evidence that our galaxies retain up to about 3 times as much of their metals compared with dwarfs, consistent with metal retention being related to galaxy mass.
Combining the above conclusions, our LSBG sample is generally consistent with the luminosity-metallicity, mass-metallicity, and gas mass fraction-metallicity plots of normal spiral galaxies and other galaxy types (e.g., dIrrs, starbursts). Thus, we speculate that LSBGs are probably not fundamentally different from normal galaxies, but rather are just the continuous extension in surface brightness to the low end.
- Abazajian et al. (2009) Abazajian K. N., Adelman-McCarthy J. K., Agueros M. A. et al. 2009, ApJS, 182, 543
- Bell et al. (2003) Bell E. F., McIntosh D. H., Katz N., Weinberg M. D. 2003, ApJS, 149, 289
- Bizyaev et al. (2014) Bizyaev D. V., Kautsch S. J., Mosenkov A. V. et al. 2014, ApJ, 787, 24
- Bothun et al. (1987) Bothun G. D., Impey C. D., Malin D. F., Mould J. R. 1987, AJ, 94, 23
- Bothun et al. (1997) Bothun G. D. et al. 1997, PASP, 109, 745
- Brodie Huchra (1991) Brodie J. P., Huchra J. P. 1991, ApJ, 379, 157
- Bruzual Charlot (2003) Bruzual G., Charlot S. 2003, MNRAS, 344, 1000
- Burkholder et al. (2001) Burkholder V., Impey C., Sprayberry D. 2001, AJ, 122, 2318
- Calzetti et al. (2000) Calzetti D., Armus L., Bohlin R. C. et al. 2000, ApJ, 533, 682
- Cardelli et al. (1989) Cardelli J. A., Clayton G. C., Mathis J. S. 1989, ApJ, 345, 245
- Ceccarelli et al. (2012) Ceccarelli L., Herrera-Camus R., Lambas D. G., Galaz G., Padilla N. D. 2012, MNRAS, 426, 6
- Chang et al. (2010) Chang R. X., Hou J. L., Shen S. Y., Shu C. G. 2010, ApJ, 722, 380
- Cid Fernandes et al. (2005) Cid Fernandes R., Mateus A., Sodre L., Stasinska G., Gomes J. M. 2005, MNRAS, 358, 363
- Cole et al. (2001) Cole S., Norberg P., Baugh C. M. et al. 2001, MNRAS, 326, 255
- de Grijs (1998) de Grijs R. 1998, MNRAS, 299, 595
- de Blok (1998) de Blok W. J. G., McGaugh S. S. 1998, ApJ, 508, 132
- de Blok (1998) de Blok W. J. G., van der Hulst J. M. 1998, AA, 336, 49
- de Naray (2004) de Naray R., McGaugh S., de Blok W. 2004, MNRAS, 355, 887
- Das et al. (2009) Das M. et al. 2009, ApJ, 693, 1300
- Denicolo et al. (2002) Denicolo G., Terlevich R., Terlevich E., 2002, MNRAS, 330, 69
- Disney (1976) Disney M. J. 1976, Nature, 263, 573
- Du et al. (2015) Du W., Wu H., Lam M.-I. et al. 2015, AJ, 149, 199
- Edmunds (1990) Edmunds M. G. 1990, MNRAS, 246, 678
- Edmunds (1999) Edmunds M. G. 1999, in Davies J. I., Impey C., Phillips S., ed., ASP Conf. Ser. Vol.170, The Low Surface Brightness Universe. Astron. Soc. Pac., San Francisco, p. 383
- Faber (1973) Faber S. M. 1973, ApJ, 179, 731
- Filippenko (1982) Filippenko A. V. 1982, PASP, 94, 715
- Galaz et al. (2011) Galaz G. et al. 2011, ApJ, 728, 74
- Garnett et al. (1987) Garnett D. R., Shield G. A. 1987, ApJ, 317, 82
- Garnett (2002) Garnett R. D. 2002, ApJ, 581, 1019
- Grevesse et al. (1996) Grevesse N., Noels A., Sauval A. J. 1996, in ASP Conf. Ser. 99, Cosmic Abundances: Proceedings of the 6th Annual October Astrophysics Conference, ed. S. S. Holt G. Sonneborn (San Francisco: ASP), 117
- Haberzettl et al. (2007) Haberzettl L. et al. 2007, AA, 471, 787
- Haynes et al. (2011) Haynes M. P., Giovanelli R., Martin A. M. et al. 2011, AJ, 142, 170
- Impey Bothun (1997) Impey C. Bothun G. 1997, ARAA, 35, 267
- Impey et al. (2001) Impey C., Burkholder V., Sprayberry D. 2001, AJ, 122, 2341
- Kauffmann et al. (2003) Kauffmann G., Heckman T. M., White S. D. M. et al. 2003, MNRAS, 314, 33
- Kelson (2003) Kelson D. D. 2003, PASP, 155, 688
- Kewley Dopita (2002) Kewley L. J. & Dopita M. A. 2002, ApJS, 142, 35
- Kochanek et al. (2001) Kochanek C. S., Pahre M. A., Falco E. E. et al. 2001, ApJ, 560, 566
- Lam et al. (2015) Lam M. I., Wu H., Yang M. et al. 2015, MNRAS, 446, 4291
- Lamareille et al. (2004) Lamareille F., Mouhcine M., Contini T., Lewis I., Maddox S. 2004, MNRAS, 350, 396
- Lee et al. (2003) Lee H., McCall M. L.,Kingsburgh R. L., Ross R., Stevenson C. C. 2003, AJ, 125, 146
- Lequeux et al. (1979) Lequeux J., Peimbert M., Rayo J. F., Serrano A., Torres-Peimbert S. 1979, AA, 80, 155
- Liang et al. (2010) Liang Y. C. et al. 2010, MNRAS, 409, 213
- Matthews Gallagher (1997) Matthews L. D. Gallagher J. S. AJ, 114, 1899
- Matthews Wood (2001) Matthews L. D. Wood K. 2001, ApJ, 548, 150
- McGaugh (1991) McGaugh S. S., 1991, ApJ, 380, 140
- McGaugh (1994) McGaugh S. S. 1994, ApJ, 426, 135
- McGaugh (1996) McGaugh S. S. 1996, MNRAS, 280, 337
- McGaugh de Blok (1997) McGaugh S. S., de Blok W. J. G. 1997, ApJ, 481, 689
- McGaugh Bothun (1994) McGaugh S. S. Bothun G. D. 1994, AJ, 107, 530
- Melbourne Salzer (2002) Melbourne J., Salzer J. J. 2002, AJ,123, 2302
- Miller Stone (1993) Miller J. S. & Stone R. P. S. 1993, Lick Obs. Tech. Rep. 66 Santa Cruz: Lick Obs.
- Minchin et al. (2004) Minchin R. F. et al. 2004, MNRAS, 355, 1303
- Morelli et al. (2012) Morelli L., Corsini E. M., Pizzella A. et al., 2012, MNRAS, 423, 962
- Oke & Gunn (1982) Oke J. B. & Gunn J. E. 1982, PASP, 94, 586
- O’Neil et al. (2000) O’Neil K. et al. 2000, AJ, 119, 136
- Osterbrock (1989) Osterbrock D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (Mill Valley, CA: University Science Books)
- Pagel et al. (1979) Pagel B. E. J. et al. 1979, MNRAS, 189, 95
- Pettini Pagel (2004) Pettini M. Pagel B. E. J., 2004, MNRAS, 348, L95
- Schombert et al. (1990) Schombert J. M., Bothun G. D., Impey C. D., Mundy L. G. 1990, AJ, 100, 1523
- Schombert et al. (2001) Schombert J., McGaugh S. S., Eder J. A. 2001, AJ, 121, 2420
- Schombert et al. (1992) Schombert J. et al. 1992, AJ, 103, 1107
- Schombert et al. (2013) Schombert J. et al. 2013, AJ, 146, 41
- Skillman et al. (1989) Skillman E. D., Kennicutt R. C., Hodge P. W. 1989, ApJ, 347, 875
- Smith et al. (2002) Smith J. A., Tucker D. L., Kent S., Richmond M. W. et al. 2002, AJ, 123, 2121
- Trachternach et al. (2006) Trachternach C., Bomans D. J., Haberzettl L., Dettmar R. J. 2006, AA, 458, 341
- Tremonti et al. (2004) Tremonti C. A., Heckman T. M., Kauffmann G. et al. 2004, ApJ, 613, 898
- van Dokkum et al. (2001) van Dokkum P. G. 2001, PASP, 113, 1420
- van der Hulst et al. (1993) van der Hulst J. M. et al. 1993, AJ, 106, 548
- van der Kruit (1981) van der Kruit P. C. 1981, AA, 99, 298
- van der Kruit (1988) van der Kruit P. C. 1988, AA, 192, 117
- van der Kruit et al. (2001) van der Kruit P. C., Kregel M., Freeman C. K. 2001, AA, 379, 374
- Wen et al. (2013) Wen X. Q., Wu H., Zhu Y. N. et al. 2013, MNRAS, 433, 2946
- Wu et al. (2005) Wu H., Cao C., Hao C. N. et al. 2005, ApJ, 632, 79
- Zaritsky et al. (1994) Zaritsky D., Kennicutt Jr. R. C., Huchra J. P. 1994, ApJ, 420, 87
- Zhu et al. (2010) Zhu Y. N., Wu H., Li H. N., Cao C. 2010, RAA, 10, 329