Abstract

Disney et al. (2008) have found a striking correlation among global parameters of HI-selected galaxies and concluded that this is in conflict with the CDM model. Considering the importance of the issue, we reinvestigate the problem using the principal component analysis on a fivefold larger sample and additional near-infrared data. We use databases from the Arecibo Legacy Fast Arecibo -band Feed Array Survey for the gas properties, the Sloan Digital Sky Survey for the optical properties, and the Two Micron All Sky Survey for the near-infrared properties. We confirm that the parameters are indeed correlated where a single physical parameter can explain 83% of the variations. When color ( ) is included, the first component still dominates but it develops a second principal component. In addition, the near-infrared color ( ) shows an obvious second principal component that might provide evidence of the complex old star formation. Based on our data, we suggest that it is premature to pronounce the failure of the CDM model and it motivates more theoretical work.

1. Introduction

One way to understand our universe is to gain insights into the structure of galaxies. For one thing, it helps to reveal the role of dark matter in their formation and dynamics. The cosmological model that consists of a cosmological constant and the cold dark matter in addition to the ordinary baryon matter and radiation ( CDM) has been able to successfully explain the evolution of the cosmic structures especially at large scales. By measuring CMB fluctuations COBE [4, 5], WMAP [69], Type Ia supernovae [10], and gravitational lensing [11], this model of cosmology has by now been established as the standard model of cosmology. Among its various implications, it suggests a hierarchical, bottom-up history of structure formation that evolves from small fluctuations to galaxies, clusters, and eventually superclusters. On the other side, the success of CDM has not been as clearcut at small scales. There still exist several inconsistencies between CDM and the observations at small scales. For example, the simulations based on CDM have revealed more number of galactic satellites [1215] and less number of disk galaxies [16] than what have been observed. Besides, the degree of emptiness in the voids is also inconsistent between theory and observation [17]. While the CDM can explain the galaxy rotational curves at large radii [18], the relatively higher density at the galactic core than that predicted by the CDM, that is, the so-called cusp-core problem, is still unresolved [19].

The hierarchical galaxy formation scenario has been actively investigated previously (e.g., [2022]), and the tight correlation between selected galactic parameters has been studied in the past [2330]. Moreover, the overall correlation among all major galactic variables has also been thoroughly investigated recently by Disney et al. ([31]; hereafter D08) and Garcia-Appadoo et al. ([32]; hereafter G09). Based on 195 galaxies, the two studies found remarkably strong correlations between six galactic variables: the 90%-light radius ( ), the 50%-light radius ( ), the HI mass ( ), the dynamical mass ( ), the luminosity ( ), and the color. G09 found strong correlations among the galactic properties, including a common dynamical mass-to-light ratio within the optical radii, a correlation between surface brightness and luminosity, and a common HI surface density. D08 further showed that all the parameters can be included in one correlation and suggested that this is in conflict with the notion of the hierarchical structure formation scenario and the CDM model.

Considering the importance of this issue, we set out to reinvestigate the analysis made by D08 and G09. We note that their claim actually involves two separate subissues. First, there is the issue of whether the galactic dynamics is indeed controlled by a single parameter. Second, even if this is true, there still remains the issue of whether such a fact necessarily concludes the failure of the hierarchical structure formation scenario and the cold dark matter model. In our view the positive conclusion of the former does not necessarily lead to the affirmative conclusion of the latter. In this paper we demonstrate tentative evidence that there may exist more than one principal component among the global parameters of galaxies with regard to the first issue. As for the second issue, we argue about the importance of the nongravitational baryon physics in galactic structure formation, which renders the naive extrapolation of the hierarchical structure formation scenario from cosmic to galactic scales questionable.

In order to scrutinize the correction issue, we performed a similar analysis on global parameters of galaxies with a significantly larger database and two additional parameters, and , based on the infrared band. We include a total of 1022 galaxies from the Arecibo Legacy Fast Arecibo -band Feed Array Survey (ALFALFA; [3336]) for the atomic gas properties, and the Sloan Digital Sky Survey (SDSS; [37]) for the optical properties. Recently, Toribio et al. published two papers ([38, 39]; hereafter T11) with these two surveys. In T11, they assembled three samples to analyse the data. They argued that HI emission provides the most reliable way to determine the morphologies and conclude that color and HImass of gas-rich galaxies cannot be very closed. In addition, they also found the correlation between mass, radius, luminosity and ( - ) by principal component analysis (PCA; [40]) for one of their subsamples. Among our 1022 galaxies, 479 of them have also been detected by the Two Micron All Sky Survey (2MASS; [41]). We use these galaxies to study their near-infrared properties, which contain more stellar information.

Near-infrared studies of HI-selected galaxies had been attempted by various groups (e.g., G09; [42, 43]). Our motivation of adding the near-infrared data in our analysis is the following. The optical emission is sensitive to young stars. The near-infrared emission, on the other hand, is less affected by the young stars and is therefore a better tracer of the total stellar mass, which dominates the baryonic matter at galactic scales. We believe that the inclusion of the infrared data would provide us additional and independent information on the baryonic mass assembly history of the galaxies.

Through the PCA, we confirm that except the color, all other observables, from HI, optical to near-infrared bands, are highly correlated and dominated by a single parameter. This is true both in the optical and in the near-infrared bands and this confirms the results in D08. In addition, we also see a second component from color, especially in our near-infrared data. Because near-infrared color provides the information of integrated star formation history, it may be an evidence for complex formation history, whereas a valid structure formation theory needs to explain this observation.

The organization of this paper is as follows. We describe the data and sources in Section 2. Then several variables are adopted and applied to statistical analysis in Section 3. In Section 4, we summarize and discuss our results.

2. Data

Our samples are the blind 21 cm survey from ALFALFA. This selects gas rich galaxies which also contain low luminosity and low surface brightness galaxies in higher proportion than those in an optical selection. The optical data for this study come from the SDSS DR7, which covers 12,000 deg2 for imaging and provides spectra of 930,000 galaxies. Here we briefly describe how we select SDSS counterparts to the ALFALFA sources, and we refer to G09 for more detailed discussion about identification.

The Arecibo Telescope has a beam size of at 21 cm. Since the majority of the HI detections have S/N [3436], we adopted a conservative searching radius of . We found 1233 SDSS galaxies that appear to be detected by ALFALFA. We then excluded Virgo galaxies, because their neutral hydrogen is known to undergo strong environmental impact (e.g., [44]). We also excluded galaxies whose half-light radii are too small ( ) comparing to the SDSS resolution. These small half-light radii either arise from misidentifications (from stars) or would result in large uncertainties. We are left with a large sample of 1022 SDSS nearby galaxies, which have distances smaller than 254  . Among these galaxies, 889 have magnitudes that are <18 and 120 have –20. The cumulative number counts of SDSS galaxies (e.g., [45]) are ~60 at and ~450 at –20. Given our search radius, we therefore expect at most 1.3 misidentifications in our 889 galaxies and additional 1.3 misidentifications in our 120 galaxies. These number are sufficiently small and misidentified galaxies should not impact our analyses.

The total stellar masses of galaxies are more directly reflected by the near-infrared observations. We added the data from 2MASS to our samples. 2MASS provides , , and -band observations of the entire sky as well as a point source catalog and an extended source catalog. We use the 2MASS All-Sky Extended Source Catalog (XSC) to find the galaxies in the HI samples. To understand the quality of the identification, we first compared the 2MASS and SDSS coordinates of the 2MASS-detected galaxies. We found that more than 90% of them have offsets between 2MASS and SDSS that are well within their half-light radii. We visually inspected all galaxies with large offsets of and found small number of cases that are likely misidentifications as well as ongoing mergers. We excluded these galaxies from our samples. Because 2MASS is shallower than SDSS, we are left with 481 reliably identified galaxies in the near-infrared. To match the same aperture of the color for SDSS and 2MASS catalog, 479 galaxies are left in our near-infrared subsamples.

From the ALFALFA-released catalog 1, 2, and 3, we obtained 1796 HI data, out of which 1265 galaxies could be found in the SDSS DR7 database. There are 32 galaxies within this set that are too faint in the optical to have reliable magnitudes and luminosities. Hence we finally used the remaining 1022 galaxies in our analyses. We have also analyzed the 195 galaxies of D08 and G09 from HIPASS [4649] and have obtained similar results. However, since the definitions of the observational variables, such as those of the rotational velocity, are not entirely consistent between the two data sets, we only report on the results from ALFALFA. These 1022 galaxies can be regarded as a blind HI-selected sample. We deduced from the data six variables, which are (half-light radius in units of pc), ( -light radius in units of pc), (luminosity in band in solar units), (HI mass in solar units), (dynamical mass in solar units), and color ( - ).

The variables and represent the radii in the Petrosian system [45, 50, 51]. In SDSS, the parameters are petroR50 and petroR90, respectively. Because the Petrosian system is based on circular objects, we corrected the radii with the major-to-minor axis ratios, which are the parameters deVAB_r or expAB_r in SDSS. To do this, we follow the result in [52] and [48, 49]. The authors fitted the corrections from circular to elliptical apertures as functions of major-to-minor axis ratios. We directly adopted their formulas for our corrections. By comparing the likelihoods of the de Vaucouleur and the exponential models, we chose the one with the larger likelihood between deVAB_r and expAB_r. was derived from the Petrosian system and calculated from the Petrosian magnitude, petroMag_r. In order to have the same aperture for 2MASS, we use Petrosian color, which is from petroMag_g and petroMag_i. The variable is calculated by , in which is the rotational velocity from the HI spectra and corrected for the inclination with the major-to-minor axis ratio as what we did for and . is acquired directly from the ALFALFA database and it is derived from the HI flux. This estimation is based on the assumption that the masses from optical radius are proportional to dark matter halos.

To make sure that the masses from the HI measurement can describe the dynamical mass, we compare them with the masses based on the stellar velocity dispersion, which is , as in [13]. Here, is the velocity dispersion, is the effective radius, , and is the Sérsic index. In SDSS, the velocity dispersion can be calculated by [53], where , is from the best fitting circularized Sérsic profile, and is the SDSS measurement within the fibers. The Sérsic data are from NYU Value-Added Galaxy catalog (VAGC [54]). After matching, we have 320 of our 1022 main galaxy sample and 164 of our 479 2MASS subsample.

Figure 1 is the comparison between the HI and the stellar dynamical masses. There is an apparent sequence, indicating that the two masses trace similar dynamics for most of the galaxies. On the other hand, the HI dynamical masses are on-average larger than the stellar dynamical masses. This ratio is indicated by the solid line in Figure 1. This is not too surprising, since neutral gas (and ) can trace the dynamical mass to a larger distance in a dark matter halo. In our subsequent analyses, we adopted HI dynamical masses. This is consistent with the work of D08/G09 and provides a larger sample here.

Because there is high degeneracy between , , and bands, we only chose band to represent the near-infrared data. Therefore, we acquired (half-light radius in band in units of pc) and (luminosity in band in solar units) of the galaxies that overlap in the ALFALFA, SDSS, and 2MASS catalogs to gain a better insight into the stars of the galaxies. More specifically, is derived from the magnitude in band, which is the parameter j_m_ext in the 2MASS database. This magnitude is based on the extrapolation of the fit to the surface brightness profile. And is the integrated half-flux radius of band, which is the parameter j_r_eff in the 2MASS database. We adopted - for the color of the near-infrared subsample. The aperture of petroMag_i in SDSS is twice of the Petrosian radius, petroRad_i. Thus we interpolate the -band magnitude in 2MASS at the same aperture as SDSS. In our 481 2MASS subsamples, 479 of them have sufficient data at different apertures to interpolate the 2MASS magnitude. Therefore our final near-infrared subsample contains 479 galaxies.

3. Methods and Results

Our sample of 1022 galaxies is not only larger than that in D08 and G09 but also covers broader ranges of size, luminosity, and mass (Figure 2). Despite that our sample is still dominated by galaxies, the minimum value of in our sample is much smaller than that of D08 and G09. As for the HI and the dynamical masses, although the median values of the two samples are similar, our sample contains a substantial amount of lower mass galaxies. Our sample also covers a broader range of the - color. Because of the larger sampling and the wider range of the galactic properties, our sample is in general more representative. However since 2MASS is shallower, our 2MASS subsample of 479 galaxies is relatively speaking less representative than that from SDSS. Even so, our 2MASS subsample is still substantially larger than that of D08 and G09.

It is important to investigate whether our HI-selected sample is biased against early type galaxies, since these galaxies are usually gas poor. To do so, we identify spheroidal and disk galaxies in our sample based on the morphology with a method similar to that in [55]. In the SDSS database, there are de Vaucouleur and exponential models for each galaxy. By comparing the likelihood and the fractions of the two models for our 1022 galaxies, we found that 804 galaxies are disk like and 218 galaxies are spheroidal like.

In Figure 3, we show a color-luminosity diagram for our 1022 galaxies. The spheroidal galaxies are in general more luminous and redder than the disk galaxies. This is consistent with what we expect for elliptical and spiral galaxies. Most importantly, in the color-luminosity space, the spheroidal galaxies are redder than the blue cloud although they do not yet form a complete red sequence. Our sample thus appears to include a fair number of red and elliptical galaxies. Although the bias against extremely gas-poor galaxies can be hardly avoided here, fortunately, we found no major difference between these two types in our subsequent studies. We therefore believe that the omission of extremely gas-poor galaxies should not have caused major systematic bias in our analysis.

Relations among the key parameters can be inferred from the 1022 galaxies in both ALFALFA and SDSS. For instance, it is found that the half-light radius is proportional to the 90%-light radius ( ; [56]); the -band luminosity is proportional to the cubic power of the half-light radius ( ; [57]); the HI mass is proportional to the square of the half-light radius ( ; [43, 58]); finally, the dynamical mass is proportional to the -band luminosity ( ; [25]). We found that all the correlations are evident even after including the near-infrared variables, except the color (Figures 4 and 5).

As a whole, the correlations between color and other variables are much weaker than other correlations. We tested various combinations of colors and found that - gives a larger PCA correlation coefficient than other colors (e.g., - , adopted by D08). The reason could be the larger wavelength difference between and . Among the three 2MASS bands, the result based on is somewhat better in the PCA, possibly because of the better signal-to-noise ratio. Hence, we adopted - for the color when we analyzed the 479 galaxies in 2MASS. Nevertheless, our analyses show that all the correlation coefficients are smaller than 0.7, indicating that they are not so highly correlated with other parameters. Intuitionally, more luminous galaxies tend to be redder because their colors are dominated by older stars. In fact, the color is more complex than any other variable because of the bias introduced by the very luminous young stars.

We conducted PCA to find correlations among the variables. PCA typically produces a series of new variables called the principal components, namely, PC1, PC2, and so forth. The correlations between these principal components and the original variables then reveal the general correlations between the particular variable and the others. In our case we found that the first principal component, PC1, is highly correlated with the six observational variables. We notice that while the color is less correlated to PC1, possibly because of the recent star formation, it is much more correlated with the second principal component, PC2. In addition to the investigation into the diagram of correlations, the eigenvalues of the correlation matrices of the original variables give quantitative information of the degree of correlations. For the 1022 galaxies based on SDSS, the eigenvalues of PC1 through PC6 are 4.29, 0.92, 0.39, 0.20, 0.18, and 0.02 (Figure 6), respectively, where the maximum possible value is 6. Based on common PCA criteria, eigenvalues larger than 1 are considered significant. We therefore only plot PC1 and PC2. Next we conducted PCA without color, and we found the eigenvalues of PC1 through PC5 to be 4.15, 0.41, 0.22, 0.20, and 0.02 (Figure 7), respectively, where the maximum possible value is 5.

The aforementioned observations confirm the finding of D08. All the observed parameters are tightly correlated with PC1. The eigenvalue of PC1 indicates that it can explain 83% of the variance in the data (when color is not included). Color itself forms a second principal component. This might be explained by the fact that the optical color tends to be strongly affected by recent star formation activities and thus carries extra information that is unrelated to the global formation history of the galaxies. We will come back to the issue of color in Section 4. D08 claimed that the strong dominance of PC1 implies a single physical parameter to govern the structure of galaxies. However, this may be an oversimplified view of galaxies due to limited observational parameters. To test this, we extend the PCA to the near-infrared.

We included the 2MASS -band radius and the luminosity to investigate the role of stars, which dominate the baryonic matter in galaxies. For the 479 galaxies detected by 2MASS, we found the eigenvalues of PC1 through PC8 to be 5.40, 1.01, 0.61, 0.45, 0.32, 0.16, 0.03, and 0.02 (Figure 8), respectively, where the maximum possible value is 8. Conducting PCA without color again, we found that the eigenvalues of PC1 through PC7 are 5.35, 0.62, 0.45, 0.31, 0.19, 0.05, and 0.02 (Figure 9), respectively, where the maximum possible value is 7.

The overall trends in the previous near-infrared PCA are similar to those in the optical PCA. When color is not included, PC1 dominates and can explain 76% of the variance in the data. The importance of PC2 slightly increases from 8% in the optical case to 9% in the near-infrared case here. When color is included, it forms another principal component by itself. These again confirm the observations of D08 that only one physical parameter governs the dynamics of galaxies.

A subtle but surprising difference between the optical and near-infrared PCAs is the behavior of color. In the optical PCA in Figure 6, although the - color forms a second principal component, it still weakly correlates with other parameters and it is part of the first principal component. On the other hand, in the near-infrared PCA in Figure 8, the - color almost does not involve in the first principal component and itself forms a second component that is almost independent of other parameters. A potential issue here is the combination of SDSS and 2MASS . It is the reason we choose consistent apertures in SDSS and 2MASS to calculate the color. To test whether other minor differences between the two photometric system hamper the correlation, we replaced the - color with the pure 2MASS - color, and we found consistent results. In addition, the median errors in and for the 1022 SDSS galaxies are 0.015 and 0.018, respectively, and the median errors in and for the 479 2MASS galaxies are 0.011 and 0.075, respectively. These translate to typical color errors of 0.023 in - and 0.076 in - . Both values are significantly smaller than the color dynamical ranges shown in Figures 4 and 5, meaning that the distribution of the data is not dominated by measurement errors. We therefore believe that the lack of correlation between the near-infrared color and other parameters is real. We will discuss more on this result in Section 4.

4. Discussion and Conclusion

For the 195 galaxies analyzed in D08 and G09, the correlations among the parameters , , , , and are obvious. By selecting significant 5 times larger, 1022 overlapping samples from SDSS and ALFALFA, we also performed the PCA and confirmed through the high eigenvalues of the correlation matrices that the correlations are similarly strong. It follows that the radius, luminosity, HI mass, and dynamical mass of those chosen galaxies are tightly correlated. D08 shows that the basic parameters of galaxies are highly correlated and that there exists only one dominant principal component. From our studies, we know that this is true in both optical and near-infrared bands. Based on the hierarchical galaxy formation assumption, the dark matter halo formation has been well studied through the merger tree process (e.g., [20, 21, 5964]). D08 believes that this scenario may not be consistent with the simple relation between all basic galactic parameters because the processes of merger would break the original galactic structures. Although our analysis has confirmed D08’s finding, we do not find it compelling or sufficient reason to reject CDM. On the other hand, the dynamical mass-to-light ratio enclosed by the optical radius that strongly anticorrelates with surface brightness has been studied [6567]. It is not shown in our data. It may be caused by the selection effects from galaxies and insufficient dynamic range in surface brightness with SDSS and 2MASS data. Nevertheless, the strong correlations between parameter in D09 and our gas-selected samples are obvious in both optical and near-infrared data.

However, the color appears to be much less correlated, as T11 mentioned. This indeed complicates the situation and may be explained by a more sophisticated theory that would include, for example, the influence of recent star formation and very luminous young stars. Indeed these studies have already been pursued by many authors (e.g., [55, 59, 6871]. D08 suggests that the optical color ( - in their case) consists of two components: the systematic component that correlates with other parameters (and therefore involves in the first principal component), and the rogue component that is more or less random (and forms the second principal component). It is tempting to assume that the systematic component comes from the established stellar populations and is related to the global formation history of the galaxies, while the random component is related to the ongoing star formation activities and can be short-lived events in the formation history. Surprisingly, our near-infrared analyses suggest a different story. In our near-infrared PCA, the color is even more uncorrelated with other parameters, comparing to the optical color.

One may believe that optical colors like - can be more strongly affected by recent star formation, whereas the - color may be a better tracer of integrated star formation history because it is less affected by ongoing star formation and dust extinction. If this is correct, the result that - is more related to other dynamical properties suggests that the recent star formation history is more controlled by the current dynamical structure galaxies. On the other hand, our result that - is unrelated to other dynamical properties suggests that the formation history of the old stellar components is a more chaotic process. In other words, the stronger second component may be an indication of complex formation and evolution history of galaxies. If this is the correct interpretation of our data, then the evidence of complex merger history predicted by the CDM hierarchical formation model had been hidden in the near-infrared colors and was not revealed by D08. However, we believe that this has to be tested by more observations and simulations, and it is premature to suggest the failure of the CDM model.

To sum up, we believe that the 5 times larger sample size strengthens the claim in D08 that there is one dominant parameter in galactic structure. In addition, our near-infrared analyses also provide additional insight into the color of galaxies. The tight correlations and the uncorrelated color part between the parameters in the optical and near-infrared data provide potentially powerful observational constraints on the hierarchical structure formation theory and any other cosmology models.

Acknowledgments

This work is supported by Taiwan National Science Council under Project no. NSC97-2112-M-002-026-MY3 (P. Chen and Y.-Y. Chang), NSC98-2112-M-001-003-MY2 and NSC99-2112-M-001-012-MY3 (W.- H. Wang and Y.-Y. Chang), and by US Department of Energy under Contract no. DE-AC03- 76SF00515. The authors also thank the support of the National Center for Theoretical Sciences of Taiwan. They are grateful to L. Lin, H. Hirashita, M. Disney, P. Kroupa, and H.-W. Rix for useful comments and discussion. The Arecibo Observatory is part of the National Astronomy and Ionosphere Center, which is operated by Cornell University under a cooperative agreement with the National Science Foundation. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US 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 website 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. This paper makes use of data products from the Two Micron All Sky Survey, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation.