Issue |
A&A
Volume 511, February 2010
|
|
---|---|---|
Article Number | A61 | |
Number of page(s) | 15 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/200913384 | |
Published online | 10 March 2010 |
Calibration of star formation rate tracers for short- and long-lived star formation episodes
H. Otí-Floranes1,2 - J. M. Mas-Hesse1
1 - Centro de Astrobiología - LAEX (CSIC-INTA), PO Box 78,
28691 Villanueva de la Cañada, Spain
2 -
Dpto. de Física Moderna, Facultad de Ciencias, Universidad de Cantabria, 39005
Santander, Spain
Received 30 September 2009 / Accepted 30 November 2009
Abstract
Context. To derive the history of star formation in the
Universe a set of calibrated star formation rate tracers at different
wavelengths is required. The calibration has to consistently take into
account the effects of extinction, star formation regime (short or
long-lived) and the evolutionary state to avoid biases at different
redshift ranges.
Aims. We use evolutionary synthesis models optimized for intense
episodes of star formation to compute a consistent calibration of the
most usual star formation rate tracers at different energy ranges, from
X-ray to radio luminosities.
Methods. We have computed the predicted evolution of the
different estimators taking into account nearly-instantaneous and
continuous star formation regimes and the effect of interstellar
extinction (attenuation at high energies, thermal reradiation in the
far infrared). We have also considered the effect of metallicity on the
calibration of the different estimators.
Results. A consistent calibration of a complete set of star
formation rate tracers is presented, computed for the most usual
star-forming regions conditions in terms of evolutionary state, star
formation regime, interstellar extinction and initial mass function. We
discuss the validity of the different tracers in different star
formation scenarios and compare our predictions with previous
calibrations of general use.
Conclusions. In order to measure the intensity of star formation
episodes we should distinguish between nearly-instantaneous and
continuous star formation regimes. While the star formation strength ()
should be used for the former, the more common star formation rate (
yr-1)
is only valid for episodes forming stars at a constant rate during
extended periods of time. Moreover, even for extended star formation
episodes, the evolutionary state should be taken into account, since
most SFR tracers stabilize only after 100 Myr of evolution.
Key words: stars: formation - galaxies: evolution - galaxies: fundamental parameters - galaxies: starburst - galaxies: luminosity function, mass function - galaxies: stellar content
1 Introduction
The availability of a continuously increasing observational data set
on galaxies at higher and higher redshifts allows the study of the
history of star formation at cosmic scales, approaching already
(presently up to )
the ages of reionization when the first population of massive stars
started to ignite. The star formation rate (SFR) of these galaxies at
different ages is derived from the so-called star formation rate
tracers, spatially integrated parameters at different wavelengths that
allow the determine of the strength of the star formation episode by
comparison with the predictions of evolutionary synthesis models. Kennicutt (1998)
discussed the uncertainties inherent to the use of optical broad-band
colors, and presented an updated calibration of the UV continuum,
recombination lines (or number of ionizing photons emitted per unit
time), forbidden lines and total far infrared emission that has become
de facto the reference standard calibration of these estimators.
While widely used, it is often neglected that the SFR tracers calibration proposed by Kennicutt (1998)
was computed for a very specific star formation history, i.e., a star
formation episode running at a nearly constant rate during a long
enough period of time (>100 Myr), so that an equilibrium is
reached between the number of massive stars dying and igniting. Under
these conditions most of the parameters considered reach an equilibrium
state with almost no evolution with time. This continuous burst scenario is certainly valid for many galaxies, especially large spirals in which massive star formation takes place at spatially integrated
nearly constant rates during long periods of time (though at different
locations within their spiral arms). But star formation seems to
proceed in a more bursty regime in other galaxies, showing
extremely high present-day star formation rates, indicating that they
might be experiencing a nearly instantaneous (i.e., extended
over few million years) massive star formation episode. Under these
circumstances, the derivation of star formation rates by comparison
with the predictions by continuous bursts might be erroneous. Moreover,
even the concept of a star formation rate, measured as yr-1,
might be misleading in these cases. The intensity of these massive star
formation episodes should be better parameterized by the star formation strength (SFS), measured directly in units of
,
indicating the total amount of mass having been transformed into stars.
On the other hand, as noted by Rosa-González et al. (2002), different studies have indicated the significant role played by dust in the estimates of star formation rates. Extinction by dust significantly alters the integrated multiwavelength spectrum escaping from a star-forming region, weakening the far-UV continuum, but also boosting the emission in the far infrared range. It is critical therefore to take into account the effects of extinction by dust in a consistent way, since the calibration of a multiwavelength set of SFR and SFS tracers will be a strong function of the dust abundance.
In this paper we present an updated calibration of a complete set of star formation rate tracers at different wavelengths, from X-rays to radio, computed in a consistent way with state-of-the-art evolutionary synthesis models considering different star formation histories, and taking into account the effects of dust extinction at all wavelengths. In Sect. 2 we summarize the properties of the tracers considered in this work, and in Sect. 3 we describe the evolutionary synthesis models we have used. In Sect. 4 we present the calibration of the different estimators, discuss their sensitivity to the properties of the star formation scenario and compare our predictions with other evolutionary codes. Finally, Sect. 5 summarizes our main conclusions.
2 Star formation rate tracers
We have selected a large sample of SFR and SFS tracers covering a broad wavelength range, from X-rays to radio, and associated to different physical processes directly linked to the strength of the star formation episode. We want to remark that most of these tracers may be contaminated by emission not related to the present episode of star formation, but associated to an existing AGN, underlying older stellar populations, evolved but hot, low-mass stars, etc. During this work we will consider that the different star formation rate tracers have been corrected by the user before applying the calibrations. Since these corrections will not always be accurate or even possible, the SFR values so derived have to be taken with the corresponding caution. The reliability of the estimates will be larger the more SFR tracers can be used simultaneously for a given object, since the contaminating effects have different strengths at different energy ranges.
2.1 Far infrared luminosity (LFIR)
Most radiation emitted by the young, massive stars is absorbed by
the dust particles surrounding the star formation region, principally
in the UV. The heated dust reemits this energy in the far infrared
(FIR) range (which we will consider as
), especially within the
domain. The fraction of energy emitted by these stars in the UV/optical
which is absorbed by dust and reemitted as FIR radiation is close to 1
for values of E(B-V) above 0.5 (Mas-Hesse & Kunth 1991). Hence, the FIR luminosity is in general a good estimator of the bolometric luminosity of the massive stars, and thus SFR(FIR) is one of the most used and reliable SFR estimators. Kennicutt (1998) derived a calibration of
using population synthesis models of extended bursts in their
asymptotic phase (when the rate of massive star formation and death is
balanced), assuming the optically thick case, but stating that the
expression strictly applies to bursts younger than 100 Myr. This
calibration would be altered if corrections for the dust heating from
old stars in quiescent galaxies were made.
Besides the total FIR luminosity, some other calibrations make use of
the SED (spectral energy distribution) of the burst at certain FIR
wavelengths. Alonso-Herrero et al. (2006) derived
an empirical, non-linear calibration for
by combining
observations of Luminous InfraRed Galaxies (LIRGs), Ultraluminous
InfraRed Galaxies (ULIRGs), normal galaxies and HII regions within
M 51, and the
values derived from the SFR
calibration from Kennicutt (1998) (assuming case B recombination). Calzetti et al. (2007) found a similar relation based on luminosity surface density values of HII knots from 33 nearby galaxies.
can become contaminated by the presence of an AGN, especially when
dealing with spatially-integrated luminosities of distant, unresolved
galaxies. Depending on the redshift of the galaxy, the width of the
Balmer or Lyman emission lines should be checked to exclude AGN
dominated objects. In galaxies dominated by massive star-forming
episodes, the SFR values derived from
and from
,
computed with the same E(B-V) levels, should be consistent.
2.2 Ionizing power (NLyc) and emission lines
Massive stars are conspicuous sources of ionizing photons, and therefore their number
could trace the burst intensity directly. However, practically all
ionizing photons are absorbed by gas and dust within the nebular
region. Recombination of the hydrogen atoms ionized by this radiation
and free electrons produces intense emission lines, such as H
or H
,
commonly used as SFR estimators since their luminosity is proportional to
under basic assumptions (Osterbrock 1989). Kennicutt (1998) provided an expression for SFR
based on the calibrations by Kennicutt et al. (1994) and Madau et al. (1998), and used it to derive a calibration of SFR
assuming case B recombination and
K. A somewhat different expression was derived by Rosa-González et al. (2002).
They studied a sample of HII, starburst and blue compact galaxies,
in order to obtain empirical SFR calibrations based on
,
and
.
Rosa-González et al. (2002) obtained an SFR
coefficient 40% lower than the value from Kennicutt (1998) since it was based on the observed
luminosity, not corrected for Balmer absorption due to the stellar population.
H
has the disadvantage of falling outside the optical spectrum when studying medium-redshift sources (
). Studying a sample of 412 local star-forming galaxies, Moustakas et al. (2006) showed that, although weaker and more affected by extinction, H
line can also be used to estimate SFR. Applying the
calibration by Kennicutt (1998), they calculated the SFR values using both
and
,
assuming case B recombination (i.e.
). They obtained very similar results, with a 40% scatter.
The luminosity of the Ly
line can also be used as a star formation rate estimator, since it is in principle a linear function of
,
with the advantage that it is visible in the optical range for galaxies at a redshift of z > 2.0. Nevertheless, Ly
photons
are affected by resonant scattering in neutral Hydrogen, strongly
affecting the intensity of the line, as we will discuss later.
2.3 Continuum luminosity
Most radiation from young, massive stars is emitted in the
UV range (912 to around 3000 Å) and could be used as a
reliable SFR and SFS estimator since the contamination from older
stellar populations is low or even negligible below 1800 Å.
However, attenuation correction plays an important role here, since
UV photons are severly affected by extinction. Kennicutt (1998) uses the calibration from Madau et al. (1998) in order to obtain an expression for SFR as a function of ,
assuming an extended star formation process with an age larger than 100 Myr and solar metallicity. According to Kennicutt (1998), the expression can be used in the wide range of
1500-2500 Å, where the spectrum of a stellar population with a
Salpeter initial mass function (IMF) is rather flat (as a function of
frequency).
Fitting GALEX and SDSS (Sloan Digital Sky Survey) data of nearly 50 000 galaxies with SEDs constructed with the population synthesis models of Bruzual & Charlot (2003), Salim et al. (2007) derived an expression for SFR(UV) which yields SFR values lower by 30% than those obtained with the Kennicutt (1998) calibration. They claim that the reasons for this difference are the low metallicity of their sample (Z=0.016), the different star formation histories of the objects of the sample and some intrinsic differences between Bruzual & Charlot (2003) and Madau et al. (1998) models.
Some authors (Hirashita et al. 2003; Iglesias-Páramo et al. 2006; Buat et al. 1999) have opted for merging in a single expression the observed values of UV and FIR luminosities. The latter component would represent the main bulge emission of massive stars, while the former would correct for the radiation which eventually escapes from the burst before it is absorbed by dust. This way, prior assumptions about extinction would not be necessary, but both the FIR and the UV luminosities would be needed, so that the method is not easily applicable for high-redshift galaxies.
Kennicutt (1998) concluded that broad-band
luminosities in the optical are poor SFR estimators, since the optical
continuum is contributed by stars at very different evolutionary states
and affected by potentially different degrees of extinction. This is
especially true for large spiral galaxies with a complex history of
star formation. Even the UV continuum might be contaminated by
blue horizontal-branch stars in those objects dominated by an old
underlying stellar population, specially in the case of early-type
galaxies. Nevertheless, continuum luminosities in the optical-near
infrared bands could be reliable estimators for compact objects whose
continuum is dominated by the present burst of star formation. We have
calibrated SFR and SFS with the luminosities at 1500 Å,
2000 Å, 3500 Å (U), 4400 Å (B), 5500 Å (V) and 2.2 m (K).
2.4 X-ray luminosity (LX)
Mechanical energy released by the massive stars' stellar winds and
supernova explosions heats the surrounding gas, originating a diffuse
X-ray emission which peaks in the soft X-ray range. This diffuse
emission adds to the X-ray point source radiation from X-ray binaries
and supernova remnants (SNR). Whereas the emission from X-ray binaries
dominates the hard X-rays range (2-10 keV), emission from diffuse
hot gas drives the radiation in the soft range (0.2-2 keV) (Cerviño et al. 2002b).
Many authors have claimed recently that the X-ray emission from
star-forming regions should be a direct function of the burst
intensity, and several SFR estimators based on different components and
energy ranges of this radiation have been published. Analyzing the
total X-ray emission from local starburst galaxies extracted from the Ho et al. (1997) atlas, which had been observed by ASCA and/or BeppoSAX, Ranalli et al. (2003) found both soft and hard X-rays SFR linear expressions using FIR as a proxy, based on the SFR(FIR) relation by Kennicutt (1998). A similar study was performed by Tüllmann et al. (2006) using XMM-Newton and Chandra data, but in this case obvious point sources were removed. On the other hand, Grimm et al. (2003)
used the hard emission from X-ray binaries of Chandra-resolved, nearby
late-type/starburst galaxies, to derive a relation which becomes
non-linear at low SFR values. Rosa-González et al. (2009) have recently analyzed XMM-Newton observations of a sample of 14 star-forming galaxies with high SFR (0.2-160 yr-1),
as derived from radio, FIR and UV tracers, confirming that the SFR
derived from soft X-rays is comparable to that obtained from
luminosities.
Using evolutionary synthesis models, Mas-Hesse et al. (2008)
proved the importance of the evolutionary state of the burst when
trying to ascertain the SFR value via the soft X-ray luminosity, since
it can increase between half and one order of magnitude within
20 Myr in the models with constant SFR. They derived two
expressions, both for young (10 Myr) and more evolved extended star formation bursts (
30 Myr).
They also provided a star formation strength expression for nearly
instantaneous star-forming episodes, yielding the initial, total
stellar mass of gas transformed into stars. We have included the
calibration of
in this paper for completeness. We refer to Mas-Hesse et al. (2008) for a deeper discussion.
2.5 Radio luminosity (Lrad)
A very tight correlation between
and radio emission has been found over the past decades, common to a
great variety of objects, such as starburst galaxies, normal spiral
galaxies, blue compact dwarfs (BCDs), E/S0 galaxies, irregular
galaxies, HII regions, etc. (Mas-Hesse 1992; Bell 2003; Condon 1992),
which proves that radio luminosity might be a reliable star formation
rate estimator. Radio luminosity in star-forming regions is composed of
both a thermal (
)
and a non-thermal (
)
component. The former appears as a byproduct of free-free interactions,
i.e. Bremmstrahlung, and free-bound transitions between the
constituents of the ionized gas and its value can be expressed in terms
of
(Rubin 1968),
while the non-thermal emission is synchrotron in nature, emitted by
electrons after being accelerated by supernovae explosions, and it is
usually assumed to be proportional to the rate at which SN explode (Ulvestad 1982). Both contributions have different spectral index values, which are
=-0.1 for
(Rubin 1968) and
=-(0.8-0.9) for
(see Mas-Hesse 1992 and references therein). SFR calibrations based on both
and
can be found in Condon (1992), assuming
for the former (Condon & Yin 1990) and an electronic temperature
K and no dust absorption for the latter. Applying the calibration of SFR(FIR) from Kennicutt (1998)
to a sample with objects of very diverse natures (spiral galaxies,
starburst galaxies, BCDs, irregular galaxies, etc), and after
considering the tight correlation between infrared and radio emissions,
Bell (2003) obtained a non-linear expression for SFR(
) using
as a proxy.
3 Evolutionary synthesis models
In order to calibrate the different SFR tracers in a consistent way, we
have computed their expected values under different scenarios using the
evolutionary population synthesis models of Cerviño et al. (2002b)
(hereafter CMHK02 models), Leitherer et al. (1999) (Starburst99, hereafter SB99 models
) and Schaerer (2003,2002) (SC02 models
).
CMHK02 models, which are based on the models by Arnault et al. (1989), Mas-Hesse & Kunth (1991) and Cerviño & Mas-Hesse (1994),
compute the evolution of a young population of massive stars which are
formed at the same time (instantaneous bursts, IB, also referred to in
the literature as single stellar populations - SSP) or which form at a
constant rate during an extended period of time (extended bursts, EB).
The different observables are calculated for the first 30 Myr
after the onset of the burst. They are scaled to the mass of gas
transformed into stars at the start of the burst in IB models (measured
in ), and to the mass of gas transfomed into stars per unit time for EB models (
yr-1). The initial mass function (IMF) of the stellar population is defined by a power-law within the mass range of 2-120
with the slopes
.
The models include different metallicities,
.
For this work we have considered a nominal scenario with
for the IMF (i.e. a Salpeter IMF), and solar metallicity
.
The presence of binary systems has not been considered. Chemical
evolution is not treated self-consistently, since all stellar
generations in EB models have the same metallicity. This should not be
a problem, since the effect becomes important at large ages (
1 Gyr), as shown by Fioc & Rocca-Volmerange (1997).
SB99 synthesis models (Leitherer et al. 1999) are based on the original models by Leitherer & Heckman (1995), but have been continuously updated with new evolutionary tracks and stellar atmospheres. Several evolutionary tracks with different metallicity values can be used, and for consistency with CMHK02 models, we opted for the Geneva standard mass-loss tracks. SB99 allows us to compute the observables up to 1 Gyr after the onset of the burst as well as the use of a step-IMF and different mass limits. As in CMHK02 models, both EB and IB star formation regimes are considered, and metallicity-evolution is not taken into account. We have used SB99 to compute the predicted values of the different parameters for EB episodes that have formed stars during more than 200 Myr at a constant rate, since the models by CMHK02 are available only for the first 30 Myr. We have checked that where they overlap in time, the predictions by both sets of models are perfectly consistent. For the IB star formation regime we have preferred to use CMHK02 models, since they are based on Monte Carlo simulations of the IMF and consider therefore in a more realistic way the stochastical nature of the massive stellar populations. Moreover, since they provide directly the expected X-ray luminosity, as discussed below, the other SFR tracers should be computed with the same models for consistency.
Finally, SC02 models have been used for completeness to compute
and
for population III (Pop. III) stellar populations, which could
correspond to the scenario valid for star-forming galaxies at very high
redshift, where metallicity is extremely low (
).
While most SFR tracers calibrations are based on the predictions by
evolutionary models which are computed for extended star formation
processes where gas is transformed into stars at a nearly constant rate
during long periods of time (hundreds to thousands Myr), we want
to stress that star formation seems to proceed in the form of short,
nearly instantaneous bursts in many galaxies. It is generally difficult
to distinguish a short, but young starburst from an evolved, extended
star formation process, since most observable parameters are related to
just the most massive and therefore young stars in the region.
Nevertheless, a long-lasting star formation process produces an
accumulation of medium-low mass stars with long lifetimes, which can
dominate the optical continuum. As an effect, these kind of episodes
show systematically lower equivalent width values of
and Wolf-Rayet star bumps at around 4686 Å (hereafter EW(WRbump)) than young instantaneous bursts (Cerviño & Mas-Hesse 1994).
![]() |
Figure 1: Evolution of EW(WRbump) predicted by CMHK02 for IB (solid lines) and EB (dashed lines) models both for Z=0.020 ( top) and Z=0.008 ( bottom), together with the histogram of the EW(WRbump) data of the sample from Brinchmann et al. (2008). |
Open with DEXTER |
Figure 1 shows that continuous star formation models in the asymptotic phase cannot reproduce the high wing of the distribution of EW(WRbump) values measured by Brinchmann et al. (2008) on their sample of 570 star-forming galaxies. Instantaneous bursts, on the other hand, of course predict the high EW(WRbump) range (see also Mas-Hesse & Kunth 1999). The sample by Brinchmann et al. (2008) is clearly biased towards star-forming galaxies with large Wolf-Rayet bump equivalent widths, i.e., experiencing short star formation episodes, but other authors have found evidences of short-lived bursts. Pellerin & Robert (2007), for example, find that EB models can not reproduce the FUSE far-UV data of the 24 starburst galaxies of their sample, whereas there is a much better agreement with the predictions by IB models.
We conclude therefore that star formation proceeds in the form of
nearly instantaneous, short bursts in a significant fraction of
objects. In the case of very large, non-resolved galaxies, the
superposition of individually instantaneous bursts at different
evolutionary states can mimic the properties of an extended star
formation process at a nearly constant rate, while large objects at
high redshift could sustain intense star formation processes during
truly extended periods of time. When deriving the intensity of the
different star formation episodes in the history of a galaxy, its star
formation regime has to be carefully taken into account. In the case of
a nearly instantaneous burst, the ionizing power decreases very fast
with time, and becomes negligible after the first 8-10 Myr (Mas-Hesse & Kunth 1991). Blindly applying an SFR estimator related to
,
calibrated for an extended star formation regime, can lead to
discrepancies of several orders of magnitude depending on the
evolutionary state. Moreover, even the concept of a star formation rate
as such becomes completely misleading in these cases, since formation
of massive stars might have ended already a few million years ago. The
concept of a star formation strength, SFS, measured as the total mass
of gas transformed into stars during the burst, is much better adapted
to this scenario. Unfortunately, deriving the SFS by comparing the
observations with the predictions of synthesis models requires a
preliminary estimate of the evolutionary state of the burst, and this
can be obtained only when parameters like the equivalent width of
can be measured. We have therefore calibrated our set of SFR and SFS
tracers for both star formation regimes, extended and instantaneous
respectively, and for different evolutionary states of the burst.
Formation of massive stars is an inherently stochastic process, so that
when a relatively low number of stars are formed, we can not guarantee
that the whole IMF will be filled. Cerviño et al. (2002a) showed that stochastic effects become small enough only for starbursts that have transfomed more than 105-6
of gas into stars. The star formation estimators should therefore not
be applied to smaller star-forming regions, where the stochastic
effects could significantly affect the integrated emission.
In the sections below we describe the way the different estimators have been computed. More details can be found in Mas-Hesse & Kunth (1991), Cerviño & Mas-Hesse (1994) and Cerviño et al. (2002b).
3.1 Far infrared luminosity (LFIR)
The FIR emission predicted by CMHK02 models is calculated with the
assumption that interstellar dust remains in thermal equilibrium. All
energy absorbed, mainly stellar UV continuum, is therefore
reradiated in the FIR range (
). Calzetti et al. (2000) noticed that this parametrization is higher by up to 75% than the value obtained using the expression by Helou et al. (1988), which models the IR output within
via the IRAS fluxes at
and
,
assuming a single dust component with T=20-80 K. The Galactic extinction law by Cardelli et al. (1989) with RV=3.1, is applied to the synthetic spectral energy distributions calculated by the models to derive
.
The extinction is parameterized by the value of the color excess E(B-V). Since
saturates for E(B-V) >0.5, the value computed for E(B-V) = 1.0 should be taken as the maximum FIR emission expected in the starburst (Mas-Hesse et al. 2008). We slightly modified the extinction law, assuming that photons with a wavelength of 912 Å
Å suffered the same extinction
.
In this way, we both a) take into account photons with 912 Å
Å, which would otherwise be ignored, and b) use a more realistic
value, bearing in mind that the attenuation law does not seem to increase in this range (as Eq. (5) in Cardelli et al. 1989 would imply), but to flatten, as concluded by Mezger et al. (1982). The effect of including the photons with 912 Å
Å in the FIR calculation is that
is increased by 5%. Similarly, the computed
changes by up to 7% due to the modifications performed in the extinction law in the range of 912 Å
Å, but this is only detectable for low E(B-V) values. A fraction 1-f=0.3
of ionizing photons, irrespective of their energy, is assumed to be
absorbed by dust before they can ionize any atom, as derived by Mezger (1978) and Degioia-Eastwood (1992) and recommended by Belfort et al. (1987).
The hypothetical presence of totally obscured stars was not considered.
Since no assumption is made about the temperature of the dust, the
models do not yield the infrared spectrum, but only the total energy
absorbed by dust, which would be reemitted in the FIR range.
The same prescriptions as those considered in CMHK02 were taken
with the SB99 models to model the FIR emission of the burst. The
extinction law by Cardelli et al. (1989),
with the minor changes explained above, was applied to the predicted
SEDs, and 30% of the ionizing radiation was assumed to be converted
into FIR emission. The estimation of
using the SB99 models is
7%
higher than predicted by CMHK02 for the same conditions, apparently due
to the differences in the stellar atmospheres used, especially in the
Lyman continuum.
3.2 Ionizing power (NLyc) and emission lines
As discussed above, we considered that a fraction 1-f=0.3 of the ionizing photons emitted by the stars is directly absorbed by dust, so that

The intensity of the emission lines has been derived from the predicted
value, assuming Case B conditions (
K,
cm-3) and the following relations (Storey & Hummer 1995):



We have computed the intensities of the emission lines for various E(B-V) values, but we have not considered the effects of radiation transfer on Ly






3.3 Continuum luminosity
We have computed the luminosity evolution of the burst's continuum
emission at several wavelengths within the UV, optical and IR spectrum:
1500 Å, 2000 Å, 3500 Å (U), 4400 Å (B), 5500 Å (V) and 2.2 m (K).
We included both the stellar and nebular continuum components and have
considered the effect of interstellar extinction as parameterized by E(B-V). The most massive stars dominate the UV emission of star-forming regions and therefore drive the evolution of L1500 and L2000.
On the other hand, the contribution of less massive stars can become
important at longer wavelengths. In EB models, medium/low-mass stars
accumulate due to their longer lifetimes and become the dominant
contributors to L4400, L5500 and
L22 200 after
100 Myr of evolution.
3.4 X-ray luminosity (L X)
The X-ray luminosity associated to a massive star formation episode is directly related to the amount of mechanical energy released into the interstellar medium (ISM) by stellar winds and supernova explosions. In order to calculate the mechanical energy relased by supernova explosions, both the CMHK02 and SB99 models assume that each SN injects 1051 erg into the ISM, but CMHK02 models subtract the energy emitted as X-ray radiation by SNR afterwards from the total energy of the supernova. See Cerviño et al. (2002b) for a more detailed description. Massive early-type stars, and especially Wolf-Rayet stars, lose a significant fraction of their mass in the form of strong winds. These winds end up interacting with the interstellar gas, which becomes extremely heated. As a result, soft X-ray emission is produced. Both models use the expressions of the gas terminal velocity from Leitherer et al. (1992) together with their own prediction of mass loss to calculate the mechanical energy released by the winds.
We have computed the evolution of the mechanical energy released by the burst per unit time,
.
This energy is converted into X-ray luminosity assuming an efficiency factor
.
As discussed by Mas-Hesse et al. (2008) the soft X-ray luminosity of star-forming galaxies can be reproduced by our synthesis models assuming
in the range of 1-10%. We assumed
=0.05
for the calibrations presented in this work. The total soft X-ray
luminosity is then obtained by adding the contribution by the SNR
present in the star-forming region, still active as X-ray emitters.
3.5 Radio luminosity (Lrad)
Both radio components, thermal and non-thermal, were calculated at 1.4 GHz. For the thermal emission
we followed the prescriptions by Lequeux et al. (1981). This thermal component of the radio emission is directly proportional to the ionizing power of the starburst, with a slope
=-0.1:

We assumed an electronic temperature


The models by Mas-Hesse & Kunth (1991) underestimated the non-thermal emission associated to a starburst region. Mas-Hesse (1992) showed that for the ratio
/
to
remain rather constant during the evolution of the burst, so that the
observational well-known correlation between radio and
FIR emission could be explained, the total non-thermal emission
observed should be
10 times
larger than computed. It was argued that interactions between
individual supernova remnants would yield a larger radio emission than
observed in isolated SNR in our Galaxy, as the ones used for the
standard calibration. The same phenomenon was found by Condon & Yin (1990)
when they applied their Eq. (7) to our Galaxy, stating that the
disagreement arises because the effect of accelerated electrons after
yr is ignored. They obtained their Eq. (8), which seems to overcome this problem. We used the expression

where






We have finally calibrated both SFS and SFR based on the total radio luminosity at 1.4 GHz
=
+
.
4 Results and discussion
Table 1: SFS estimators which do not depend on E(B-V).
Table 2: SFS estimators which do depend on E(B-V) for an IB model at age = 4 Myr.
Table 3: SFS estimators which do depend on E(B-V) for an IB model at age = 5 Myr.
Table 4: SFS estimators which do depend on E(B-V) for an IB model at age = 6 Myr.
Table 5: SFR estimators which do not depend on E(B-V).
Table 6: SFR estimators which do depend on E(B-V) for an EB model at age = 10 Myr.
Table 7: SFR estimators which do depend on E(B-V) for an EB model at age = 30 Myr.
Table 8: SFR estimators which do depend on E(B-V) for an EB model at age = 250 Myr.
Table 9:
SFS estimators based on
and L1500 for Pop. III.
Table 10:
SFR estimators based on
and L1500 for Pop. III.
In Tables 1 to 8 we present the calibration of the different estimators computed with our models. We list the values of the factor
by which the different observables must be multiplied to obtain either SFS or SFR, i.e. SFS
,
SFR
.
As discussed above, Mas-Hesse et al. (2008)
showed that it is very important to take into account the evolutionary
state of the burst when estimating its strength, especially in the case
of IB regimes, but also in some cases for EB episodes. For IB models we
have calibrated the estimators for ages of 4, 5 and 6 Myr, which
is the typical range of ages measured for star-forming galaxies (Pellerin & Robert 2007; Mas-Hesse & Kunth 1999). The ionizing power decreases very rapidly, so that the emission lines are barely detectable after 7 Myr: Mas-Hesse & Kunth (1999) predict
values below 2 Å after this age for solar metallicity starbursts.
On the other hand, in an EB regime stars are replaced by new massive
stars as they die, eventually reaching a steady state after a few tens
of Myr. For this reason, we studied a wider range of ages: a
non-evolved continuous episode with 10 Myr, a young burst after
30 Myr of evolution, and a star formation process already in the
steady state, at 250 Myr. Not all parameters remain constant when
the steady state is reached, since medium/low-mass stars accumulate and
contribute more and more to the optical - IR stellar continuum. As a
result the calibration of some tracers is a function of the
evolutionary state as well, though the effect is much weaker than for
IB regimes. For
and hence for
and
,
the asymptotic state is reached very soon, so that the same
value is displayed for the three points at 10 Myr, 30 Myr and 250 Myr.
For the parameters which are affected by interstellar extinction we computed the calibration of the different estimators for different E(B-V) values within the range E(B-V) = 0.0 - 1.0 under the assumption of the Galactic extinction law by Cardelli et al. (1989).
As mentioned above, our calibrations have been computed for a Salpeter initial mass function with limits at 2 and 120 .
To ease the comparison with other calibrations we list the correction factors in Tables 11 and 12 by which values from Tables 1-4, and 5-8, respectively, must be multiplied to convert our calibrations to other usual IMF mass ranges of
M=0.1-100
and M=1-100
.
Similarly, in Tables 13 and 14 we list the conversion factors which would be required for two limiting IMF slopes, with
and 3. These correction factors are the values
,
where
and
are the new mass limits,
is the new IMF slope and L the magnitude. For instance, the SFR(
) calibration for a burst that has been creating stars with a constant rate for 30 Myr, and assuming an IMF with mass range
M=0.1-100
and
,
is SFR(
)
(luminosities are measured in erg s-1 Hz-1 and SFR in
yr-1).
We calculated the values of this correction factor using the SB99
models, since these models allow the IMF to be defined by the user.
This correction between different IMFs depends on the age of the burst,
on the parameter studied and on the star formation history assumed in
the model and has two different origins: a) when considering
different mass limits, both SFR and SFS measure the mass converted into
stars in a different mass range, and b) if
or
changes, so does the amount and eventual type of massive stars that
contribute to the emission, which can be significantly affected.
Emission by low-mass stars does not contribute dramatically to any of
the computed parameters at the ages considered, and therefore any
change in
only exerts a variation of type a) in the calibration, as long as it
does not increment unphysically. Also, since stars with masses above
100
have already died in IB models at the ages considered in the calibrations, only correction of type a) is needed (if
is not modified), which can be performed analytically by the calculation of the ratio of the integrated IMFs (Schaerer 2003; Wilkins et al. 2008),
and which does not depend on the evolutionary state. On the contrary,
in EB models there are always stars more massive than 100
,
hence a proper correction for both a) and b) using SB99 as explained is needed. As can be checked in Tables 11 and 12, and also by a comparison of these values with those from Schaerer (2003),
values of correction factors obtained analytically agree with those
obtained for EB models using SB99 typically within 5%, with the
difference being higher only for
(
10-15%) due to the fact that it is solely dominated by the most massive stars of the burst. Since SB99 models do not yield
,
it was not possible to follow the same procedure in this case as for
the other magnitudes, and therefore only the analytical value for
is given in Table 12.
The effect of the age of the burst in the correction of SFR estimators when different mass limits are assumed is very low (8% at most) for the ages considered, hence we decided to leave only the asymptotical value at 250 Myr in Table 12. Similarly, the effect of E(B-V) value in the SFR(FIR) correction factor is typically below 5%, therefore we opted for showing only the correction obtained for E(B-V) =0.3, which lies somewhat in the middle of the range studied, where
has not attained the saturation described in Sect. 3.1. Since
and
are proportional to
,
only the correction for the latter was included. We have applied these
corrections when we compared our calibrations with the results of Kennicutt (1998) and Salim et al. (2007), which were computed for a Salpeter IMF, but with mass limits of 0.1-100
.
Condon (1992) assumed mass limits of 5-100
,
but an IMF slope
,
therefore
L(2-120,2.35)/L(5-100,2.50) was calculated explicitly in order to compare his radio calibrations with our predictions.
Table 11:
IMF correction factors for all SFS estimators, to be applied on values from Tables 1-4 when IMF mass limits 0.1-100
or 1-100
are assumed. See Sect. 4 for a more thorough discussion.
In the sections below we discuss the properties and applicability of each of the SFR and SFS estimators, and we compare our results with previous calibrations.
4.1 Far infrared luminosity (LFIR)
Table 12:
IMF correction factors for SFR estimators, to be applied on values from Tables 5-8 when IMF mass limits of 0.1-100
or 1-100
are considered. See Sect. 4 for a more thorough discussion.
![]() |
Figure 2:
Evolution of
|
Open with DEXTER |
Evolution of
is shown in Fig. 2 for both the CMHK02 IB and SB99 EB models and different metallicities
,
computed for E(B-V)= 1.0. Since
in star-forming regions is directly linked to the UV-optical emission
produced by the most massive stars of the bursts, it shows a peak in IB
models at
3 Myr, followed by a steep decline afterwards. On the other hand, when a constant SFR is assumed,
starts to stabilize after the first 20 Myr of evolution, when an
equilibrium between formation and destruction of massive stars is
reached. Nevertheless, since medium/low-mass stars also contribute to
,
this parameter increases slowly but steadily as long as the star formation episode is active. As discussed by Mas-Hesse et al. (2008),
saturates for E(B-V) >0.5 since the attenuation is already so high that most of the continuum emission is absorbed. We have calibrated
as SFR and SFS estimator in Tables 1 to 8 for various E(B-V) values. If E(B-V)= 1.0 is assumed,
becomes indeed a good estimate of the bolometric luminosity of the star formation process.
Table 13:
IMF correction factors for SFS estimators, to be applied on values from Tables 1-4 when IMF slope ,
3.
Table 14:
IMF correction factors for SFR estimators, to be applied on values from Tables 5-8 when IMF slope ,
3.
The evolution of
with time is sensitive to the metallicity, since low-metallicity stars evolve more slowly. As a result,
is 15-30% larger for Z=0.008
than for solar metallicity models at ages 4-6 Myr. In EB models,
the effect is to accumulate more massive stars at a given time, so that
becomes
8-11% larger for Z=0.008 once the burst is older than 30 Myr. On the other hand, during the initial phases of evolution,
is essentially independent of the metallicity.
Together with our predictions we show in Fig. 2 the
value derived from the expression by Kennicutt (1998), which was IMF-corrected as explained above to be consistent with the IMF we have used. The adapted Kennicutt (1998) value becomes
=43.90, which agrees quite well with our predictions for solar metallicity at ages above
100 Myr. On the other hand, the discrepancy becomes larger for less evolved bursts, differing by
21% at 30 Myr, and by much more for the younger the episode. Although Kennicutt (1998) affirms that his calibration applies only to bursts younger than 100 Myr, it can be seen in Fig. 2 that it could be used accurately for solar metallicity extended bursts at 50-250 Myr.
4.2 Ionizing power (NLyc) and emission lines
![]() |
Figure 3:
Evolution of
|
Open with DEXTER |
![]() |
Figure 4:
Evolution of
|
Open with DEXTER |
We show in Fig. 3 the evolution of
for IB and EB models at different metallicities, as predicted by CMHK02
models. The emission of ionizing photons is driven exclusively by the
most massive stars in the burst. Since these stars have short
lifetimes, the ionizing power of instantaneous bursts decreases
drastically with time after 3 Myr, when the most massive stars
begin to explode as supernovae. For the same reason, their number
reaches an equilibrium value quite rapidly in models with a constant
star formation rate, so that
reaches an asymptotic value in just
8 Myr in EB models. This steady state of
lasts hundreds of Myr in EB models, but we show in the figure just
the first 20 Myr of evolution of the burst to remark its rapid
stabilization. As we discussed, we considered that a fraction 1-f=0.3
of the ionizing photons are absorbed by dust, and here we show their
number which can eventually ionize the atoms in the surrounding gas.
When considering
,
we find that the slower evolution of high-mass stars translates into a
delay in the emission of ionizing photons in IB models, as can be
observed in Fig. 3. The effect is to increase
significantly, by 0.4 dex at 4 Myr and 0.7 dex at 6 Myr for Z=0.001. In EB models the accumulation of massive stars originates
values 0.2 dex higher for Z=0.001.
We include in Fig. 3 the
value by Kennicutt (1998), adapted to our IMF limits. This
value is larger than our prediction for solar metallicity models by 50%, which is apparently due to the fact that Kennicutt (1998) did not apply any f correction. If we add 1-f=0.3 correction to the Kennicutt (1998)
value, both predictions are within 6% in the asymptotic phase, as can be observed in Fig. 3. We insist that, as discussed by Mas-Hesse & Kunth (1991),
it should be more realistic to assume that a fraction around 30% of
ionizing photons are absorbed by dust, even in relatively dust-free
environments, and that they therefore do not contribute to the
ionization process.
In Tables 1 to 8 we list the predicted
values for IB and EB models at different evolutionary states and the predicted
and Ly
calibrations for different values of interstellar extinction, computed as discussed in Sect. 3.2. In Tables 9 and 10 we list the corresponding
predictions for starbursts dominated by Pop. III stars. As discussed above, 1-f=0.0 has been assumed in this case.
We wish to remark that there may be significant differences between the
extinction affecting the nebular lines and the stellar continuum. Calzetti et al. (2000) estimated that on average the color excess (E(B-V))
derived from the stellar continuum was a factor of 0.4 lower
than the value derived from the emission lines. The analysis of
NGC 4214, a well-resolved starburst galaxy, by Maiz-Apellaniz et al. (1998)
showed that the dust appeared concentrated at the boundaries of the
ionized region in its central star-forming region and affected mainly
the nebular emission lines, while the stellar continuum itself was
located in a region relatively free of dust and gas. But on the other
hand, no spatial decoupling between massive stars, ionized gas and dust
was found on other star-forming regions of the same galaxy, which
indicates that the specific geometry of each star-forming knot drives
the decoupling of continuum and emission line extinctions. Erb et al. (2006) concluded from the analysis of a large sample (114 objects) of star-forming galaxies at
that the average SFR values derived from both
and
,
applying the same E(B-V)
for the extinction correction, were completely consistent (indeed,
almost identical). Nevertheless, the scatter was significant
(0.3 dex), which again indicates that the geometry of individual
objects plays an important role. This effect has to be taken into
account if the SFR is derived simultaneously from tracers based on both
the emission lines intensity and continuum level. Finally, we want to
stress that the strength of the extinction assumed in our computations
results from the combination of the listed E(B-V) values and the Cardelli et al. (1989) law. To get the same extinction with different laws, correspondingly different E(B-V) values have to be used (for example, the same extinction is obtained at
with E(B-V) = 0.3 and the Cardelli et al. (1989) law as with E(B-V) = 0.23 and the Calzetti et al. (2000) parameterization).
![]() |
Figure 5:
Evolution of
|
Open with DEXTER |
![]() |
Figure 6:
Evolution ( from top to bottom) of L1500, L2000 and L3500 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER |
4.3 Radio emission
The evolution of
and
is presented in Figs. 4 and 5, respectively. Since the thermal radio component is assumed to be proportional to
,
its evolution directly follows that of
,
as discussed above. But we still show the plots in Fig. 4
to allow for a direct comparison between the intensities of the thermal
and non-thermal components at different evolutionary states of the
burst. We have also plotted on Figs. 4 and 5 the values predicted by Condon (1992) for EB models after they were adapted to our IMF.
When constant star formation is assumed,
is observed to increase two orders of magnitude within
35 Myr
after the SN start to explode, reaching then a steady state with almost
no sensitivity to the metal content of the star-forming region. The
predictions for IB models, on the other hand, are more sensitive to
metallicity, especially at ages between 3 and 12 Myr, when the
most massive stars explode (their lifetime and evolution depend on
metallicity).
In Tables 1 to 8 we list the calibration of
,
which includes both the thermal and non-thermal components at any time.
It can be seen in the figures that the total radio emission at 1.4 GHz
is dominated by the non-thermal component as soon as supernova
explosions take place in the star-forming region. Only during the first
2-3 Myr in IB bursts can the radio luminosity be dominated by
thermal emission.
Although the calibration derived by Condon (1992) is consistent with our predictions, Fig. 4 shows that the
value is 12% higher than predicted by our models. However, when using the
value from Kennicutt (1983) as Condon (1992) did, with the prescriptions we have explained in Sect. 3.2, it can be checked that this difference falls to 5%. But we wish to note that Condon (1992) ignores the dust absorption of ionizing photons, but assumes a
30%
lower ionizing power of the burst than predicted by the SB99 models.
Both effects seem to cancel each other, which renders both
predictions consistent.
![]() |
Figure 7:
Evolution ( from top to bottom) of L4400 and L5500 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER |
![]() |
Figure 8:
Evolution of
L22 200 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER |
4.4 Ultraviolet and U-band continuum
UV emission is dominated by the most massive stars in the burst (via stellar and nebular continuum radiation), and therefore constitutes a direct tracer of the presence of young, massive stars, since it is free of contamination by the continuum emission of low-mass stars. Nevertheless, interstellar extinction is very strong in this energy range, and therefore its correction plays a key role in the use of the UV continuum luminosity as a tracer of star formation. In Fig. 6 we show the evolution of the luminosity at 1500, 2000 and 3500 Å for IB models as predicted by CMHK02, and for EB models by SB99, for Z=0.020 and Z=0.008. No extinction has been applied in these figures.
In IB models, the UV continuum follows the well-known trend, with a peak of the emission at 3 Myr,
and a severe decrease when the most massive stars begin to explode. The
slower evolution at low metallicity is also evident in the figures. In
the case of EB models, the UV continuum asymptotically stabilizes
after the first
50 Myr of evolution. The stabilization is slower for L3500, since it is contributed by stars of a lower mass than L1500.
We show also in Fig. 6 the L1500 values by Kennicutt (1998) and Salim et al. (2007), after they are adapted to our IMF mass limits. Both calibrations are close to our predictions, with a perfect match with the one by Kennicutt (1998) at ages around 100 Myr for solar metallicity models. The predictions by Salim et al. (2007) are somewhat overestimated with respect to both the predictions by Kennicutt (1998) and by us. But these authors used the Bruzual & Charlot (2003) synthesis models in order to fit the photometric data points of their sample, assuming an exponentially-declining continuous star formation history, over which random-in-time starbursts were superimposed. The star formation histories are so different between these models that a direct comparison is not straightforward.
Kennicutt (1998) argued that his calibration
is applicable to bursts which have been forming stars for 100 Myr
or longer. Our results show that these calibrations can indeed be used
accurately for ages above Myr, since the stronger source of error will be the correction of the interstellar extinction.
In Tables 9 and 10 we list the corresponding L1500 predictions for starbursts dominated by Pop. III stars. As discussed above, no extinction has been assumed in this case, since no dust particles are expected to be present.
![]() |
Figure 9: Evolution of L22 200 predicted by SB99 for EB models as in Fig. 8. |
Open with DEXTER |
4.5 B-, V- and K-band continuum
![]() |
Figure 10:
Evolution of
|
Open with DEXTER |
The evolution of L22 200 is strongly affected by the formation and evolution of red supergiant stars (RSG), which are very dependent on metallicity. We plot in Fig. 8 the evolution of L22 200 as predicted by CMHK02 for IB models and by SB99 for EB models for three values of the metallicity. The formation of RSG is favoured in high metallicity environments, dominating completely L22 200 in Z=0.020 IB models in the period of 5-15 Myr. On the other hand, their contribution for Z=0.001 becomes negligible (see Cerviño & Mas-Hesse 1994 for a more detailed discussion).
In extended episodes the population of RSG stars tends to stabilize at high metallicities. Nevertheless, the accumulation of low-mass stars, whose contribution to the K-band is important, drives a continuous increase of L22 200, which is specially evident for low metallicities where the contribution of RSG is negligible. In Fig. 9 we zoom on the first 20 Myr of evolution, to show the effect of metallicity on L22 200 early evolution.
4.6 X-ray luminosity (L X)
At the early stages of the burst, the injection of mechanical energy
into the ISM is dominated by stellar winds, which tend to be stronger
as the most massive stars evolve. Around 3 Myr after the onset of
the burst, massive stars end their lifetimes, and this component begins
to decrease monotonically. Then the first SN start to explode, becoming
the dominant sources of mechanical energy release in IB models after
3-4 Myr. In EB models both components tend to be similar once the
stellar population has attained a steady state, around 30-40 Myr
after the onset of the burst. This can be seen in Fig. 10, where we have plotted the evolution of
as predicted by SB99 and CMHK02 for EB and IB models respectively, and
for different metallicities. The strongest release of mechanical energy
in IB models occurs at
4 Myr,
when the first SN appear and add their contribution to the winds
dominated by the Wolf-Rayet stars. The release of mechanical energy
stabilizes in IB models after
6 Myr of evolution, since the supernova rate decreases very slowy for the next 20 Myr for a Salpeter IMF (see Cerviño & Mas-Hesse 1994). In EB models the evolution of
stabilizes completely 30-40 Myr after the onset of the burst,
depending on metallicity, when both the population of the most massive
stars, responsible for the winds and SN, and medium-mass stars (M > 6-8
), originator of SN as well, have reached an equilibrium.
is very sensitive to metallicity as long as it is dominated by stellar
winds, since the mass-loss rate varies largely with the metal abundance
of stellar atmospheres. This is especially evident in the first
3 Myr of evolution of IB models, as shown in Fig. 10, but the effect can also be seen on the lowest metallicity EB models.
As discussed above we have assumed
= 0.05 to calibrate
,
i.e., we consider that 5% of the mechanical power that is released is
reconverted into soft X-ray emission. We want to stress that while
seems to be in the range
= 0.01-0.1 for star-forming galaxies, it could attain higher or lower
values in individual objects, depending on the details of the
interaction between the accelerated winds and the surrounding
interstellar medium. In Table 1 we list only the
)
calibration for an IB model at 5 Myr. As discussed by Mas-Hesse et al. (2008) there are large uncertainties in
at these ages, therefore the value listed should be considered as the
best average for the 4-6 Myr period. The original SFR calibration
by Mas-Hesse et al. (2008)
did not provide a value for a burst age of 250 Myr, but as we have
explained, a steady state is expected at 30-40 Myr in EB models,
so we assume the same estimator value for 250 Myr as for
30 Myr, as listed in Table 5.
5 Conclusions
We have computed a self-consistent calibration of the most usual star formation rate tracers using evolutionary synthesis models optimized for both short-lived and extended star formation processes. We have taken into account the effects of the evolutionary state, the star formation regime and the effect of interstellar extinction on the different estimators. We have also included the prediction for population III starbursts. The main results can be summarized as follows:
- The intensity of star formation episodes should be measured as a star formation rate (SFR), in units of
yr-1, when star formation proceeds at nearly constant rates for long periods of time (tens to hundreds of Myr), and as a star formation strength (SFS), in units of
, when the episode is nearly instantaneous or short-lived (few Myr). Deriving an SFR value for an instantaneous burst would not have any physical meaning. Moreover, the evolution is so rapid in these conditions that the results would be misleading and even erroneous.
- We have provided calibrations for SFS estimators at the ages of 4-6 Myr. The ionizing power of these bursts decays so rapidly that they would show just very weak emission lines at older ages, and would not be easily recognized as massive starbursts.
- The evolutionary state should also be taken into account even for extended episodes of star formation, whenever possible. In general most of the SFR tracers have stabilized after 100 Myr of evolution at a constant rate.
- Some of the most used SFR and SFS tracers are strongly
affected by interstellar extinction. We have provided their calibration
as a function of E(B-V) assuming the Galactic extinction law by Cardelli et al. (1989), to allow for a more precise result if the extinction can be somehow estimated. SFR(FIR) has also been computed as a function of the extinction. For E(B-V) values above 0.5
can be considered as a good estimate of the bolometric luminosity of the burst.
- SFR and SFS tracers allow to derive the amount of gas
transformed into stars, but they are calibrated for a given initial
mass function. The choice of different usual IMF mass limits within the
range 0.1-120
, and IMF slopes
can lead to variations in the derived rates of up to a factor 4. A comparison of results by different authors is therefore meaningful only after a normalization to a common IMF.
- Our calibration of star formation estimators at different energy ranges (from X-rays to radio) can also be used inversely, (taking care of the different time constants and properties of the different tracers): for example, once the SFR or SFS is estimated from the UV or optical continuum luminosity, the expected emission at other ranges or emission lines can be directly computed, corrected for the effects of interstellar extinction in a consistent way. Nevertheless, the inverse use of the SFR and SFS calibrators to derive the expected luminosity at other energy ranges has to be taken with caution, since they would predict only the luminosity associated to the star formation episode, which could be just a fraction of the total luminosity of the galaxy at certain wavelength ranges dominated by an old underlying stellar population.
J.M.M.H. and H.O.F. are partially funded by Spanish MICINN grants CSD2006-00070 (CONSOLIDER GTC), AYA 2007-67965 and AYA2008-03467/ESP. H.O.F. is funded by Spanish FPI grant BES-2006-13489. We thank members from the Estallidos team for giving us the idea of creating the website and for comments on its first versions, and Raúl Gutiérrez Sánchez and Carlos Rodrigo Blanco for their help in its implementation. We also thank Miguel Cerviño for his helpful comments about the CMHK02 models. We want to acknowledge the use of the Starburst99 models. We are also very grateful to D. Schaerer for providing his calculations for population III models.
References
- Alonso-Herrero, A., Alonso-Herrero, A., Rieke, G. H., et al. 2006, ApJ, 650, 835 [NASA ADS] [CrossRef] [Google Scholar]
- Arnault, P., Kunth, D., & Schild, H. 1989, A&A, 224, 73 [NASA ADS] [Google Scholar]
- Atek, H., Kunth, D., Schaerer, D., et al. 2009, A&A, 506, L1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Belfort, P., Mochkovitch, R., & Dennefeld, M. 1987, A&A, 176, 1 [NASA ADS] [Google Scholar]
- Bell, E. F. 2003, ApJ, 586, 794 [Google Scholar]
- Berkhuijsen, E. M. 1984, A&A, 140, 431 [NASA ADS] [Google Scholar]
- Brinchmann, J., Kunth, D., & Durret, F. 2008, A&A, 485, 657 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bruzual, G., & Charlot, S. 2003, MNRAS, 344, 1000 [NASA ADS] [CrossRef] [Google Scholar]
- Buat, V., Donas, J., Milliard, B., & Xu, C. 1999, A&A, 352, 371 [NASA ADS] [Google Scholar]
- Calzetti, D., Armus, L., Bohlin, R. C., et al. 2000, ApJ, 533, 682 [NASA ADS] [CrossRef] [Google Scholar]
- Calzetti, D., Kennicutt, R. C., Engelbracht, C. W., et al. 2007, ApJ, 666, 870 [NASA ADS] [CrossRef] [Google Scholar]
- Cardelli, J. A., Clayton, G. C., & Mathis, J. S. 1989, ApJ, 345, 245 [NASA ADS] [CrossRef] [Google Scholar]
- Cerviño, M., & Mas-Hesse, J. M. 1994, A&A, 284, 749 [NASA ADS] [Google Scholar]
- Cerviño, M., Valls-Gabaud, D., Luridiana, V., & Mas-Hesse, J. M. 2002a, A&A, 381, 51 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Cerviño, M., Mas-Hesse, J. M., & Kunth, D. 2002b, A&A, 392, 19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Condon, J. J. 1992, ARA&A, 30, 575 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
- Condon, J. J., & Yin Q. F., 1990, ApJ, 357, 97 [NASA ADS] [CrossRef] [Google Scholar]
- Degioia-Eastwood, K. 1992, ApJ, 397, 542 [NASA ADS] [CrossRef] [Google Scholar]
- Erb, D. K., Steidel, C. C., Shapley, A. E., et al. 2006, ApJ, 647, 128 [NASA ADS] [CrossRef] [Google Scholar]
- Fioc, M., & Rocca-Volmerange, B. 1997, A&A, 326, 950 [NASA ADS] [Google Scholar]
- Grimm, H.-J., Gilfanov, M., & Sunyaev, R. 2003, MNRAS, 339, 793 [NASA ADS] [CrossRef] [Google Scholar]
- Hayes, M., et al. 2009, in preparation [Google Scholar]
- Helou, G., Khan, I. R., Malek, L., & Boehmer, L. 1988, ApJS, 68, 151 [NASA ADS] [CrossRef] [Google Scholar]
- Hirashita, H., Buat, V., & Inoue, A. K. 2003, A&A, 410, 83 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ho, L. C., Filippenko, A. V., & Sargent, W. L. W. 1997, ApJS, 112, 315 [NASA ADS] [CrossRef] [Google Scholar]
- Iglesias-Páramo, J., Buat, V., Donas, J., Boselli, A., & Milliard, B. 2004, A&A, 419, 109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kennicutt, R. C., Jr. 1983, ApJ, 272, 54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kennicutt R. C., Jr. 1998, ARA&A, 36, 189 [Google Scholar]
- Kennicutt R. C., Jr., Tamblyn, P., & Congdon, C. E. 1994, ApJ, 435, 22 [Google Scholar]
- Kunth, D., Mas-Hesse, J. M., Terlevich, E., et al. 1998, A&A, 334, 11 [NASA ADS] [Google Scholar]
- Leitherer, C., & Heckman, T. M. 1995, ApJS, 96, 9 [NASA ADS] [CrossRef] [Google Scholar]
- Leitherer, C., Robert, C., & Drissen, L. 1992, ApJ, 401, 596 [NASA ADS] [CrossRef] [Google Scholar]
- Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 [NASA ADS] [CrossRef] [Google Scholar]
- Lequeux, J., Maucherat-Joubert, M., Deharveng, J. M., & Kunth, D. 1981, A&A, 103, 305 [NASA ADS] [Google Scholar]
- Madau, P., Pozzetti, L., & Dickinson, M. 1998, ApJ, 498, 106 [NASA ADS] [CrossRef] [Google Scholar]
- Maiz-Apellaniz, J., Mas-Hesse, J. M., Muñoz-Tuñón, C., Vilchez, J. M., & Castañeda, H. O. 1998, A&A, 329, 409 [NASA ADS] [Google Scholar]
- Mas-Hesse, J. M. 1992, A&A, 253, 49 [NASA ADS] [Google Scholar]
- Mas-Hesse, J. M., & Kunth, D. 1991, A&AS, 88, 399 [NASA ADS] [Google Scholar]
- Mas-Hesse, J. M., & Kunth, D. 1999, A&A, 349, 765 [NASA ADS] [Google Scholar]
- Mas-Hesse, J. M., Kunth, D., Tenorio-Tagle, G., et al. 2003, ApJ, 598, 858 [NASA ADS] [CrossRef] [Google Scholar]
- Mas-Hesse, J. M., Otí-Floranes, H., & Cerviño, M. 2008, A&A, 483, 71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mezger, P. O. 1978, A&A, 70, 565 [NASA ADS] [Google Scholar]
- Mezger, P. G., Mathis, J. S., & Panagia, N. 1982, A&A, 105, 372 [NASA ADS] [Google Scholar]
- Moustakas, J., Kennicutt R. C., Jr., & Tremonti, C. A. 2006, ApJ, 642, 775 [NASA ADS] [CrossRef] [Google Scholar]
- Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and Active Galactic Nuclei (California: University Science Books) [Google Scholar]
- Östlin, G., Hayes, M., Kunth, D., et al. 2009, AJ, 138, 923 [NASA ADS] [CrossRef] [Google Scholar]
- Pellerin, A., & Robert, C. 2007, MNRAS, 381, 228 [NASA ADS] [CrossRef] [Google Scholar]
- Ranalli, P., Comastri, A., & Setti, G. 2003, A&A, 399, 39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Rosa-González, D., Terlevich, E., & Terlevich, R. 2002, MNRAS, 332, 283 [NASA ADS] [CrossRef] [Google Scholar]
- Rosa-González, D., Terlevich, E., Jiménez Bailón, E., et al. 2009, MNRAS, 399, 487 [NASA ADS] [CrossRef] [Google Scholar]
- Rubin, R. H. 1968, ApJ, 154, 391 [NASA ADS] [CrossRef] [Google Scholar]
- Salim, S., Rich, R. M., Charlot, S., et al. 2007, ApJS, 173, 267 [NASA ADS] [CrossRef] [Google Scholar]
- Schaerer, D. 2002, A&A, 382, 28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Schaerer, D. 2003, A&A, 397, 527 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Storey, P. J., & Hummer, D. G. 1995, MNRAS, 272, 41 [NASA ADS] [CrossRef] [Google Scholar]
- Tammann, G. E. 1982, Supernovae: A Survey of Current Research, ed. M. J. Rees, & R. J. Stonehan (Dordrecht: Reidel), 371 [Google Scholar]
- Tüllmann, R., Breitschwerdt, D., Rossa, J., et al. 2006, A&A, 457, 779 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ulvestad, J. S. 1982, ApJ, 259, 96 [NASA ADS] [CrossRef] [Google Scholar]
- Verhamme, A., Schaerer, D., & Maselli, A. 2006, A&A, 460, 397 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Wilkins, S. W., Hopkins, A. M., Trentham, N., & Tojeiro, R. 2008, MNRAS, 391, 363 [NASA ADS] [CrossRef] [Google Scholar]
Footnotes
- ... models
- Downloadable from http://laeff.cab.inta-csic.es/users/mcs/
- ... models
- Downloadable from http://www.stsci.edu/science/starburst99/
- ... models
- Downloadable from http://obswww.unige.ch/sfr/sfr_tls/pop32/
All Tables
Table 1: SFS estimators which do not depend on E(B-V).
Table 2: SFS estimators which do depend on E(B-V) for an IB model at age = 4 Myr.
Table 3: SFS estimators which do depend on E(B-V) for an IB model at age = 5 Myr.
Table 4: SFS estimators which do depend on E(B-V) for an IB model at age = 6 Myr.
Table 5: SFR estimators which do not depend on E(B-V).
Table 6: SFR estimators which do depend on E(B-V) for an EB model at age = 10 Myr.
Table 7: SFR estimators which do depend on E(B-V) for an EB model at age = 30 Myr.
Table 8: SFR estimators which do depend on E(B-V) for an EB model at age = 250 Myr.
Table 9:
SFS estimators based on
and L1500 for Pop. III.
Table 10:
SFR estimators based on
and L1500 for Pop. III.
Table 11:
IMF correction factors for all SFS estimators, to be applied on values from Tables 1-4 when IMF mass limits 0.1-100
or 1-100
are assumed. See Sect. 4 for a more thorough discussion.
Table 12:
IMF correction factors for SFR estimators, to be applied on values from Tables 5-8 when IMF mass limits of 0.1-100
or 1-100
are considered. See Sect. 4 for a more thorough discussion.
Table 13:
IMF correction factors for SFS estimators, to be applied on values from Tables 1-4 when IMF slope ,
3.
Table 14:
IMF correction factors for SFR estimators, to be applied on values from Tables 5-8 when IMF slope ,
3.
All Figures
![]() |
Figure 1: Evolution of EW(WRbump) predicted by CMHK02 for IB (solid lines) and EB (dashed lines) models both for Z=0.020 ( top) and Z=0.008 ( bottom), together with the histogram of the EW(WRbump) data of the sample from Brinchmann et al. (2008). |
Open with DEXTER | |
In the text |
![]() |
Figure 2:
Evolution of
|
Open with DEXTER | |
In the text |
![]() |
Figure 3:
Evolution of
|
Open with DEXTER | |
In the text |
![]() |
Figure 4:
Evolution of
|
Open with DEXTER | |
In the text |
![]() |
Figure 5:
Evolution of
|
Open with DEXTER | |
In the text |
![]() |
Figure 6:
Evolution ( from top to bottom) of L1500, L2000 and L3500 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER | |
In the text |
![]() |
Figure 7:
Evolution ( from top to bottom) of L4400 and L5500 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER | |
In the text |
![]() |
Figure 8:
Evolution of
L22 200 predicted by CMHK02 for IB models ( left, normalized to
|
Open with DEXTER | |
In the text |
![]() |
Figure 9: Evolution of L22 200 predicted by SB99 for EB models as in Fig. 8. |
Open with DEXTER | |
In the text |
![]() |
Figure 10:
Evolution of
|
Open with DEXTER | |
In the text |
Copyright ESO 2010
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.