Issue 
A&A
Volume 599, March 2017



Article Number  A76  
Number of page(s)  13  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201629728  
Published online  02 March 2017 
The L−σ relation for HII galaxies in green
^{1} European Southern Observatory, Av. Alonso de Cordova 3107, Santiago, Chile
email: jmelnick@eso.org
^{2} Observatorio Nacional, Rua José Cristino 77, 20921400 Rio de Janeiro, Brasil
^{3} Cavendish Laboratory, University of Cambridge, 19 J. J. Thomson Ave, Cambridge CB3 0HE, UK
^{4} Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
^{5} Instituto Nacional de Astrofísica Óptica y Electrónica, AP 51 y 216, 72000 Puebla, Mexico
^{6} Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
^{7} Institute for Astronomy of the University of Hawaii, 2680 Woodlawn Drive, 96822 Honolulu, HI, USA
^{8} Physics Department, Aristotle University of Thessaloniki, 54124 Thessaloniki, Greece
^{9} Academy of Athens, Research center for Astronomy and Applied Mathematics, Soranou Efesiou 4, 11527 Athens, Greece
Received: 16 September 2016
Accepted: 6 December 2016
The correlation between emissionline luminosity (L) and profilewidth (σ) for HII galaxies provides a powerful method to measure the distances to galaxies over a wide range of redshifts. In this paper, we use SDSS spectrophotometry to explore the systematics of the correlation using the [OIII]5007 lines instead of Hα or Hβ to measure luminosities and line widths. We also examine possible systematic effects involved in measuring the profilewidths and the luminosities through different apertures. We find that the green L−σ relation, defined using [OIII]5007 luminosities, is significantly more sensitive than Hβ to the effects of age and the physical conditions of the nebulae, which more than offsets the advantage of the higher strength of the [OIII]5007 lines. We then explore the possibility of mixing [OIII]5007 profilewidths with SDSS Hβ luminosities using the Hubble constant H_{0} to quantify the possible systematic effects. We find the mixed L(Hβ) − σ_{[OIII]} relation to be at least as powerful as the canonical L−σ relation as a distance estimator, and we show that evolutionary corrections do not change the slope and the scatter of the correlation and, therefore, do not bias the L−σ distance indicator at high redshifts. Locally, however, the luminosities of the giant HII regions that provide the zeropoint calibrators are sensitive to evolutionary corrections and may bias the Hubble constant if their mean ages, as measured by the equivalent widths of Hβ, are significantly different from the mean age of the HII galaxies. Using a small sample of 16 adhoc zero point calibrators we obtain a value of H_{0} = 66.4^{+ 5.0}_{4.5}km s^{1} Mpc^{1} for the Hubble constant, which is fully consistent with the best modern determinations, and which is not biased by evolutionary corrections.
Key words: distance scale / dark energy / galaxies: starburst / galaxies: kinematics and dynamics
© ESO, 2017
1. Introduction
The relation between the integrated emissionline luminosity (L) and the velocitywidth (σ) of the emissionline profiles for HII galaxies provides a powerful method for measuring the distances to these compact strong emissionline objects, and this, out to redshifts of high cosmological relevance (z> 2; Chávez et al. 2016, and references therein).
In principle, the underlying astrophysics of the L−σ relation should be surprisingly simple: HII galaxies are powered by young starburst clusters that are a few million years old, so both the number of ionising photons and the turbulence of the gas are ultimately controlled by the total mass of the star cluster and gas complex.
A first indication that this naive scenario may be incomplete is that the scatter in the L−σ relation is substantially larger than the observational errors. Since the starburst clusters that ionise the nebulae evolve very rapidly in the first few million years, age is the natural culprit to explain the scatter. This effect has been well understood since the 1980s (e.g. Terlevich & Melnick 1981) and was analysed in detail by Chávez et al. (2014) who show that luminosity evolution, as parametrised by the equivalent widths of Hβ,does not explain the scatter in the relation.
Instead, Chávez et al. (2014) show that including the radius of the galaxy, R, as a second parameter significantly reduces the scatter as would be expected if HII galaxies are virialized systems. However, as pointed out by Terlevich et al. (2004), HII galaxies have complex starformation histories, and many are also powered by more than one young starburst, so the physical meaning of the observed radii is difficult to interpret.
In fact, a significant fraction of the objects in the sample compiled by Chávez et al. (2014) present multiple peaks in their emissionline profiles. Because of thermal broadening, these multiple profiles are often difficult to discern in the Balmer lines, but are clearly seen in the [OIII] lines, notably the green line at λ5007 Å, which in these objects have strengths comparable to Hα and are significantly stronger that any other line in the spectrum. Moreover, for the highest redshift objects (z> 2.5), [OIII]5007 is often the only line that can be measured with reasonable exposure times with stateoftheart instrumentation on 8−10 m class telescopes.
For these reasons, it is appealing to study the properties of the L−σ relation as seen in the green light of the [OIII]5007 emission line, which is the aim of this paper. In a parallel paper (FernandezArenas et al., in prep.) we present the analysis of our most recent application of the socalled canonical L−σ relation (i.e. using the Balmer lines) to measure the Hubble parameter (H_{0}) using a new sample of 36 giant HII regions in nearby galaxies as zeropoint calibrators. In the present analysis, we assume Ω_{Λ} = 0.71 and a flat Universe. H_{0} is, of course, left as a free parameter, but when required as a prior we adopt a value of H_{0} = 71 km s^{1} Mpc^{1}.
2. The data
We have combined data sets from our previous works, to which we have added new measurements as described below.
2.1. Bordalo and Telles
Bordalo & Telles (2011) observed a sample of 120 local HII galaxies (D< 500 Mpc) and we refer to the original paper for full details of the sample selection and the observational procedures. For the present purpose we recall that their velocity dispersions were determined with the FEROS fibrefed echelle spectrograph on the 1.52 m and 2.2 m telescopes on La Silla through a 2.7′′ entrance aperture. The fiber was positioned on the brightest knot of the starforming region (known as the kinematic core) that dominates the internal motions (Telles et al. 2001; Bordalo et al. 2009). The resolution of the spectra was ℛ = 48 000, which corresponds to an instrumental velocity dispersion σ_{inst} = 2.50 ± 0.20km s^{1}.
The spectrophotometry was obtained with the Boller and Chivens spectrograph on the 1.52 m telescope using a 2′′ slit (Kehrig et al. 2004). This information is relevant because we combine samples observed with different telescopes and different aperture sizes, and we must be aware of possible systematic effects. The present analysis is restricted to the best 70 objects from the sample of Bordalo & Telles (2011), as defined by V. Bordalo in his Ph.D. Thesis, using diagnostic diagrams to eliminate outliers, which we will refer to as the BT11 sample.
As we show below, the importance of including the BT11 sample is that it overlaps significantly in both luminosity and velocity dispersion with the giant HII regions used as zeropoint calibrators.
2.2. Chavez et al.
Chávez et al. (2014) studied a new sample of 128 HII galaxies: 122 selected from the SDSS on the basis of equivalent width of Hβ (EW(Hβ) > 50 Å) and redshift z (0.02 ≤ z ≤ 0.2), the latter chosen to minimize the effects of local peculiar velocities in particular induced by the Local Supercluster and the Great Attractor, as well as six galaxies in the southern hemisphere taken from the catalogue of Terlevich et al. (1991).
The emissionline profiles were obtained with two different spectrographs: HDS on Subaru and UVES on the VLT, while the total emissionline fluxes and equivalent widths were measured from lowdispersion spectrophotometry obtained with the Mexican 2 m telescopes at San Pedro Martir (Baja California) and Cananea (Sonora) through apertures larger than 8′′. The HDS observations were taken through a 4′′ slit, which nominally corresponds to an instrumental velocity dispersion σ_{inst} = 12.3 km s^{1}, although the instrumental profile is flattopped and not Gaussian. We restrict our analysis to the subsample of 107 galaxies, which is defined by Chávez et al. (2014) as their best data set, and which we refer to as the Chavez14 sample.
As shown in Fig. 1, the diameters of most objects in the Chavez14 sample are comparable to 4′′, so the instrumental resolution in the HDS spectra is determined not by the entrance slit but by the surface brightness profiles of the objects, which are very difficult to quantify using the available data. Chávez et al. (2014) corrected their observations for instrumental broadening using a sort of equivalent SDSS Petrossian diameter, but since this procedure can introduce a distancedependent error in the velocity dispersions, we chose to use a different approach, as discussed below.
The UVES observations were obtained through a 2′′ slit corresponding to a nominal instrumental resolution of σ_{inst} = 4.65km s^{1}, but again the instrumental profiles are also boxshaped. Although the seeing during the observations was substantially better than 2′′, the sizes of most objects are larger than the slit width, so no elaborate procedure was required for the instrumental corrections as was the case for the HDS data. We note, however, that both for UVES and HDS the flattopped instrumental profiles imply that the lineprofiles have nonGaussian cores that can be clearly appreciated in the residuals of the Gaussian fits (Chávez et al. 2014 and Appendix A).
Following the earlier work of Melnick et al. (1987), Chávez et al. (2014) measured their photometric fluxes through wide apertures. However, as shown in Fig. 1, the objects in the Chavez14 sample are quite compact, so the use of wide slits may introduce systematic aperture effects. Therefore, to combine the Chávez et al. (2014) and the Bordalo & Telles (2011) samples, it seems safer to match as closely as possible the photometric and spectroscopic apertures. Thus, for the Chavez2014 sample here, we use the spectrophotometry from SDSS taken through 3′′diameter fibre apertures.
Fig. 1
Distribution of SDSS Petrossian radii measured in the u′band. 
Their benchmark sample (S3) of 107 galaxies, defined by Chávez et al. (2014) as objects with log (σ) < 1.8km s^{1}, contains 104 objects from the SDSS and three southern galaxies for which no SDSS data is available. We downloaded the SDSS spectra and remeasured all relevant line intensities, in particular [OII], [OIII], and the Balmer lines, using a dedicated pipeline that rebins the spectra to zero redshift and corrects the fluxes for foreground extinction using redshifts and extinctions provided by in the DR9 release of the SDSS. We find that, even using the highquality SDSS spectra available for these objects, the ratios of Hγ to Hβ are much noisier than Hα to Hβ, so we have only used the latter to compute the internal extinction. We thus chose not to correct the Hβ fluxes for underlying absorption, as was done by Chávez et al. (2014) using the socalled Q method. Our final catalog of Chavez14/SDSS objects is listed in Table 1 and the BT11 objects are listed in Table 2.
Objects from Chávez et al. (2014).
2.3. Giant HII regions
Several authors (e.g. Hippelein 1986; Bordalo & Telles 2011) have noted that for HII galaxies (HIIGs) and giant HII regions (GHIIRs), the Balmer lines are systematically broader by about 2 km s^{1} than the [OIII] lines. Since the reasons for this difference are not well understood (see e.g. Hippelein 1986), for the purposes of calibrating the L−σ relation using GHIIRs it seems safer to rely on the measured [OIII]5007 linewidths rather than blindly applying an empirical correction^{1}.
Unfortunately, however, the highresolution spectroscopic observations of both Chávez et al. (2012) and FernandezArenas et al. (in prep.) used narrowband orderseparating filters to isolate Hα, so we must rely on the scanning FabryPerot observations of Hippelein (1986), which provide widths for both Hα and [OIII]5007 for 21 GHIIRs, for which spectrophotometry is available from the former authors. The FabryPerot data have the additional advantage that they are collected through wide apertures of diameters comparable to the spectrophotometric observations of Chávez et al. (2012) and FernandezArenas et al. (in prep.).
The GHIIRs in NGC 6822 have EW(Hβ) lower than our threshold of 46 Å (see below), while those of NGC 4256 lack equivalent width measurements in the literature, so we have excluded these four objects from the sample. We also eliminated one of the HII regions of M101, NGC 5447, because it is a complex of several HII regions that are elongated along a spiral arm, and not a bonafide single giant HII region. This leaves us with a sample of 16 zero point calibrators.
2.4. Line profiles of HII galaxies
As mentioned, most of the objects in our sample present some degree of nongaussianity, probably owing to multiple objects superimposed along the line of sight, to large expanding bubbles (which are actually clearly visible in HST images of local HII galaxies), or both. In fact, 3D kinematical studies show that the brightest knots in HII galaxies tend to have Gaussian profiles, while the more extended, low surface brightness regions have significantly broader profiles (see e.g. Moiseev et al. 2015, for a comprehensive discussion and references). Thus, the integrated profiles are expected to have one or more dominant Gaussian cores and extended nonGaussian wings.
Bordalo & Telles (2011) deal with this problem by eliminating galaxies with clear multiple profiles on the basis of an analysis of the profile moments, while Chávez et al. (2014) resort to fitting multiple Gaussians. Examples of the range of typical [OIII]5007 lineprofiles among the galaxies in our sample are presented in Appendix A.
It seems risky, however, to eliminate objects that have multiple profiles because multiplicity will sometimes appear as elevated velocity dispersions and sometimes as elevated luminosities depending on the relative radial velocities of the starbursts that are superposed along the line of sight. Therefore, to avoid introducing subjective biases in our analysis of the L−σ relation, it seems safer to pay the price of including all the objects in the sample: a larger scatter.
Still, we eliminated two galaxies with flagrant multiple profiles: J013258 and J142342 shown in Fig. A.1. The third “pathological” object in the figure, J212043, was not in the original Chavez14 sample. Thus, our final sample of Chavez14 with SDSS data − hereafter the Chavez14/SDSS sample − contains 102 galaxies.
To deal with nonGaussian profiles, we chose to use GaussHermite functions that are commonly used for other problems in astronomy (e.g. Riffel 2010, and references therein). Thus, we fitted the original FEROS, HDS, and UVES [OIII]5007 lines with GaussHermite functions using the package PAN as implemented in IDL. Details of this can be found in Westmoquette et al. (2007).
To minimise possible contamination by underlying old stellar populations, we restricted the sample to objects with equivalent widths, EW(Hβ) > 46 Å, following the rationale explained by Chávez et al. (2014)^{2}. This reduces the sample to 32 objects from which we further eliminated two more: IIZw40, which is very extended, and UM559, which has an anomalous EW(Hβ) of 535 Å. This reduces our final sample of BT11 galaxies to 30 objects.
Objects from Bordalo & Telles (2011).
2.5. Errors
2.5.1. Lineprofile widths
Our treatment of the FEROS data is very similar to that of Bordalo & Telles (2011; having been performed by the same person!), so we retain the original errors that were evaluated from measurements of different lines of hydrogen and oxygen as a function of the S/N as an internal test, and by measuring different spectra of the same object as an external test.
The same is not possible for the Chavez14 data because the [OIII]5007 lines have significantly higher S/N than Hβ so we need to propagate the errors including the corrections for instrumental (σ_{ins}) and thermal (σ_{th}) broadening. The thermal corrections are small for the oxygen lines, so can safely adopt a constant electron temperature of T_{e} = 12 100 K ± 1440 K for all galaxies (that corresponds to the mean and standard deviation of the Chavez14/SDSS sample as measured by our pipeline).
For the instrumental profiles, we adopt: σ_{FEROS} = 2.5 ± 0.2km s^{1}; σ_{UVES} = 4.65 ± 2.32km s^{1}; and σ_{HDS} = 6.15 ± 3.07km s^{1}. This is equivalent to assuming that, on average, the objects in the Chavez14 sample have effective diameters of r ≃ 2′′ ± 1′′. We do not have the same problem with the FEROS data because the spatial information is efficiently scrambled by the fibers. From the standard formula for the velocity dispersion, , the corresponding errors are calculated in the usual way by the square root of the sum of the variances in each variable.
We note that our final velocity dispersion errors are significantly smaller than those of Chávez et al. (2014) for Hβ. However, the [OIII]5007 lines are, on average, five times stronger than Hβ and the errors become rather comparable if we multiply our values by .
2.5.2. Luminosities
For our Chavez14/SDSS sample, the errors in luminosity are dominated by the errors in the fluxes and the extinction corrections. As mentioned before, we only used Hα/Hβ for the extinction corrections, whereas the errors in the fluxes were determined from the SDSS spectra using the statistical weights for each pixel given by the SDSS pipeline. For the BT11 data, we use their indication that their fluxes, published by Kehrig et al. (2004), have typical errors of 5%, which we propagate to the luminosities using the dispersion in the Balmer decrements from Hα to Hδ weighted by the corresponding line intensities to estimate the errors in the extinction. As for the Chavez14 sample, we adopted the 30 Dor/SMC extinction law of Gordon et al. (2003), thus,
(1)where F(Hβ) is the observed flux, A_{β} is the extinction at 4861 Å, and δ denotes the measurement errors in each variable.
We recall that, while the Chavez14 sample was selected to minimise distance errors owing to local perturbations of the Hubble flow, so we can safely ignore distance errors, the Bordalo & Telles (2011) sample contains many very local objects. As such, we are probably underestimating the errors in distance and hence luminosity.
3. Results
3.1. The L−σ relation in green
We begin by examining the L−σ relation for the [OIII]5007 line, which we will call the green L−σ relation as opposed to the socalled canonical relation, where we use Hβ for both the luminosities and the velocity dispersions. The green relation is presented in Fig. 2. The three lines correspond to the direct (L vs. σ) and inverse (σ vs. L) standard leastsquares (LSQ) fits, both ignoring the observational errors, and to a maximumlikelihood (ML) fit that considers the observational errors in both variables. For the latter we used an iterative algorithm that was originally proposed by Leon Lucy (priv. comm.) and described in detail in Appendix C. The relevant parameters for the three fits are listed in Table D.1.
Throughout this paper we use the ML fit to define the slope of the L−σ relation. We note that, as expected since the errors in the independent variable (σ) are substantially larger than the errors in the dependent variable (L), the slope of the ML fit is very similar to the slope of the inverse LSQ fit.
Fig. 2
Green L−σ relation for our 102 HII galaxies with SDSS [OIII]5007 luminosities and GaussHermite velocity dispersions. The black solid and dashed lines represent direct (L vs. σ) and inverse (σ vs. L) standard Least Squares (LSQ) fits to the data. The orange line shows a maximumlikelihood (ML) fit that takes the observational errors into account as described in Appendix B. The relevant parameters of the three fits are listed in Table D.1. The 30 galaxies from the BT11 sample are not plotted because we lack [OIII]5007 luminosities. 
The slope of the green L−σ relation is the same as that of the canonical relation (4.65 ± 0.14; Chávez et al. 2014), but steeper than the slope that we obtain using the SDSS Hβ luminosities (see below). The scatter of the green relation, however, is significantly larger than for Hβ. Chávez et al. (2014) studied the sources of scatter in the L−σ relation using multiparameter fits and concluded that the radius of the HII galaxies, measured on SDSS u′band images, seems to be the dominant second parameter, as would be expected for virialized systems, and that age, metallicity, and excitation also contribute to the scatter.
The radius R also appears to be the dominant second parameter for the green relation. Including log R as a second parameter reduces the rms scatter to 0.22 (compared to 0.20 for Hβ)^{3}. The other two parameters that contribute significantly to the scatter of the green relation are age, as parametrised by the equivalent widths of Hβ, EW(Hβ), and excitation, [OII]/[OIII].
Thus, fitting a linear equation of the form, (2)using standard LSQ techniques (i.e. ignoring the observational errors) yields (3)Introducing the radius as a second parameter in the distance indicator is not practical because it not only weakens the distance sensitivity, but mostly because radii are difficult to measure even at low redshifts. Introducing the other two parameters is also not practical because both require measuring lines ([OII]3727 and Hβ) that are substantially weaker than [OIII]5007 thus completely erasing the advantage of using [OIII] (using Hβ luminosities has the additional advantage that the evolutionary signal (c_{3}) is significantly weaker than that of L([OIII]) ). We conclude that because of its significantly larger scatter and sensitivity to additional parameters, the green L−σ relation is not as useful as a distance indicator as the canonical relation. In the next section we examine the possibility of combining [OIII]5007 lineprofile widths with Hβ luminosities.
3.2. Back to Hβ luminosities: the mixed L−σ relation
Figure 3 shows SDSS Hβ luminosities, L(Hβ), plotted as a function of velocity dispersion from our GaussHermite fits to the [OIII]5007 line profiles. The lines show the same three fits described in the previous section. We recall that the luminosities are obtained using our new measurements of the fluxes from the SDSS spectra, and therefore are not necessarily the same as those given by Chávez et al. (2014).
Fig. 3
Mixed L−σ relation for the Chavez/SDSS sample of 102 HII galaxies. As in the previous figure, the lines represent the direct, inverse LSQ, and ML fits and the slope correspond to the value from the ML fit. We have not included the 30 galaxies from our BT11 sample to facilitate comparison with the green relation. The relevant parameters of the three fits are given in Table D.1. 
The slope of the mixed correlation for the Chavez14/SDSS sample is somewhat flatter than that of the green relation, but the scatter is significantly lower. The evolutionary signal from Eq. (2) is significantly smaller than for [OIII]5007 luminosities: c_{3} = 0.0027 ± 0.0004, indicating that, as expected, the [OIII]5007 luminosities evolve faster than L(Hβ) as the clusters age. Curiously, the effect of excitation ([OII]/[OIII]), c_{4} = 0.25 ± 0.04, is stronger for L(Hβ) than for L([OIII]).
4. The Hubble constant
The purpose of this paper is not to measure the Hubble constant (in particular because we only have a restricted sample of giant HII regions with [OIII] velocity dispersions), but rather to explore the systematics of the method. The determination of the Hubble constant using the L−σ relation was discussed in Chávez et al. (2012) and is revisited in a separate paper (FernandezArenas et al., in prep.).
However, since we actually know the value of H_{0} pretty accurately (Riess et al. 2016; FernandezArenas et al., in prep.), we can use the value of H_{0} that results from our relations as a parameter to quantify the systematics of the L−σ relation. We thus use the slope (α) of the L−σ relation together with the data for our sample of 16 giant HII regions with accurate Cepheid distances to calibrate the zero point (Z_{p}) of the distance indicator as follows (4)where L_{GHR,i} is the Hβ luminosity of each giant HII region and σ_{GHR,i} the corresponding velocity dispersion. The statistical weights W_{i} are calculated as (5)We stress that, although the observational errors are not correlated, the covariance between L and σ must be considered in the propagation of errors as described in Appendix B. Thus, from Eq. (B.5) here ⟨ log σ_{HIIG} ⟩ is the average velocity dispersion of the sample of HII galaxies used to define the slope of the relation.
The calibrated L−σ relation is then: log L(Hβ) = αlog σ + Z_{p} and to calculate the Hubble constant we minimise the function, (6)where μ_{i} is the logarithmic distance to each HII galaxy calculated from the distance indicator and the reddeningcorrected Hβ fluxes F(Hβ) as (7)and μ_{h,i} is the distance from the redshifts (z_{i}), using the standard cosmological equations for the luminositydistance, (8)where the independent variable h− the Hubble constant − is embedded in the equation for the luminosity distance, which we calculated using the analytical approximation of Wickramasinghe & Ukwatta (2010), assuming a flat universe and Ω_{Λ} = 0.71. This choice has a negligible effect in the distance moduli at the redshifts of our objects.
The best value of h is then calculated minimising χ^{2} with statistical weights calculated as (9)Figure 4 repeats Fig. 3 but this time including the 30 objects from our BT11 sample, and the 16 giant HII region from our Hippelein/Arenas sample for the determination of the zero point, and the corresponding value for the Hubble constant and χ^{2} curve.
Fig. 4
Mixed L−σ relation for the combined Chavez14/SDSS (black squares) and BT11 sample (red triangles) of 132 HII galaxies. The lines show the LSQ and ML fits as in the previous figures and the relevant parameters are given in Table D.1. The green dots correspond to the 16 giant HII regions of the Hippelein/Arenas sample used to calibrate the zero point of the relation and thus establish the value of H_{0}. The inset shows the χ^{2} curve (actually ) and the value of H_{0} determined from a cubic fit to the curve. 
We emphasise that we do not have a large enough sample of calibrator giant HII regions with [OIII] velocity dispersions, and therefore that our values for H_{0} are only intended as a way to quantify the possible systematic effects involved in calibrating the L−σ relation as a distance indicator.
5. Systematics
5.1. Aperture matching
From our initial studies in the 1980s, our conventional wisdom has been to measure the velocity dispersions through relatively narrow slits to preserve the spectral resolution, and the luminosities through wide apertures to include all the flux. The underlying assumption was that the turbulence is isotropic and the internal extinction modest so that even through a narrow slit we still sample the full turbulent cascade. This is a pretty good assumption for single objects, but 50% of our objects have complex profiles, and half of these show multiple peaks. Thus, even narrow slits encompass more than one starburst along the line of sight (see Fig. A.1), so using different entrance apertures for luminosities and velocity dispersions may introduce systematic effects.
For this reason here we chose to use the SDSS spectrophotometry, that matches rather well the apertures used for the velocity dispersions (2′′−4′′), to measure luminosities. The fact that the galaxies in our high EW subsample of BT11 objects, observed through a 2′′ slit, are not systematically shifted relative to the Chavez14/SDSS sample despite being significantly more local, and that the scatter of the L−σ relation in Chávez et al. (2014) – using wide apertures – is almost the same as that of our mixed relation, indicates that the SDSS fluxes do not miss a significant fraction of the light.
Spatially resolved kinematic studies of HII galaxies (e.g.. Moiseev & Lozinskaya 2012) indicate that the integrated profile widths of these objects should, in principle, become broader as a function of aperture size. However, in practice the integrated profiles are dominated by the brightest cores of the galaxies, while the extended low surfacebrightness regions appear as nonGaussian wings that are effectively taken into account by our GaussHermitte functions.In fact, Lagos et al. (2007) show, using narrowband imaging of some of the galaxies in the BT11 sample, that even in these local objects the luminosities increase only marginally for apertures larger than 2′′−3′′, while the equivalent widths remain substantially unchanged.
5.2. Age effects
The ionising fluxes of starburst clusters decrease steeply after about 3 Myr, and thus their emissionline luminosities. The objects in our sample, on the other hand, span a range of ages, as estimated from the equivalent widths of the Balmer lines (Leitherer et al. 1999) from 2.5 to 5 Myr so, in principle, we must correct the Hβ luminosities for evolution. Figure 5 presents the evolution of L(Hβ). This is as a function of equivalent width EW(Hβ) computed for single starbursts using STARBURST99 (Leitherer et al. 1999) with Geneva isochrones without rotation and a standard Kroupa IMF, and for a range of abundances spanning the values observed for our HIIGs.
Fig. 5
Hβ luminosity evolution of single starburst clusters computed from Starburst99 models with parameters indicated in the legends. The lines show the fit for 45 Å <EW(Hβ) ≤ 175 Å and EW(Hβ) > 175 Å. 
For all abundances, the run of luminosity with equivalent width can be approximated reasonably well by two straight lines of slopes:
over the range of equivalent widths spanned by our objects.
Luminosity evolution, therefore, should leave a strong imprint in the scatter of the L−σ relation. However, the turbulence of the gas may also evolve as the ionising clusters fade although we lack a full understanding of the origin of the supersonic gas motions observed in GHIIRs and HIIGs. The gravitational potential of the ionising cluster contributes a large fraction of the energy, but detailed studies of nearby objects, such as 30 Doradus in the LMC and NGC 604 in M 33, show that gravity is only part of the story, and that hydrodynamic effects − winds, bubbles, filaments − also play an important role (e.g. Medina Tanco et al. 1997; Melnick et al. 1999).
Unfortunately, however, EW(Hβ) is not a precise age indicator for HII galaxies (Terlevich et al. 2004). IR imaging by Lagos et al. (2011) shows that HIIGs harbour many old and intermediateage clusters and, in fact, the majority of the galaxies in our sample show clear absorption components in Hδ, which is the telltale sign of a significant intermediateage stellar population. This is confirmed by population synthesis models of the broadband spectral energy distributions (SED) of 238 SDSS HII galaxies that span more than two decades in wavelength (from the FUV at 0.15 μm to the midIR at 22 μm; Telles & Melnick, in prep.).
Fig. 6
Mixed Lσ relation corrected for evolution at a fiducial “age” of ⟨ EW ⟩ = 133 Å (corresponding to about 3.5 Myr for the Geneva isochrones). The remaining details of this figure are similar to those of Fig. 4, except that here the mean value of EW(Hβ) is shown for each subsample and for the complete sample of HII galaxies. The relevant parameters of the fits are listed in Table D.1. 
On the other hand, at least in principle, GHIIRs should not be affected by this problem and, in fact, correcting their observed luminosities for evolution using the models shown in Fig. 5 leads to a significant reduction in the scatter of the L−σ plane. Therefore, unless we correct both HIIGs and GHIIRs to the same fiducial age, we introduce a systematic bias in the calibration of the zero point and thus in the value of the Hubble constant.
The mean equivalent width of our sample of 16 calibrators is ⟨ EW ⟩ = 133 Å, compared to ⟨ EW ⟩ = 95 Å for our 132 HII galaxies. Preliminary results from our population synthesis models (Telles & Melnick, in prep.) indicate that, on average, the observed EW(Hβ) must be corrected by between 30 Å and 50 Å (depending on the details of the modelling) to account for the contribution of the underlying intermediateage stellar populations to the continuum at Hβ. Thus, our best estimate for the average equivalent width of the starburst component(s) of our sample is between 125 Å and 145 Å.
Figure 6 shows the luminosities corrected to a fiducial age of ⟨ EW ⟩ = 133 Å (corresponding to ~3.5 Myr for the Geneva models). The figure shows that the slope, the zeropoint (i.e. H_{0}), and the scatter of the L−σ relation remain practically unchanged by the evolutionary corrections (see Table D.1 for a comparison of the different slopes). However, if we correct the luminosities to the lower value of the mean corrected equivalent width of our HII galaxies (⟨ EW ⟩ = 125 Å), we obtain while, for ⟨ EW ⟩ = 145 Å, we find , on average not significantly different from our uncorrected value of , shown in Fig. 4. The smaller errors of the corrected values are due to the reduced scatter of the zeropoint calibrators (GHIIRs) resulting from the agecorrected luminosities, although these corrections introduce an additional uncertainty owing to the correction of EW(Hβ) for the contribution of the intermediateage populations in the HII galaxies.
Clearly, therefore, unless we can estimate accurately the ages of the starburst components of HII galaxies, we cannot refine our estimation of H_{0} to better than 4−5 km s^{1} Mpc^{1}.
6. Conclusions
We have studied the systematics of the L−σ relation for HII galaxies, and investigated the possibility of using the green lines of [OIII] at 5007 Å to define both the velocity dispersion (σ) and the luminosity (L) of the correlation. Our results can be summarised as follows:

The scatter of the L−σ relation using only [OIII]5007 is substantially larger than for L(Hβ). This is mostly because the [OIII] luminosities of the starburst that power HII galaxies evolve faster than the Hβ luminosities as the ionising stars age. Thus, while the advantages of using [OIII]5007 for measuring σ are strong, the same is not the case for the luminosities. The green L−σ relation is not as good as the canonical L−σ relation as a distance indicator;

Combining Hβ luminosities with [OIII] velocity dispersions provides the best compromise; the scatter of this mixed L−σ relation is somewhat smaller than that of the canonical relation using Hβ for both parameters, with the additional advantage that, in HII galaxies, [OIII]5007 is significantly stronger than Hβ and much less affected by thermal broadening;

Our mixed L−σ relation, established using GaussHermite fits to the emissionline profiles and SDSS spectrophotometry, is at least as good, if not better, than the canonical relation as a distance indicator that uses wideaperture photometry and a completely different technique to fit the lineprofiles. This demonstrates that the emissionline luminosities of HII galaxies are dominated by the compact cores of these objects, and that the influence of the diffuse extended components is negligible;

The scatter of the mixed L−σ relation is substantially larger than the observational errors. This is partly due to the presence of one or more additional parameters in the relation. We confirmed the results of Chávez et al. (2014), i.e. the radius, as measured on SDSS u′band images, appears to be a second parameter of the correlation. We also found that starburst age, as measured by the equivalentwidth of Hβ EW(Hβ), has no influence on the scatter of the relation, contrary to what is expected from stellarevolutionary models. However, EW(Hβ) is only an approximate ageindicator for HII galaxies that contain massive intermediateage populations, which contribute significantly to the continuum at Hβ;

Using an ad hoc sample of 16 giant HII regions to determine the zero point, we used the calibrated mixed L−σ relation for a sample of 132 HII galaxies to measure the Hubble constant. Our value agrees within the errors with the most recent determinations in the literature. We showed that these errors can be reduced by correcting the Hβ luminosities of the zeropoint calibrators (giant HII regions) for evolutionary dimming. This, however, requires that similar corrections be applied to the HII galaxies, which could only be performed approximately. Our preliminary results suggest that our value for H_{0} is probably not biased by age corrections, but this needs to be confirmed by detailed populationssynthesis modelling.
The main result from this investigation is that the L−σ relation is robust to a number of potentially serious systematic biases, particularly aperturematching and nonGaussian line profiles. While the scatter of the relation is larger than that of SNIa, for the purposes of measuring the Hubble constant, HII galaxies are subject to different and independent systematics, while for measuring the expansion history of the Universe, HII galaxies can be observed over a wider range of redshifts. The full power of the L−σ relation to measure the Hubble constant will be presented in a forthcoming paper (FernandezArenas et al., in prep.).
Acknowledgments
J.M. acknowledges support from a CNPq Ciência sem Fronteiras grant at the Observatorio Nacional in Rio de Janeiro, and the hospitality of ON as a PVE visitor. E.T. and J.M. acknowledge the hospitality of INAOE in Puebla, during which parts of this research were conducted. R.C. is grateful to the Mexican research council (CONACYT) for supporting this research under Grant 263561. Funding for the SDSS and SDSSII has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/. The SDSS is managed by the Astrophysical Research Consortium for the Participating Institutions. The Participating Institutions are the American Museum of Natural History, Astrophysical Institute Potsdam, University of Basel, University of Cambridge, Case Western Reserve University, University of Chicago, Drexel University, Fermilab, the Institute for Advanced Study, the Japan Participation Group, Johns Hopkins University, the Joint Institute for Nuclear Astrophysics, the Kavli Institute for Particle Astrophysics and Cosmology, the Korean Scientist Group, the Chinese Academy of Sciences (LAMOST), Los Alamos National Laboratory, the MaxPlanckInstitute for Astronomy (MPIA), the MaxPlanckInstitute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.
References
 Bordalo, V., & Telles, E. 2011, ApJ, 735, 52 [NASA ADS] [CrossRef] [Google Scholar]
 Bordalo, V., Plana, H., & Telles, E. 2009, ApJ, 696, 1668 [NASA ADS] [CrossRef] [Google Scholar]
 Chávez, R., Terlevich, E., Terlevich, R., et al. 2012, MNRAS, 425, L56 [NASA ADS] [CrossRef] [Google Scholar]
 Chávez, R., Terlevich, R., Terlevich, E., et al. 2014, MNRAS, 442, 3565 [NASA ADS] [CrossRef] [Google Scholar]
 Chávez, R., Plionis, M., Basilakos, S., et al. 2016, MNRAS, 462, 2431 [NASA ADS] [CrossRef] [Google Scholar]
 Deming, W. 1964, Statistical Adjustment of Data, 3rd edn. (Dover), 291 [Google Scholar]
 Gordon, K. D., Clayton, G. C., Misselt, K. A., Landolt, A. U., & Wolff, M. J. 2003, ApJ, 594, 279 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hippelein, H. H. 1986, A&A, 160, 374 [NASA ADS] [Google Scholar]
 Kehrig, C., Telles, E., & Cuisinier, F. 2004, AJ, 128, 1141 [NASA ADS] [CrossRef] [Google Scholar]
 Lagos, P., Telles, E., & Melnick, J. 2007, A&A, 476, 89 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Lagos, P., Telles, E., NigocheNetro, A., & Carrasco, E. R. 2011, AJ, 142, 162 [NASA ADS] [CrossRef] [Google Scholar]
 Leitherer, C., Schaerer, D., Goldader, J. D., et al. 1999, ApJS, 123, 3 [NASA ADS] [CrossRef] [Google Scholar]
 McDonald, J., & Thompson, J. 1992, Am. J. Phys., 60, 66 [NASA ADS] [CrossRef] [Google Scholar]
 Medina Tanco, G. A., Sabalisck, N., JatencoPereira, V., & Opher, R. 1997, ApJ, 487, 163 [NASA ADS] [CrossRef] [Google Scholar]
 Melnick, J., Moles, M., Terlevich, R., & GarciaPelayo, J.M. 1987, MNRAS, 226, 849 [NASA ADS] [CrossRef] [Google Scholar]
 Melnick, J., TenorioTagle, G., & Terlevich, R. 1999, MNRAS, 302, 677 [NASA ADS] [CrossRef] [Google Scholar]
 Moiseev, A. V., & Lozinskaya, T. A. 2012, MNRAS, 423, 1831 [NASA ADS] [CrossRef] [Google Scholar]
 Moiseev, A. V., Tikhonov, A. V., & Klypin, A. 2015, MNRAS, 449, 3568 [NASA ADS] [CrossRef] [Google Scholar]
 Riess, A. G., Macri, L. M., Hoffmann, S. L., et al. 2016, ApJ, 826, 56 [NASA ADS] [CrossRef] [Google Scholar]
 Riffel, R. A. 2010, Ap&SS, 327, 239 [NASA ADS] [CrossRef] [Google Scholar]
 Telles, E., MuñozTuñón, C., & TenorioTagle, G. 2001, ApJ, 548, 671 [NASA ADS] [CrossRef] [Google Scholar]
 Terlevich, R., & Melnick, J. 1981, MNRAS, 195, 839 [NASA ADS] [CrossRef] [Google Scholar]
 Terlevich, R., Melnick, J., Masegosa, J., Moles, M., & Copetti, M. V. F. 1991, A&AS, 91, 285 [NASA ADS] [Google Scholar]
 Terlevich, R., Silich, S., RosaGonzález, D., & Terlevich, E. 2004, MNRAS, 348, 1191 [NASA ADS] [CrossRef] [Google Scholar]
 Tremaine, S., Gebhardt, K., Bender, R., et al. 2002, ApJ, 574, 740 [NASA ADS] [CrossRef] [Google Scholar]
 Westmoquette, M. S., Exter, K. M., Smith, L. J., & Gallagher, J. S. 2007, MNRAS, 381, 894 [NASA ADS] [CrossRef] [Google Scholar]
 Wickramasinghe, T., & Ukwatta, T. N. 2010, MNRAS, 406, 548 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Profiles
Figure A.1 illustrates the zoo of emissionline shapes encountered in our sample of 132 HII galaxies. The solid lines show the GaussHermite profilefits to the data, and the lower panel in each figure shows the residuals from the fit. The majority of the objects in our sample show some kind of profile asymmetry, be it extended wings or multiple components. Some cases, like J142342 in the lowerright corner of the figure, show two components of very similar strength, plus a third of slightly lower intensity. In this case we are integrating three bursts of comparable ages along the line of sight.
We have verified that, while in some cases the residuals in the L−σ relation are correlated with complex profiles, this is not usually the case. Therefore, in our study we have made no attempt to select objects on the basis of the profile shapes. This clearly increases the scatter, but eliminates any subjective bias involved in the selection of acceptable profiles.
Fig. A.1
Examples of the different emissionline profile shapes in our sample of HII galaxies, together with the GaussHermite fits. The lower panel shows the residuals from the fit. 
The figure also illustrates the usefulness of using [OIII]5007 instead of Hβ. The double structure of J002425 that is clearly discerned in [OIII]5007, is almost impossible to detect in Hβ due to the larger thermal broadening. It is also interesting to notice that the prototypical nearby HII Galaxy IIZw40, whose emissionline luminosity is dominated by a very compact core, shows a clear profile asymmetry even through the 2.7′′ FEROS entrance aperture.
Appendix B: Propagation of errors
When we use an empirical relation, such as L−σ, to predict a value of the dependent variable we must include the covariance of the two variables when we propagate the errors (Deming 1964). In general, the error in the prediction of a linear correlation of the form (B.1)at a given value of x = x_{m}, δY_{m}, when the parameters a and b are determined using leastsquares techniques, can be expressed as (Deming 1964) (B.2)where and are the variances of the two variables, ⟨ x ⟩ is the mean value of x, N is the number of data points, and r^{2} is the correlation coefficient defined as (B.3)For our particular application it is convenient to recall that for standard leastsquares, the errors in the coefficients are given by (Deming 1964), (B.4)so the error in the value predicted for a measurement x_{m} can also be expressed as, (B.5)In order to minimise the effect of the covariance it is useful to replace the independent variable x by x′ = x− ⟨ x ⟩ so that ⟨ x′ ⟩ = 0. This change of variables does not affect the error in the slope (δb), but does change the error in the zero point δa. In cases like ours where the zero point Z_{p} is determined from a set of independent calibrators (in our case giant HII regions), the covariance error is given by, (B.6)When the data are subject to experimental errors in both variables (δx_{m},δy_{m}), the error in the prediction becomes, (B.7)As we will discuss below, the standard leastsquares solution is in general not unbiased when the independent variable is subject to error. In such cases Eq. (B.7) is a very good approximation, but is not exact.
Appendix C: The maximum likelihood method
There is a substantial body of literature dealing with the question of how to fit straight lines to data affected by observational errors in both the dependent and the independent variables (see e.g. McDonald & Thompson 1992; Tremaine et al. 2002). For historical reasons, however, in this paper we have used the method originally suggested to us by Leon Lucy back in 1985, which is the one we used in the first papers dealing with the calibration of the L−σ relation as a distance indicator (Melnick et al. 1987). In fact, here we actually used the same FORTRAN code used in our early publications.
The central assumption of the method is that the errorfree variables (ξ_{i},η_{i}) fit exactly a straight line of the form, (C.1)The actual observations (x_{i},y_{i}) are subject to errors with standard deviations (σ_{xi},σ_{yi}), so the likelihood L of the observational sample is, (C.2)or in logarithmic form, (C.3)where and and C is a constant. Taking partial derivatives with respect to m,b, and ξ we derive the equations, (C.4)(C.5)(C.6)which we solve by iteration. First guess the value of m and b, usually by standard leastsquares fits, and with these values estimate ξ_{i} from Eq. (C.4). Then obtain new values for m and b using Eqs. (C.5) and (C.6) and iterate the procedure until it converges.
To estimate the errors in the coefficients we make use of the basic formula for the variance of a quantity θ, (C.7)where is the MaximumLikelihood estimate of the variable θ.
Appendix D: Parameters for the different fits
The relevant parameters (slopes and rms) for the different linear fits employed in the text are presented in Table D.1 to facilitate comparison between the different methods, and between the different parametrizations of the L−σ relation as described in the text.
The zero points for each fit are arbitrary and are not used in the determination of the Hubble constant, which is done using Giant HII Regions as zeropoint calibrators.
Fit parameters.
All Tables
All Figures
Fig. 1
Distribution of SDSS Petrossian radii measured in the u′band. 

In the text 
Fig. 2
Green L−σ relation for our 102 HII galaxies with SDSS [OIII]5007 luminosities and GaussHermite velocity dispersions. The black solid and dashed lines represent direct (L vs. σ) and inverse (σ vs. L) standard Least Squares (LSQ) fits to the data. The orange line shows a maximumlikelihood (ML) fit that takes the observational errors into account as described in Appendix B. The relevant parameters of the three fits are listed in Table D.1. The 30 galaxies from the BT11 sample are not plotted because we lack [OIII]5007 luminosities. 

In the text 
Fig. 3
Mixed L−σ relation for the Chavez/SDSS sample of 102 HII galaxies. As in the previous figure, the lines represent the direct, inverse LSQ, and ML fits and the slope correspond to the value from the ML fit. We have not included the 30 galaxies from our BT11 sample to facilitate comparison with the green relation. The relevant parameters of the three fits are given in Table D.1. 

In the text 
Fig. 4
Mixed L−σ relation for the combined Chavez14/SDSS (black squares) and BT11 sample (red triangles) of 132 HII galaxies. The lines show the LSQ and ML fits as in the previous figures and the relevant parameters are given in Table D.1. The green dots correspond to the 16 giant HII regions of the Hippelein/Arenas sample used to calibrate the zero point of the relation and thus establish the value of H_{0}. The inset shows the χ^{2} curve (actually ) and the value of H_{0} determined from a cubic fit to the curve. 

In the text 
Fig. 5
Hβ luminosity evolution of single starburst clusters computed from Starburst99 models with parameters indicated in the legends. The lines show the fit for 45 Å <EW(Hβ) ≤ 175 Å and EW(Hβ) > 175 Å. 

In the text 
Fig. 6
Mixed Lσ relation corrected for evolution at a fiducial “age” of ⟨ EW ⟩ = 133 Å (corresponding to about 3.5 Myr for the Geneva isochrones). The remaining details of this figure are similar to those of Fig. 4, except that here the mean value of EW(Hβ) is shown for each subsample and for the complete sample of HII galaxies. The relevant parameters of the fits are listed in Table D.1. 

In the text 
Fig. A.1
Examples of the different emissionline profile shapes in our sample of HII galaxies, together with the GaussHermite fits. The lower panel shows the residuals from the fit. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext 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 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.