abstract: We measure star formation rates (SFRs) of $\approx 50,000$ optically-selected galaxies in the local universe ($z\approx0.1$), spanning a range from gas-rich dwarfs to massive ellipticals. We obtain dust-corrected SFRs by fitting the *GALEX* (ultraviolet) and SDSS (optical) photometry to a library of population synthesis models that include dust attenuation. For star-forming galaxies, our UV-based SFRs compare remarkably well with those derived from SDSS-measured emission lines (primarily H$\alpha$). Systematic deviations from perfect agreement between these two methods is shown to be due to differences in the dust attenuation estimates. In contrast to measurements based on H$\alpha$, UV provides reliable SFRs for galaxies with weak or no H$\alpha$ emission, and where H$\alpha$ is contaminated with an emission from an AGN (1/2 of the sample). We use full-SED SFRs to calibrate a simple prescription that uses *GALEX* far-UV magnitude and the UV slope to produce good dust-corrected SFRs for normal star-forming galaxies. The specific SFR (SFR normalized by stellar mass) is considered as a function of stellar mass for (1) star-forming galaxies with no AGN, (2) those hosting an AGN, and for (3) galaxies without H$\alpha$ emission (the latter two groups forming the bulk of the optical red sequence). We find that the three have distinct star formation histories, with AGN lying *intermediate* between the star-forming and the quiescent galaxies. Normal star forming galaxies (without an AGN) lie on a relatively narrow linear sequence. Remarkably, galaxies hosting a strong AGN appear to represent the *massive continuation* of this sequence. On the other hand, weak AGN, while also massive, have lower SFR, sometimes extending to the realm of quiescent galaxies. We propose an evolutionary sequence for massive galaxies that smoothly connects normal star-forming galaxies to quiescent (red sequence) galaxies via strong and weak AGN. We confirm that some galaxies with no H$\alpha$ emission show signs of star formation in the UV. We derive a UV-based cosmic star formation density at $z=0.1$ with significantly smaller total error than previous measurements.
author: Samir Salim, R. Michael Rich, Stéphane Charlot, Jarle Brinchmann, Benjamin D. Johnson, David Schiminovich, Mark Seibert, Ryan Mallery, Timothy M. Heckman, Karl Forster, Peter G. Friedman, D. Christopher Martin, Patrick Morrissey, Susan G. Neff, Todd Small, Ted K. Wyder, Luciana Bianchi, Jose Donas, Young-Wook Lee, Barry F. Madore, Bruno Milliard, Alex S. Szalay, Barry Y. Welsh, Sukyoung K. Yi
title: UV Star Formation Rates in the Local Universe

# Introduction

Current studies of galaxies are characterized by two major features: availability of large samples of objects (galaxy surveys), and the utilization of the multiwavelength approach. Such studies extend from the nearby galaxies to those close to the epoch of the formation of the first galaxies. A rate at which a galaxy forms stars is one of the more important properties in studying galaxy evolution. The multiwavelength approach allows us to employ a suite of star formation (SF) indicators—from X-rays to radio wavelengths. Major effort in recent years was made to provide reliable calibrations for different SF estimators and to understand their differences, advantages and limitations. The most straightforward evaluation is achieved by comparing two or more star formation indicators for the *same* set of objects. While new galaxy surveys provide large statistical samples with which one can attempt such studies, the sample of galaxies for which more than a single star formation indicator can be applied is not necessarily large. Moreover, various SF indicators are often applied for different redshift regimes. Global cosmic star formation history is therefore the result of studies that employ different SF indicators.

Among the most frequently used star formation indicators are the UV continuum (usually at $\lambda <2000$ Å), nebular recombination lines (primarily H$\alpha$, but also \[OII\]), far-IR dust emission, and the synchrotron radio continuum at 21 cm . A comprehensive review of most of these methods was presented in , together with simple formulae for the conversion of the true flux into a star formation rate (SFR). Recently, the use of other SF indicators has been explored, such as the X-ray continuum (e.g., ), or the luminosity of PAH features in the mid-IR . Star formation indicators are only as good as the assumptions that connect a certain observed luminosity to the actual star formation rate, and therefore all are sensitive to various systematic uncertainties. One factor has proved more frustrating to account for than the others—the effect of dust obscuration on the UV and H$\alpha$ star formation rates. In contrast to these two dust-sensitive indicators, bolometric IR luminosity and the radio luminosity are often considered to be "true" SF indicators. While this is certainly an oversimplification, such perception is bolstered by the very tight correlation of the IR and the radio luminosities, at least as observed for the normal galaxies in the local universe . However, even if they were perfect, IR or radio methods cannot make UV or H$\alpha$ methods obsolete for many practical reasons.

There have been a number of studies that compared the UV and H$\alpha$ star formation rates—either as the observed luminosities, or by applying various schemes to correct for the dust attenuation.[^1] These comparisons were often carried out with respect to some dust-free SF indicator. used a sample of several hundred objects with H$\alpha$ and UV measurements (with the $U$-band photometry serving as a proxy for the UV), and compared them to the far-IR measurements compiled by . They find that in order to reconcile H$\alpha$ and far-IR luminosities one needs to apply a dust attenuation that is not fixed, but rather depends on the SFR itself. Actually, this dependence is a consequence of the relationship between the dust attenuation and the far-IR luminosity found by . use an extinction curve to extrapolate the dust attenuations obtained for H$\alpha$ into the UV regime, but find that such correction fails to bring UV luminosities to agree with the far-IR, i.e., the simple application of the same attenuation mechanism to both the H$\alpha$ and the UV emission does not appear to be correct. used actual UV observations of 50 nearby galaxies in two ultraviolet bands obtained with the Ultraviolet Imaging Telescope (UIT), and compared them to H$\alpha$ luminosities from the narrow-band imaging. They determine UV attenuation for 13 galaxies using the correlation with the UV slope , and measure H$\alpha$ attenuation for 21 objects using the thermal radio continuum fluxes. They find that both can reach high levels even for normal galaxies ($\sim 4$ mag for H$\alpha$, and $\geq 5$ mag for far-UV), and lend support to previous notions that the attenuation increases with SFR. confirmed that a better agreement between the UV observations from the FOCA balloon mission , and the fiber H$\alpha$ spectra is achieved when attenuation corrections are taken to be luminosity-dependent. They also suggest that one perhaps cannot use simple extinction curve scalings to convert H$\alpha$ attenuations into UV attenuations. The breakthrough in resolving this problem came with the introduction of the two-component dust attenuation model of . This model was indeed motivated by the need to produce a consistent model for dust attenuation affecting H$\alpha$ and UV continuum photons. It postulates the existence of short-lived (10 Myr) giant molecular clouds that affect photons producing the H$\alpha$ line. On the other hand, the attenuation of the UV continuum, having timescales longer than the lifetime of giant molecular clouds, is predominantly produced by the diffuse ISM (after the molecular clouds have dispersed), at levels that are typically 3 times lower (for a given wavelength) than those in the molecular clouds. In this paper we will be using the model, thus testing it for the first time on a large scale.

In addition to systematic trends, we should mention that some previous studies were finding that the measurement errors are smaller than the observed scatter. offer the explanation for this scatter as arising from the differing timescales over which H$\alpha$ and UV SFRs are sensitive, so that UV could in some cases (especially in low-mass galaxies) detect a short starburst that is no longer observable in H$\alpha$.

Sometimes, the systematics are present in the observations themselves. showed that even the relatively reliable determination of H$\alpha$ attenuation can be affected by the systematics when not correcting for the underlying stellar absorption in Balmer emission lines. Systematics can also arise with galaxy samples selected at different wavelengths .

The main obstacle in obtaining UV measurements for a large number of galaxies in the local universe is that they need to be made outside of Earth's atmosphere. For this reason in 2003 NASA launched the *Galaxy Evolution Explorer* (*GALEX*, ). *GALEX* is currently conducting the first ever survey of the UV sky. The imaging (in two UV bands) is executed in several modes—from a shallow all-sky survey, to the ultra-deep fields. In this paper we are using measurements obtained in the medium-deep survey that is designed to image regions of the sky covered by the Sloan Digital Sky Survey (SDSS). Thus we obtain a large sample of galaxies with both UV and optical photometry, as well as spectroscopic redshifts. *GALEX* and SDSS data, and the resulting sample, are described in §§ and . Dust-corrected star formation rates (and some other physical properties such as the stellar mass) are obtained by comparing the observed colors to stellar population synthesis models to which the dust attenuation has been applied (§ ). While we use full UV to optical SED, the SFRs are essentially constrained by the UV . The large size of our sample permits a robust statistical analysis. We compare our UV-based SFRs to the results of the major study of , who use SDSS spectra to derive H$\alpha$-based SFRs for $\sim 10^5$ galaxies (§ ). They employ an intricate scheme to correct for the fiber aperture effects. Thus, our study also serves as a check on the reliability of their methodology. Finally, we discuss star formation histories of different classes of galaxies (§) and derive a UV-based estimate of the global star formation density at $z=0.1$ (§).

# Data

## *GALEX* data

Technical aspects of *GALEX* telescope, detectors and data products are presented in . Here we give a summary of relevant information. *GALEX* surveys the sky in either the imaging or the grism spectroscopy mode. It simultaneously produces a far-UV (FUV) and a near-UV (NUV) image having a circular field of view of $1.25 \degr$ diameter. FUV and NUV filters are broadband, with effective wavelengths of $1528$ Å and $2271$ Å, respectively. We will denote the magnitudes measured in these photometric bands as *FUV* and *NUV*. A single field of view imaged by *GALEX* is called a *tile*. A tile can consist of one or more *visits*, i.e., individual exposures. *GALEX* surveys the sky in several imaging modes, which differ in the exposure time per tile. In this study we use the Medium Imaging Survey (MIS), which is designed to maximize the coverage of the sky that is included in the SDSS. Typical exposure times in MIS are 1500 s, yielding limiting magnitudes of $FUV=NUV=22.7$ mag (AB system is used throughout). *GALEX* cannot point in the vicinity of bright sources (usually stars), which inevitably produces some gaps in the coverage.

*GALEX* data used in this paper come from the MIS portion of the Internal data release 1.1 (IR1.1), which is an expanded version of the first public GALEX release (GR1) of MIS (in both IR1.1 and GR1 the same pipeline, version 4, is used to reduce the data and produce source catalogs). Dataset consists of 705 *GALEX* tiles. Because of the anomaly with the FUV detector, 98 tiles lack FUV images. The 705 tiles cover 797 sq.deg. of the sky. Each *GALEX* field is restricted to $0.6\degr$ radius, close to the maximum field of view.[^2] Because of the overlap between *GALEX* tiles, the total unique area is 741 sq.deg. Source catalog for each tile was produced using the SExtractor software .

In our study we use FUV and NUV fluxes measured in Kron elliptical apertures. We recalculate flux errors because they were incorrect in the pipeline version 4 reductions (thus also affecting GR1, cf. ). Kron magnitudes should represent a good measure of a galaxy's total flux. If both the FUV and the NUV detections are present, we use the aperture (size, shape and position) derived from the NUV image to measure the FUV light. We go from NUV to FUV because we are generally more sensitive in the NUV. We add zero-point calibration errors of 0.052 and 0.026 magnitudes to $FUV$ and $NUV$ respectively to account for systematic inaccuracies. The calibration errors were estimated by analyzing the repeat imaging of a calibration stellar object, and confirmed by a large number of repeat observations of the same field . If a UV detection is present in only one band, we still measure the formal flux and its error in the other band (using the aperture defined in the detected band). In some cases, the optical source is not detected at all by *GALEX*, in which case we note the sky background at the position of the SDSS source and compute the flux error. In this calculation we use galaxy sizes that have been derived (using a calibration based on objects detected by *GALEX*) from SDSS Petrosian radii in $r$-band. Finally, in some cases the FUV cannot be used because the FUV image does not exist (15% of all galaxies).

## SDSS data

In addition to the ultraviolet data from *GALEX*, we use optical data from SDSS Data Release 4 (DR4, ). SDSS is providing five-band broadband photometry ($ugriz$ bands), and the spectroscopic follow-up of most galaxies with $r\leq 17.77$ (main galaxy spectroscopic survey, ). In addition to the official data products from the SDSS collaboration, we use value-added galaxy catalogs produced by MPA/JHU SDSS team. These catalogs include reprocessed SDSS spectroscopic data, and some derived galaxy physical parameters that are based on the spectroscopic data.

SDSS catalog lists magnitudes measured in several different ways. We use MODELMAG magnitudes, which are the measurements of choice for relative fluxes (i.e., provide stable colors), while still capturing most of the total light. We apply slight adjustments ($-0.04$ and 0.02 mag) to $u$ and $z$ magnitudes respectively to bring them closer to the actual AB system . We convert SDSS magnitudes and errors into fluxes using the transformations given in . In addition to catalog photometric errors, we add 0.01 mag of uncorrelated calibration error to each of the bands (see §), and, specifically for the $u$ band, we add an additional color-dependent error due to the red leak ($\sigma_{u, {\rm
RL}} = 0.0865(r-i)^2+0.0679(r-i)$, based on ). In rare cases an SDSS magnitude would be missing and such band has to be excluded from the analysis. We further exclude from analysis individual SDSS magnitudes that are nominally fainter than $ugriz=25$. These are invariably spurious, regardless of the listed photometry error. Finally, we exclude magnitudes with errors larger than $10\times$ the typical error in that band (we find typical errors to be 0.07, 0.009, 0.007, 0.007 and 0.017 mag for $u$, $g$, $r$, $i$, and $z$ bands respectively). While large, we find these errors to be significantly underestimated. To summarize, except in cases listed above when we exclude certain individual flux points, our input data consists of 7-band photometry and spectroscopic redshifts from SDSS.

## *GALEX*/SDSS matched catalog

Of 741 sq.deg. of *GALEX* unique imaging, 645 sq.deg. overlaps with SDSS DR4 spectroscopic area, thus defining the solid angle of our sample.[^3] We estimate this area by counting the SDSS galaxies that: (1) have spectra, (2) have a dereddened magnitude $14.5<r_{\rm Petro}<17.5$ (the faint end is taken to be comfortably brighter than the spatially variable spectroscopic limit) and (3) lie at redshifts $0.005<z<0.22$. We count such galaxies both in our survey, and in the full DR4, whose spectroscopic area is known. The ratio of the two counts gives us the size of our survey area. This estimate should be good to within 1 sq.deg. There are 67,883 objects from the SDSS DR4 spectroscopic survey (not restricted to the main galaxy survey) that lie in this area, and are spectroscopically classified as galaxies. For each of these objects we search for a match in the *GALEX* source catalog (which already combines FUV and NUV detections) within $4''$. Our analysis of SDSS point sources with a match in *GALEX* indicates that *GALEX* positions have a random error of $0\farcs 8$ in either R.A. or declination (becoming somewhat larger at the edges of the field). In addition, there are overall tile-to-tile offsets between *GALEX* and SDSS coordinate systems of several tenths of an arcsecond. In any case, astrometric uncertainties are significantly smaller than our matching radius (see also ). A genuine match can be missed if the centroid of the optical light is offset by more than $4''$ compared to the centroid of UV light. We expect such cases to be quite rare, since at the mean redshift of the sample ($z= 0.104$) this offset would translate into a 7 kpc separation. A problem with matching in general is that what is considered to be a single object in one catalog can be resolved into multiple objects in another catalog, whether they are indeed separate objects (blending), or actually belong to the same system (shredding). This problem is more pronounced when combining catalogs with different resolutions, as is the case here (4 to $5''$ for *GALEX* vs. $1\farcs4$ for SDSS). However, in our particular sample this problem is not severe. SDSS galaxies with spectra are relatively bright objects, and if more than one *GALEX* object is found within the search radius our procedure was to take the one that is positionally the closest. Since we are dealing with high-latitude fields where the density of foreground stars is not that high, this simple matching procedure produces reliable results. Besides, we do have mechanisms of identifying the majority of incorrect matches at the later stage, when we perform an SED fitting to the combined *GALEX*/SDSS photometry (see §). Since we combine photometry from two different catalogs there is a concern of photometric zero-point mismatch. We explore such a possibility in § . In our final matched catalog we eliminate duplicate *GALEX* observations (stemming from overlaps or repeat observations) by keeping those that lies closer to the center of *GALEX* field of view.

# The sample

## Sample selection

We initially define our sample by applying magnitude and redshift cuts to galaxies with SDSS spectroscopy (note that objects spectroscopically classified as QSOs are excluded). We require the dereddened Petrosian magnitude to be in the $14.5<r_{\rm
Petro}\leq17.77$ range, where the faint limit is the nominal limit of the SDSS main galaxy spectroscopy survey (see also § ), while the bright limit is chosen to avoid objects with saturated SDSS photometry. We require redshifts to lie within $0.005<z\leq 0.22$ range. The lower limit is chosen to eliminate galaxies that may deviate from the Hubble flow, i.e., whose redshift distances could be unreliable. Redshift distribution in SDSS main galaxy sample peaks sharply around $z=0.1$ with few galaxies beyond our upper redshift limit. These magnitude and redshift limits are identical to those used in , to which we will be comparing many of our results. We will refer to the above redshift range as representing the "local universe".

There are 49,346 galaxies that meet the conditions defining our initial sample. Good quality SED fitting (see §) is obtained for 98% of them (48,295 galaxies). In the rest of the paper we will use only this *latter* sample and refer to it as "all" galaxies. Note that since we retain objects regardless of whether they were detected by *GALEX* (as long as they fall within the area covered by *GALEX*), our sample is only *optically* ($r$-band) selected.

Throughout the paper we will use $\Omega_m=0.3$, $\Omega_\Lambda=0.7$, $H_0= 70\, {\rm km\, s^{-1}\, Mpc^{-1}}$ (i.e., $h_{70}=1$) cosmology.

## Emission line diagnostics and sample classification

This work focuses on the physical properties of galaxies, such as their star formation rates (SFRs) and stellar masses. Before we start analyzing the sample based on these *derived* properties, we would like to characterize it in terms of its *observable* quantities.

Throughout this paper we will rely on optical emission lines to classify galaxies in our sample. This classification, based on the BPT diagram , plots the flux ratio of \[OIII\]$\lambda$`<!-- -->`{=html}5007Å and H$\beta$ lines against the flux ratio of \[NII\]$\lambda$`<!-- -->`{=html}6584Å and H$\alpha$ lines. In this paper we fully adopt the BPT classification of (hereafter B04), which is based on emission-line strengths corrected for the underlying stellar absorption (see their Fig. 1). The position of a galaxy in the BPT diagram indicates the nature of its ionizing source. The majority of galaxies in the diagram fall within one of the two branches. One branch lies mostly above the maximum line ratios expected from star formation. These galaxies must have some ionizing source in addition to young stars. Most researchers attribute this emission to an active galactic nucleus (AGN). Specifically, the emission is associated with a narrow-line AGN (a LINER or a Seyfert 2), also known as Type 2 AGN . Note that any Type 1 AGN spectrum (QSO or Seyfert 1) would have broad lines, and would be classified by the SDSS pipeline as a QSO and thus not included in this sample. Following , B04 classify galaxies that lie on the lower portion of the AGN branch, but with line ratios allowed by star formation, as "composite" (star forming/AGN) galaxies. For galaxies that lie on the star forming branch we expect very little emission line flux to come from an AGN . B04 required a S/N ratio in all four lines to be $>3$ in order to apply the BPT diagram classification. However, in some cases a secure AGN classification is possible even when only \[NII\]6584 and H$\alpha$ have S/N$>3$, i.e., in cases when this ratio is larger then the one allowed by SF. B04 calls this class a "low-S/N AGN (LINER)". Following B04, we will study this class together with the AGN. B04 also introduce the category of "low-S/N SF" galaxies. These are the galaxies that do not fall in any of the previous categories because their lines have low S/N, yet they have an H$\alpha$ detection with S/N$>2$. While their lines are not strong enough for secure placement on the BPT diagram, if we nevertheless do so, we find that many object lie in the high-mass end part of the SF sequence, as well as in the lower portion of the AGN branch. Therefore, these objects represent a heterogeneous class. Finally, there are galaxies without detectable lines, thus precluding the classification in the above scheme. B04 call this group "unclassifiable". We will call them "No H$\alpha$" class. We find that in $\sim 2\%$ of galaxies in this class the H$\alpha$ non-detection is due to some artifact in the spectrum or line-extraction pipeline. We exclude these galaxies from this class (but not from the whole sample). The fraction of galaxies in our sample belonging to different classes is as follows: Star forming (SF) – 27%, Low S/N SF – 19%, Composite (SF/AGN) – 8%, AGN – 12% and No H$\alpha$ – 33%.

## UV-to-optical color-magnitude diagram

Color-magnitude diagram (CMD) is a powerful tool in assessing the basic properties of a sample of galaxies. Historically, the study of CMDs was preceded by the studies of a relationship between color and morphology. Basic morphological segregation into disk-like spiral galaxies and spheroidal elliptical galaxies was established by . Afterward (e.g., ) it was realized that spiral (late type) galaxies have bluer colors than the ellipticals (early types). Optical CMDs were constructed for cluster early type galaxies (e.g., ), where they were found to form a narrow sequence (the so called red sequence). Field spiral galaxies also displayed the color-magnitude relationship , albeit with larger scatter. The bimodal nature of the field galaxy CMD became much more apparent with recent large scale surveys, in particular SDSS . However, unlike the traditional optical CMDs, a CMD in which the color is constructed from an *ultraviolet* and an optical magnitude has a particular diagnostic power. By virtue of contrasting the current (or recent) star formation as indicated by the UV light (modulo attenuation) to the total past star formation as indicated by the optical light, the UV to optical color is a good proxy of a galaxy's SF history (e.g., ). On the other axis, one plots absolute optical magnitude, which is related (modulo variations in the optical mass-to-light ratio) to galaxy's current stellar mass. We present the observed CMD of our sample in Figure . In this and many subsequent figures, the individual data points have been converted into a greyscale density scatter plot, in such a way that the shade of gray is proportional to the number of objects occupying a given "pixel". Such representation is desirable when dealing with large samples where it is easy to saturate a traditional scatter plot. Unless noted otherwise, the figures are constructed from raw, unweighted counts. We will include volume corrections later, when appropriate. Both the color and the absolute magnitude have been K-corrected to $z=0$ rest-frame bands (§). The upper left panel shows the CMD of all galaxies with an NUV detection (85% of the total sample; see Table for a breakdown of UV detection rates per class). A striking feature of a UV to optical CMD is the pronounced bimodality of blue and red galaxies. Blue galaxies form a well-defined sequence extending to faint luminosities. The red sequence is somewhat more narrow than the blue sequence (note that this becomes evident only after K corrections have been applied), and extends to intrinsically more luminous galaxies than the blue sequence. This is related to well-known fact that the most massive elliptical galaxies are more massive than the most massive spirals . The two sequences are separated by $\sim 3$ mag. If the two sequences are modeled as gaussians, one finds that there is an excess of galaxies in the gap . This is not the case in classical optical CMDs . We will refer to the gap region and its population as the "green valley" . A detailed quantitative study of the *GALEX* UV to optical CMD is presented in .

```latex
\begin{deluxetable}{lrrrr}
\tablewidth{0pt}
\tablecaption{Average errors of ``UV'' and ``{H$\alpha$}''-based
  star formation rates\tablenotemark{a}}
\tablehead{ Class & Number & UV detected\tablenotemark{b}  &
``UV'' & ``H$\alpha$'' \\
 & & &  $\langle \sigma({\rm log~SFR}) \rangle$ &
 $\langle \sigma({\rm log~SFR}) \rangle $ }
\startdata
All\tablenotemark{c}  &  48295  &  86\%  &  0.38  &  0.43 \\
SF              &  12901  &  99\%  &  0.20  &  0.29 \\
SF (low S/N)    &   9060  &  93\%  &  0.30  &  0.39 \\
Comp            &   3966  &  96\%  &  0.28  &  0.40 \\
AGN             &   5827  &  90\%  &  0.41  &  0.49 \\
No {H$\alpha$}&  16159  &  68\%  &  0.60  &  0.54 \\
\enddata

\tablenotetext{a}{UV-based SFRs are averaged over 100 Myr.}
\tablenotetext{b}{Either FUV or NUV detection. Of objects detected in
  FUV, 98\% are also detected in NUV.}
\tablenotetext{c}{Includes 382 objects for which classification was
  not possible, see \S \ref{ssec:bpt}.}
\label{tab:errors}
\end{deluxetable}
```

In the subsequent panels in Figure we display the CMDs of various classes of galaxies as defined in §. Each panel has the greyscale normalized to the number of galaxies in the given class. For better reference with respect to the full sample, in each subsequent panel we repeat the contour containing 90% of *all* galaxies. Not surprisingly, the star-forming galaxies (SF) occupy the blue portion of the CMD. Part of the width of the sequence is due to the intrinsic dust attenuation. The CMD of composite galaxies (showing both signatures of SF and AGN) are shown in the upper right panel. Their $NUV-r$ colors are offset to the red compared to those of the "pure" SF galaxies. Galaxies with narrow line AGN (lower left panel) occupy the regions of the red sequence *and* of the green valley (). Moreover, most galaxies with intermediate colors are either AGN or composites. Star formation histories of AGN will be discussed more extensively in § . Next, we have low-S/N SF galaxies (middle lower panel), which indeed mostly lie on the blue, star-forming sequence (preferentially its more luminous part). However, there is a tail of red galaxies which are probably contaminated by AGN as observed in § . Finally, we have galaxies with no detectable H$\alpha$ (lower right panel). As expected, these galaxies form the bulk of the red sequence. However, there is a tail of galaxies of this class extending towards the blue colors. Since the spectral classification used here is nuclear (central $3''$) it is possible that some of these galaxies are classified as not having H$\alpha$ because of the dominant bulge, while the relatively faint star forming disk is giving a galaxy an overall blue color. However, as we discuss in § , this is not the case for most of them. The CMDs presented in Figure require a detection in NUV. The remaining 15% of our sample with no NUV detection falls mostly in the red sequence (as evidenced from their $u-r$ colors). These galaxies are too faint to be detected in the UV.

# Obtaining galaxy properties with SED fitting

## SEDs of model galaxies

Spectral energy distribution (SED) fitting is becoming a widely used technique for deriving galaxy properties. It was pioneered in works of and . To first order it consists of comparing the observed SED to a set of model or template SEDs, and searching for the best match. It is assumed that the physical properties of model or template galaxies are known, and that one can use this knowledge to deduce the properties of an observed galaxy.

In this study, we compare the observed SED to a large number of model SEDs constructed from population synthesis code. To construct the models we first choose input model properties by randomly selecting values from prior distributions defined in the following way. Formation of a galaxy is uniformly distributed between 0.1 Gyr and the age of the universe. Global metallicity is uniformly distributed between 0.1 and 2 $Z_{\odot}$. Star formation histories are not single stellar populations, but the combination of an exponentially declining continuous star formation, SFR$(t)\propto\exp(-\gamma t)$, with $0\leq\gamma\leq1\,$Gyr$^{-1}$ uniformly distributed over that range, and the random *starbursts* superimposed on the continuous SF. Bursts are parametrized to have a duration uniformly distributed in the 30–300 Myr range, with a strength such that the mass produced in a burst is between 0.03 and 4 times the mass produced in the continuous SF over the present lifetime of a galaxy (distributed uniformly in log). Such parametrization conforms with most observational studies. Finally, bursts are randomly produced so that the probability that a given galaxy has undergone at least one burst over a 2 Gyr period is 50%. The above parameters define how the population synthesis for each model will be carried out. Next, each of the above model SEDs is subjected to dust attenuation parametrized according to the two-component model of . For the $V$-band optical depth we randomly chose a value from a distribution $0<\tau_{V}<6$ (it peaks at $\tau_V=1.2$ and has most values in the $0<\tau_{V}<2$ range). The prior distribution for $\tau_V$ is empirical, and comes from Balmer decrements in SDSS spectra. The choice of $\tau_V$ prior distribution will be discussed in §. Fraction of the optical depth that affects stellar populations older than 10 Myr (i.e., most of the UV continuum flux) is given by factor $\mu$, with values ranging from 0.1 to 1, peaking around $\mu=0.3$. Altogether, we produce 100,000 model spectra spanning a range of SF histories, metallicities and dust attenuations. Note that since we pick input parameters randomly, we do not call our set of models a grid, which would suggests a set of points with equal spacing in parameter space.

Finally, we convolve the resulting model spectra with *GALEX* and SDSS bandpasses at five redshifts equally spaced in the \[0.05, 0.25\] interval, producing the libraries of model broad-band photometry. In each library we keep only models that have an age smaller than the age of the universe at that redshift. This effectively reduces the number of models from 95,000 at $z=0.05$ to 78,000 at $z=0.25$. We also add effects of the intervening intergalactic extinction according to . Our final libraries list model photometry as well as a number of galaxy parameters associated to that model (such as the SFR averaged over several timescales, stellar masses, dust attenuation parameters, etc.)

## SED fitting

In our study, the observed SED is constructed from *GALEX* and SDSS broadband photometric fluxes (the broad-band SED). For an observed galaxy at some redshift we select the model library with the closest redshift. We step through a library one model at the time. Model flux points will have some arbitrary zero point, and in order to see how well an observed flux compares to the model flux we first need to find a factor $a$ that minimizes the $\chi^2$ between the observed and the model points. In other words, for model $i$, we need to minimize the following expression:

$$\chi^2_i = \sum_X \left(\frac{F_{{\rm obs},X}-a_i F_{{\rm mod}_i,X}}
{\sigma(F_{{\rm obs},X})}\right)^2 \label{eqn:chi2}.$$

Here, the summation is over $X=(FUV,NUV,u,g,r,i,z)$, the observed flux points are $F_{\rm obs}$, and their errors $\sigma(F_{\rm obs})$. Flux points of model $i$ are $F_{{\rm
mod}_i}$. By taking the derivative of Equation with respect to $a$, and equating it with zero, we find the scale factor $a^*$ that best matches the observed and the model flux points:

$$a^*_i = \frac{\sum_X \frac{F_{{\rm obs},X} F_{{\rm mod}_i,X}}
{\sigma^2(F_{{\rm obs},X})}} {\sum_X \left(\frac{F_{{\rm
mod}_i,X}}{\sigma(F_{{\rm obs},X})}\right)^2}.$$

Because we include non-detections that have meaningless fluxes (small or even negative values), we use the observed fluxes *without* correcting for Galactic reddening, and instead apply the reddening to model fluxes, which are always positive. This allows us to treat non-detections like any other flux point with a known flux error. Reddening corrections for SDSS magnitudes come from SDSS database, and those for *GALEX* from relations in (Their Eq. 2 for $NUV$).

The goodness of the fit between an observed galaxy and model $i$ is then obtained by inserting $a^*_i$ in place of $a_i$ in Equation . Note that only a single parameter is being fit—the scale factor $a$. We emphasize that the word "parameter" when fitting is discussed should not be confused with (in principle) arbitrary number of *galaxy* parameters (i.e., properties) that *correspond* to each model SED. Therefore, since we have seven photometric points, there are six degrees of freedom. For a given galaxy we evaluate $\chi^2_i$ for *each* model as outlined above. Now, each of these $\chi^2_i$ values will determine the relative *weight* $w_i$ for the given model $i$ (and therefore the weight of the galaxy parameters associated with that model) as $w_i=
\exp(-\chi^2_i/2)$. We then build a probability distribution function (PDF) for every galaxy parameter of interest by compounding weights at the corresponding parameter value. We repeat the procedure for each of the $\sim 10^5$ models in the given redshift library. We normalize the final PDFs and note the parameter values corresponding to the 2.5, 16, 50 (median), 84 and 97.5 percentiles of the cumulative PDF. We use the average of the PDF as our nominal estimate of the parameter value. Since in the general case the PDF is not symmetric, medians, averages and modes will differ. For most relevant parameters in this study the differences between the median and the average are small. We do not consider using a mode, because due to the discreteness of the model parameter space, the mode can be stochastically offset from the bulk of the PDF, and is sensitive to PDF binning. For similar reasons we do not use a *single* model with the best $\chi^2$ (i.e., the maximum likelihood model) as representative of parameter values, but we do keep track of the best $\chi^2$ values in order to evaluate the overall quality of the fitting for a given galaxy (§ ). Additionally, for certain parameters we preserve more detailed information on the shape of the PDF (in the case of the specific SFR and stellar mass, we actually keep their mutual two dimensional PDFs). In certain cases we will discuss the errors of the derived parameters. We use 1/4 of the 2.5-97.5 percentile range as a proxy for what would have been a $1 \sigma$ error in the gaussian distribution. We will refer to these quantities as "formal" errors.

We prefer the above-described Bayesian approach (e.g., ) to often-used maximum likelihood, since it allows one to transparently test the dependence on the prior, i.e., how well the data constrain certain parameters given that the shape of the resulting PDF can to some extent depend on the distribution of the given parameter in the model library (especially when the observational constraints are poor). Another advantage is that while the maximum-likelihood model might indeed correspond to a high probability density, in the case of extensive models such as ours, it might actually correspond to a negligible probability "mass", and might therefore lead to misleading parameter estimation (see for instance the discussion in ). Note also that our technique properly accounts for the degeneracies inherent in some galaxy properties. Imagine that we have two galaxy parameters whose values can be picked in such a way as to produce *identical* model SEDs. In such case our fitting will give *equal* probabilities to all of the degenerate models, and the resulting PDFs for each of these two parameters will therefore reflect the entire range that these parameters can take (i.e., there will be nothing to constrain the PDFs). If we then characterize the errors of those parameters using the width of the PDF, we will obtain very large values, i.e., we will know that a certain parameter is poorly constrained. This, in effect, is in contrast to fitting using a *single* best-fit model, which would pick one of the degenerate models as the best-fitting (and the corresponding parameter value), without giving a researcher an idea that there is an underlying degeneracy making the obtained value to be quite uncertain. Finally, we want to emphasize that the widths of the PDFs (and also their typical values) have a meaning only if the observational errors have been estimated reasonably correctly, and if the overall fit (which can be well characterized by the $\chi^2$ of the best-fitting model) is reasonably good (§).

## Quality of a fit

For our initial sample of 49,346 galaxies we obtain model parameters in the way described in §. For each galaxy we also note the $\chi^2$ of the *single* best-fitting model, $\chi^2_{\rm
best}$. This value serves as an indicator of the overall quality of the fitting for a given galaxy. Since we include galaxies close to the edge of *GALEX* field of view where we know that image PSFs get distorted, we would first like to check if this affects the photometry, and therefore the quality of SED fitting. Thus we look at the average $\chi^2_{\rm best}$ as a function of distance from the *GALEX* image center. We find no correlation with the distance or degradation close to the edge.

In the ideal case, $\chi^2_{\rm best}$ values would form a theoretical $\chi^2$ distribution. The average of such distribution should equal the degrees of freedom (in our case 6). Before looking at the $\chi^2$ distribution we notice that $\chi^2_{\rm best}$ values are systematically higher for blue galaxies compared to the red. We believe that this difference arises from the presence of emission lines in blue galaxies that effectively increase the discrepancy between the observed SED and the model (line-free) continuum. Therefore, we start by comparing the observed and the theoretical $\chi^2$ distribution of red galaxies. Initially, the observed distribution was shifted towards the values that were too low. This was indicative of photometry errors being overestimated. We arrive at the $\chi^2$ distribution that was well-matched to the theoretical one when we adopt 0.01 mag calibration errors in each of the 5 SDSS bands.[^4] While now the bulk of the red galaxies follows a $\chi^2$-like distribution, there is also present a tail of high $\chi^2_{\rm best}$ values, in numbers well above any expected distribution. One can imagine that the main reason for the presence of $\chi^2$ outliers is that the data do not reflect the actual SED for some reason. We will investigate possible causes momentarily. The high-$\chi^2$ tail of red galaxies contains some 2% of objects. We decide to exclude these galaxies from further analysis (although we do account for them statistically when making any kind of space density calculation). Assuming that the reasons driving extreme $\chi^2$ values are color independent, we also clip 2% of the blue-galaxy $\chi^2_{\rm best}$ distribution tail. This leaves us with 48,295 galaxies—our final sample.

While removing the outliers is necessary to have results not be biased by unreliable values, it is instructive to understand what causes some $\chi^2$ to assume such high values. We thus examine 22 galaxies with extreme $\chi^2_{\rm best}$ values, distributed evenly in the $-1<NUV-r<7$ range (thus sampling galaxies with various SF histories). In 10 cases we find obvious problems with cataloged photometry. This is most often the case with SDSS (9 objects), where the visual inspection of broadband SEDs and comparison with spectra indicates that some flux points are large outliers. Next, in 7 cases we are dealing with either close galaxies being blended in *GALEX* (4), or with SDSS shredding galaxies into smaller photometric objects (2). For one object it is unclear whether we have blending or shredding. Further, in one case an object shows a clear QSO spectrum, i.e., the continuum is affected by the AGN, yet it is spectroscopically (mis)classified as a galaxy. Finally, in four cases it appears that the broad-band SEDs (and spectra) are genuinely unusual. While their number is too small to affect the present study, this type of outliers deserves a special scrutiny. It needs to be explored if there exists any set of model input parameters (SFH, dust attenuation, metallicity) beyond our already wide range, which can produce a model SED matching these observations. For our full sample we also expect on the order of 10 cases where a foreground star of similar brightness is superimposed with a galaxy, and falls within the fiber, thus producing a "composite" galaxy/star spectrum which would fail the fitting.

Since our model libraries are calculated at some fixed redshifts, there is some error associated with the library redshift not perfectly matching the galaxy redshift. We thus contrast the $\chi^2_{\rm best}$ values of galaxies having the maximum redshift offset (close to 0.025) to those with no offset. We find a systematic increase in $\chi^2_{\rm
best}$ of 10%, which is too low to affect the overall $\chi^2$ distribution.

Finally, we want to assess the effects of the "resolution" of our model parameters, i.e, we can ask how extensive should a model library be that the resolution is not an issue. To that extent, we artificially degrade our model libraries by running full SED fitting using every 2nd, every 4th and every 8th model SED. We then compare average $\chi^2_{\rm best}$ values from each run. Overall, the quality of the fits are surprisingly stable. The run with every 2nd model produces only a 4% increase in $\chi^2_{\rm best}$ values, while the run using every 16th model is 16% worse compared to the full library. This implies that increasing the extent of the library (for the given range of input parameters) would not bring significant improvements. This, of course, is not a general conclusion, but rather states that with the precision of our data the current number of models is adequate, i.e., having a finer resolution of input model parameters would create model SEDs that are degenerate from the observational point of view. In addition to this, the models themselves are not likely to be of such an accuracy to warrant more extensive model library.

## K corrections and volume corrections

For deriving the galaxy properties from the SED fitting, we are not required to know the K corrections. This is because for a galaxy at a given redshift we compare the observed magnitudes to model magnitudes with bandpasses that have been shifted to the matching galaxy's redshift (i.e., the redshift of the model library). Therefore, the "K corrections" are present in the SED fitting implicitly. However, for some applications, such as constructing a CMD (§), or obtaining $V_{\rm max}$ weights (§§,), we need to know the K corrections explicitly. The process of SED modeling allows us in principle to derive the K corrections alongside galaxy parameters. It even has the advantage of allowing one to construct a PDF for any K-corrected magnitude, which thus yields an estimate of a K correction *error*. We find this approach to be especially useful in K correcting higher redshift ($z\sim1$) samples (Salim et al., in prep.) In the current case, given that our redshift range is covered by only 5 model libraries, we choose to calculate the K corrections using the publicly available code KCORRECT V4_1_4[^5] . The code allows *GALEX* magnitudes to be used together with SDSS magnitudes to constrain the SED fit from which K corrections are derived. We obtain K corrections using the same combined *GALEX*/SDSS photometry used for the SED fitting, except that $<3\sigma$ UV detections are not used. We derive K corrections without a priori assuming an evolution of the luminosity function. We will explicitly correct for evolution where appropriate.

In cases when we require a volume-corrected sample, we assign a weight to each galaxy according to its $V_{\rm max}$ value—the volume in which a galaxy would be visible taking into account redshift and apparent magnitude limits, and the solid angle of a survey. While our faint limit ($r=17.77$) is also the nominal limit of the SDSS main spectroscopic survey, we take into account that the actual spectroscopic magnitude limit varies from one spectroscopic plate to another (i.e., that it can be brighter than $r=17.77$). Thus, for each object we modify the volume by the spectroscopic completeness (usually around 0.9). We take spectroscopic limits and completeness values from NYU-VAGC catalog . We caution that in order to calculate $V_{\rm max}$ correctly, it has to be done iteratively, because the K correction depends on redshift limits that determine $V_{\rm max}$.

# Comparison of "UV" and "H$\alpha$" Star Formation Rates

## "UV" SFR estimates

We derive our SFRs from photometric constraints that extend from the UV to the $z$-band. However, it is the UV luminosity that constrains the SFR the most. Therefore, in the remainder of the paper we will refer to SFRs derived from the SED fitting as "UV", where quotation marks stand to remind us that it is not only UV information that contributed to these SFR estimates. In our model SED libraries we report "current" star formation rates averaged over several time intervals, most notably over the most recent 10 Myr, 100 Myr and 1 Gyr. By comparing the formal errors of SFR estimates (derived from 95 percentile range of each galaxy's SFR PDF), we find that the best-constrained SFR for the overall sample is the one averaged over the last 100 Myr (Table ), which is the timescale for UV bright stars. For the star-forming class specifically, we find that the SFRs have slightly smaller formal errors over 1 Gyr timescales than over 100 Myr. On the other hand, SFR estimates on timescale of 10 Myr for SF galaxies (a timescale that matches H$\alpha$ SFR) have drastically larger formal errors, as expected since we do not have information to constrain them. In the remainder of the paper we will use SFR estimates averaged over 100 Myr, which represent a good compromise between the quality and a timescale that is not too long compared to that of the H$\alpha$ SFRs.

## "H$\alpha$" SFR estimates

B04 have developed a method of deriving SFRs from SDSS spectra that is primarily based on the intensity of an H$\alpha$ line. This represents the largest currently available sample of precise H$\alpha$-based SFRs. Obtaining accurate SFRs from SDSS spectra is far from straightforward. For the purposes of this paper we briefly describe the procedures employed in B04. The essence of the B04 approach was to model SDSS spectra by first removing the absorption line spectrum using a combination of burst models from , and then to model the emission lines using the models which combine the galaxy evolution models with the CLOUDY photoionization code . Also, the same two-component dust attenuation prescription of was used by B04 to model attenuation of the emission lines, and by us to model broad-band continuum attenuation. While conventionally the attenuation of H$\alpha$ flux is determined by comparing the observed H$\alpha$ to H$\beta$ ratio (Balmer decrement) to the theoretical one, the attenuation estimate in B04 is constrained using many emission lines (although it is dominated by H$\alpha$/H$\beta$). Using a suite of models, B04 apply a Bayesian approach to produce the probability distributions for each of the four parameters (gas metallicity, ionization parameter, dust attenuation, and dust-to-metal ratio). B04 thus simultaneously produce attenuation estimate as well as the attenuation-corrected SFR within SDSS $3''$ fiber. SFR within the fiber is constrained by many emission lines, with the greatest weight carried by H$\alpha$. This is one reason why we will denote B04 SFR and attenuation estimates by "H$\alpha$", with quotation marks again serving to indicate that not only H$\alpha$ was involved in those estimates. The above procedure is directly applicable *only* to galaxies for which the ionizing source is predominantly star formation, and not, for example, an AGN activity. B04 therefore calculate *fiber* SFRs directly from emission lines only for galaxies classified in the BPT diagram as star forming and low-S/N star forming (§). For other classes of galaxies (AGN, composite, galaxies with no H$\alpha$) for which either there is no emission line detection or the lines are contaminated by a non-SF ionizing source, B04 use the relation between the D4000 spectral index and the specific SFR (SFR normalized by stellar mass) that has been calibrated using the *star-forming* galaxies. Therefore, for these galaxies, emission-line SFRs are used indirectly. This is another reason why we denote B04 SFRs as "H$\alpha$" with quotation marks. We will return to this calibration in § . The above procedure gives only the SFRs within the $3''$ aperture of SDSS spectroscopic fibers, i.e, the *fiber* SFRs. In order to obtain *total* SFRs, an aperture correction is required. The total SFR is the sum of the SFR within the fiber and outside of it. To estimate SFR outside of the fiber, galaxies are first divided into a color-color grid, based on their *fiber* colors. Then, for each color cell, a distribution of $i$ band-luminosity normalized SFRs (SFR/$L_i$) is constructed by adding up all galaxies of SF or low-S/N SF class that fall within the given color cell. Then, the color of a galaxy *outside* of fiber is used to select the appropriate color cell, and the SFR/$L_i$ in that cell is used with the $i$-band luminosity to derive the SFR outside of fiber. This is added to previously found SFR within the fiber to obtain the total SFR. In other words, the aperture correction procedure is based on two assumptions. First, that SFRs within fibers have the same dependency on color as the SFR outside of it, and second, that this dependency, calibrated with star-forming galaxies, applies to other galaxy classes. The aperture correction factors range from close to one (no correction) to around a hundred-fold. On average they are 0.9 dex for the entire sample and 0.6 dex for the star-forming class.

## Comparison of "UV" and "H$\alpha$" SFRs for all galaxies

The comparison between SF indicators serves the obvious purpose of providing a better understanding of each indicator, but also ensures the mutual check on the reliability of the techniques used to produce the SFR estimates—in this case, the SFRs based on UV (this study) and on H$\alpha$ (from B04).

We begin by comparing the formal errors of the two methods. Note that B04 used a sample based on the expanded SDSS DR1. Here, when we refer to B04 SFRs we mean the B04 calculations applied to SDSS DR4, and available as a part of MPA/JHU SDSS value-added catalog. In Table we show average formal errors in SFR for the entire sample, and for each galaxy class. We see that the formal errors of the two SF estimates happen to be comparable.

To facilitate the comparison, we first convert B04 SFRs that were calculated for IMF to IMF used in this study (both with 0.1-100 $M_{\odot}$ limits). From models we find the conversion factor of 1.06 (in the sense that Kroupa IMF SFRs are slightly higher). In Figure we plot B04 "H$\alpha$" SFR estimates against the "UV" for all galaxies in the sample (throughout the paper SFR always means dust-corrected SFR). Note that the "UV" SFRs exhibit 5 orders of magnitude of a dynamic range— from all but negligible rates to those approaching 100 ${\rm M_\odot yr^{-1}}$. We notice that a large number of galaxies (those with higher "UV" SFRs) has comparable SFRs. However, there is a plume of galaxies for which B04 SFRs are up to 2 orders of magnitude higher than the "UV" ones. This systematic discrepancy cannot be explained by random errors of the two methods. Recall that B04 do not derive SFRs directly for classes other than star-forming (and low S/N star-forming), but instead use relations calibrated on SF class sample. Thus, in Figure we break the comparison into galaxy classes (omitting for the moment the SF class). For the low-S/N star-forming class (upper left panel), for which B04 SFR estimates do come from the emission lines, the comparison is generally good for the majority of galaxies (apart from the overall offset that makes "H$\alpha$" SFRs higher). For the other three classes in Figure , B04 SFRs rely on calibrations made using the SF galaxies. Composite (SF/AGN) galaxies compare relatively well, but the comparison becomes significantly worse for AGN, where we can see both the significant offset for galaxies with higher SFRs, and the prominent plume for those with lower "UV" SFRs. Finally, for galaxies for which no H$\alpha$ detections were possible, the discrepancy is very large for almost all galaxies. It is this class that contributes the most to the plume seen in the comparison of all galaxies (Figure ). Given that the galaxies without H$\alpha$ detections belong almost exclusively in the red sequence (Figure ), the low SFRs rates estimated by the SED fitting ("UV") appear more realistic (see also §). Also, such high SFRs as estimated by B04 would result in a high fraction of No H$\alpha$ galaxies to have a detectable H$\alpha$ emission, which is obviously not the case. In general, it appears that the largest discrepancies between the "H$\alpha$" and "UV" rates occur for galaxies with low SFRs, especially when B04 estimate SFRs based on calibrations that employed SF galaxies.

B04 SFRs are the sum of SFRs within and outside of fiber. To trace the source of the SFR discrepancies, we revisit the calibrations used to derive B04 SFRs for non-SF classes. First, within the fiber, B04 base the SFRs for non-SF classes on the relation between the specific SFR (in the fiber) and the D4000 index calibrated using the SF and low-S/N SF galaxies. We show this calibration in both panels of Figure as the dashed line (this can be compared to their Figure 11). Then, B04 use D4000 to determine the specific SFR (and therefore the SFR itself) of galaxies of other classes. In the upper panel of Figure , we compare B04 relationship to our ("UV") relationships for the star forming (bold solid line) and the low-S/N star-forming (thin solid line) galaxies. Note that in principle we cannot compare the B04 specific SFR inside of fiber to the total specific SFRs obtained from the UV, so the comparison is more for illustrative purposes. At D4000 $<1.7$ either SF or low-S/N SF lines from "UV" compares well to B04. However, notice that there are *no* galaxies classified as star forming beyond D4000 = 1.9. Therefore, from that point on, B04 relationship is based only on low-S/N SF class which are most likely just AGN contaminants (see § ). In contrast to B04, the "UV" relationship for low S/N SF continues to decline as D4000 increases. In the lower panel of Figure we show the same B04 calibration against "UV" rates for the AGN (thin solid line) and galaxies with no H$\alpha$ detections (bold solid line). In these cases, the relationship for these two types as indicated by the "UV" is entirely different from that of the star forming galaxies, even at lower D4000 values. For clarity we omit showing the composite class. Altogether, this leads to a conclusion that the relationship between the specific SFR and the D4000 index is not unique for different classes of galaxies. We tested if this is an artifact of the fact that D4000 comes from the central $3''$, while the specific SFR for "UV" is integral. We find that the relationships are qualitatively the same for $z<0.07$ and for $z>0.12$ samples for which D4000 probes different physical sizes. Instead, the most likely explanation is that since D4000 and specific SFR are sensitive to SF activity on different timescales, the galaxies with differing SF histories will have different relationships. To much lesser extent, some of the differences could be due to metallicity .

Given the level of aperture corrections (a factor of several), the analysis given above does not provide the explanation for the large discrepancies we see when comparing the *total* SFRs. Since the SFR outside of fiber dominates in B04 total SFR, we will now focus on the calibration B04 use to determine it. As outlined in § B04 aperture correction relies on the calibration of the luminosity-normalized SFR as a function of $g-r$ and $r-i$ color. Using our UV-based SFRs we can test one of the two assumptions behind B04 aperture corrections—whether the SFR/$L$ calibration against color holds for different classes of galaxies. In Figure we plot SFR/$L$ against two colors that B04 use to determine SFRs outside of fiber. The bold line represents the relation from galaxies classified as star-forming. This is the basis for B04 calibration which they apply to other classes as well. The dashed and the thin solid lines show the relations for AGN and galaxies with no H$\alpha$. We see (left panel) that as the $g-r$ color increases, the discrepancy between SF and other classes rises, and reaches some 0.8 dex for AGN and 1.3 dex for No H$\alpha$. Similar levels of difference are present against $r-i$ (right panel). The level of these differences matches the discrepancies (plumes) in the comparison of B04 and "UV" SFRs in Figure . We can interpret these differences as reflecting the fact that the colors used for B04 calibration trace old populations, while the current SFR traces young populations, and it is not surprising that the two will differ for different classes of galaxies. Also note that the discrepancies of specific SFRs of different classes of galaxies are qualitatively similar in the case of D4000 and in the case of $g-r$ color, which is not surprising since they are both sensitive to population age on the similar timescales, and are similarly not affected by the dust. To conclude, UV-derived SFR is subject to fewer limitations, so it can be applied to a more diverse types of normal galaxies.

## Comparison of "UV" and "H$\alpha$" SFRs for star-forming galaxies

To avoid the problems indicated in the previous section, we need to compare the SFR estimates from the "H$\alpha$" and the "UV" methods in galaxies where H$\alpha$ is well-detected (S/N $> 3$) and arises predominantly from star formation, i.e., to galaxies classified as star-forming (SF). In Figure we show the B04 "H$\alpha$" star formation rates plotted against our "UV". We see that the comparison, spanning some three orders of magnitude in SFR, is exceptionally good. Formal error bars are comparable, with the "UV" being somewhat smaller (see Table ). The scatter (standard deviation of the difference) of the two measurements is 0.50 dex. When $3\sigma$ outliers are excluded, the scatter is reduced to 0.36 dex. This is very well matched by the sum (in quadrature) of the formal errors of the two methods (0.35 dex), confirming that the two measures are predominantly independent. There is an average offset between the two SFRs of only 0.06 dex in the sense that "H$\alpha$" SFR is higher (which reduces to 0.02 dex, i.e., 5%, when $3\sigma$ outliers are excluded).

However, Figure alone can potentially hide some systematic trends between the two SFR estimates. Where can the differences in the SFR estimates arise from? Deriving SFRs for either methods incorporates several steps: (a) obtaining full "H$\alpha$" or "UV" observed luminosity (which in the case of "H$\alpha$" involves aperture corrections), (b) correcting the observed luminosity for dust attenuation, and (c) converting the luminosity into a SFR. While both methods perform these three steps simultaneously, we can still design tests that would differentiate between. Given in some cases the large aperture correction applied by B04, we first plot in Figure the difference in SFR estimates (SFR residuals) with respect to the level of aperture correction applied to "H$\alpha$" rates. We see no correlation. We also check, but do not find a correlation between the residuals and the apparent sizes of galaxies (not shown). Those are good indications that the aperture corrections of B04 are quite reliable for galaxies classified as star-forming. Finally, we find that the SFR residuals exist when plotted against the stellar mass of a galaxy (shown in Figure ). We have a change of SFR residual of 0.38 dex over the $8.5<\log\, M_* <11$ range, with $1\sigma$ dispersion of 0.33–0.41 dex around the running mean. To double check if the residuals are in any way connected with B04 aperture corrections, we limit the sample to farther, and therefore on average smaller galaxies ($z>0.12$). We still find that the residuals correlate with stellar mass.

Before ruling out step (a) as being responsible for the mass-dependent residuals, we perform two additional checks. First, imagine there is an offset in zero points between *GALEX* and SDSS magnitudes (after all, they come from different surveys and are not measured in identical apertures). E.g., if *GALEX* fluxes were systematically overestimated, this could cause "UV" SFRs to be overestimated as well, possibly in such a way to preferentially boost blue, low-mass galaxies (as the trend in Figure would suggest). To test this, we perform full SED-fitting runs in which we make $FUV$ and $NUV$ fainter by 0.1, 0.2 and 0.4 mag. We first check the quality of these new fits by comparing the distribution of $\chi^2_{\rm best}$ values with the nominal (no magnitude offset) run. Runs with 0.2 and 0.4 mag offsets produce evidently inferior fits, already suggesting that any potential offset cannot be that large. The run with 0.1 mag offset, however, appears as good as the nominal run, with an even slightly smaller average $\chi^2_{\rm
best}$. However, comparing the "UV" SFRs produced with offset *GALEX* magnitudes to "H$\alpha$" SFRs, we find that while slightly flattening the slope of the residuals with respect to mass, the trend is not eliminated. In fact, the scatter of the residuals around the mean is larger when modified "UV" SFRs are used instead of nominal. We also look at the quality of fits with *GALEX* magnitudes offset in the opposite direction (making them brighter), but such fits are evidently inferior.

The second test concerns SDSS magnitudes. Namely, B04 use CMODEL SDSS magnitudes to transform fiber SFRs into total SFRs, while we use MODEL SDSS magnitudes when performing SED fitting.[^6] Therefore, we perform another SED fitting using CMODEL magnitudes instead. First, we notice that the quality of fitting with CMODEL magnitudes is noticeably inferior, stressing their inadequacy for producing reliable color estimates. Also, the use of CMODEL magnitudes does not remove the trend of SFR residuals, and the scatter of the residuals becomes larger than in the nominal run, especially at higher masses. With this test we exhaust the possibilities that discrepancies arise in step (a) above.

We now move onto step (b), i.e., correcting the observed flux for dust attenuation. Note that both B04 and this study use prescription for dust attenuation, the aim of which was to produce *consistent* treatment for H$\alpha$, and UV continuum attenuation. Is the presence of the SFR residuals an indication that this dust attenuation model does not produce fully consistent answers? In our SED fitting we keep track of estimates on $\tau_V$—the dust opacity in rest-frame $V$ band. We can thus compare our $\tau_V$ values with those obtained by B04 (which we denote $\tau_V$("UV") and $\tau_V$("H$\alpha$"), respectively). While we use the same model to constrain attenuations, B04 obtain them from the emission lines (to first order from the Balmer decrement), while we obtain them from the broad-band SED (to first order from the UV spectral slope, $\beta$, or equivalently, the UV color). Also, B04 measurements are restricted to the fraction of the galaxy inside of the fiber, which can in principle produce some systematic differences. Direct one-to-one comparison, shown in Figure , shows a rough agreement, albeit with a large scatter due to relatively large errors of both estimates (average formal errors shown as the error bar). Looking instead at the distributions of the two $\tau_V$ estimates (not shown), we find them to be very similar on the whole, with the $\tau_V$("UV") slightly offset towards the larger values with respect to $\tau_V$("H$\alpha$"). Since we find SFR residuals to correlate with the mass, it is more instructive to check whether such a trend exists in the difference of $\tau_V$ estimates from "H$\alpha$" and "UV". In Figure we plot $\tau_V$ residuals against the stellar mass. Despite the large scatter, we find a systematic trend. The gradient is 0.47 (equivalent to 0.52 mag) over the $8.5<\log\, M_* <11$ range. For lower masses, "UV" $\tau_V$ attenuation is higher than "H$\alpha$", which implies that the observed UV flux is corrected more than $\tau_V$("H$\alpha$") would require, resulting in "UV" SFRs being higher than SFR("H$\alpha$"). The situation is reversed at high masses. This trend matches the sense of the trend of SFR residuals.

In order to establish whether the difference in $\tau_V$ estimates is the dominant cause of the trend of SFR residuals, we will try to modify one of the two SFR estimates (e.g., the "UV" SFR) based on the difference in $\tau_V$ values. Note that this is an approximate technique, since our SED SFRs are not "corrected" for dust by the application of a single number, but rather by a complex application of prescription on many different populations that constitute the model SED. Nevertheless, for the purposes of this test we would assume that the attenuation in the FUV drives the correction of the "UV" SFR estimate. We therefore modify it with:

$${\rm SFR(``UV'')_{mod}} = {\rm SFR(``UV'')} \times 10^{0.4\Delta
  A_{FUV}},
\label{eqn:sfrmod}$$

where $\Delta A_{FUV}$ is the difference between the attenuation in the rest-frame FUV as implied by $\tau_V$("H$\alpha$") and that determined from the SED fitting. We look at the parameters from the SED fitting to calibrate the relationship between the attenuation in FUV and the $V$-band opacity, and find the following relationship:

$$\Delta A_{FUV} = \mu(5.4-0.84 \tau_V) \Delta \tau_V,
\label{eqn:deltaafuv}$$

where $\Delta \tau_V= \tau_V\rm{(``H\alpha'')} -
\tau_V\rm{(``UV'')}$, while $\tau_V$ and $\mu$ come from the SED fitting (i.e., they are "UV" measurements). The parameter $\mu$ is the coefficient in the model that determines the fraction of the total opacity affecting the diffuse ISM. We modify SFR("UV") according to Equations and , and in Figure re-plot the difference in SFRs against the mass. The trend of the SFR residuals with mass has disappeared, and the difference stays within $\pm 10\%$ across the full mass range. The dispersion around the running mean is 20% larger than with the non-modified SFR residuals, which is to be expected due to the approximate nature of the modification applied. Therefore, we conclude that the "UV" continuum and the "H$\alpha$" in essence produce identical answers for the star formation rate, the apparent difference stemming from different estimates of the attenuation that the two methods provide. With this in mind, we can also rule out other possibilities for the cause of the original differences in the SFRs, such as the reliability of converting the luminosity into a SFR (point (c) above). Also, these results show that, contrary to some notions, models cannot be too much off in the UV.

Although we can now account for the SFR differences, we would still like to understand the origin of attenuation differences that cause them. Related to this, we should try to determine which attenuation estimate (emission-line or the UV continuum) is more accurate. At this point we need to emphasize that our SED fitting has a non-flat prior distribution of $\tau_V$ (described in §), while B04 uses a flat prior, with attenuation taking values in the $0.01<\tau_V<4$ range. While the non-flat prior used in our study is a reasonable assumption, and the one that helps constrain the parameter space to physically realistic values, we would like to check whether the differences in priors induce any systematic effects. To that effect we first perform an additional SED fitting run, with the model libraries created with a flat attenuation prior, taking values in the $0<\tau_V<3$ range. At lower ($\tau_V<0.7$) values, we find the distribution of the resulting $\tau_V$ values to be similar to the original distribution, while at higher values it displaces the original peak, located at $\tau_V\approx 1.3$, to a broader distribution peaking at $\tau_V\approx 2.0$. The latter is the obvious consequence of assigning equal probability to models with high attenuations. In any case, we end up with a $\tau_V$ distribution that differs more with respect to the B04 $\tau_V$("H$\alpha$") than did the original $\tau_V$("UV"). Nevertheless, we proceed and check the relationship between the new SFR residuals and the stellar mass. The trend observed in the original relationship is still present, although it is somewhat weaker at low masses. Anyhow, it appears that we cannot force the two attenuation estimates to reach an exact agreement by a simple modification of the prior distribution used in SED fitting.

Next we will try to get a sense of which $\tau_V$ estimate is more realistic. Since the direct comparison with an external independent measurement is not readily available, we will try to evaluate the differences using internal relations. For star-forming galaxies we expect attenuation to be correlated with the stellar mass. In Figure we show this relation for our "UV" estimates (upper panel) and B04 "H$\alpha$" estimates (lower panel). It is apparent that the latter defines a more tight relationship. Also, $\tau_V$("H$\alpha$") spans a larger range of values than $\tau_V$("UV") (1.4 vs. 0.8 over $8.5<\log\, M_* <11$). Finally,, the dispersion around a running mean is some 20% smaller for $\tau_V$("H$\alpha$"). This all indicates that B04 $\tau_V$ estimates are probably more accurate. While the two $\tau_V$ estimates agree well on the whole (Figure ), the "UV" estimate appears more poorly constrained, which may cause it to more often take values away from the extrema, suppressing the span in relationship with the mass. As a result, the attenuation estimates for low-mass galaxies will be overestimated, and for high-mass galaxies underestimated. It is difficult to directly test this explanation. Namely, if we restrict the analysis only to galaxies with small formal errors in $\tau_V$("UV") (or alternatively small errors in $FUV-NUV$), we at the same time bias the sample to smaller values of $\tau_V$("UV"). Alternative explanation is that the $\propto \lambda^{-0.7}$ extinction curve deviates from the true extinction curve in a way that would be mass-dependent (or some quantity related to mass, such as the metallicity). A possibility that the extinction curve is mass-dependent is raised in . We would be able to acquire further insight by comparing both attenuation estimates to some external measure, such as the IR excess (IRX). have recently obtained IRX attenuations from *Spitzer* MIPS and *GALEX* observations of SDSS galaxies, therefore such sample can be used to include UV (i.e., SED) derived attenuations. Also, the SED-based attenuation estimates can possibly be improved by adding the near-IR data, such as the $JHK_s$ photometry from 2MASS survey , which would improve attenuation estimates by placing stronger constraints on the stellar metallicity.

# Obtaining the UV SFRs without SED modeling

In this work we are using sophisticated SED modeling to obtain UV-based star formation rates. In many applications such a detailed approach is not practical, or not even possible. In such cases one would like to obtain a reasonably good estimate of a dust-corrected star formation by applying some simple transformations to the UV photometry. While such methods have been used in many previous studies, here we will use the results of the detailed SED analysis to *calibrate* such simplified models, with a special emphasis on users of UV data obtained with *GALEX*. Since for the foreseeable future *GALEX* will remain the only facility capable of observing a large number of galaxies in the UV, it is not without justification to treat its filters as defining some standard photometric bands in the UV domain. Researchers who study high-$z$ galaxies in rest-frame UV can calibrate their blueshifted filters against *GALEX* FUV and NUV response curves.[^7]

Obtaining a UV-based SFR consists of K-correcting the UV magnitudes, estimating the dust attenuation of the FUV flux, and converting the dust-corrected FUV luminosity into a SFR. We assume that the user will correct their data for the Galactic extinction, and then apply some standard K-correction procedure to obtain rest-frame $FUV$ and $NUV$ magnitudes. Of course, obtaining reliable K corrections requires optical or near-IR photometry. To calibrate the relations in this section, we use KCORRECT V4_1_4 (§).

We start by using the results of the full SED fitting to calibrate the well-known correlation between the attenuation and the UV spectral slope (e.g., ). This is often refereed to as IRX-$\beta$ relation, although strictly speaking IRX indicates IR excess, which is *correlated* with UV attenuation. IRX-$\beta$ relationship for *normal* galaxies, such as those in our sample, has been previously studied with *GALEX* data by . In Figure we show FUV attenuations ($A_{FUV}$) obtained from our SED modeling, for galaxies classified as star-forming, plotted against the rest-frame UV color, which is linearly correlated to the UV spectral slope. We see that the majority of objects lies along the ridge. This confirms that there exists an "IRX-$\beta$" relationship for normal galaxies, not just for starburst galaxies as usually assumed. To quantify the relationship for our sample, and thus allow $A_{FUV}$ to be determined from *GALEX* observations of normal galaxies, we fit a linear function to running medians (with value from each 0.05 mag bin weighted equally). Using the medians is necessary in order to avoid the fit to be affected by numerous outliers. Also, we find that past some red UV color the attenuation does not seem to rise, so we adopt a constant value. In any case, $A_{FUV}$ exhibits a large scatter for such red colors.

$$A_{FUV} = \left\{\begin{array}{ll}
           3.32\, ^0(FUV-NUV)+0.22, &  ^0(FUV-NUV)<0.95 \\
       3.37, &  ^0(FUV-NUV)\geq0.95, \end{array} \right.
\label{eqn:afuv1}$$

where $A_{FUV}$ is in magnitudes, and superscript 0 designates rest-frame colors. This relation is plotted in Figure as a solid line. In many applications, especially at higher redshift, spectroscopic classification may not be available. In those cases one can select star-forming galaxies based on their color. For blue-sequence galaxies ($^0(NUV-r)<4$, and without applying any class selection) one should use a slightly modified relation:

$$A_{FUV} = \left\{\begin{array}{ll}
           2.99\, ^0(FUV-NUV)+0.27, &  ^0(FUV-NUV)<0.90 \\
       2.96, &  ^0(FUV-NUV)\geq0.90, \end{array} \right.
\label{eqn:afuv2}$$

shown in Figure as a dashed line. Note that these relationships are optimized for *GALEX* bandpasses *and* for normal star-forming galaxies. We overplot relation as a dotted line (we used transformation to obtain $FUV-NUV$ from the UV spectral slope). Apparently, this relation, constructed from a sample of starburst galaxies, is well above the majority of galaxies in our sample. and , using smaller samples (RC3 and cluster galaxies, respectively) also find that relation overpredicts FUV attenuation. derive a relation that exactly bisects our SF fit and the relation. Also note that these relations pertain to the optically selected sample (since practically all SF galaxies are detected in UV bands, there is no additional UV selection). In general, samples selected at different wavelengths will have different IRX-$\beta$ relations, as demonstrated by , using *GALEX* NUV and *IRAS* 60$\mu$m selected samples.

For a source with known redshift, the dust-corrected, rest-frame $FUV$ is converted into a FUV luminosity. The final step requires converting the luminosity into SFR. Following B04 notation, we define:

$$\eta^{0}_{FUV} = L^{0}_{FUV}/{\rm SFR(``UV'')},$$

to be the inverse conversion factor between a dust-corrected rest-frame FUV luminosity (in erg s$^{-1}$ Hz$^{-1}$) and the SED-derived SFR. The conversion factors comes from the stellar population modeling that was used to perform the SED fitting. Like the equivalent conversion factor for H$\alpha$ luminosity, $\eta^{0}_{FUV}$ is sensitive to metallicity, albeit more weakly so. For our sample (the metallicity of which is on average $0.8 Z_{\odot}$), the median conversion factor is:

$$\log \eta^{0}_{FUV} = 28.165.\label{eqn:conv1}$$

The above factor is given for IMF, which is used in this paper. The often-used conversion factor given by for the IMF is:

$${\rm SFR}(M_{\odot} {\rm yr^{-1}}) = 1.4 \times 10^{-28} L_\nu ({\rm
  erg^{-1} s^{-1} Hz^{-1}}).\label{eqn:rob}$$

We find using models that the appropriate transformation factor between UV-derived SFRs that assume Chabrier and Salpeter IMFs is 1.58 (both with 0.1-100 $M_{\odot}$ limits)[^8]. Thus our "empirical" conversion factor (Equation ) for Salpeter IMF becomes:

$${\rm SFR}(M_{\odot} {\rm yr^{-1}}) = 1.08 \times 10^{-28} L^{0}_{FUV}({\rm
  erg^{-1} s^{-1} Hz^{-1}}). \label{eqn:conv2}$$

This implies that the conversion factor (Equation ) is 30% higher. We verify that the effect of the difference in bandpasses ( factor was given for the 1500–2800 Å range, while the *GALEX* FUV filter spans 1300–1800 Å) is completely negligible. The actual reasons for the difference are two-fold. First, even when obtaining the *model* conversion factor using and the *same* assumptions as given by (i.e., solar metallicity and a constant star formation history), we still find conversion to be 15% higher (i.e., the corresponding $\eta$ is lower than one produced by ). Further, we have that our sample has average metallicity somewhat lower than the solar. This accounts for another 5% difference. Finally, the remaining 10% difference stems from the fact that our sample has a variety of star formation histories, giving $\eta$ that is on average higher than that for constant SF history. To summarize, for optically selected samples similar to ours, we suggest using "empirical" conversion given in Equation or .

Finally, we combine the above ingredients and "predict" the SFR from K-corrected $FUV$ and $NUV$ magnitudes. In Figure we compare the SFR obtained using this simple prescription to the original SFR("UV"). We observe a good agreement for the majority of objects (the greyscale scatter plot). In addition, for each bin in SFR("UV") we show a 16–84 percentile range of the predicted SFR values, i.e., their conditional distribution. This allows us to track the trend in the regions where there are few points to show on the greyscale scatter plot. We see that the predicted SFR starts to deviate at both the low and the high SFR ends, by up to 1 dex. Good agreement is obtained in the $-1<\log\, {\rm SFR} < 1$ range, i.e., where the majority of galaxies lie. The reason for the deviations at low and high-SFR end stem from the inability to estimate $A_{FUV}$ reliably only from $FUV$ and $NUV$ for galaxies having extreme SF. No simple modifications of Equations and are able to reproduce $A_{FUV}$ estimate of galaxies with extreme SFRs without destroying the agreement for the majority of normal galaxies. In $-1<\log\,{\rm SFR} < 1$ range, the 16–84 percentile range is 0.35 dex wide, indicating that the predicted values have an equivalent $1\sigma$ scatter of 0.17 dex. Therefore, the total error of the predicted SFRs in this range is 0.26 dex (with the error of SFR("UV") itself added in quadrature). This is a small penalty considering the ease of the method. Note that the prescription given here is optimized to produce good results for the large majority of galaxies, and that a result for any given galaxy should be used with caution.

# Star formation in the local universe

## Star formation history and stellar mass

It is now commonly accepted that the star formation history of a given galaxy depends strongly on its mass (for early work based on UV see ), and to some extent its environment. To first order, we can characterize the star formation history using the specific star formation rate, which is defined as the ratio of the current SFR to the current stellar mass. Thus, higher values of the specific SFR indicate that a larger fraction of stars was formed recently. The unit of specific SFR is inverse time, and it is thus often referred to as the galaxy build-up time. However, this should not be confused with an actual galaxy age, since the specific SFR tells us only how long it would have taken to build a galaxy assuming it had a current SFR throughout its lifetime, and assuming that the current stellar mass is close to the total mass formed in stars during a galaxy's lifetime (the two are not the same because of gas recycling).

For each galaxy in our sample we construct a two-dimensional probability distribution function in $(\log\,({\rm SFR}/M_*),\log\,
M_*)$ plane as outlined in §. We use full 2-D PDF instead of collapsing it into a single value in order to better represent low-SFR galaxies whose PDFs are often very wide (i.e., not well constrained) in specific SFR, and in which cases the single value, such as the average of the marginalized PDF, would tend to artificially concentrate in the middle of the PDF. In Figure we show the co-added 2-D PDFs in 0.1/0.05 dex wide bins. In the upper panel of Figure , as in all previous figures in the paper, we show unweighted values, i.e., each galaxy in our sample contributes equally to the greyscale density. Given our redshift and magnitude limits, we adopt the $r$-band absolute magnitude limit of $M_r=-14.75$ (which we verify by constructing a luminosity function). This luminosity limit leads to the completeness in stellar mass of $\log\, M_* \approx 7.5$, which is below the plotted mass range. The dashed line is shown for reference and represents a constant SFR of $1\, M_\odot \rm {yr^{-1}}$. We see that lower-mass galaxies appear to be confined in a relatively narrow sequence which declines as the mass increases. At masses above $10^{10}\, M_{\odot}$, some galaxies assume a much lower value of the specific SFR, i.e., the sequence turns over. In the lower panel, each galaxy's PDF is weighted by its $V_{\rm max}$ value, i.e., the sample is volume-corrected. Since the volume-corrected values span a very wide range of values, we now show logarithms of number densities, scaled so that the full greyscale range displays densities from maximum (black) to 2 dex below the maximum value (lightest shade of gray). Each horizontal slice in effect represents a mass function of galaxies of a given specific SFR. It now becomes apparent that the sequence seen in the unweighted plot extends to lower masses, with an increasing number density of galaxies. Also, the mass distributions at low specific SFR extends to lower masses than implied by the non-corrected plot. Considering the size of the sample, the accuracy of SFRs, and the sensitivity to low levels of SF, this figure represents the best available assessment of the star formation history in the local universe. The typical error in stellar mass is only 0.07 dex, and 0.40 dex in the specific SFR (note that the errors in specific SFR vary by a factor of three across the sample). We complement the above figure with Figure (upper panel), where we display *conditional* probability distribution in each 0.05 dex wide mass bin, which is equivalent to normalizing each column to its maximum value. Such representation allows us to see what specific SFR dominates at a given mass, and it also displays more clearly the behavior at masses where there are too few objects to show on the standard plot. At masses below $10^{10}\,
M_{\odot}$, high specific SFRs dominate. At $M_*=10^{10}\, M_{\odot}$, more galaxies assume values below the dominant high-specific SFR "sequence", down to very low levels. At $M_*=2\times10^{11}\,
M_{\odot}$, close to where the high-specific SFR "sequence" terminates, the galaxies with *low* specific SFRs begin to dominate. Above $\log\, M_*=11.8$ the plot is dominated by noise from very few objects with such high mass. This figure compares well to Figure 24 in B04.

## The star-forming sequence

We now focus on star formation histories of various classes of galaxies. In Figure we show specific SFR against the stellar mass for the star-forming, AGN (together with SF/AGN composites), and the class without H$\alpha$ detection. For each class the upper panels show nominal, unweighted data, while the lower panels are volume-corrected. Dashed line has the same meaning as in Figure . First we notice that the three classes occupy relatively distinct portions of the parameter space. This is especially pronounced in the unweighted plots. Thus the three classes appear to have had quite different star formation histories. This was to some extent indicated in the CMDs (Figure ), but is more striking now.

The star forming class (SF), which forms a blue sequence in the CMD (Figure ), dominates the high values of the specific SFR. The sequence that was obvious in Figure stands completely isolated here. In the weighted plot we see an ever-increasing density towards the lower mass, which reflects a rising low-mass slope of the mass function of late type (blue) galaxies. While the typical SFR *rises* from $0.1\, M_\odot {\rm
yr^{-1}}$ at the low-mass end to $10\, M_\odot {\rm yr^{-1}}$ at $10^{11}\, M_{\odot}$ (consistent with the assumption that the more massive galaxies contain more gas and therefore have higher SFRs), the specific SFR actually declines by a factor of 10. The tightness of the "star forming" sequence (note that the formal error of the specific SFR is only 0.20 dex) represents an important indication that the galaxy's mass regulates the overall star formation history (see also ). It appears reasonable to assume that in the absence of an event that may disturb galaxy's gas reservoir, a galaxy would "naturally" sit on this tight sequence. This appears contrary to the suggestions that the red (non-SF) sequence forms by simple gas exhaustion, since we would then have a much wider range of (specific) SFRs at the given mass, for galaxies observed in the various stages of gas exhaustion. Note that for masses below $\log\,
M_* = 9.5$ we have an additional tail of high specific SFRs (best seen in volume-weighted plot), which is indicative of *starbursts*. We will return to the width of the SF sequence in §.

## Specific SFRs and stellar masses of galaxies with an AGN

We now turn to galaxies that lie on the AGN branch of the BPT diagram. These represent both galaxies classified as SF/AGN composites, and "pure" AGN (those lying in the region of the BPT diagram where non-SF ionizing source is necessary). We have already seen that a significant fraction of AGN hosts populates the green valley of the UV-to-optical CMD (Figure )—the intermediate region between the blue and the red sequences. (However, note that in the optical CMD they would lie predominantly in the red sequence). A possible reason for the presence of these galaxies in the green valley is that they are entering the stage in which their star formation is being shut off, and the colors are becoming red as a result. , using the *GALEX* data, has demonstrated that if the SF of galaxies lying in the green valley is indeed being quenched, then their color evolution implies a mass flux (transfer from blue to red) that roughly matches the mass increase required to populate the red sequence between the epochs $z\sim 1$ and $z \sim 0$ .

In the middle panels of Figure we present specific SFR against mass for AGN and the SF/AGN composites combined. The typical error of the specific SFR is 0.35 dex. Galaxies with AGN populate the region *between* the SF galaxies with no AGN and the galaxies with no H$\alpha$ (adjacent panels), in a more striking manner than even in the UV-to-optical color-magnitude diagram. We note that these AGN are of Type 2 (Seyfert 2s and LINERs), which are not expected to affect the stellar continuum, especially not in the integrated galaxy light (, and reference therein). In galaxies having a luminous Type 2 AGN, limit the AGN continuum contribution to less than a few percent in $r$ band. Assuming a power-law AGN continuum of the form $\propto
\nu^{-1.5}$, this limit corresponds to less than $\sim 15\%$ contribution in the UV. Another indication that the UV contribution from an AGN is low in the UV comes from the fact that the quality of our SED fitting for galaxies with AGN is nearly the same as for the similar galaxies with no AGN. Also, explore *GALEX* imaging of a more nearby sample of AGN hosts and find that their UV emission is quite extended, thus unlikely to emanate from the central AGN.

The fact that most AGN have intermediate star formation histories, and that this region is *not occupied* by galaxies unless they have an AGN, is remarkable in the light of recent suggestions that the AGN feedback may be regulating the star formation of massive galaxies . Still, we would like to know if there is a connection of galaxies hosting AGN with normal star-forming galaxies with no AGN. In Figure , we plot SF class (thin solid contours), and separate AGN into "pure" AGN (dashed contours) and SF/AGN composites (bold solid contours). Composites form the lower part of the AGN branch in the BPT diagram, and the "pure" AGN the upper part. Contours encompass 10, 30, 50 and 70% of each class's objects (i.e., their PDF densities), with no $V_{\rm max}$ weighting. We notice that composites "bridge" the normal SF population and that of "pure" AGN, i.e., they have values of both the mass and the specific SFR that are intermediate between the two classes. AGN, in addition, have a tail extending to very low values of the specific SFR. Another, more physically motivated way of splitting the AGN hosts is based on the strength (i.e., the accretion rate) of the AGN itself, as indicated by its \[OIII\]$\lambda$`<!-- -->`{=html}5007Å luminosity. We use dust-corrected \[OIII\]$\lambda$`<!-- -->`{=html}5007Å luminosities, and divide AGN hosts into strong and weak AGN, using their demarcation of $L$\[OIII\] = $10^7
L_{\odot}$. We apply this division on all galaxies lying on the AGN branch of the BPT diagram (we refer the reader to for a visual presentation of where AGN of various strengths lie on the BPT diagram). We show in Figure the relation between the specific SFR and mass for SF galaxies (thin solid contours), strong (bold solid contours) and weak (dashed contours) AGN. Apparently, strong AGN lie on the massive *continuation* of the SF sequence, with SFRs as *high* as those of the massive galaxies of the SF class (i.e., the offset from the dashed line is similar). Higher masses with similar SFRs makes their specific SFRs take more intermediate values. Weak AGN fall at lower star formation rates relative to the strong AGN, with a tail toward very low star formation rates that extends well into the domain of the quiescent galaxies (see below). , using the optical spectroscopy, found that strong AGN have younger populations than weak AGN, which we qualitatively confirm, and provide a closer connection between strong AGN and normal SF galaxies.

Here, we propose an evolutionary progression for massive galaxies: from SF galaxies with no AGN, to strong AGN that are more massive yet with similar SFRs, to weak AGN that have masses similar to strong AGN but have lower (more suppressed) SFRs. Thus, the weak AGN would represent a population with fading star formation. We emphasize that while most optical and structural properties of the majority of AGN would place them in the optical *red* sequence, we find that they represent a smooth, massive extension of the star forming sequence defined by optically *blue* galaxies. This picture is supplemented by the results of . They confirm that the majority of blue $NUV-r$ galaxies with high velocity dispersion host AGN, and suggest that strong AGN require extended gas-rich disks to provide a reservoir of gas that is fueled into an AGN. Previously, without the UV data, suggested that many type 2 AGN are post-starburst systems. Such a possibility, by itself, allows at least two different evolutionary scenarios. One in which a normal red-sequence galaxy experiences a burst (say, through a minor gas-rich merger) which also fuels an AGN, and another, in which the AGN is originally concurrent with normal SF. Our results, showing a smooth sequence that begins with SF galaxies without AGN and extends at its massive end to AGN hosts with different levels of star formation, appear more compatible with the second scenario.

While neither nor our study directly proves that AGN are actually responsible for the subsequent quenching of the star formation, they do provide more evidence for the connection between the AGN and the fading star formation in massive galaxies. Such picture is qualitatively compatible with recent theoretical work based on $N$-body simulations of galaxy evolution and semi-analytic treatment of feedback. Specifically, the "radio"-mode AGN feedback, proposed by , demonstrates that the galaxies situated in dark halos above some mass start to suppress the accretion of gas, and at the same time suppress further star formation. These are called "radio" AGN because they can have black hole accretion rates well below Eddington, as is the case in our sample. Also, this mechanism is efficient in very massive galaxies. From the weighted plot of AGN (lower middle panel of Figure ) we find that the mass function of AGN peaks at around $4\times 10^{10} M_{\odot}$, and has an abrupt turn-over at lower masses. This was noted by who also demonstrate that the low-mass cutoff is not likely to be a selection effect; if a low-mass galaxy contained anything but the very weakest AGN it would be classified as an AGN in the BPT diagram. The connection that *GALEX* has provided between the SF galaxies and the AGN feedback is quite reassuring, and evidently further work is required to complete this picture.

## Specific SFRs and stellar masses of galaxies with no H$\alpha$ emission

Finally, we concentrate on galaxies for which no H$\alpha$ emission was detected in SDSS spectra. These galaxies have red *optical* colors, which puts them in the "red sequence" category. These galaxies are often referred to as "old, red and dead", because of the presumably old populations and the lack of current SF. However, in the UV to optical colors (Figure ) we saw that some of these galaxies have relatively blue colors (which indicates that one needs to be careful when using the term such as "red sequence", since it apparently depends on the color used). In the right-hand panels of Figure we show specific SFRs of galaxies with no H$\alpha$. In the unweighted plot they occupy a region of the parameter space distinct from the star-forming and AGN galaxies, with specific SFRs typically lower than for those two classes. The volume-corrected plot still puts the bulk of them lower than the other two classes, but the overlap with the low specific-SFR tail of the AGN is more prominent. Also, the mass extent is now displayed accurately, with the mass function turning over and flattening at lower masses.

If we expect most red-sequence galaxies to have no star formation at all, how do we interpret the (specific) SFRs derived here? First, we note that the specific SFRs (and also just the SFRs) are more poorly constrained than for any other class (average formal error in specific SFR is 0.63 dex). This should be interpreted in a way that the observed SEDs are in many cases compatible with a very wide range of SFRs, mostly very low ones. Let us note that our model SF histories do not contain truncated models, but only exponentially declining models (with added starbursts), which can never reach exact zero SFR. The dynamic range of SFR in models is determined by the range of exponential timescales, which in our case have a minimal e-folding time of 1 Gyr. This produces a lower limit of $\log\, ({\rm
SFR}/L)\approx -14$, which is roughly equivalent to the lower limit in the specific SFR shown in Figure . While many galaxies in this class are compatible with this lower limit, they do not cluster close to it, because in the absence of stronger constraints, higher SFRs will also be allowed (i.e., to firmly constrain that the SFR is, say, smaller than $0.001\, M_{\odot} {\rm
yr^{-1}}$ one needs to measure a very low amount of UV light with a very high precision). However, there are clear cases in which a galaxy classified as having no H$\alpha$ has a relatively well constrained SFR. For some of them it is as high as $1\, M_\odot {\rm yr^{-1}}$ (crossing the dashed line). What is the morphology of galaxies with high (specific) SFRs, yet no H$\alpha$ emission? There are 837 galaxies with $\log\, M_*>10$ and $\log\, ({\rm SFR}/M_*)>-11$. We inspect SDSS images of these galaxies. Most appear as normal early-type galaxies (with several obvious E+A types), which was already indicated in who find that some morphologically and spectroscopically early-type galaxies show blue UV-to-optical colors. However, one quarter shows some sign of structure or a disturbed light profile. This fraction would very likely be higher in deeper images or images with higher spatial resolution. One quarter of galaxies with structure show traces of spiral arms. In most cases these spiral arms are faint, and possibly represent cases where the star formation was recently shut down (leading to the lack of emission lines), yet the residual UV light and the spiral structure remain. Such galaxies have been recognized in galaxy clusters as "passive spirals" , and are believed to represent a transition between a spiral and an S0 galaxy. In the inspected sample these galaxies are found at all redshifts, indicating that the reason H$\alpha$ is not detected is most likely not due to the aperture size being smaller than the size of the galaxy from which the integral flux is measured (especially for more nearby galaxies). Many of the galaxies with structure exhibit rings—from very faint ones to bright blue ones. In some of those cases SF could be occurring outside of the spectroscopic fiber, leading to non-detection of H$\alpha$.

One can be concerned that our models attribute the UV light from old populations (the UV "upturn") to young stars, thus implying an artificially high star formation rate. models do include post-AGB phases, but not the extreme blue horizontal branch, which is believed (by no means unanimously) to be responsible for the UV upturn in early-type galaxies . However, explore nearby "pure" elliptical galaxies and find that both the UV and the optical colors that they observe are fully reproduced in the models that we use in this work, and that if the UV excess from the actual star-formation is present, it is discernible from the combination of UV and optical colors. Therefore, our models do appear to account for UV light from old populations and therefore some of the objects with no H$\alpha$ emission are quite likely to have had actual recent SF. This conclusion is also reached in other studies based on *GALEX* data , where the blue UV-to-optical colors of morphologically selected early-type galaxies are interpreted as signs of recent star formation. It is quite conceivable that these galaxies are in a post-starburst phase, so that the H$\alpha$ emission is no longer present, while the traces of star formation are still visible in the UV. In evolutionary terms, as in the case of AGN, a question remains whether this recent SF is the sign of a previously normal SF that is now quenched, or is it the result of a small SF episode in an already quiescent galaxy. In the latter case, the starburst could be due to a minor merger with a gas-rich galaxy. Most likely both processes are present, and the forthcoming work will address the relative importance of each.

## Relationships between specific SFR and mass

One would like to quantify the trends of the specific SFR trends with respect to mass in order to facilitate easier comparison of our optically-selected, local-universe sample to studies at other redshifts, or samples selected in a different manner (keeping in mind that a conversion may be required for a different choice of IMF or $H_0$).

If one can identify "pure" star forming galaxies (with no AGN contribution) using the BPT diagram, then the star forming sequence is best described with a linear fit:

$$\log\, ({\rm SFR}/M_*) = -0.35(\log\, M_*-10)-9.83.
\label{eqn:ssfr_mass_sf}$$

In many instances one does not have all the requisite spectral lines to make an exact classification, but instead relies on color to separate star-forming galaxies. Adopting a color cut of $NUV-r=4$, which passes in between the blue and the red sequence, and taking galaxies on the blue side, we encompass all "pure" star forming galaxies, but also many SF/AGN composites (see Figure . In that case, the star-forming sequence is better described with a piecewise linear fit:

$$\log\, ({\rm SFR}/M_*) = \left\{\begin{array}{ll}
                       -0.17(\log\, M_*-10)-9.65, & \log\, M_*>9.4 \\
                       -0.53(\log\, M_*-10)-9.87, & \log\, M_*\leq 9.4.
\end{array} \right.
\label{eqn:ssfr_mass_blue}$$

This relation differs from Equation mostly because the high-mass end is somewhat lower—a consequence of including SF/AGN composites which have lower specific SFRs than pure SF galaxies of the same mass. In both cases we have performed least-square fitting to the *conditional* volume-corrected composite PDFs (analogous to the upper panel of Figure ), binned into running modes in the $8.1<\log\, M_*<11.7$ range for SF galaxies, and $8.0<\log\, M_*<11.6$ range for blue galaxies.

The trend of the specific SFR with respect to mass for the *entire* sample is best observed in Figure . We find that the parametrization having the form similar to luminosity function provides an excellent overall fit to this trend. Recently, have used Schechter functions to describe the upper *envelopes* of their SFR/$M_*$ vs. $M_*$ data points, but here we find that it represents a good description of the *main* trend. We use the Schechter function of the form:

$${\rm SFR}/M_* = ({\rm SFR}/M_*)_0\, 10^{(\alpha+1)(\log\, M_*-\log\, M_0)}
\exp[-10^{(\log\, M_*-\log\, M_0)}],
\label{eqn:schechter1}$$

where $({\rm SFR}/M_*)_0$, $M_0$, and $\alpha$ describe the normalization, characteristic stellar mass, and the low-mass slope, respectively. Taking modes of the conditional distribution in the $8.0<\log\, M_*<11.9$ range, we find the following best-fitting parameters:

$$\begin{array}{l}
({\rm SFR}/M_*)_0 = 5.96\times10^{-11}\,{\rm yr^{-1}}, \\
\log\, M_0 = 11.03, \\
\alpha = -1.35.
\end{array}
\label{eqn:schechter2}$$

This fit, and the data points that were used to construct it are shown in the lower panel of Figure .

In the paradigm of "downsizing", at higher redshifts the SF was more active in galaxies with larger mass. This should be reflected in the relative flattening of the SF sequence at higher redshifts. While most current studies observe an overall shift towards higher specific SFRs, the slope does not appear to change much (e.g., ). Our quantification of the specific SFR should enable more direct comparison with the local universe. In many cases, simple overplotting of Equations and to other sets of data would be indicative (unless data are color-selected to include only the blue sequence, in which case we recommend using Equation ).

One often wishes to determine the average specific SFR of their sample. The average specific SFR can be defined in several ways, but the most used definition is one in which the total SFR in a volume-limited sample is simply divided by the total stellar mass. We call this the *cosmic* specific SFR:

$${\rm SFR}/M_*({\rm Cosmic}) = \frac{\sum ({\rm ^{0.1}SFR}/V_{\rm max})}
{\sum ({M_*}/V_{\rm max})},$$

where ${\rm ^{0.1}SFR}$ is the SFR corrected to $z=0.1$ (see §). From our sample we find ${\rm SFR}/M_*({\rm
Cosmic}) = 6.8\times 10^{-11}\, {\rm yr^{-1}}$, which is in very good agreement with the equivalent measure in B04. How much is the estimate of the cosmic specific SFR dependent on the limiting mass of the sample? In Figure we plot cosmic specific SFR as it would have been determined had we had a mass limit (and assuming the same magnitude limits). Right-hand $y$-axis shows a fraction of the "actual" (no mass limit) cosmic ${\rm SFR}/M_*$ at the given mass limit.

Finally, we ask with what precision can one use the relations to "predict" specific SFR from stellar mass? For blue galaxies, we find the average scatter around the mean sequence to be 0.5 dex. This result is only slightly affected after taking out in quadrature the formal measurement error of 0.2 dex. Therefore, the sequence has an intrinsic width, and is likely to depend on other factors besides the mass. The width itself is a function of mass. It declines at the rate of $-0.11$ dex$^{-1}$, i.e., the scatter of the specific SFR at $\log
M_* =8.0$ is 65% larger than at $\log\, M_* =10.5$. This is an indication that the star formation becomes less stochastic as the mass increases, which is what one expects considering that massive galaxies have a larger number of star-forming regions than dwarf galaxies.

# Global star formation density

In this section we use our SFRs to determine the global star formation density at $z=0.1$, the mean redshift of our sample. In order to take into account that the sample has a spread of redshifts, and therefore probes a range of cosmic star formation history, in the remainder of the section each galaxy's SFR estimate is "corrected " to the mean redshift of the survey, where the uncertainties in the correction are minimized. To make the correction, one assumes a cosmic star formation density evolution of the form $(1+z)^{\beta}$, where B04 take $\beta=3$. Therefore, the corrected SFR for a galaxy lying at redshift $z$ is ${\rm ^{0.1}SFR} = [1.1/(1+z)]^{\beta}\, {\rm SFR}$.

Since we will be comparing our results to those obtained by B04, we first describe their methodology and apply it to our data. B04 obtain star formation density ($\rho_{\rm SFR}$), and its formal (random) error, by Monte Carlo method applied to SFR PDFs in the following way. First, 100 different bootstrap samples are constructed from the original sample, using drawings with replacements. Bootstrap samples have the same number of objects as the original sample. A random number is then used to select a value from the cumulative PDFs for each galaxy in the bootstrap sample. The random sampling of the PDF is repeated 30 times. Each SFR PDF is weighted by $1/V_{\rm max}$ and summed to get $\rho_{\rm SFR}$. From the ensemble of such samplings the mode and the 16% and 84% confidence levels are determined for the entire sample, and for each galaxy class. The results of this exercise are presented in B04 Table 3. The random error of the total $\rho_{\rm SFR}$ is found to be around 0.018 in units of $10^{-2}\,
h_{70}\, M_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}$. For our sample, we perform only the multiple probing of each galaxy's cumulative PDF, without bootstrapping the sample. Performing 100 Monte Carlo summations, we find a random error of the total $\rho_{\rm SFR}$ to be 0.017, in the same units. Therefore, the two studies have comparable formal random errors, which is not surprising given the fact that the formal errors of the SFRs are comparable. In any case, these random errors are smaller than some systematic uncertainties that we will explore later.

```latex
\begin{deluxetable}{lrrrrr}
\tablewidth{0pt}
\tablecaption{Comparison of SFR densities at $z=0.1$ in this study and
  in B04, derived using compatible methodology}
\tablehead{ Class & $\rho_{\rm SFR}$ & \% &
  $\rho_{\rm SFR}$ & \%  & fractional \\
 & This paper & of total & B04 & of total & change }

\startdata
All             &  1.953  &  100.0\%  &  1.980  &  100.0\%  &  0.99\\
SF              &  1.395  &   71.4\%  &  1.013  &   51.2\%  &  1.38\\
SF (low S/N)    &  0.296  &   15.2\%  &  0.455  &   23.0\%  &  0.65\\
Comp            &  0.147  &    7.5\%  &  0.207  &   10.5\%  &  0.71\\
AGN             &  0.076  &    3.9\%  &  0.084  &    4.2\%  &  0.90\\
No {H$\alpha$}&  0.038  &    2.0\%  &  0.221  &   11.2\%  &  0.18\\
\enddata

\tablecomments{For this comparison we adopt $\beta=3$ SFR density
  evolution. SFR density is given in units of $10^{-2} h_{70}
  M_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3}$}
\label{tab:sfrd_comp}
\end{deluxetable}
```

How do the actual values of SFR density compare? First, we notice that the B04 method is essentially equivalent to the summation of the medians of the PDFs since the average value of a random number that probes the cumulative PDF is 0.5. In Table we compare the values obtained from the summation of the medians of our SFR PDFs, to the average of the 16 and 84 percentiles of B04 values (which should be close to the median). Our estimate includes the correction for galaxies excluded from the sample because of the poor SED fits (§). These corrections are 1.6–2.8%, depending on the class. Looking at the first row, we see that our estimate of the total SFR density is very close to that of B04. However, this is a mere coincidence when we look at how the densities compare by class. Our SFR density for the SF class is 38% higher; however, for all other classes the estimate is lower—between 10% less for AGN to 5.7 times less for galaxies without H$\alpha$. The causes for these differences can be traced from the analysis of SFRs in §. Consequently, the contribution of each class to the total SFR density has changed as well. SF class now contributes some 70%, while it was around 50% in B04. This comes at the expense of other classes (except for the AGN), and most significantly from much lower contribution of galaxies without H$\alpha$ (2%, vs. 11% in B04).

We next describe the calculation of our SFR density, which somewhat differs from B04. Direct summation that was used to find SFR density in B04 assumes that the losses from the incompleteness of the sample are insignificant. B04 estimates them to be between 1 and 2 percent. We fit a Schechter fit to the number density of galaxies of various classes and integrate the luminosity density down to our completeness limit of $M_r=-14.75$, and to a limit which is 10 magnitudes fainter. We compare the two total luminosity densities and find that 1.1% of luminosity density is below our limit for SF (and low-S/N SF) galaxies. For other classes, the loss is practically zero. Therefore, for our SFR density estimate we perform the direct summation by taking galaxies $M_r<-14.75$, and then correcting the SFR density. This result differs insignificantly from one in which we would sum all galaxies with no magnitude limit. Our SFR density estimates are obtained by summing the *averages* of the SFR estimates. We prefer summing averages over medians because it is equivalent to summing an entire PDF. Averages produce a result that is 5% higher than the medians. Also, we adopt the exponent of the SFR density evolution $\beta=2.5\pm 0.7$, determined by from *GALEX*-VVDS measurements.

```latex
\begin{deluxetable}{lrr}
\tablewidth{0pt}
\tablecaption{SFR densities at $z=0.1$ determined in this study}
\tablehead{ Class & $\rho_{\rm SFR}$ & \% of total }

\startdata
All             &  1.828  &  100.0\%  \\
SF              &  1.328  &   72.6\%  \\
SF (low S/N)    &  0.277  &   15.2\%  \\
Comp            &  0.136  &    7.5\%  \\
AGN             &  0.067  &    3.6\%  \\
No {H$\alpha$}&  0.021  &    1.1\%  \\
\enddata

\tablecomments{We adopt $\beta=2.5$ SFR density evolution. SFR density is
  given in units of $10^{-2} h_{70} M_{\odot} {\rm yr}^{-1} {\rm
  Mpc}^{-3}$. Random error of SFR density for 'All' galaxies is 0.017 in
  same units. See text for the discussion of systematic errors.}
\label{tab:sfrd}
\end{deluxetable}
```

The results, and the breakdown by class, are presented in Table , and shown in Figure . These are our nominal SFR densities. Galaxies for which no classification was possible (see §) are included in the SF class. They constitute 0.8% of the total SFR densuty. Altogether, 88% of star formation occurs in "normal" SF galaxies, and another 11% in galaxies that host an AGN. Up to 1% occurs in galaxies that are quiescent based on the optical spectroscopy.

The uncertainty in the evolution exponent leads to the systematic uncertainty in the total SFR density of 2%. By far the largest uncertainty stems from the systematic uncertainties in the determination of the dust attenuation, and hence the SFR of SF galaxies (§). We estimate the degree of this error by simply replacing the SFRs of SF galaxies by those obtained by B04. This error is preferentially in the direction that the "H$\alpha$" estimate becomes 11% higher for the SF class, leading to the total density that is 8% higher. We also take into account that our SFR estimates are averaged over $10^8$ yr. Using $10^7$ yr instead, we obtain a SFR density that is 1% lower, which we will treat as the systematic error. Finally, another small source of systematic error comes from the error in the estimate of survey area of 0.3%. Adding up in quadrature random and systematic errors in the positive and in the negative directions, we arrive at our total, dust-corrected SFR density estimate at $z=0.1$:

$$\rho_{\rm SFR} = 1.828^{+0.148}_{-0.039} \times 10^{-2} h_{70}
M_{\odot} {\rm yr}^{-1} {\rm Mpc}^{-3},$$

which assumes the Chabrier IMF. For a compilation of recent measurements of the local $\rho_{\rm SFR}$ please refer to . Our $\rho_{\rm SFR}$ agrees very well ($<1 \sigma$) with four of the six estimates listed in their Table 1, but has significantly smaller error bars. Note that except for B04, all other estimates are based on samples with fewer than 300 galaxies, and many apply fixed dust correction for all galaxies. We can find the average FUV attenuation by "de-correcting" our SFRs by adding our FUV attenuation estimates to each galaxy. Such uncorrected SFR density is smaller by a factor of 6.34. In other words, the average FUV attenuation for the *volume-complete* sample is $A_{FUV}=2.01$ mag.

We next explore the distribution of global SFR density across the galaxies having different properties. In Figure we show the distribution of SFR density against the rest-frame $NUV-r$ color. The distribution sharply peaks at $^{0}(NUV-r)=2.2$. While the number density distribution is strikingly bimodal in this color, the SFR distribution is very even. The blue edge is more abrupt then the red side, which shows a tail extending to red-sequence colors, in part due to the intrinsic dust reddening. We have previously used an $NUV-r=4$ cut, which goes across the green valley, to conveniently divide the blue and the red sequences. From Figure we see that the large fraction (94%) of star formation occurs blue of this cut. The first bin sums SFR density in galaxies lacking $>3\sigma$ measurement in NUV (2% of total).

UV to optical color is easy to measure, but gives only a very rough idea of a galaxy's SF history. Therefore, in Figure we plot the SF density distribution against the specific SFR. The picture is somewhat different than from the distribution against the color. First, the high specific-SFR end does not fall as evenly as did the blue $NUV-r$ color end; instead, we see an excess of galaxies with $\log\, ({\rm SFR}/M_*) \approx
-8.3$. These are the bursty, low-mass galaxies mentioned in § . The reason that $NUV-r$ color does not distinguish this population is because it "saturates" for galaxies with $\log\,
({\rm SFR}/M_*) \approx -9$, so that all galaxies with higher specific SFR will have similar $NUV-r$ colors. FUV (included in the specific SFR) is needed to break the degeneracy. On the side of the low specific SFR (equivalent to redder $NUV-r$ colors), the decline is more steep than for the colors. The difference arises from the fact that the specific SFR, unlike $NUV-r$ color, does take out the effects of the intrinsic dust attenuation.

Finally, in Figure we show the distribution of SFR density across the galaxies with different stellar masses. The fraction of SF occurring in galaxies with increasing mass rises almost linearly. The peak of the SF today occurs in galaxies having a mass of $\approx 4\times 10^{10} M_{\odot}$. This mass is close to the "transitional" mass of , which is effectively a mass at which blue and red mass functions cross, i.e., the mass above which the red galaxies become more numerous. In the paradigm of downsizing, we would expect this distribution to move towards the higher masses at higher redshifts. This, however, is often difficult to establish since the completeness limits of many high-redshift surveys are uncomfortably close to the peak of the distribution whose position one tries to establish. In Figure we show the cumulative SFR density as a function of stellar mass, summing up from the *highest* mass galaxies. This figure shows how much SF is missed (in the local universe) by being limited to a given stellar mass (and assuming the same magnitude limits). Thus, a limit of $10^{11} M_{\odot}$ captures only 6% of the total SFR density. In order to reach 50%, one needs to probe galaxies more massive than $10^{10} M_{\odot}$. Today, 1/2 of the total SFR density comes from galaxies in the mass range $9.3<\log\, M_*<10.6$.

# Conclusions

We have considered the measurement of the star formation rate using a sample of 50,000 optically selected galaxies with *GALEX* and SDSS photometry and SDSS redshifts. We have given a detailed description of our approach to SED fitting, including the use of a large library of population models that include a physically motivated treatment of dust attenuation. These give dust-corrected SFRs ("UV" SFRs) of galaxies in the local universe. We show that there are no inherent problems that would prevent UV from serving as a reliable SFR indicator, despite past concerns over population synthesis and dust attenuation models. We recite some of the highlights from our study.

-   25 min imaging from *GALEX* in two UV bands can be used to derive SFRs with a formal accuracy of 0.2 dex for $r<18$ mag star-forming galaxies with known redshifts.

-   For galaxies for which the dominant ionizing source are young stars, (as determined using $3''$ fiber measurements and the BPT diagram) the SFRs obtained through SED fitting ("UV" SFRs) agree exceptionally well with dust-corrected SFRs obtained from modeling of the nebular lines—essentially H$\alpha$, corrected using a Balmer decrement ("H$\alpha$" SFRs), as long as the same dust model is applied to both methods. The comparison presented here is at a scale that is two to three orders of magnitude larger than the comparison of these SF indicators done in the past.

-   The dust attenuation model that we apply to "UV" SFRs and that apply to "H$\alpha$" SFRs succeeds in producing UV and H$\alpha$ dust attenuation estimates that lead to compatible SFRs.

-   Minor differences between the "UV" and "H$\alpha$" SFRs exist at the extremes of galaxy stellar mass distributions (equivalent to the extremes of the dust attenuations). They arise either from somewhat weaker constraints from the SED-derived dust attenuation, or from the mass-dependent deviations of the assumed $\propto \lambda^{-0.7}$ extinction law in the UV regime.

-   If the systematics in the attenuation estimates are accounted for, the "UV" and "H$\alpha$" SFRs agree within 10% across the entire range of galaxy masses.

-   We calibrate the relations that allow SFRs to be determined directly from K-corrected $FUV$ and $NUV$ magnitudes. They involve an FUV attenuation estimate from the UV color (equivalent to the UV slope) that we calibrated for normal star-forming galaxies, and find it to be less steep than the relation. To obtain SFRs we find that one requires a conversion factor between the FUV luminosity and SFR that is 30% lower than value. Such "simple" UV SFRs compare well with those obtained from the full SED fitting in the 0.1–10 $M_{\odot} {\rm yr^{-1}}$ range, i.e., for normal star-forming galaxies.

-   For some 50% of the optically selected galaxies from SDSS, the emission lines are either too weak (red sequence galaxies), or they are contaminated by narrow-line AGN, which precludes obtaining SFRs from H$\alpha$ directly. In those cases use calibrations between SFR and optical colors from SF galaxies, and apply them to other classes of galaxies. We find that such calibrations are not reliable and lead to large overestimation of SFRs. UV, on the other hand, provides more secure SFR estimates.

-   We present stellar-mass dependent star formation histories for various classes of optically selected SDSS galaxies, from the starbursting dwarf galaxies with specific SFRs (SFR/$M_*$) of $10^{-8}
    {\rm yr^{-1}}$ to essentially dead, red-sequence galaxies with 6 orders of magnitude lower specific SFR. We also present a volume-corrected version of these relationships which give a better account of the SF at low masses. We quantify these relations in order to facilitate comparison with similar results obtained at higher redshifts.

-   Star forming galaxies with no AGN form a very well defined linear star-forming sequence in specific SFR vs. mass plot, with the intrinsic $1\, \sigma$ width of 0.5 dex in log SFR/$M_*$. The sequence is 65% wider at the low mass end ($\sim 10^8 M_{\odot}$) than at high masses, consistent with a more stochastic (bursty) SF history for gas-rich dwarfs. The star-forming sequence terminates at $3\times
    10^{11} M_{\odot}$.

-   Galaxies whose ionizing flux (as measured in $3''$ fiber) comes primarily from young stars (star formation), those for which there is a narrow-line AGN component (mostly Seyfert 2s and LINERs), and those which are spectroscopically quiescent (no detection of H$\alpha$ lines), all occupy distinct positions in the specific SFR vs. stellar mass diagram, with the AGN hosts lying in between SF galaxies with no AGN, and the quiescent galaxies. Star forming galaxies with no AGN *do not* occupy the intermediate regions of the specific SFR—those are occupied exclusively by galaxies with AGN.

-   If we take all galaxies forming the AGN branch in the BPT diagram (again from fiber measurements) and divide them by their AGN strength (as indicated by \[OIII\]$\lambda$`<!-- -->`{=html}5007Å luminosity), we find that strong AGN represent the massive *continuation* of the SF sequence defined by galaxies with no AGN. Strong AGN have similar SFRs as the most massive galaxies with no SF, but are more massive, driving their specific SFRs to intermediate values. Weak AGN have similar masses as strong AGN, but fall at lower star formation rates relative to strong AGN, with a tail toward very low star formation rates that extends well into the domain of the optically red galaxies with no H$\alpha$. While AGN appear to be a massive continuation of the blue, star-forming sequence, their optical colors would place most of them onto the *red* sequence. We suggest an evolutionary sequence from massive star forming to quiescent galaxies via strong and weak AGN that bridge them. This scenario supports a picture of AGN feedback suppressing star formation in galaxies situated in very massive dark matter halos, the so-called radio mode AGN . Further study is needed.

-   While mostly inactive, some galaxies with no detectable H$\alpha$ in their fiber measurements exhibit SFRs as high as $\sim 1 M_{\odot}
    {\rm yr^{-1}}$. The optical colors of these galaxies are predominantly red. We find that many (1/4) of such galaxies display disturbed light profiles, sometimes in the form of fading spiral arms with no current SF detectable in H$\alpha$. In many cases we see ring-like structures. Some of these are likely post-starburst systems. We find that up to 1% of UV SFR density comes from galaxies with no H$\alpha$ lines.

-   We determine total SFR density in the local universe ($z=0.1$) of $1.828^{+0.148}_{-0.039} \times 10^{-2} h_{70} M_{\odot} {\rm yr}^{-1}
    {\rm Mpc}^{-3}$ (for Chabrier IMF), which is in good agreement with most recent estimates, yet has a smaller error (which includes random and most systematic errors).

-   Large fraction (95%) of SF takes place in blue sequence galaxies ($NUV-r<4$). SFR density increases in galaxies of increasing mass, and peaks in galaxies with masses $\approx 4\times 10^{10}
    M_{\odot}$—close to the "transition" mass of (i.e., the mass at which the mass functions of red (old) and blue galaxies intersect). Galaxies in the mass range $9.3<\log\, M_*<10.6$ account for 1/2 of the total SFR density.

We demonstrate the ability of *GALEX* UV observations of the galaxies in the local universe to constrain SFRs of a statistically large number of galaxies with a wide dynamic range of SF activity. Besides characterizing the galaxies in the local universe, this study provides a better reference point for high-redshift investigations. This and the forthcoming studies based on *GALEX* will continue to provide a better understanding of the galaxy evolution, specifically the origin of the galaxy bimodality, the inter-connection between the SF and AGN activity, and the general star formation history of galaxies spanning a range from gas-rich dwarfs to massive ellipticals.

We thank Janice C. Lee for insightful comments and helpful discussions. We thank the referee for numerous valuable suggestions. We also thank Alessandro Boselli and Michael Blanton. *GALEX* (*Galaxy Evolution Explorer*) is a NASA Small Explorer, launched in April 2003. We gratefully acknowledge NASA's support for construction, operation, and science analysis for the GALEX mission, developed in cooperation with the Centre National d'Etudes Spatiales of France and the Korean Ministry of Science and Technology. 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 U.S. 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. This research has made use of NASA's Astrophysics Data System.

*Facilities:*

[^1]: In this paper we will use the term "attenuation", rather than the more commonly used term "extinction", thus emphasizing the complex processes of absorption and scattering in a galaxy, rather than the dimming of light along a single line of sight.

[^2]: Many other studies with *GALEX* restrict analysis to 0.5 or $0.55\degr$ radius field of view. This is because artifacts are more common near the detector edge, and the PSF becomes distorted. Since in this study we match SDSS objects to *GALEX* sources, the chances of a match with an artifact are small. Also, while the PSF at the edges is distorted, the total *GALEX* flux is not affected (see also §).

[^3]: Matching with the current SDSS data release DR5 (which uses the same pipeline as DR4) would not significantly increase the overlap, since current *GALEX* pointings mostly follow the footprint of SDSS DR1 and DR2.

[^4]: The "official" SDSS calibration errors are given as 2% for $r$ band, and between 2% and 3% for each of the four *colors*. However, the individual *magnitude* errors (which figure in Equation ) are correlated, which is why we find them to be smaller than 2%. Our results stay practically the same even if errors are increased, but the $\chi^2$ distribution then deviates from the expected one.

[^5]: Available from <http://cosmo.nyu.edu/blanton/kcorrect>.

[^6]: SDSS MODEL magnitudes (defined as either the exponential or de Vaucouleurs magnitude, depending on which profile better describes a galaxy) are preferred magnitudes for SED fitting as they preserve relative fluxes (colors) better. CMODEL magnitudes (defined as a composite of exponential and de Vaucouleurs magnitude) should provide a good measure of a total galaxy light.

[^7]: *GALEX* filter response curves can be obtained from
    <http://galexgi.gsfc.nasa.gov/tools/Resolution_Response>.

[^8]: Note that the conversion between the IMFs is not generally the same for H$\alpha$-derived SFR, UV SFR, or stellar mass.