A&A 490, 571-581 (2008)
DOI: 10.1051/0004-6361:200810604

The relationship between gas content and star formation rate in spiral galaxies. Comparing the local field with the Virgo cluster

M. Fumagalli1,2 - G. Gavazzi1

1 - Istituto di Fisica G. Occhialini, Universita' di Milano-Bicocca, Piazza della scienza 3, Milano, Italy
2 - Department of Astronomy and Astrophysics, University of California, Santa Cruz, CA 95064, USA

Received 15 July 2008 / Accepted 1 August 2008

Context. Despite many studies of star formation in spiral galaxies, a complete and coherent understanding of the physical processes that regulate the birth of stars has not yet been achieved, nor has unanymous consensus been reached, despite the many attempts, on the effects of the environment on the star formation in galaxy members of rich clusters.
Aims. We focus on the local and global Schmidt law and we investigate how cluster galaxies have their star formation activity perturbed.
Methods. We collect multifrequency imaging for a sample of spiral galaxies, members of the Virgo cluster and of the local field; we compute the surface density profiles for the young and for the bulk of the stellar components, for the molecular and for the atomic gas.
Results. Our analysis shows that the bulk of the star formation correlates with the molecular gas, but the atomic gas is important or even crucial in supporting the star formation activity in the outer part of the disks. Moreover, we show that cluster members that suffer from a moderate HI removal have their molecular component and their SFR quenched, while highly perturbed galaxies show an additional truncation in their star forming disks.
Conclusions. Our results are consistent with a model in which the atomic hydrogen is the fundamental fuel for the star formation, either directly or indirectly through the molecular phase; therefore galaxies whose HI reservoirs have been depleted suffer from starvation or even from truncation of their star formation activity.

Key words: Galaxy: evolution - galaxies: ISM - galaxies: spiral - galaxies: clusters: individual: Virgo - stars: formation

1 Introduction

A satisfactory understanding of the physical processes governing the formation of giant molecular clouds (GMC) from primeval atomic hydrogen (HI) and their instability leading to the formation of stars is still far from achieved. This is not surprising due to the complex mix of hydrodynamical and gravitational physics involved in the process of star formation. Beside the theoretical difficulties, what surprisingly limits our present understanding of this issue is the lack of a robust phenomenology of star formation in nearby galaxies. Despite the many recipes that have been proposed[*], the simplest and most commonly used parametrization is the Schmidt law (Schmidt 1959) that establishes a correlation between the star formation rate density  $\rho_{\rm sfr}$ and the gas density  $\rho_{\rm gas}$:

\rho_{\rm sfr}=a\rho_{\rm gas}^\alpha;
\end{displaymath} (1)

derived for our Galaxy, the validity of this law has been tested also for external galaxies (e.g. Kennicutt 1983), for which a similar relation, usually written in terms of surface densities, holds:

\Sigma_{\rm sfr}=A\Sigma_{\rm gas}^n.
\end{displaymath} (2)

Many studies of the global Schmidt law, i.e. the relationship between the disk-averaged star formation rate (SFR) and the total gas content, have been carried out on large samples of galaxies. A milestone paper on this issue is by Kennicutt (1998), who studied the connection between the integrated SFR surface density and the integrated gas surface density in a sample of normal and starburst galaxies. Using integrated quantities, he found that the same Schmidt law with an index n=1.4 holds for both normal and starburst galaxies, covering a range of 5 dex in the gas density and over 6 dex in the SFR; he also found that this relation is mainly driven by the correlation between the SFR and the HI, while there is a poorer correlation between the SFR and the H2, significant only when massive ( $L_B>10^{10}~L_\odot$) objects are considered. Owing to the fact that star formation takes place in molecular regions on parsec scales, this result was somewhat unexpected, but Kennicutt argued that the poor knowledge of the CO-to-H2 conversion factor might be responsible for the scatter in the H2/SFR correlation. Using seven galaxies from the BIMA survey (Helfer et al. 2003), Wong & Blitz (2002) studied the applicability of the Schmidt law on local scales, i.e. through a comparison between the SFR and the gas surface density profiles. They found that a Schmidt law between the SFR and the gas content holds also on local scales and, even if the molecular gas content alone does not always control the SFR, they concluded that the correlation found for the total gas is entirely driven by the molecular component. Consistent results were found by Boissier et al. (2003), who used low resolution CO data but included an X factor that varies as a function of the metallicity. Another open issue is to asses the influence of the environment on the star formation. Galaxies in clusters suffer from a variety of environmental perturbations, mainly due to gravitational interactions between galaxies themselves or between a galaxy and the cluster potential well (see the review by Boselli & Gavazzi 2006) or to various hydrodynamical interactions between the ISM and the intergalactic medium (IGM). Tidal interactions between galaxies or between individual galaxies and the cluster potential (Byrd & Valtonen 1990; Merritt 1984) produce both the removal of the outer and looser components and the gas infall towards the galaxy center; harassment (Moore et al. 1996) can induce sinking of gas towards the galaxy center and can shape the stellar profiles; the ram pressure stripping (Gunn & Gott 1972) or the viscous stripping (Nulsen 1982) can effectively remove the gas components, mainly from the outer part of the disks. Furthermore, a combination of gravitational and hydrodynamic effects can cause galaxy starvation or strangulation (Bekki et al. 2002; Larson et al. 1980). According to starvation, the gas that feeds the star formation in the local universe probably comes from the infall of an extended gas reservoir, therefore the effect of removing the outer galaxy halo would be that of preventing further gas infall; as a consequence, the star formation exhausts the available gas, quenching further activity. The literature about the classification and the effects produced by the environment on the evolution of cluster galaxies is very rich (e.g. Boselli & Gavazzi 2006), but, despite many studies on nearby and distant clusters, it is not clear what the dominant processes are, nor how the star formation activity is affected (Koopmann & Kenney 2004, and reference therein). Unanimous consent seems only to hold on the statement that environmental processes are effective at removing the outer ISM and therefore affecting the HI component (Giovanelli & Haynes 1985; Cayatte et al. 1994), leaving the molecular component, well-bounded inside the galaxy potential well, undisturbed (Boselli et al. 2002b; Kenney & Young 1989). In order to study the environmental effects on star formation, we collect both objects that are members of the Virgo cluster and local field galaxies. Since the different components involved in the star formation are observable along a broad stretch of the electromagnetic spectrum, we adopt a multifrequency analysis. To quantify the atomic hydrogen we collect observations at 21 cm; the molecular hydrogen content is estimated indirectly via CO observations at 2 mm; we quantify the stellar mass with observations at the near-infrared (NIR) or visible bands and we study the presence of new stars indirectly, observing the hydrogen recombination line (H${\alpha }$) at 6563 Å. The sample is presented in Sect. 2, together with the description of the data reduction procedures; the analysis is given in Sect. 3, while we discuss our results and conclusions in Sects. 4 and 5.

2 Data

2.1 The sample

The sample of galaxies chosen for the present analysis is selected primarily according to the availability of high sensitivity, extensive CO mapping, such as in the Nobeyama CO atlas of nearby spiral galaxies by Kuno et al. (2007). Secondly we require that reliable H${\alpha }$ imaging exists in the literature. 28 giant spiral galaxies match these criteria. For all of them we were able to find additional imaging material in the red continuum (i or H band) and for 25/28 we found sensitive HI mapping. Table 1 lists some observational parameters and references for the selected galaxies. The sample is well balanced between isolated objects and members of the Virgo cluster, as stressed by the local density parameter (see Col. 11). One characteristic of our sample is that it is composed of giant spirals with a typical H-band luminosity of $10^{10}{-}10^{11}~L_\odot$. No BCDs nor Irr/dIrr are included. The fact that galaxies in the Nobeyama atlas have not been selected according to a CO emission criterion (Kuno et al. 2007, Sect. 2) protects against a possible bias towards CO rich galaxies. One could question why we preferred the Nobeyama atlas to higher resolution interferometric CO material, such as the BIMA Survey of Nearby Galaxies (Helfer et al. 2003) that provides images for 44 nearby spiral galaxies. This is mainly to ensure that our flux measurements are unaffected by severe missing flux problems on scales larger than 20-30'' at 115 GHz. Even though spectra from the central galaxy regions have been used to recover the missing flux or single dish observations provide the zero-spacing data, the algorithms used to recover the extended flux are somewhat unreliable (see Helfer et al. 2003, Appendix A).

Table 1: Observational parameters. Columns 1, 2: NGC, VCC or UGC designations. Columns 3, 4: coordinates (J2000) from GOLDMine and HyperLeda. Column 5: recessional velocities in km s-1 from GOLDMine and HyperLeda. Column 6: distances in Mpc, assigned following the subcluster membership criteria of Gavazzi et al. (1999) or from the Nearby Galaxies Catalogue by Tully (1994). Column 7: galaxy inclinations in deg from GOLDMine and HyperLeda. Column 8: morphological types from GOLDMine and HyperLeda. Column 9: H-band luminosity in $L_\odot $ from GOLDMine or computed from H magnitude from SIMBAD. Column 10: radius in arcsec at the B-band 25th isophote from HyperLeda. Column 11: density in Mpc-3 for galaxies brighter than -16 mag in the vicinity of the galaxy considered, determined on a 3D-grid at 0.5 Mpc spacing (Tully 1994). Column 12: reference to H${\alpha }$ imaging: a = Boselli & Gavazzi (2002); b = Koopmann et al. (2001); c = Boselli et al. (2002a); d = Marcum et al. (2001); e = Knapen et al. (2004); f = Meurer et al. (2006); g = Kennicutt et al. (2003); h = James et al. (2004). Column 13: reference to metallicity: i = Skillman et al. (1996); j = GOLDMine; k = Kennicutt et al. (2003); l = van Zee et al. (1998); m = Zaritsky et al. (1994). Column 14: reference to HI mapping: n = THINGS; o = VIVA; p = Cayatte et al. (1990); q = Verheijen & Sancisi (2001). Column 15: reference to NIR/Visual imaging: r = SDSS; s = GOLDMine; t = Kuchinski et al. (2000); u = Mulchaey et al. (1997).

2.2 H${\alpha }$ data

We compute the star formation rate from the intensity of the H${\alpha }$ emission line, following the method of Kennicutt (1998,1983). H${\alpha }$ images for the VCC galaxies are retrieved from GOLDMine (Gavazzi et al. 2003), while for the remaining galaxies we retrieved images from NED[*]. Images from GOLDMine are calibrated following Boselli & Gavazzi (2002) and Boselli et al. (2002a), i.e. using the photometric zero-point included in the image headers. For the images collected from NED instead, we derive the photometric calibration using published fluxes (Young et al. 1996; Romanishin 1990; Kennicutt & Kent 1983)[*]. Before converting the H${\alpha }$ flux into a star formation rate, we correct for the [NII] emission; for almost all the galaxies we apply a correction derived from drift-scan spectra that accounts for the mean [NII] emission from the entire disk. Where individual correction factors are not available, we apply the standard value $\rm [NII]/H\alpha = 0.53$ (Kennicutt 1992). Since we consider [NII]/H${\alpha }$ factors averaged over the entire disk, the correction underestimates the [NII] emission from the bulges of galaxies that host active galactic nuclei (AGN), where nuclear [NII] can be even stronger than H${\alpha }$. This is particularly severe in our sample composed of giant spirals, since almost all massive galaxies are known to host an AGN (Decarli et al. 2007). In addition, even if the nuclear emission was corrected by the suitable factor, it would not be possible to disentangle the H${\alpha }$ flux due to star formation from the one due to the AGN ionization. Therefore, by inspecting both the H${\alpha }$ imaging and the available spectra, we identify a nuclear region for each galaxy, and in the following analysis we do not include the emission from these regions. As for the dust absorption, since the H${\alpha }$/H$\beta$ or the H${\alpha }$/Br$_\gamma$ ratios are unknown for the majority of our galaxies, we prefer to apply a standard dust extinction A=1.1 mag (Kennicutt 1983) to other more sophisticated, yet unreliable corrections. From the corrected H${\alpha }$ luminosity we compute the star formation rate using the conversion (Kennicutt 1998)

{\it SFR}\; [M_\odot~{\rm yr}^{-1}]=7.94 \times 10^{-42}~L({\rm H}\alpha)~[{\rm erg/s}],
\end{displaymath} (3)

based on a Salpeter IMF ( $\alpha=2.35$) over $0.1{-}100~M_{\odot}$. Furthermore we assume that a fraction f=0.57 of Lyman continuum photons ionizes the hydrogen atoms inside the star forming regions, where dust is mixed with gas (Hirashita et al. 2003). On the contrary, we neglect the correction for the escape fraction of continuum Lyman photons, since it has been shown that it is less than 6% in starburst galaxies and even negligible for normal galaxies (Heckman et al. 2001). We estimate that the typical flux uncertainty is 0.05 mag for VCC galaxies; 0.1 mag for fluxes from Kennicutt & Kent (1983), 5% for those from Young et al. (1996) and up to 20% for Romanishin (1990). These errors produce a final uncertainty on the star formation rate which ranges from 0.04 dex to 0.09 dex, with a mean value of 0.05 dex. An additional uncertainty of 0.2 dex affects all data due to the uncertainty in the dust correction.

2.3 H2 data

In the absence of a direct tracer of the H2 molecule, the molecular hydrogen abundance can be inferred from CO observations. We collect flux calibrated maps of the $^{12}{\rm CO}(J{:}1{-}0)$ transition at 115 GHz from the Nobeyama CO atlas of nearby spiral galaxies taken with the 45-m radiotelescope by Kuno et al. (2007). Despite their spatial resolution of 15 arcsec (at 115 GHz), these observations have been preferred to higher resolution interferometric data because of their higher sensitivity to the large-scale emission from galaxy disks. In fact the typical rms noise of 40-100 mK combined with a velocity resolution of 5 km s-1, gives a sensitivity of $N({\rm H}_2)=5{-}10$ $\times $ $10^{19}~{\rm cm}^{-1}$. The most challenging task related with the CO observations, both from a theoretical and an observational point of view, is the determination of the CO-to-H2 conversion factor (X in $\rm cm^{-2}/(K~km~s^{-1})$). Many attempts to determine a reliable CO-to-H2 conversion factor can be found in the literature; nevertheless, there is not yet unanimous consent on the dependence of X on parameters such as density, temperature or even the geometry of the clumps that constitute the molecular regions (Wall 2007). Due to these difficulties, a constant conversion factor X, as derived for our Galaxy, has been often assumed, probably underestimating the H2 in low luminosity galaxies. In addition, it is expected that the CO-to-H2 conversion strongly depends on dust and metallicity that regulate the H2 formation and photodissociation or the CO heating. Therefore, we avoid a single conversion factor for our sample, but we choose individual X as a function of the metallicity, following the empirical calibration by Boselli et al. (2002b), made on fourteen independent measurements of CO and molecular hydrogen

\log X=-1.01\pm0.14 \times \left(12+\log({\rm O/H})\right) + 29.28\pm0.20.
\end{displaymath} (4)

With this equation we compute the conversion factor for each galaxies (See Col. 13 in Table 1 for references to metallicity data). For seven galaxies without published $12+\log({\rm O/H})$ we derive the metallicity using (Boselli et al. 2002b)

12+\log({\rm O/H})=0.37\pm0.04 \times \log L_{\rm H}~[L_\odot]+4.98\pm 0.15,
\end{displaymath} (5)

calibrated on a sample of nearby and Virgo cluster galaxies. The mean discrepancy between the metallicity directly measured from spectra and derived using (5) is 0.2 dex in our sample, not much greater than the uncertainties on the measurements themselves and significantly lower than the statistical uncertainty on the X conversion factor derived using (4), that is on average 1.3 dex. One should keep in mind that other unknown dependencies, beside the one on metallicity, would increase the uncertainty on X. Since spiral galaxies exhibit a radial metallicity gradient, when available in the literature, we account also for them to compute X as a function of the galaxy radius. We find that X spans 0.4 dex on average in individual galaxies and in some cases up to 1 dex. With the X conversion factor, we compute the molecular hydrogen column density using:

N({\rm H}_2)~[M_\odot~{\rm cm}^{-2}]=2\times1.36\times m_{\rm p}\times X \times I_{\rm co},
\end{displaymath} (6)

where 1.36 is the helium correction factor, $m_{\rm p}$ is the proton mass and $I_{\rm co}$ is the CO flux in K km s-1.

2.4 HI data

To study the distribution of the HI mass, we collect images at 21 cm from various sources[*]. For seven Virgo galaxies, we obtained moment 0 maps from VIVA, The VLA Imaging of Virgo in Atomic gas (Chung et al. 2008, in preparation; Chung et al. 2007); VIVA images have a spatial resolution of 15-17 arcsec and they reach a typical sensitivity of 1.5 $\times $ 1019 cm-2. For nine field galaxies, instead, we obtain HI images from THINGS, The HI Nearby Galaxy Survey (Walter et al. 2008); THINGS images have a spatial resolution of 11 arcsec and a sensitivity of 1.2 $\times $ 1019 cm-2. For the remaining galaxies, we collect HI data from the literature, either from Cayatte et al. (1990) or from Verheijen & Sancisi (2001). These maps have a lower spatial resolution, spanning from 40 to 70 arcsec. The sensitivity ranges from 1 $\times $ 1019 cm-2 to 1 $\times $ 1020 cm-2. For three galaxies (NGC 1068, NGC 4212 and NGC 6951) we cannot find HI maps, but only a measurement of their integrated HI mass. From the flux calibrated images, we compute the neutral hydrogen surface density using

N({\rm HI})~[M_\odot~{\rm cm}^{-2}]=1.36\times m_{\rm p}\times I_{\rm HI},
\end{displaymath} (7)

where 1.36 is the correction factor for the helium and $I_{\rm HI}$ is the column density in cm-2. The total HI mass in solar units is computed using:

\log M_{{\rm HI}}~[M_\odot]=\log(2.36\times10^5\times {\rm SintP} \times D^2),
\end{displaymath} (8)

where the total HI flux SintP (Jy km s-1) is from ALFALFA, The Arecibo Legacy Fast ALFA Survey (Giovanelli et al. 2005), or from the HyperLeda database and D is the galaxy distance in Mpc.

2.5 Optical data

The stellar mass distribution is obtained from imaging in the NIR H band from GOLDMine or the i band from the Sloan Digital Sky Survey (SDSS-DR6, Adelman-McCarthy et al. 2008). Exceptions are NGC 5236 and NGC 6951 (not imaged in the SDSS) for which we collect respectively I band and K band images from NED. In Table 1 we show individual references. Due to their large angular dimension, our galaxies are often cut in two or more SDSS fields, thus we retrieve images using Montage, a toolkit for assembling FITS images over multiple frames into a mosaic. The photometric calibration for galaxies from GOLDMine relies on the zero-points provided in the image headers, obtained according to Gavazzi et al. (1996b); GOLDMine images for NGC 4192 and NGC 4569 are from The Two Micron All Sky Survey (Skrutskie et al. 2006), thus for these galaxies we adopt the 2MASS standard calibration procedure. For all images, the typical uncertainty is <0.05 mag. SDSS images are calibrated (within 1% accuracy) using the zero-points derived from a set of parameters (calibration constant, air mass and extinction coefficient) provided in the tables associated with the images. To derive the stellar mass, we apply to each galaxy a mass-to-light ratio, computed through an extrapolation along the IMF according to Bell et al. (2003):

\log (M/L)_\lambda~[M_\odot/L_\odot]=a_\lambda+b_\lambda \cdot {\rm Color},
\end{displaymath} (9)

with $a_\lambda$ and $b_\lambda$ chosen accordingly to the band used. We use B-V Color from HyperLeda in combination with NIR data and g-i for i band data. The mean dispersion on M/L is 0.1 dex.

\includegraphics[width=6.5cm,clip]{0604fi1d.ps}\end{figure} Figure 1: Example of multiband images and profiles: NGC 4051. Top-left: i band; top-right: H${\alpha }$; bottom-left: CO; bottom-right: surface brightness profiles. $\Sigma _{\rm sfr}$ (red) in $M_\odot $ pc-2 Gyr-1; $\Sigma _{\rm HI}$ (magenta), $\Sigma _{{\rm H}_2}$ (blue) and $\Sigma _{\rm star}$ (green) in $M_\odot $ pc-2. The arrow indicates r25 as in Col. 10, Table 1.
Open with DEXTER

2.6 The surface brightness profiles

For each galaxy we derive a calibrated surface brightness profile at four frequencies. All images are firstly astrometrically calibrated using GAIA, then aligned and resampled to 3 arcsec/pixel. Stars near the galaxies are masked using the IRAF task imedit and the local sky background is estimated and subtracted using marksky and skyfit from Galphot[*]. Using a Gaussian kernel we smooth all images to a spatial resolution of 15 arcsec, in common with the CO observations and most of the HI images. To derive the surface brightness profiles we use the IRAF task ellipse that fits elliptical isophotes to the galaxy images. Starting from a set of initial parameters given manually for the NIR or i band imaging, the fit maintains as free parameters the ellipse center, the ellipticity and the position angle; the ellipse semi-major axis is incremented by a fixed fraction of its value at each step of the fitting procedure and the routine halts when the surface brightness found in a given annulus equals twice the sky rms. This task provides the statistical error of the surface brightness. From the same set of ellipses we derive the surface brightness profiles for the remaining images, with the exception of the HI data from Cayatte et al. (1994) and Verheijen & Sancisi (2001) for which we interpolate the published profiles on our grid. The derived surface density profiles are corrected for projection effects using the galaxy inclination on the plane of the sky. For each galaxy we finally obtain the HI, H2 and stellar surface density profiles  $\Sigma _{\rm HI}$, $\Sigma _{{\rm H}_2}$ and $\Sigma _{\rm star}$ in $M_\odot~{\rm pc}^{-2}$; the star formation rate surface density profiles  $\Sigma _{\rm sfr}$ are in $M_\odot~{\rm pc^{-2}~Gyr^{-1}}$. Moreover, we compute the color profiles: ${\rm sfr}-{\rm sta}=\log \Sigma_{\rm sfr} - \log \Sigma_{\rm star}$, ${\rm sfr}-{\rm H}_2=\log \Sigma_{\rm sfr} - \log \Sigma_{\rm H_2}$ and ${\rm H}_2-{\rm sta}=\log \Sigma_{\rm H_2} - \log \Sigma_{\rm star}$. Figure 1 illustrates an example of the imaging material and density profiles for the galaxy NGC 4051. From top-right (anticlockwise): H${\alpha }$ net; i-band (SDSS); CO image; surface brightness profiles of stars (green), HI (magenta), SFR (red) and H2 (blue).

3 Analysis

3.1 The Schmidt law

Following Kennicutt (1998), Wong & Blitz (2002) and Boissier et al. (2003) we test the validity of the Schmidt law on both local and global scales.

3.1.1 The local Schmidt law

\includegraphics[width=7cm,clip]{0604fi2b.ps}\end{figure} Figure 2: The local Schmidt law for the entire sample, separately for the molecular ( left panel) and the total gas ( right panel). Profiles are normalized to the maximum values for each galaxy. The dashed line in the left panel represents the index n=1, in the right panel n=1.4.
Open with DEXTER

Figure 2 (left) shows the correlation between the SFR and the H2 surface brightness in the individual galaxies belonging to the studied sample; in doing this comparison, we exclude the nuclear regions as previously discussed. Since in each galaxy the profiles do not have exactly the same radial extent, we prefer to truncate the most extended profile at the radius of the shortest one, avoiding any extrapolation. Typically, the SFR and the H2 emission have a comparable spatial extent, in which on average there is a fair correlation between the SFR and the molecular hydrogen surface density. In fact 80% of our sample shows a correlation index r>0.8. Only a few galaxies show a weak correlation, mainly due to peculiarities such as prominent bars empty of star formation (NGC 4303), starburst rings (NGC 4321) or significant gas displacements due to known severe environmental perturbations (NGC 4548). For the correlated galaxies, we perform a least square fit; despite the overall validity of the Schmidt law, 8/20 objects cannot be fitted using a single power law or the obtained rms is greater than 0.15 dex, arbitrarily assumed as a threshold. Inspection of the left panel of Fig. 2 shows that the faintest points in the SFR/H2 relation that belong to the outer parts of the disks significantly deviate from the correlation driven by the points associated with the inner part of the galaxy disks. For the remaining 60% of the correlated galaxies, we obtain a reliable fit using a single power law; the mean index obtained is n=0.93 $\pm$ 0.11, in good agreement with the finding of Wong & Blitz (2002) (n=0.8 $\pm$ 0.3). Despite this overall agreement, there are slightly differences in the individual indexes n for the five galaxies we have in common, mainly due to the fact that we exclude the innermost galaxy regions and that Wong & Blitz choose a conversion factor X that does not vary as a function of the radius nor from one galaxy to another.

We checked that there is no evident correlation between the star formation rate and the atomic hydrogen, especially in the inner part of the galaxy disks where the H${\alpha }$ emission is often strong but usually there is a low abundance of HI. To study the role of the total gas in the Schmidt law we compute the surface density profiles  $\Sigma_{\rm HI+H_2}$, combining the atomic and the molecular hydrogen profiles without any extrapolation. In Fig. 2 (right panel) we plot the Schmidt law for the total gas surface density on the same scale used for the molecular hydrogen alone (left panel). Figure 2 shows that the addition of the atomic hydrogen makes the correlation between the SFR and the gas better than considering the H2 alone; the most striking feature is that the faintest points in the local Schmidt law that deviate the most when considering the molecular hydrogen alone, now follow a single power law. Furthermore, we notice that the dispersion in the inner part of the disk is on average reduced, perhaps because by adding the atomic component we are taking into account the HI produced by the photodissociation in regions that have a high star formation activity. Inspecting individual galaxies, we find that all but two galaxies with r>0.8 can be fitted using a single power law with an rms < 0.15, i.e. the 60% of the entire sample can be described with a single Schmidt law. Using a least square fit, we derive a mean index n=1.39 $\pm$ 0.10, in good agreement with the value obtained by Kennicutt using integrated quantities. Despite confirming the correlation between the SFR and the H2that was suggested by Wong & Blitz, the important difference between our results and theirs is that the addition of HI improves this correlation, significantly reducing the number of galaxies for which the Schmidt law fails in the outer part of the disks. Therefore we cannot agree with the conclusions of Wong & Blitz that the correlation for the total gas is entirely driven by the molecular component. A possible cause for this disagreement might be their use of CO imaging from the BIMA survey. Inspecting the H2 radial profiles for the five galaxies we have in common, we detect molecular hydrogen emission that is up to two times more extended than in the BIMA observations. The analysis of Wong & Blitz seems more limited to the inner part of the galaxy disks, for which we agree that the correlation is almost entirely driven by the molecular component. Due to the superior sensitivity of the Nobeyama mapping to large-scale emission, we can push this analysis to the outer parts, where we find that the HI is the most important quantity that drives the correlation between the star formation rate and the total gas.

\includegraphics[width=7cm,clip]{0604fi3b.ps}\end{figure} Figure 3: The Schmidt law for integrated quantities, separately for the molecular ( left panel) and the total hydrogen ( right panel). Filled circles are for our measurements, while empty circles are data from Kennicutt (1998), for the galaxies we have in common. Uncertainties are 0.2 dex for SFR and up to 1 dex for H2 and for the total gas.
Open with DEXTER

3.1.2 The global Schmidt law

After having studied the local Schmidt law, we focus on its applicability while considering integrated quantities; therefore, we compute for each galaxy the total SFR and the gas surface density, by integrating along the profiles. In Fig. 3 we show the correlation between the SFR and the H2 (left panel) and the correlation between the SFR and the HI+H2 (right panel), for our measurements (filled circles) and for the galaxies we have in common, using data from Kennicutt (1998) (open circles). Uncertainties are not drawn, but as previously discussed the SFR has a typical error of 0.2 dex while the uncertainty for the H2 may be as large as 1 dex; the uncertainty for the total gas is entirely driven by that of the molecular hydrogen. We see that there is a fairly good agreement between our values and the ones from Kennicutt, confirming the validity of our calibration procedures, in spite of the slightly different corrections. Considering the left panel, we find a good correlation (r=0.76) between the SFR density and the molecular hydrogen density; therefore, for our sample composed of massive galaxies, we can assume that even on a global scale the Schmidt law is a suitable parametrization for the relationship between the SFR and the molecular gas. Adding the HI (right panel), we can see the same good correlation (r=0.77) between the SFR density and the total hydrogen density; therefore, considering a sample of massive galaxies, we find that the Schmidt law for the total gas does not differ significantly from the one for H2 alone, thus the bulk of the star formation that takes place in the inner part of the disks correlates with the molecular hydrogen surface density, in agreement with what is found on local scales. This finding supports the hypothesis that the poor correlation with the H2 found by Kennicutt is driven by the inclusion of low luminosity galaxies for which it is difficult to obtain reliable CO observations and for which a constant X conversion factor does not hold.

3.2 The environmental effects

In the previous section we have shown that the atomic hydrogen seems to play some role in the star formation process, especially at the edge of galaxy disks; moreover, the atomic hydrogen, built up during the cosmological nucleosynthesis, is the primeval constituent of the molecular hydrogen and therefore provides the essential fuel for the star formation, even if indirectly. Due to this fact, we now focus on the possible perturbations to the star formation activity that might be induced by the removal of the neutral hydrogen in the environment of a rich cluster of galaxies. The most suitable parameter to quantify the HI removal in a galaxy is the HI deficiency, computed according to Haynes & Giovanelli (1984) as the logarithmic difference between $M_{\rm HI}$ of a reference sample of isolated galaxies and $M_{\rm HI}$ actually observed in individual objects: ${\rm def}_{\rm HI}= \log M_{\rm HI~ref} - \log M_{\rm HI~obs}$. $\log M_{\rm HI~ref}$ has been found to be linearly correlated to the galaxy linear diameter d as $\log M_{\rm HI~ref}=a+b \log d$, where d (kpc) is determined at the $25{\rm th}$ B-band isophote, and a and b are weak functions of the Hubble type. Here we adopt the updated coefficients from Solanes et al. (1996). In the following analysis we choose to consider a galaxy HI deficient if the ratio between the measured and the expected HI mass is >2.5, i.e. if ${\rm def}_{\rm HI}\geq 0.4$. Therefore, our sample is composed of 54% non deficient galaxies that we consider unperturbed and 46% deficient galaxies that suffer from significant HI removal.

3.2.1 Global abundances

\includegraphics[width=7cm,clip]{0604fi4b.ps}\end{figure} Figure 4: The relationship between ${\rm def}_{\rm HI}$ and color residuals $({\rm H}_2-{\rm sta})_R$ ( right panel) and $({\rm sfr-sta})_R$ ( left panel). Correlation coefficients are respectively r=0.60 and r=0.65. Uncertainties are discussed in the text.
Open with DEXTER

\includegraphics[width=5.4cm,clip]{0604fi5c.ps}\end{figure} Figure 5: Comparison between the mean $\rm sfr-sta$ profile for deficient (dashed lines) and non deficient (continuous lines) galaxies; from left to right, deficiency threshold is set to be ${\rm def}_{\rm HI}\geq 0.4$, ${\rm def}_{\rm HI}\geq 0.65$ and ${\rm def}_{\rm HI}\geq 0.9$. Non deficient galaxies are always considered for ${\rm def}_{\rm HI}<0.4$.
Open with DEXTER

Firstly we focus on the relationship between the global SFR, H2 mass and stellar mass as a function of the deficiency. Therefore, for each galaxy we compute a mean color along the radial profile between three pairs of quantities: $\rm sfr-sta$, ${\rm sfr}-{\rm H}_2$ and ${\rm H}_2-{\rm sta}$, averaged over the entire disk. Since the color profiles typically exhibit a radial gradient that in some cases can be significant, considering a mean value for each galaxy may sometimes be a too crude assumption. We quantify this effect calculating the standard deviations $\sigma$ and we obtain:

\begin{eqnarray*}&& 0.07\leq\sigma\leq0.40~{\rm with~a~mean~of}~\bar\sigma=0.20~...
...h~a~mean~of}~\bar\sigma=0.31~{\rm dex~for}~{\rm H}_2-{\rm sta};

these deviations, plotted in Fig. 4, provide an estimate of the uncertainty in addition to the errors in the profile calibrations. To clarify the physical meaning of these colors we notice that $\rm sfr-sta$ is similar to the H${\alpha }$ equivalent width. In fact the $\rm sfr-sta$ color and the H${\alpha }$ equivalent width (EW) (retrieved from GOLDMine or collected from Romanishin 1990; Koopmann & Kenney 2004; Kennicutt & Kent 1983) are found to be tightly correlated (r=0.84). Similarly the ${\rm H}_2-{\rm sta}$ color quantifies the H2 emission, normalized to the stellar continuum; finally ${\rm sfr}-{\rm H}_2$ can be considered a tracer of the efficiency at which the molecular gas is converted into new stars. It is known that galaxies obey empirical scaling laws according to their dynamical mass and therefore to their H-band luminosity, owing to the tight correlation between these two quantities (Gavazzi et al. 1996a). As the luminosity increases, the HI mass to luminosity ratio  $M_{\rm HI}$/$L_{\rm H}$ (Gavazzi et al. 2008), the H2 mass to luminosity ratio  $M_{\rm H_2}$/$L_{\rm H}$ (Boselli et al. 1997) and the H${\alpha }$ equivalent width (Gavazzi et al. 1996a) decrease; in other words, the smallest galaxies are, relative to their mass, the richest in gas and the most active in star formation. Before studying if the above three sets of colors display a residual dependence on ${\rm def}_{\rm HI}$ independently of their mass, we need to correct for the primary luminosity dependence. Unperturbed galaxies in our sample ( ${\rm def}_{\rm HI}<0.4$) obey the following relationships, derived using a least square fit:

{\rm H}_2-{\rm sta}=-0.25\times \log L_{\rm H} + 1.55~({\rm rms=0.30~dex})
\end{displaymath} (10)

{\rm sfr-sta}=-0.26\times \log L_{\rm H} + 1.95~({\rm rms=0.32~dex}).
\end{displaymath} (11)

In Fig. 4 we plot, as a function of the HI deficiency, the residuals $({\rm H}_2-{\rm sta})_R$ (left panel) and $({\rm sfr-sta})_R$ (right panel) computed using the previous equations. Inspecting the left panel, we can see a significant correlation (r=0.65) between the $({\rm sfr-sta})_R$ and the deficiency; since environmental processes do not significantly affect the stellar distribution in massive galaxies, the meaning of this correlation is that HI removal causes an increasing reduction in the SFR with increasing deficiency. In the right panel a weaker but still significant correlation (r=0.60) exists between $({\rm H}_2-{\rm sta})_R$ and the ${\rm def}_{\rm HI}$, showing that molecular hydrogen, even if well bounded in the galaxy potential well, can be perturbed by the HI removal, at least indirectly. Conversely we find (not shown) that ${\rm sfr}-{\rm H}_2$ and ${\rm def}_{\rm HI}$ are uncorrelated quantities (r=0.20), suggesting that the efficiency of molecular gas to star conversion does not depend significantly on the amount of HI removal.

3.2.2 Radial behavior

\includegraphics[width=5.3cm,clip]{0604fi6c.ps}\end{figure} Figure 6: Comparison between the mean ${\rm H}_2-{\rm sta}$ profile for deficient (dashed lines) and non deficient (continuous lines) galaxies; from left to right, deficiency threshold is set to be ${\rm def}_{\rm HI}\geq 0.4$, ${\rm def}_{\rm HI}\geq 0.65$ and ${\rm def}_{\rm HI}\geq 0.9$. Non deficient galaxies are always considered for ${\rm def}_{\rm HI}<0.4$.
Open with DEXTER

After having shown that both the SFR and the molecular abundance are reduced in deficient galaxies on a global scale, a new question arises: is this reduction mainly due to a global quenching or it does take place primarily in the external disk, producing a truncation in the star forming disk? This issue has been addressed, among others, in a series of studies by Koopmann et al. They find that Virgo galaxies have their H${\alpha }$ disk truncated, while the case of starvation that produces a global reduction in star formation across the entire disk is rare (Koopmann & Kenney 2004). To assess this result in our sample and to verify if similar conclusions hold also for the molecular gas, we study the radial behavior of the color profiles. We normalize the radius of each galaxy using r25 from Table 1 and we compute the mean colors in bins of r/r25, averaging over all galaxies in each bin. In Fig. 5 we show a comparison between the mean $\rm sfr-sta$ profile for deficient (dashed lines) and non deficient (continuous lines) galaxies; from left to right, the deficiency threshold increases, i.e. in the left panel deficient galaxies are considered with $\rm def_{HI}\geq0.4$, in the middle those with $\rm def_{HI}\geq0.65$ and in the right panel the ones with $\rm def_{HI}\geq0.9$. Conversely, non deficient galaxies have $\rm def_{HI}<0.4$ in all panels. The first panel shows a gap of 0.4/0.5 dex between the $\rm sfr-sta$ color for non deficient and deficient galaxies, suggesting that the SFR in deficient galaxies is on average reduced. Increasing the deficiency threshold (middle panel) the gap increases up to 0.6 dex, while for the highest deficient galaxies $\rm def_{HI}\geq0.9$ (right panel) the gap increases up to 0.8 dex. Moreover, there are no galaxies with H${\alpha }$ emission beyond r/r25=0.8; therefore, for extremely deficient galaxies, besides the global quenching of the star formation rate, there is an evident truncation effect. This result suggest that the correlation found between the global  $\rm sfr-sta$ and the $\rm def_{HI}$ is mainly driven by quenching, while truncation is important only for highly deficient galaxies. For $\rm H_2-sta$, an effect similar to the one described for the $\rm sfr-sta$ is present in Fig. 6. In fact, in the left panel the gap between non deficient and $\rm def_{HI}\geq0.4$ galaxies is 0.4/0.6 dex, it increases to 0.8 dex when $\rm def_{HI}\geq0.65$ (middle panel); for $\rm def_{HI}\geq0.9$ (right panel) the gap increases significantly with increasing radii and there is no emission beyond r/r25=1. Once again, this finding is consistent with the fact that the correlation found between the global  $\rm H_2-sta$ and the $\rm def_{HI}$ is mainly driven by quenching, while truncation is important for highly deficient galaxies only. Conversely we find no significant difference between non deficient and deficient galaxies as far as the radial variation of the $\rm sfr-H_2$ efficiency is concerned, being on average within $\pm$0.2 dex. The existence of a significant correlation between the extent of the HI disk and the $\rm def_{HI}$ parameter is a well known result (e.g. Cayatte et al. 1994) implying that environmental effects, and in particular ram pressure stripping, perturb the ISM distribution outside-in, as also found in the simulations (Abadi et al. 1999). In this work we consider the HI isophotal radius, taken arbitrarily at 5 $\times $ 1019 cm-2, as a suitable measure of the HI disk. This threshold well matches the sensitivity of most of our HI data, except for NGC 3521, 4192, 5055, 5236 and 5457 where the HI column density stays above the threshold at all radii, but their exponential decay can be confidently extrapolated outside; only for NGC 4402 and NGC 4579, whose HI data are limited to the very inner regions, would any extrapolation be meaningless and we exclude them. In Fig. 7 we show the tight correlation (r=0.81) between the isophotal HI radii normalized to r25 and $\rm def_{HI}$, allowing us to quantify that at very large deficiencies ( $\rm def_{HI}\geq0.9$) it occurs that the extent of the HI disk becomes smaller than the optical radius. Is there something fundamental in the transition between galaxies that have their HI disk smaller or bigger than the optical radius? Does this transition also govern the difference in the SFR and H2 radial profiles between extremely deficient and normal galaxies? To test these hypotheses, we further investigate if the outer parts of the SFR and H2 radial profiles exhibit a different behavior as a function of the deficiency. We derive for each galaxy the radial gradient by fitting the color profiles ( $\rm sfr-sta$ and $\rm H_2-sta$) in the outermost regions ( $r/r_{25}\geq0.50$). The fit is performed with a maximum likelihood method, assuming a constant uncertainty along the profiles and using an exponential, a reasonable model for the $\rm H_2-sta$ since both the molecular and the stellar component typically behave as exponential disks. The issue is more complex for the $\rm sfr-sta$ since the H${\alpha }$ emission tends to be rather erratic on small scales; however, after a 15 arcsec smoothing even the H${\alpha }$ profiles can be fitted with an exponential function. The results of this analysis are summarized in Fig. 8, where we plot against $\rm def_{HI}$ the mean radial gradients for the outer part of the $\rm sfr-sta$ (top panel) and $\rm H_2-sta$ (bottom panel) color profiles; the mean values are derived averaging in four bins of deficiency (0-0.4, 0.4-0.65, 0.65-0.90, >0.90) and the error bars refer to the uncertainties on these averages. The first panel of Fig. 8 shows that the mean gradients of the $\rm sfr-sta$ color profiles are similar for non deficient and moderately deficient galaxies (up to $\rm def_{HI}=0.90$), while when considering highly deficient galaxies the discrepancy becomes significant at $1\sigma$. Finding steeper $\rm sfr-sta$ profiles for highly deficient galaxies implies that when the $\rm def_{HI}>0.9$ and thus the HI disk is reduced inside the optical radius, the SFR sharply decreases in the outer part of the disk, consistently with the truncation effect. The second panel shows that the discrepancy between the mean gradients of the $\rm H_2-sta$ color profiles in non deficient or moderately deficient galaxies and extremely deficient galaxies is less sharp. In fact, due to the great dispersion in the 0.65-0.90 and >0.90 bins, the difference between highly deficient galaxies and moderate deficient galaxies is less marked; however, it is evident that the steepest  $\rm H_2-sta$ are found in galaxies with $\rm def_{HI}>0.9$, once again supporting the idea that also the H2 profiles sharply decrease in the outer part of the highly deficient disks, suggesting that the truncation may also affect the H2 distribution.

\par\includegraphics[width=7cm,clip]{0604fi7.ps}\end{figure} Figure 7: The relationship between $\rm def_{HI}$ and the HI isophotal radius at 5 $\times $ 1019 cm-2, normalized to the r25. Correlation coefficient is r=0.81. The dashed lines emphasize that $r_{\rm HI}> r_{25}$ for $\rm def_{HI}>0.9$.
Open with DEXTER

\par\includegraphics[width=6.8cm,clip]{0604fi8b.ps}\end{figure} Figure 8: Mean radial gradients for the outer part of the $\rm sfr-sta$ ( top panel) and $\rm H_2-sta$ ( bottom panel) colour profiles, in four bins of deficiency (0-0.4, 0.4-0.65, 0.65-0.90, >0.90). Error bars are the uncertainties on the mean values in each bin.
Open with DEXTER

4 Discussion

4.1 The Schmidt law

The study of the local Schmidt law addresses one of the main open questions about the star formation activity in spiral galaxies: is the process of star formation related to the atomic or the molecular component? Clearly, the formation of stars takes place inside the GMCs, in cold and dense molecular regions on a scale of few parsecs; however, while studying the molecular distribution in our sample at a resolution of 15'', our observations are sensitive to the diffuse component on typical scales of 1 kpc; this molecular medium can sometimes be not self gravitating and therefore it may not take part to the star formation activity (Elmegreen 2002). From our analysis of the local Schmidt law we find that there is a good correlation between the SFR and the H2, that eventually becomes poorer in the outer part of the disks (see Fig. 2, left). Since the bulk of the star formation activity takes place in the inner disks, it can be concluded that the Schmidt law for the molecular component is a good model of the SFR in galaxies. Therefore it seems plausible that the molecular medium on the kpc scale traces the self gravitating molecular clouds which form at much smaller scales and sustain the star formation. However we found that the Schmidt law computed for the total gas, i.e. including the atomic hydrogen, is significantly improved, particularly in the outer disks (see Fig. 2, right). This result supports the theoretical argument by Elmegreen (2002) that some clouds can be in a molecular state without being strongly self-gravitating and, on the contrary, there might be instabilities in the HI component capable of triggering the star formation via direct collapse of HI; therefore it is reasonable to accept that the star formation rate scales with the total column density of gas in all forms, including a possible contribution from photodissociation of H2 molecules. We thus disagree with the conclusion of Wong & Blitz (2002) who state that the correlation in the Schmidt law is driven entirely by the molecular component. The fact that we find a mean index n=1.4 for the local Schmidt law, similar to the one originally proposed by Kennicutt using integrated quantities, seems to imply that this relation is scale invariant, i.e. it is able to describe the process of star formation on scales ranging from 1 kpc to the entire galaxy disk.

4.2 Environmental effects

The HI is the primary fuel for the star formation in galaxies, both directly and indirectly. We have shown that at the disk periphery the HI might directly contribute to the star formation by collapsing into molecular clouds and indirectly by providing the constituent of the diffuse H2 medium, composed of molecules which form on dust grains inside the stellar potential wells (Blitz & Rosolowsky 2004). Therefore, it is not surprising that the HI removal caused by the environmental processes may affect the star formation activity in cluster galaxies; we have shown that HI deficient galaxies have their SFR and H2 abundances reduced with respect to normal galaxies. Given the fact that we find no significant dependence of the SF efficiency on the environmental conditions, the observed trend of decreasing SFR and H2 density in deficient galaxies is consistent with models of galaxy starvation, as proposed by Larson et al. (1980) and elaborated by Bekki et al. (2002), according to which the effect of reducing the HI halo halts the replenishment of H2, leading to a quenching of the star formation rate. This model implies a continuous HI inflow from outside, as observed in some nearby galaxies, although at a rate that seems insufficient to sustain the current rate of star formation (Sancisi et al. 2008). When galaxies become extremely HI deficient, i.e. when the hostile environment has produced the ablation of their neutral gas up to a factor of ten and the HI disk reduces its size within the optical radius, superimposed on starvation, significant truncation of the star formation appears to exist, i.e. the SF disk is sharply eroded outside-in. The disagreement with Koopmann & Kenney (2004), who maintain that the reduction in the global star formation in the Virgo cluster is caused primarily by the spatial truncation of the star-forming disks, and not due to starvation, is only suggested. Their sample is dominated by extremely deficient spirals (among their deficient galaxies, 65% have $\rm def_{HI}\geq0.9$ and 85% have $\rm def_{HI}\geq0.8$). Beside the agreement with previous studies on the SFR in clusters, our analysis points out for the first time that the H2 content might be affected by environmental processes, even if indirectly as a consequence of HI depletion.

5 Conclusions

With the aim of exploring the relations between the process of star formation and the physical properties of the ISM in various galaxy environments, we collected state-of-the-art imaging material and maps for 28 massive spiral galaxies belonging to the Virgo cluster and to the local field. The observational material includes: images of the stellar continuum (taken in the red/NIR bands); H${\alpha }$ images of the young (<4 $\times $ 106 yrs) stars; radio maps at 2 mm from the recent Nobeyama CO atlas of nearby spiral galaxies that combines good sensitivity to large-scale CO emission with a 15 arcsec spatial resolution, and sensitive radio maps at 21 cm from the ongoing, yet unpublished, VIVA and THINGS surveys carried out at the VLA. Physical parameters have been derived homogeneously and carefully for the individual galaxies applying corrections based on measured quantities, rather than average values. We hope that the effort we put in the analysis of the data reduces the statistical limitations that arise from the paucity of imaging material currently available.

The present analysis indicates that the bulk of the star formation in spiral galaxies is supported by a diffuse molecular medium which forms through the conversion of the atomic hydrogen, due to the pressure exerted by the stellar potential. However, at the edge of the star forming disks, the HI plays a more important role than previously believed in directly sustaining the star formation activity. When environmental processes cause significant removal of the outer part of the HI disks, galaxies suffer from starvation; the replenishment of the atomic gas inside the optical disks is reduced, leading to a depletion of the molecular component that is consumed during the star formation activity, thus causing the quenching of the star formation activity itself. When the HI removal is so severe that the HI disk shrinks inside the optical radius, there is also a truncation of both the molecular and the star forming disk.

We thank A. Leroy, E. Brinks, A. Chung, J. van Gorkom and J. Kenney for their kind permission to use HI maps from THINGS and VIVA, prior to publication. We thank A. Boselli, L. Cortese, C. Bonfanti, P. Franzetti and B. Devecchi for useful hints and discussions. This research has made use of the GOLDMine Database. This research 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. We acknowledge the usage of the HyperLeda database (http://leda.univ-lyon1.fr). This research has made use of the SIMBAD database, operated at CDS, Strasbourg, France. This research made use of Montage, funded by the National Aeronautics and Space Administration's Earth Science Technology Office, Computational Technologies Project, under Cooperative Agreement Number NCC5-626 between NASA and the California Institute of Technology. The code is maintained by the NASA/IPAC Infrared Science Archive. IRAF is the Image Analysis and Reduction Facility made available to the astronomical community by the National Optical Astronomy Observatories, which are operated by AURA, Inc., under contract with the US National Science Foundation. STSDAS is distributed by the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc., under NASA contract NAS 5-26555. Funding for the Sloan Digital Sky Survey (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, and 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 (ARC) 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, The University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, The 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.



Copyright ESO 2008