A&A 487, 147-160 (2008)
DOI: 10.1051/0004-6361:200809706
O. Vaduvescu1 - M. L. McCall2
1 - Instituto de Astronomía, Universidad Católica del Norte,
Avenida Angamos 0610, Antofagasta, Chile
2 -
York University, Department of Physics and Astronomy,
4700 Keele Street, M3J 1P3, Toronto, ON, Canada
Received 3 March 2008 / Accepted 19 May 2008
Abstract
Aims. Near-IR photometry of dwarf irregular galaxies (dIs) has been acquired to expand the sample of objects suitable for investigating the fundamental plane (FP) discovered previously by us. Data for blue compact dwarfs (BCDs), dwarf ellipticals (dEs), and dwarf spheroidal galaxies (dSphs) are amalgamated to evaluate how closely these classes of galaxies are related to dIs.
Methods. Surface brightness profiles for dIs are modeled using a hyperbolic secant function (sech), which has been shown to sample the old galaxy component observed in the NIR (). Also, profiles for BCDs are modeled using a sech function, but to this is added a Gaussian burst. Kinematics for dIs and BCDs are constrained by the width W20 of integrated HI line profiles, and motions within dEs and dSphs are judged from stellar velocity dispersions.
Results. The existence of a FP for dIs is confirmed. Omitting extreme deviants, the scatter of dIs about the plane is 0.43 mag, somewhat higher than expected from observational errors alone. Corrections of line widths for tilt do not reduce the residuals, so motions must be predominantly random. Corrections of surface brightnesses for tilt have no effect either. BCDs and dEs lie precisely on the FP defined by dIs, suggesting strong underlying connections among dIs, BCDs, and dEs.
Key words: galaxies: dwarf - galaxies: irregular - galaxies: photometry - galaxies: structure - infrared: galaxies - cosmology: distance scale
Since the discovery of the Faber-Jackson relation for elliptical galaxies (Faber & Jackson 1976) and the Tully-Fisher relation for spiral galaxies (Tully & Fisher 1977), many authors have looked for inter-connections among the physical parameters of massive galaxies. Djorgovski & Davis (1987) were the first to recognize the so-called ``fundamental plane'' of elliptical galaxies, which links kinematics to length and brightness scales. During the last few years, fundamental planes have been identified for other massive entities such as spiral galaxies (Han et al. 2001), early type galaxies (Reda et al. 2005), radio galaxies (Bettoni et al. 2003), AGNs (Woo et al. 2004; Wang et al. 2006), quasars (Hamilton et al. 2006), and black holes (Kording et al. 2006; Merloni et al. 2003). Although Kormendy (1985) did a pioneering study of pair-wise correlations of parameters characterizing dwarfs, more sophisticated relationships have been elusive. What is clear is that the Tully-Fisher relation breaks down at low luminosities, probably because motions are no longer dominated by rotation.
Recently, deep imaging in the near-infrared and a breakthrough in the
characterization of light profiles have led to the recognition of a fundamental
plane (FP) for dwarf irregular galaxies (dIs). Vaduvescu et al. (2005
- V05 from now on) found that dI surface brightness profiles can be modeled
using a sech function. They found a strong correlation among the sech absolute
magnitude measured in
(
), the central surface brightness of the
sech function (
), and the neutral hydrogen line width (W20).
Furthermore, Vaduvescu et al. (2006) discovered that light
profiles for blue compact dwarf galaxies (BCDs) can be modeled as a Gaussian
burst on top of a sech profile, making ready separation of young and old
components possible. Thus, the FP for dIs provides an avenue to establish with
greater clarity how closely BCDs (and other classes of dwarfs) may be related
to dIs.
In this paper, we continue our study of the FP for dIs by expanding the sample of objects with good near-infrared photometry. We also make use of improved distances published recently in the literature. To complete the dwarf picture, we also address in the FP context two other classes of dwarf galaxies, namely BCDs and dwarf ellipticals (dEs).
In Sect. 2 we list criteria for selecting the sample of dwarfs presented in this paper. Our previous observations which defined the FP initially are reviewed in Sect. 3. Data from 2MASS of value to this study is discussed in Sect. 4. New observations acquired at CFHT and CTIO are presented in Sect. 5. Homogenized distances for our galaxies are derived in Sect. 6. The FP for dIs is examined in Sect. 7, and its relation to other dwarf classes is discussed in Sect. 8. In Sect. 9 we present conclusions.
To investigate the FP for dIs, good distances are required. Consequently, our sample of dIs is restricted to objects in the Local Volume. The Local Volume (LV) is defined to encompass all galaxies within about 10 Mpc on a scale where the distance modulus of the LMC is 18.50 (Karachentsev et al. 2004). Cepheids (CEPH) and/or the tip of the red giant branch (TRGB) have been detected for many dIs in the LV, making possible the determination of accurate distances. For a few objects lacking such data, though, the Hubble law (HUB) or the three brightest stars (3BS) are employed to derive distances.
Further guiding the selection of dIs was size. To maximize observing efficiency, we emphasized targets which would fit into one quarter of the field of NIR detectors available at CFHT or CTIO. This made it possible to measure the sky background without moving the target off the infrared array. Data for the largest objects were extracted from 2MASS.
In this paper, we enlarge the sample of 19 dIs presented in V05 by 25. Of those, 8 objects were observed by 2MASS, 8 were observed at CFHT, and 9 were observed at CTIO.
In a previous paper (V05), we conducted an extensive examination of the
surface brightness profiles of dIs in the NIR in an effort to identify a
unifying mathematical description. Many functions were assessed, including
exponential, de Vaucouleurs, and Sersic profiles, as well as ``mixed'' and
``truncated'' profiles. We concluded that a simple hyperbolic secant function
(sech) fits the surface brightness profiles of the majority of dIs to a high
degree of accuracy (mostly within 0.1 mag).
According to this law, the fitted flux I at the distance r from the centre
along the semi-major axis can be expressed as a function of only two parameters,
namely the central surface brightness, I0, and scale length, r0:
Values of r0, I0,
,
and
for our original sample are
provided in V05. For the 2MASS, CFHT, and CTIO samples presented here, they
are given in Table 2.
Between 1997 and 2001, the Two Micron All Sky Survey (2MASS) imaged virtually
the entire sky in the NIR (J, H, and ). The survey was taken using
relatively small exposures (7.8 s co-added), employing two 1.3 m telescopes
located at Mount Hopkins Observatory in Arizona and CTIO in Chile (Skrutskie
2006). Although most dwarfs failed to be imaged by 2MASS due to their
faint surface brightness, a few bright dIs were detected. We searched the 2MASS
Atlas in order to extract all field dwarfs imaged with a S/N at
sufficiently good to make possible solid fits to the light profiles. Eight dIs
were selected for this paper. For the calibration, we employed the zero points
included in the original image headers.
Following the approach in V05, we removed first all resolved sources from the original 2MASS image. Usually, these include Galactic stars in the foreground, young stars resolved in the host dwarf, and some extragalactic sources in the background. To remove all these sources, we used KILLALL, an IRAF script built around the DAOPHOT package (Buta & McCall 1999). The output is a ``cleaned'' image showing only the old unresolved component of the galaxy and a resolved sources catalog. Surface photometry was performed on the cleaned image using the task ELLIPSE in STSDAS. Then a sech law was fitted to the unresolved component (using the task NFIT1D).
Each of the dIs studied is displayed in Fig. 1. Panels on the left
show the original 2MASS
images, while the right panels display the surface
photometry and the sech fits. Fitting parameters and magnitudes are given in
the first section of Table 2.
It is difficult to observe faint extended sources in the NIR due to the fact
that the sky is very bright and variable in a short time-frame. An additional
problem is that the background due to the detector and thermal radiation may
vary both temporally and spatially (Vaduvescu & McCall 2004). Given
these factors, we adopted a prudent strategy of observing the sky for as long a
time as the targets, shifting between target and sky at a rate fast enough to
resolve temporal variations in the background. This approach allowed us to
observe dwarfs as deep as
in
,
consistent
with our previous work. We present a log of our observations in
Table 1.
![]() |
Figure 1:
Images and surface photometry of dIs observed by 2MASS.
Left panels: ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 1: continued. |
Open with DEXTER |
Table 1: The observing log for our CFHT 2004 and CTIO 2007 observations.
On Mar. 5-8, 2004 we acquired NIR imaging at the f/8 Cassegrain focus of the 3.6 m Canada-France-Hawaii Telescope (CFHT) located atop Mauna Kea, Hawaii (Run ID: 2004-AC23). Due to a technical problem at CFHT, one additional night (Feb. 23/24, 2004) was allocated to the same proposal, being observed by the CFHT staff in service mode. From four nights of the run, data were acquired during three, one being clouded out (Mar. 5/6, 2007). During the last two clear nights, humidity was variable and going high at times. The Moon was within five days of full.
We employed the CFHT-IR camera equipped with a HgCdTe array with
pixels and the
filter. The scale was
,
yielding a
FOV.
Data were reduced in the manner described by V05, with calibrations being
determined using 2MASS stars in the fields of view of the targets (typically 5,
but as low as 1 or 2 for UGC 5764 and UGC 7408). Surface
brightness profiles and fits to them were derived as described above for the
2MASS data.
Imaging and photometry of eight targets are presented here. Another three galaxies (UGC 4483, UGC 5423, and Sextans A) were presented in V05. In the left panels of Fig. 2, we show the reduced CFHT images, while in the right panels we display the surface photometry and the sech fits. In the second section of Table 2, we list the fitting parameters and magnitudes.
On Aug. 1-3, 2007, we acquired NIR imaging using the 4 m Blanco telescope
located at Cerro Tololo Interamerican Observatory, Chile (Run ID: 2007B-0907).
We employed the ISPI camera equipped with a Hawaii array with
pixels mounted at the f/8 Cassegrain focus. The scale was
,
yielding a
FOV. For all targets, we
used only the
filter. The first night was clouded out, while the second
was clear, the sky being bright (two days from Full Moon).
Data were reduced, calibrated, and analysed in the same manner as described
above for CFHT. Typically, 10 to 15 2MASS stars were employed to calibrate each
field.
Imaging and photometry for nine targets are presented here. In the left panels of Fig. 3, we show the reduced CTIO images, while in the right panel we display the surface photometry and the sech fits. In the third section of Table 2, we give the fitting parameters and magnitudes.
Most distances adopted in this paper are based upon either Cepheid variables or the TRGB. For some galaxies, it has been necessary to use the brightest blue stars or the Hubble law to derive distances.
Following V05, all distances have been homogenized to a common distance scale
anchored to the maser distance of NGC 4258 (Herrnstein et al. 1999).
Additionally, care was taken to correct all distance indices for extinction
using an algorithm which appropriately accommodates effective wavelength shifts
(McCall 2004).
HST Key Project distances to Cepheid calibrators were appropriately revised
(see Fingerhut et al. 2007), and those for which the TRGB was observed
were used to calibrate absolute magnitudes. In the I band, the TRGB has been
found to be nearly invariant over a wide range of metallicities, thereby serving
as a powerful standard candle. Our calibration yields
(on the maser zero-point). Fifteen new dwarfs presented in
this paper have distances determined from the tip of the red giant branch.
To calculate Hubble distances (for most CFHT targets), we employed velocities
listed in NED in the Virgo and Great
Attractor flow model (following Mould et al. 2000).
Given the choice of the maser zero-point, we adopted
(Fingerhut et al. 2007).
All Galactic extinctions presented in Table 2 ()
originate from reddenings of Schlegel et al. (1998), E(B-V), calculated
by NED. The reddenings were converted to broadband extinctions using the York
Extinction Solver (YES), the foundations of which are described by McCall
(2004)
.
![]() |
Figure 2:
Images and surface photometry of dIs observed at CFHT.
Left panels: ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 2: continued. |
Open with DEXTER |
Table 2: Photometric parameters for the 2MASS sample (first group), CFHT 2004 sample (second group), and CTIO 2007 sample (third group).
Given the reddening for a target, YES computes the corresponding optical depth
at 1
(
)
using the spectral energy distribution for an elliptical
galaxy (the type of galaxy to which the reddenings of Schlegel et al.
(1998). YES then converts this value for a particular target into
broadband extinctions in desired filters using an appropriately redshifted and
extinguished spectral energy distribution for the target (in our case, Im). By
employing YES to calculate broadband extinctions, we avoid source-dependent
shifts in the effective wavelengths of broadband filters.
In a previous paper (V05), we discovered the ``fundamental plane'' (FP) of dwarf
irregular galaxies (dIs) using photometric data for 19 objects observed in the
field and having radio data from the literature. Most targets were selected to
have accurate distances derived from Cepheids and the TRGB. We defined the plane
in ,
expressing the absolute sech magnitude (
)
derived from the sech
model (i.e.,
)
as a function of the central surface brightness for the
sech model (
)
and the 20% width of the HI line profile (W20).
In this paper, we augment our study of the FP by adding the 25 new dwarfs presented in Sects. 4 and 5, which increases the size of our sample to 44 objects. From this sample, we select 34 dIs to refine the FP equation.
For most brighter dwarfs (V05 and CFHT samples), we employ line widths W20from the RC3 catalog (de Vaucouleurs 1991), while for fainter targets (CTIO 2007 sample) we take W20 from Huchtmeier et al. (2000,2001) and HIPASS (Koribalski et al. 2004). In four cases (ESO 320-14, ESO 321-14, ESO 325-11, and ESO 349-31) no W20 is available, so we calculate W20 using the scaling factor W20/W50=1.45, taking W50 from Karachentsev et al. (2004).
In Fig. 4, we plot the Tully-Fisher relation for all dwarfs presented in this paper. This figure definitively shows that the TF relation breaks down for dwarf galaxies. In Fig. 5, we plot the dwarf fundamental plane, coefficients for which have been defined by the 34 dIs for which data are most reliable (see below). The galaxies from V05 are marked by filled circles. 2MASS galaxies are shown as filled triangles, CFHT 2004 galaxies by crosses, and the CTIO 2007 galaxies by empty circles. Symbols for the objects with the most accurate distances (TRGB, CEPH) are encompassed by squares.
Of the 25 new objects, 10 deviate by more than
(about 1 mag) with
respect to the 34 galaxies defining the FP. NGC 3077 and Holmberg II
(from the 2MASS sample) have accurate distances, but they are known to have tidal
interactions in the M 81-M 82-NGC 3077 system (e.g., Walter
et al. 2006). UGC 7408 and UGC 7490 (from the CFHT sample) are
the most distant dwarfs observed, being 8 and 9.2 Mpc away, respectively (Hunter & Elmegreen
2004). Distances had to be derived from the Hubble law. Moreover, both
galaxies appear to have light profiles which look more exponential than sech,
UGC 7490 being listed by NED as a spiral (SAb). Notably, too, very few 2MASS
stars were available to calibrate the UGC 7408 data. Profiles of
ESO 223-09
and ESO 245-05 (from the CTIO sample) also appear exponential in shape. At the
bright end of the luminosity scale, NGC 2915 and NGC 4214 (from the 2MASS
sample) have accurate distances, but they appear deviant, possibly reflecting the
limitations of 2MASS photometry. At the faint end, the deviants are
ESO 320-14
and ESO 321-014 (from the CTIO sample). Deeper imaging will be necessary to pinpoint
why.
We derived the FP using the ``best'' 34 dIs (19 from the first sample, plus
15 new objects with deviations less than ). These 34 are what remained
after iteratively rejecting the most extreme deviants until a point was reached
where further rejections had little impact on the fit. The improved equation of
the FP is:
Uncertainties in central surface brightnesses are dominated by calibration errors.
Typically, we estimate the uncertainties to be under 0.1 mag. Uncertainties in the
line widths are typically 10%. As can be seen from Eq. (2), the line
width errors dominate uncertainties in the abscissae of Fig. 5. We estimate
a typical error of 0.2 mag in total, which is displayed in the upper left corner of
Fig. 5. Uncertainties in the absolute magnitudes are dominated by
errors in distances. For the TRGB method, the average uncertainties amount to
about 0.18 mag. HUB distances are uncertain by at least 12% based upon the
uncertainty in the Hubble constant (V05), but of course actual uncertainties
must be considered to be higher for such nearby objects due to peculiar motions.
Apparent sech magnitudes are coupled to central surface brightnesses, so their
uncertainties are typically around 0.1 mag. Overall, we estimate that
uncertainty in
is typically 0.2 mag, which is displayed in the upper
left corner of Fig. 5.
Considering the uncertainties in both abscissae and ordinates, it is reasonable to expect observational errors to lead to vertical scatter about the fundamental plane amounting to about 0.3 mag. This is less than the observed rms of 0.43 mag, indicating that there remain unexplained contributions to the noise.
![]() |
Figure 3:
Images and surface photometry of dIs observed at CTIO.
Left panels: ![]() ![]() ![]() |
Open with DEXTER |
![]() |
Figure 3: continued. |
Open with DEXTER |
![]() |
Figure 3: continued. |
Open with DEXTER |
To explore whether varying inclinations might be a source of scatter, we examined
how the FP changes when line widths or central surface brightnesses are corrected
for tilt. First, we estimated the inclination i for each galaxy from the ellipticity
listed in Table 2 using the prescription of
Staveley-Smith et al. (1992) for dwarfs. Then, we made the
assumption that motions are predominantly rotational, and estimated the
rotational velocity from
.
Using
instead
of W20 as the gauge of kinematics did not yield any improvement to the fit
to the FP. The result indicates that internal motions in the majority of the dwarfs
are predominantly random. This is a major advantage of the dI FP over the
TF relation
for spirals, in the sense that distances can be determined without
correcting W20 for tilt.
How
is corrected for tilt requires an assumption about the intrinsic
shape of the galaxies. Both prolate and oblate geometries were considered.
We computed the surface brightness at the orientation yielding circular isophotes
from
,
with the plus applicable to
the prolate case and the minus to the oblate. Neither approach to correcting surface
brightnesses led to any improvement in the fit to the FP.
In summary, using an enlarged sample of galaxies, we have confirmed the existence of the dI FP. However, there remains scatter which is unexplained by formal observational errors, although some of it is likely due to photometry which is too shallow. Also, a few galaxies appear to be distorted by tidal interactions and a few may be massive enough that rotational motions predominate over random ones (in which case, deviance can arise from not correcting the line width for tilt). With the aid of additional data, we plan to examine outliers in more detail in a future paper.
![]() |
Figure 3:
continued. Top panel: WLM image mosaiked from North and South frames
(North is up, East to the left; field of view
![]() ![]() |
Open with DEXTER |
![]() |
Figure 4: The Tully-Fisher relation for dIs. No corrections for inclination have been applied. The vertical scale is the same as for the fundamental plane plotted in the next figure. |
Open with DEXTER |
![]() |
Figure 5:
The fundamental plane for dIs. No corrections for inclination have
been applied. Different symbols are used to distinguish the different
sources of observations. Objects with accurate distances (TRGB, CEPH) are
emphasized with squares. Ten labeled objects lie off the plane (deviating by
more than ![]() |
Open with DEXTER |
![]() |
Figure 6: dIs, BCDs, and dEs on the dI fundamental plane. dIs were fitted using a sech law, BCDs using a sech plus Gaussian, and dEs using an exponential. BCDs and dEs lie on the FP defined by the dIs, suggesting deep physical connections between the three dwarf classes. |
Open with DEXTER |
The dI FP provides a baseline with which to assess linkages with other kinds of dwarfs, such as Blue Compact Dwarfs (BCDs) and dwarf Ellipticals (dEs). Preliminary work by Vaduvescu (2005) suggested that BCDs and dEs lie close to the dI FP. Here, we examine relationships further using our improved FP determination.
In a previous paper (Vaduvescu, Richer & McCall 2006) we presented deep NIR data for 16 BCDs located in the Virgo cluster. We showed that surface brightness profiles of BCDs can be fitted accurately (within 0.1 mag) by the sum of two components: a hyperbolic secant (sech) to model the overall profile, and a Gaussian component to trace the starburst near the centre of the galaxy. By comparing parameters of the sech fits, it was shown that BCDs appear to be physically similar to dIs. Also, their chemical properties track those of dIs (Vaduvescu et al. 2007).
To place our BCDs on the dI FP, we adopt a distance modulus of 30.62 for the Virgo Cluster (Freedman et al. 2001), which is the HST Key Project value but anchored to the maser distance for NGC 4258. For magnitudes and surface brightnesses, we adopt those of the sech component.
Figure 6 includes 16 BCDs from Virgo and 3 BCDs in the field plotted with respect to the FP defined by Eq. (2) (represented by the dotted line). BCDs are plotted as triangles. Most BCDs lie close to the FP.
Unlike dIs and BCDs, dwarf elliptical galaxies (dEs) and dwarf spheroidals (dSphs) are known to harbour old stellar populations with no traces of gas or star forming regions. A linkage between dIs and dEs was signaled originally through pair-wise comparisons of structural and kinematic parameters (Kormendy 1985). Also, dEs and dIs proved to be more closely related than dEs and giant ellipticals . However, Papaderos et al. (1996) did not find a clear evolutionary connection between dEs, dIs, and BCDs.
The dI FP offers a more refined avenue for exploring how dIs and dEs may be related. A first attempt was presented by Vaduvescu & McCall (2005) and Vaduvescu (2005) who employed 22 dEs located in the Virgo Cluster and the Local Group. The same sample is used here. Data were selected as follows:
Surface brightness profiles of dEs are peaked at small radii (i.e., a convex function), so cannot be modeled by a sech function. Combining a sech function with a Gaussian, as we did for BCDs, also does not work because of the sharpness of the cores. Instead, we adopted an exponential model to fit surface brightness profiles of dEs. An exponential was shown to fit profiles of dEs by Kormendy (1977), Faber & Lin (1983) and Binggeli et al. (1984). The exponential law has two parameters, the central surface brightness and the scale length, i.e., the same degrees of freedom as the sech law for dIs and BCDs. To fit the profiles, we restricted the exponential to the linear part of the fit, preventing the core which in most cases was a small fraction from the entire profile (around 10%). The absolute magnitude was computed from the apparent magnitude of the exponential model.
Figure 6 shows the dEs plotted with respect to the FP defined by Eq. (2) (the dashed line). Despite the approximations made to extract photometric parameters, the dEs fall on the FP defined by the dIs over a 9 mag interval in absolute magnitude.
New NIR photometry for dIs has been presented, leading to an expanded sample of
44 galaxies with which to study the fundamental plane discovered by Vaduvescu et al.
(2005). As before, it has been appropriate to define photometric parameters
by fitting a sech function to the surface brightness profiles. The FP remains clearly
defined, with the rms
absolute magnitudes for the 34 galaxies used in the fit
being 0.43 mag (as in the original paper). Correcting line widths for inclination
did not lead to any reduction in the rms, indicating that internal motions are
predominantly random. Correction of surface brightnesses for tilt led to no
improvement, too, regardless as to whether geometries were assumed to be prolate
or oblate. The results suggest that the FP can be used to derive dwarf distances
without knowing tilts.
The new fit to the dI FP has been used to evaluate linkages to BCDs and dEs. Light profiles of BCDs were fit as a combination of a Gaussian burst on top of an underlying distribution defined by a sech function. Adopted photometric parameters were derived from the shape of the sech function. With photometric parameters defined this way, BCDs fall precisely on the FP defined by dIs. dEs light profiles cannot be modeled using a sech function. Instead, we fitted their overall envelopes with an exponential, and used the parameters ot the exponential to define photometric parameters for the FP. When parametrized this way, dEs also lie on the FP defined by dIs. The outstanding overlap of the three dwarf classes suggests strong linkages among dIs, BCDs, and dEs. The simplicity of light profiles and chemistry of dIs motivates using them as a control for the evolutionary studies of other dwarf galaxies.
Acknowledgements
This publication makes use of data products from the Two Micron All Sky Survey (2MASS), which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the National Science Foundation. The paper included data based on observations obtained at the Canada-France-Hawaii Telescope (CFHT) which is operated by the National Research Council of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique of France, and the University of Hawaii. We acquired some data at Cerro Tololo Inter-American Observatory (CTIO), which is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under a cooperative agreement with the National Science Foundation. This research has made use of the NASA/IPAC Extragalactic Database (NED) which is operated by the Jet Propulsion Laboratory, California Institute of Technology, under contract with the National Aeronautics and Space Administration. This research has made use of the GOLDMine Database (Gavazzi et al. 2002). For our data reductions, we used IRAF, distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. M.L.M. thanks the Natural Sciences and Engineering Council of Canada for its continuing support. Special thanks are due to the anonymous referee who provided valuable feedback which allowed us to improve the paper.