Issue |
A&A
Volume 562, February 2014
|
|
---|---|---|
Article Number | A49 | |
Number of page(s) | 22 | |
Section | Extragalactic astronomy | |
DOI | https://doi.org/10.1051/0004-6361/201220700 | |
Published online | 06 February 2014 |
What will blue compact dwarf galaxies evolve into?⋆
1
Astronomisches Rechen-Institut, Zentrum für Astronomie der Universität
Heidelberg,
Mönchhofstraße 12-14,
69120
Heidelberg,
Germany
e-mail:
HTM@x-astro.net,TL@x-astro.net
2
Division of Astronomy, Department of Physics,
PO Box 3000, 90014 University of Oulu, Finland
3
Centro de Astrofísica and Faculdade de Ciências, Universidade do
Porto, rua das
Estrelas, 4150-762
Porto,
Portugal
Received:
6
November
2012
Accepted:
30
October
2013
We present and analyse the photometric properties of a nearly complete sample of blue compact dwarf (BCD) and irregular galaxies in the Virgo cluster from multi-band SDSS images. Our study intends to shed light on the ongoing debate of whether a structural evolution from present-day star-forming dwarf galaxies in a cluster environment into ordinary early-type dwarf galaxies is possible based on the structural properties. For this purpose, we decompose the surface brightness profiles of the BCDs into the luminosity contribution of the starburst component and that of their underlying low surface brightness (LSB) host. The latter dominates the stellar mass of the BCD. We find that the LSB-components of the Virgo BCDs are structurally compatible with the more compact half of the Virgo early-type dwarfs, except for a few extreme BCDs. Thus, after termination of starburst activity, the BCDs will presumably fade into galaxies that are structurally similar to ordinary early-type dwarfs. In contrast, the irregulars are more diffuse than the BCDs and are structurally similar to the more diffuse half of the Virgo early-type dwarfs. Therefore, the present-day Virgo irregulars are not simply non-starbursting BCDs. If starbursts in cluster BCDs are transient phenomena with a duration of ~100 Myr or less, during which the galaxies could not travel more than ~100 kpc, then a substantial number of non-starbursting counterparts of these systems must populate the same spatial volume, namely the Virgo cluster outskirts. The majority of them would have to be early-type dwarfs, based on the abundance of different galaxy types with similar colours and structural parameters to the LSB-components of the BCDs. However, most Virgo BCDs have redder LSB-host colours and a less prominent starburst than typical field BCDs, preventing a robust conclusion on possible oscillations between BCDs and early-type dwarfs.
Key words: galaxies: dwarf / galaxies: evolution / galaxies: clusters: individual: Virgo
Appendix is available in electronic form at http://www.aanda.org
© ESO, 2014
1. Introduction
Galaxy clusters like the Virgo cluster are characterised by a wide variety in galaxy morphology. This morphological variety depends on the local galaxy density, and therefore shows a clear trend with the clustercentric distance, as was first revealed by Dressler (1980) and confirmed by many subsequent studies (e.g. Binggeli et al. 1990; Jerjen 2012). At short clustercentric distances the dominant galaxy types are early-type galaxies (ETGs): elliptical (E) and lenticular (S0), dwarf elliptical (dE), dwarf lenticular (dS0), and dwarf spheroidal (dSph) galaxies, spanning a range from high to low luminosities. Traditionally, ETGs are associated with a smooth, regular appearance and no signs of star formation. However, the structural characteristics of early-type dwarfs have recently been shown to be unexpectedly complex (Janz et al. 2012, 2013). Furthermore, they comprise a particular subclass of galaxies with relatively blue cores, which indicate recent or still ongoing low-level star formation at very low star formation rates (e.g. Lisker et al. 2006b).
Studies of ETGs have shown that they have the oldest luminosity-weighted stellar ages among galaxies in clusters, resulting in red colours. On the other hand, dynamically young galaxy clusters (e.g. the Virgo galaxy cluster) also contain late-type galaxies, which in contrast to the ETGs, are located at greater clustercentric distances. At low luminosities, these comprise star-forming dwarf galaxies such as blue compact dwarfs (BCDs) and irregulars. As shown by Vilchez (1995), the Hβ equivalent width (EW) of late-type galaxies in the Virgo cluster increases with increasing distance from the center, echoeing the strong impact of the cluster environment on galaxy evolution. For this and other reasons, it has become a common hypothesis that late-type galaxies may be environmentally transformed to early-type dwarf galaxies (e.g. Kormendy 1985; Boselli & Gavazzi 2006).
Tempting questions in this respect are i) how could the morphology-density relation be explained; ii) which mechanisms are responsible for the gradual transformation of late-type galaxies into ETGs as one approaches the cluster center; and iii) are the descendants of today’s low-mass late-type galaxies structurally compatible with early-type dwarfs in the Virgo cluster? In the literature several mechanisms have been proposed as drivers of galaxy transformations within the dense cluster environment, most notably ram pressure stripping (Gunn & Gott 1972), tidal stirring (Mayer et al. 2001) and harassment (Moore et al. 1996; Mastropietro et al. 2005), all of which have the removal of gas in common (see e.g. Hensler 2012, for a recent review).
Dwarf galaxies with a low central stellar density ρ⋆, such as irregulars, are expected to be particularly prone to gas removal when they plunge into the hot intracluster medium (ICM), whereas high-ρ⋆ systems, such as BCDs (e.g. Papaderos et al. 1996b), might be able to retain some fraction of their gaseous reservoir down to a shorter clustercentric distance. But what would these dwarfish late-type galaxies in the cluster periphery look like after some billion years, once their gas has been removed, in the course of one or several passages through the dense ICM core and their ensuing long passive photometric evolution? Addressing this question is fundamental to understanding the morphological diversity of the dwarf galaxy population in clusters. Another issue of special interest is what impact the initial contact of late-type dwarfs with the cluster periphery has on their star formation activity and whether, for certain conditions, starbursts can be ignited, transforming them into BCDs. If so, then how do these cluster-BCDs differ from the main population of field BCDs in their recent star formation history and the morphological properties of their star-forming component?
BCDs are low-luminosity galaxies (MB > −18mag) with a compact optical appearance and blue integral colours (Thuan & Martin 1981). Many studies over the past decades have shown that BCDs are metal deficient with a median oxygen abundance of 12 + log (O/H) ~ 8.0 (see e.g. Kunth & Östlin 2000, for a compilation of literature data) and a low percentage of systems with a gas-phase metallicity as low as 7.0 ≲ 12 + log (O/H) ≲ 7.6 (Searle & Sargent 1972; Izotov et al. 1999; Kunth & Östlin 2000; Kniazev et al. 2004; Papaderos et al. 2008). Various lines of evidence consistently suggest that these systems (XBCDs1) have formed most of their stellar mass at late cosmic epochs. They are therefore the best nearby analogues of young low-mass galaxies in the early Universe (see e.g. the discussion in Papaderos et al. 2008). Apart from their low metallicity, BCDs also exhibit strong bursts of star formation, which are fed by a relatively large amount of gas (Thuan & Martin 1981; Staveley-Smith et al. 1992; van Zee et al. 1998). Various arguments suggest, in line with evolutionary synthesis models, that starbursts in typical field BCDs do not last longer than ~107 yr (Thuan 1991; Krüger et al. 1995; Mas-Hesse & Kunth 1999; Thornley et al. 2000) and have to be separated by long (~1 Gyr) quiescent phases. Dwarf irregulars, on the other hand, are characterised by prolonged low-level star formation that lasts 450 Myr up to 1.3 Gyr, as shown by, e.g., McQuinn et al. (2010) through colour−magnitude studies of 20 of these systems.
The seminal study by Loose & Thuan (1986b) has shown that BCDs are composed of two main stellar components. The first one comprises the region of the ongoing star formation where young OB-stars dominate, and which therefore has a low stellar mass-to-light ratio (M/L). Due to the simultaneous formation of new massive stars in the starburst and their death within a few Myr, metal-enriched galactic winds driven by multiple supernova explosions are expected to influence the chemical evolution of BCDs. However, studies have shown that in most cases starburst-driven feedback is insufficient for the expulsion of the entire ISM from a BCD (e.g. Silich & Tenorio-Tagle 1998; Ferrara & Tolstoy 2000; Tajiri & Kamaya 2002; Recchi & Hensler 2006).
The starburst component contributes, on average, ~50% of the optical emission of a BCD and, in some cases, up to 90% (Papaderos et al. 1996b; Noeske 1999; Cairós et al. 2001b; Amorín et al. 2009). Quite importantly, in several XBCDs nebular emission has been determined to be extraordinarily intense (EW ~ 1600−2000 Å) and to contribute 30–50% of their total optical luminosity (Papaderos et al. 1998, 2002), in line with theoretical predictions for young starbursts (e.g. Krüger et al. 1995).
The second component of BCDs is dominated by an old population of low-mass stars. This component is referred to in the following as the host galaxy or low-surface brightness (LSB)2 component and is characterised by a high stellar M/L. Various studies indicate that the LSB-component dominates the baryonic-, and in some cases even the virial mass of BCDs within their Holmberg radius3. Therefore, it has to have a significant influence on the gas collapse characteristics and the starburst activity in these systems (e.g. Papaderos et al. 1996a; Östlin et al. 1999; Elmegreen et al. 2012; Lelli et al. 2012; Micheva et al. 2013).
Different subtypes of BCDs according to the morphology of the starburst- and LSB-component based on Loose & Thuan (1986b).
According to Loose & Thuan (1986b,a), BCDs can be classified in four main subtypes, based on the morphology of their star-forming and LSB-component (Table 1; see Salzer et al. 1989; Telles et al. 1997; Cairós et al. 2001a and Sung et al. 2002 for alternative classification schemes.). Nuclear-elliptical (nE) and irregular-elliptical (iE) BCDs are both characterised by an extended circular or elliptical LSB-component, while BCDs with an irregular LSB-component are classified “iI”. A recent study by Zhao et al. (2013) found that iI-type BCDs are on average more gas-rich and metal-poor than nE/iE-type BCDs. Interestingly, the cometary (iI,C) subtype (see e.g. Noeske et al. 2000, for two nearby examples) is remarkably common among XBCDs (Papaderos et al. 2008), a fact pointing to a connection between gas-phase metallicity, morphology and evolutionary status. This is also indirectly suggested by the high frequency of cometary galaxies (also referred to as tadpoles) among comparatively unevolved high-z galaxies (e.g. Elmegreen & Elmegreen 2010).
Cometary morphology in field XBCDs has been proposed to result from unidirectional sequential star formation activity with a typical velocity of the sound speed in the warm ISM (~20 km s-1, Papaderos et al. 1998, 2008). This scenario is supported by stellar age gradients along the “comet’s tail” (e.g. Guseva et al. 2003). In the cluster periphery, however, cometary morphology may arise from extranuclear star formation that is triggered through the interaction with the ICM (see Hensler 2012, for a recent review), as, e.g. the impressive cases of two star-forming dwarf galaxies in Abell 1367 (Gavazzi et al. 2001) demonstrate.
The evolution of star-forming dwarf galaxies like BCDs to “red and dead” early-type galaxies is still under debate and there is no satisfying answer yet. Early studies by Thuan (1985) and Davies & Phillipps (1988) introduced the idea of a possible evolution of irregulars to BCDs in several bursts and finally, after reaching a higher metallicity and the depletion of gas, the fading to early-type dwarfs. Thuan (1985) concluded that the metallicities of BCDs and early-type dwarfs are very different and that only periods of three to ten bursts over the Hubble time could be able to produce the metallicity range of early-type dwarfs. The study of Marquart et al. (2007) on the BCD He 2−10 showed that the stars in this galaxy have random motions and show signatures of a merger. They concluded that, due to the velocity dispersion of the stars, in the future this BCD may evolve into a nucleated early-type dwarf. In comparison to this, the simulations of Bekki (2008) showed that dwarf-dwarf mergers are able to produce BCDs, but the further evolution of these BCDs into gas-free dwarfs would require additional influence to remove the extended gas discs of the simulated BCDs.
In contrast to BCDs, early-type dwarf galaxies typically show almost no evidence for ongoing star formation. However, several studies (e.g. Vigroux et al. 1984; Gu et al. 2006; Lisker et al. 2006b) found early-type dwarfs with blue colours in their central region, indicating recent or ongoing star formation, partly with post-starburst spectra that may hint at a population of recent arrivals to the cluster environment (Gavazzi et al. 2010). Among 476 early-type dwarfs in the Virgo cluster, Lisker et al. (2006b) identified 23 galaxies with blue centers and classified them as “dE(bc)”. Since the debate of the evolutionary connection of early-type dwarfs and BCDs is still ongoing, it is worth comparing the structural properties of BCDs to those of early-type dwarfs and in particular to those of dE(bc)s. Can BCDs in the periphery of the Virgo cluster evolve into objects similar to the present-day early-type dwarfs? Would this be different for field BCDs?
This paper is structured as follows: Sect. 2 defines our sample of cluster BCDs and the classification criteria used. In Sect. 3 we outline our photometric measurements, with particular emphasis on the decomposition methodology for BCDs. Our results on the structural properties of BCDs and the comparison to early-type dwarfs and irregulars are presented in Sect. 4. In Sect. 5, we discuss the results and conclude with an assessment of possible evolutionary connections between the different dwarf types.
2. Sample selection and SDSS data
Our sample is based on the Virgo Cluster Catalog (VCC) by Binggeli et al. (1985), which includes galaxies of all types within the Virgo cluster area. Due to incomplete velocity information for the VCC galaxies, the VCC includes certain and possible cluster members (updated by Binggeli et al. 1993), apart from background galaxies. Since new velocities have become available in the meantime, largely due to the Sloan Digital Sky Survey (SDSS, Adelman-McCarthy et al. 2007), the membership was revised by one of us (T. Lisker, see appendix of Weinmann et al. 2011) using the NASA/IPAC Extragalactic Database4. If a galaxy is listed as certain or possible member in the VCC, but has a velocity above 3500 km s-1, it is considered as a background galaxy. We do not change possible to certain members or vice versa.
Here we consider both certain and possible members to a magnitude limit of mB ≤ 18.0 mag, to which the VCC was found to be complete (Binggeli et al. 1985). When applying a constant distance modulus of m − M = 31.09 mag (d = 16.5 Mpc; Mei et al. 2007) to all galaxies, this corresponds to an absolute magnitude limit of MB < −13.09 mag.
VCC galaxies classified as BCD or candidate BCD.
Within these limits, the VCC contains 57 galaxies which have the term “BCD” included in their morphological type, 38 of them with BCD as the only or primary class (Table 2). Among the galaxies whose class was “?” and “pec” in the VCC, a visual inspection of the SDSS images5 suggests a reclassification of VCC 0429 and VCC 1713 to the mixed type “Im / BCD”, and of VCC 1411 to a possible “BCD?”, thus adding up to 60 galaxies, and 39 with BCD as primary class (Table 2). Since the VCC photographic plates were more sensitive in the blue, it needs to be kept in mind that this could have partly influenced the classification, as any underlying red stellar population would appear less prominent in the blue. Therefore, the galaxies with primary class BCD are chosen as the working sample for our analysis, whereas the uncertain BCD candidates are treated as a separate sample that is compared to the BCDs in Sect. 4.4. Figure 1 shows the distribution of both samples within the Virgo cluster.
Fig. 1 Distribution of Virgo cluster galaxies. Red symbols refer to objects with BCD as primary class, black symbols to uncertain BCD candidates, small black dots to early-type galaxies, green asterisks to dE(bc)s, and the black cross to M87. Open symbols denote galaxies that could not be analysed (see text for details). |
Five of the 39 BCDs (VCC 0464, 0468, 0741, 0772, and 0985) were not covered by the SDSS Data Release Five (DR5, Adelman-McCarthy et al. 2007), two were excluded because of a nearby other galaxy (VCC 0340 and 1944), and two were excluded because of a nearby bright star (VCC 1750 and 2033). Our primary working sample thus consists of 30 galaxies. Two of the 21 uncertain BCD candidates (VCC 0737 and 0848) were not covered by the SDSS DR5, leaving us with 19 such objects.
Our photometric study relies on the SDSS DR5 images in u, g, r, i and z, each with an exposure time of 54 s. Due to insufficient sky subtraction of the SDSS pipeline for nearby galaxies of large apparent size (Lisker et al. 2007), we used the sky-subtracted images provided by Lisker et al. (2007). All images were flux calibrated and corrected for Galactic extinction following Schlegel et al. (1998).
3. Surface photometry and BCD decomposition
3.1. Global photometry of late-type and early-type galaxies
The analysis of the sample of ETGs (classes E, S0, dS0, and dE, including dE(bc) ) was performed by Lisker et al. (2007) and Janz & Lisker (2008). In order to obtain global photometric charateristics for the irregulars and for the uncertain BCD candidates, we closely followed the procedure of Lisker et al. (2007), as outlined below, but using the co-added g, r, and i images (Lisker et al. 2006a) to benefit from a higher signal-to-noise ratio (S/N). While we also applied this procedure to our primary BCD working sample – yielding global characteristics that include the light from the starburst region – we then proceeded to perform a thorough decomposition of their light profile, as outlined in the following sections. We point out that we did not consider it appropriate to apply a similar decomposition to the irregulars, since these are not clear multi-component objects like the BCDs.
We first determined the circular Petrosian radius (Petrosian 1976) on the co-added image, using a visual estimate of the galaxy center6 that aims to reflect the overall light distribution, not only the inner light peak(s) that an irregular may exhibit. We then computed an elliptical shape, i.e. position angle and axis ratio, from the image moments of the enclosed light. This is used to determine a “Petrosian semi-major axis” aP on elliptical apertures. At aP, we let the IRAF task ellipse (Jedrzejewski 1987) determine the isophotal position angle, axis ratio, and centre. With these values, we performed another iteration on deriving aP as well as position angle and axis ratio b/a. In this process, contaminating objects that do not belong to the galaxy were masked.
We obtained total r magnitudes from the flux within 2aP on the r-band image, which also yielded the values of the half-light semi-major axis ahl and the average surface brightness ⟨μ⟩eff within ahl. In our figures we show the effective radius, (1)Colours were obtained from the flux within 2aP on the images of the different bands. Photometric errors were derived using Eqs. (3) and (4) in Lisker et al. (2008), based on the signal-to-noise ratio and on estimations of the uncertainties in the sky subtraction and determination of Petrosian radii.
3.2. BCD surface brightness profiles
The luminosity contribution of the starburst in field BCDs amounts, on average, to ~50% of the B band emission, and in some extreme cases reaches up to ~90% (Papaderos et al. 1996b; Noeske 1999; Salzer & Norten 1999; Amorín et al. 2009). On the other hand, the starburst component is almost negligible, in terms of its fraction by mass, as its stellar M/L is several times lower than that of the LSB-component. An adequate separation of its emission via 1D or 2D decomposition is clearly necessary for isolating the emission and studying the structural properties of the host galaxy. Indeed, the composite surface brightness profile of a BCD holds little insight into the photometric structure prior to decomposition. For example, Papaderos et al. (1996b) have pointed out that the profile of a nE-type BCD can closely resemble the steep inner profile of a giant elliptical, whereas a typical feature of iE-type BCDs is an extended plateau. Similarly, the effective radius Reff can vary by up to a factor ~3 depending on the starburst luminosity fraction (cf. Papaderos et al. 2006).
For galaxies with a smooth appearance like elliptical galaxies the easiest way to obtain surface brightness profiles is to use elliptical apertures and sum up the enclosed flux. Due to the irregular morphology of BCDs, however, such an approach is impractical, on the one hand because the choice of the “galaxy centre” is subjective and wavelength-dependent, and, on the other hand, because profiles derived in this way show, in the case of iE- and iI-BCDs, discontinuous jumps. In some cases, the latter can significantly affect the profile of the LSB-component, thereby biasing studies of the photometric structure of BCDs. 2D axisymmetric models to BCDs other than those falling into the nE class also yield systematic residuals, unless one carefully masks out and avoids fitting of the star-forming component, as was done in e.g. Amorín et al. (2009).
With these considerations in mind, we therefore applied method iv of Papaderos et al. (2002), which was also used by Noeske et al. (2006) under the name LAZY. This method, which has as input a set of co-aligned multi-band images of the same pixel scale (SDSS: 0.396 arcsec/pix) and point spread function (PSF), does not require a choice for the “centre” of a galaxy. In our study, an average full width at half maximum of the PSF of 4 pix (1.584 arcsec) was applied to smooth the input images with a Gaussian convolution kernel. From these co-aligned images a S/N weighted average image (called reference frame) is created, which is used to calculate a mask for each intensity interval ΔI within the range Imin to Imax. Pixels outside the intensity interval of I − ΔI ≤ F ≤ I are set to zero, while pixels within the interval are given full consideration. To apply LAZY, it was necessary to remove contaminating sources by replacing their flux with the median flux of the immediate environment.
In the next step, the mask of a given object and intensity interval is multiplied by the image in a given band, and the product’s total flux is computed. The corresponding photometric radius R⋆ is calculated as (2)where AI and AI − ΔI are the areas with intensities above I and I − ΔI. In the case of multiple star-forming regions the R⋆ derived in this way corresponds to the sum of their area, a concept which translates into a monotonous increase in radius with decreasing intensity threshold.
Fig. 2 Colour images from SDSS (upper panels), spectra (middle panels), as well as surface brightness and colour profiles (lower panels) of VCC 1744 (left) and VCC 0130 (right). In the case of VCC 1744 a pure exponential law was assumed to fit the LSB-component, while for VCC 0130 an inner profile flattening was assumed (black dashed lines). The horizontal bar denotes the interval in R⋆over which the LSB-component was fitted. The images in this figure and throughout the paper have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. They were obtained from the online Image List Tool of the SDSS (http://skyserver.sdss3.org/dr8/en/tools/chart/list.asp) and were slightly enhanced at low intensity, using the same correction curve for all. |
The errors of the data points in the surface brightness profiles (see e.g. Fig. 2) were calculated by estimating the Poisson photon noise and the variation of the background of each SDSS image, as detailed in Noeske et al. (2003), as well as in Papaderos et al. (1996b) and Cairós et al. (2001b).
3.3. LSB-component: exponential model
To distinguish the starburst from the LSB-component in the galaxies of our primary BCD sample, we use (g − i)-colour profiles. At smaller R⋆ the (g − i)-colour profile is due to the superposition of the contribution of the young population and the old population of the LSB-component, which results in relatively blue values. At larger R⋆ the contribution of the starburst component vanishes and the old LSB-component dominates, which results in relatively constant red colours. The radius beyond which the colour index levels off to a red, nearly constant value is referred to as transition radius Rtr (Papaderos et al. 1996b; Cairós et al. 2003; Noeske et al. 2003). We used this characteristic radius to define the minimum R⋆ for fitting an exponential law of the form (3)to the LSB-component. In units of mag/arcsec2 Eq. (3) reads as (4)where ILSB,0 and μLSB,0 are the central intensity and central surface brightness of the LSB-component, and α is its exponential scale length.
3.4. LSB-component: inner flattening
In some cases the extrapolation of the exponential fit to R⋆ = 0 arcsec exceeds the intensity observed at intermediate and small radii. This implies that the exponential law is not applicable in the central part of the LSB-component and must flatten to a central surface brightness that is fainter than the extrapolated value μLSB,0 of the fit for R⋆ ≥ Rtr. This type of outer exponential profile with a flat core was noticed by Binggeli & Cameron (1991) for early-type dwarfs, who have called them type V. Papaderos et al. (1996b) introduced a modified exponential law to approximate such profiles, motivated by two considerations: firstly, a Sérsic law with n ≲ 0.5 cannot fit type V profiles without producing systematic residuals (see Noeske et al. 2003 for a detailed discussion of this subject). Secondly, as shown by radiation transfer models by Papaderos et al. (1996b), a Sérsic law with n ≲ 0.5 implies an extended “hole” in the intrinsic luminosity density of a spherically symmetric emitter, if radiation isotropy, and a uniform M/L ratio and intrinsic extinction are assumed. As these authors considered the evacuation of the dwarf galaxy centres to be improbable, they introduced an intensity profile for which the intrinsic luminosity density increases monotonously with decreasing radius (for small flattening parameters; see below) and has a finite central value. This modified exponential fitting law (modexp) has the form (5)where (6)with the typical ratio b/q being of the order of three (Papaderos et al. 1999 and Papaderos & Östlin 2012). The parameter b is a measure of the radial extent of the central core in units of the exponential scale length α, and describes the attenuation of the modexp fit, as compared to the central intensity ILSB,0 predicted by the pure exponential fitting law. Due to the poor knowledge of the structure of the LSB-component in its central part it is not clear at the moment which parameter combination of (b, q) satisfactorily describes the original shape. Therefore, the choice of (b, q) is not a straightforward task and has to rely on plausibility considerations (cf. Noeske et al. 2003). In this study we chose (b,q) = (2.40,0.80), by which the total magnitude of the LSB-component becomes 0.48 mag fainter and μ0 is reduced by 1.747mag/arcsec2 as compared to the purely exponential LSB-component. Since this magnitude offset is applied to all filter bands, the integrated colours of the LSB-components are not affected. We apply such a flattening only when μLSB,0 exceeds the total central surface brightness μtot,0 and the spectrum shows clear signs of star formation. This is only the case for three BCDs (VCC 0074, 0130, and 0641), and is exemplified in Fig. 2 for VCC 0130.
3.5. LSB-component: outer tail
Fig. 3 Surface brightness profile of a BCD (VCC 0001). Left panel: the line corresponds to an exponential fit to the LSB-component and the red data points show the outer tail in the r-band. Right panel: same, but for the g, r and i-filters. Obvious are the outer tails of the profiles at low surface brightness levels. The horizontal bars indicate the radius at which the profile starts deviating from a pure exponential fit. |
The inspection of the surface brightness profiles of our BCD sample reveals that in many cases the profile of the LSB-component flattens out at large radii. This is illustrated in Fig. 3 for VCC 0001, and is also clearly visible in Fig. 2 for VCC 1744 for radii R⋆ > 15arcsec. We note that such a flattening has been reported in the outskirts of some BCDs, e.g. II Zw 71 and Mrk 178, where it was found to contribute no more than ~3% of the total luminosity (Papaderos et al. 2001, 2002). For some of our sample galaxies, however, already a quick estimate based on the surface brightness profile reveals a significant contribution (10% for VCC 1744, see Fig. 2). Since this leads to an increase in the derived effective radius, and since we will discuss in Sect. 4 whether BCDs are smaller in size than early-type dwarfs, we cannot neglect this systematic effect on the radius.
In the following, we will call the radius where the observed profile of the LSB-component starts deviating from a pure exponential slope tail radius Rtail (Fig. 3). The exponentially fitted LSB-component was integrated up to Rtail and was then continued with the observed data points of the tail. From this LSB+tail profile – also taking the possible inner flattening into account – we derive the Petrosian radius (Petrosian 1976; Graham & Driver 2005), where the intensity drops to a fifth of the mean intensity within that radius (Blanton et al. 2001; Yasuda et al. 2001).
LAZY, just like all surface photometry codes, is very sensitive to the quality of the sky subtraction and the objective selection and removal of fore/background sources. Especially in the faint extended outskirts of a galaxy, this can be difficult and may introduce additional uncertainties. To overcome this problem the parameters of the LSB-component were derived with LAZY only for R⋆ ≤ Rtail, and were subsequently corrected for the missing flux. The latter was estimated by integrating the SDSS imagesfrom Rtail out to two Petrosian radii, using elliptical apertures (see Sect. 3.1), masking contaminating sources, and correcting for the missing flux of these masked patches.This analysis resulted in a correction value δf for the missing flux, such that the total flux of the LSB-component becomes (7)The colours of the LSB-components were determined from the corrected total flux in the different bands. The effective radius Reff,LSB was derived in the r-band by the integration of the combined LSB profile (inner flattening, exponential section, outer tail) to the radius within which half of the total flux is enclosed. The mean effective surface brightness ⟨μ⟩eff,LSB was calculated as (8)
3.6. LSB-component: error estimation
The errors in our determination of magnitude and effective radius of the LSB component were estimated from a set of artificial BCD images, which we created based on the parameter range of the observational sample. These images consist of LSB-components with various profile shapes and a superposed starburst component with various strengths. They were analysed in the same manner as the Virgo BCDs, allowing us to compare input and output parameters. Details are given in Appendix A.
The analysis of the artificial images reveals statistical errors of 0.10 mag in Mr, 6% in Reff, and 0.16 mag/arcsec2 for ⟨μ⟩eff. For the fainter BCDs of our sample, the radii derived with our method are somewhat smaller (11%) than the true values, and the magnitudes are systematically brighter (0.15 mag). These add up to a systematically brighter effective surface brightness by 0.41 mag/arcsec2 for the fainter BCDs. In the following figures we indicate the statistical errors with representative error bars and the systematic offsets with arrows pointing from the center of these error bars to the true value.
3.7. Decomposition examples
We illustrate our decomposition with two examples in Fig. 2. VCC 1744 (left) shows a very strong star-forming region, which is off-centred. On the other hand, the star-forming regions in VCC 0130 (right) spread over the entire galaxy, resulting in an irregular LSB-component. The SDSS spectrum of VCC 1744 shows very strong emission lines, as is typical for BCDs. VCC 0130 also shows emission lines and a rising blue continuum, but these emission lines are not as strong as in VCC 1744.
The transition radius of VCC 1744 – where the colours are becoming roughly constant, or the slope of the colour profile changes – is Rtrans ≈ 4arcsec (see the colour profile). At smaller radii R⋆the starburst region starts to dominate. By extrapolating the exponential slope of VCC 1744 to R⋆ = 0arcsec one obtains a central surface brightness of μ0,LSB = 21.6mag/arcsec2, two magnitudes fainter than the value when including the starburst region.
In contrast, the profile of VCC 0130 shows that a pure exponential approximation of the LSB-component would overestimate the central surface brightness of the LSB-component. Therefore, we assume a central flattening and obtain μ0,LSB = 22.1mag/arcsec2, which by construction (see the modexp fitting function in Sect. 3.4) is 1.75 mag/arcsec2 fainter than the unmodified profile. Since the optical image of VCC 0130 does show several star-forming regions, we consider the application of an inner flattening for separating the LSB and starburst component to be justified, even though the galaxy’s colour profile is rather flat.
4. Results
Derived structural parameters of the BCDs and their LSB-components.
Fig. 4 The Mr − (g − i) − Reff − ⟨μ⟩eff-plane of BCDs and ETGs. No decomposition into starburst and LSB-component was applied. Red squares represent the primary BCD sample, blue squares denote the candidate sample. Black dots represent the sample of early-type galaxies taken from Janz & Lisker (2008, 2009), for which colours were derived from the light within 2Reff. Green asterisks are dE(bc)s from Lisker et al. (2006b). For each (vertical) interval in magnitude or surface brightness, the (horizontal) ± 2σ-region of the ETGs is indicated by the black bars. VCC 0428 and 1313 have colours ((g − i)tot = −0.049 and −0.263mag, respectively), lying outside the colour range of the plot and are therefore only indicated by an arrow. |
Fig. 5 Left: colour−magnitude diagram of the LSB-components of our primary BCD sample, compared to the Virgo early-type dwarf galaxies (labelled “dE” but including dS0s), for which colours were derived within 2Reff. Red squares correspond to values of the LSB-components, where black crosses show LSB-components with inner profile flattening. The notation “LSB(corr)” indicates that the parameters of the LSB-components are corrected for the outer tail (see Sect. 3.5). The black vectors display the change of the parameters when an inner profile flattening is applied. Green asterisks correspond to dE(bc)s and blue squares to galaxies with a uncertain morphological classification (see Sect. 4.4). Typical errors are shown on the left-hand side. Right: same diagram, but additionally plotted are the irregulars (cyan triangles) and the entire population of early-type galaxies (black dots). In both panels we indicate the 2σ-region of the entire ETG population with vertical bars. |
Fig. 6 Effective radius Reff vs. absolute magnitude (top panels) and mean effective surface brightness (bottom panels). Symbols are the same as in Fig. 5. The black error bars in the left-hand panels illustrate the statistical errors, while the arrows indicate the systematic error, as derived from analysing artificial BCD images (see the Appendix). The error bars are placed on the ordinate so that they correspond to the average surface brightness and the magnitude of the respective set of artificial BCDs. The arrows point in the direction towards which the measurements would have to be shifted to correct for the systematic error. Note that the upper arrow is almost invisible due to the very small error. |
Fig. 7 Absolute magnitude Mr vs. mean effective surface brightness ⟨μ⟩eff,r in r-band. Symbols are the same as in Fig. 5. The black error bars in the left panel illustrate the statistical errors, while the arrows indicate the systematic error, as derived from analysing artificial BCD images (see the Appendix). The error bars are placed on the abscissa so that they correspond to the magnitude of the respective set of artificial BCDs. The arrows point in the direction towards which the measurements would have to be shifted to correct for the systematic error. Note that the left arrow is almost invisible due to its very low value. |
The results of the photometric analysis are summarised in Table 3. For each BCD values for the entire galaxy and for the underlying LSB-component are given (Cols. 5−12). The membership (ms, Col. 2) to the Virgo cluster was adapted from the VCC and later updates (see Sect. 2). The BCD-classification according to Table 1 is given in Col. 3. Column 4 indicates whether a flattening towards smaller radii was applied (see Sect. 3 for more details). Column 14 shows the total absolute magnitude of the entire BCD when applying the distance (Col. 13) given by the GOLDMine database7 (Gavazzi et al. 2003). The differences between a distance of 16.5 Mpc and the GOLDMine distance are discussed in detail in Sect. 5.1.
The structural and colour properties of the BCDs are illustrated in Fig. 4 for the entire galaxies, and in Figs. 5 through 7 for only the LSB-components. As a result of the decomposition, the LSB-components are less compact and redder than the entire galaxies.
4.1. Colour−magnitude diagrams
Figure 5 shows the colour−magnitude diagram (CMD)8 of the LSB-components (red squares) of our BCD sample. Galaxies with a profile flattening toward smaller radii are indicated with a black cross and the corresponding change in the magnitude is described by a vector. As mentioned in Sect. 3.4 the profile flattening parameters are uncertain; therefore, the vector is to be regarded as an aid to the eye, pointing to the locus of the diagram where the true values are expected to be.
Blue squares show the results for galaxies with an uncertain morphological classification (see also Sect. 4.4). Additionally, we plot with black dots the early-type galaxies of the Virgo cluster (Janz & Lisker 2009). According to the VCC these can be divided into the early-type dwarf (dE/dS0) and giant (E/S0) galaxies, based mainly on surface brightness and central concentration of light (Sandage & Binggeli 1984; Binggeli et al. 1985). In Figs. 5 through 7 the early-type dwarfs are shown in the left panels, while all early types are shown in the right panels. The horizontal 2σ-region of all early types is indicated in magnitude or surface brightness intervals by the vertical bars. The green asterisks correspond to dE(bc)s, i.e. early-type dwarfs with a blue core in the centre, taken from Lisker et al. (2006b).
On average the LSB-components of the BCDs are still bluer than the ETG population. The LSB-components of the BCDs show a large interval in colour with an average colour and 2σ-scatter of ⟨(g − i)⟩LSB = 0.64 ± 0.45mag, in contrast to the early-type dwarfs with ⟨g − i⟩dE = 0.88 ± 0.19mag.
There is a very blue BCD (VCC 1313) at (g − i)LSB ≈ 0.06mag. Such an extremely blue (g − i)LSB-colour can naturally arise from extended nebular emission, as is the case for the XBCDs SBS 0335-052 E (Papaderos et al. 1998) and I Zw 18 (Papaderos et al. 2002). Extended nebular emission could diminish morphological asymmetries; the structure of VCC 1313 thus needs to be interpreted with some caution. In distant, poorly resolved starburst galaxies, extended emission can easily be mistaken for a stellar disc (Papaderos et al. 2002; Papaderos & Östlin 2012).
From Fig. 5 it can be seen that the dE(bc)s are slightly offset from all early-type dwarfs on the CMD plot, showing bluer colours with an average value of ⟨g − i⟩ = 0.81mag. The majority of them is brighter than the LSB-components of the BCDs. While focusing on dE(bc)s, Lisker et al. (2006b) also analysed an additional sample of galaxies, which includes some BCDs of our study presented here (marked with “*” in Table 3). Due to the sample selection, all BCDs in the sample of Lisker et al. (2006b) show regular elliptical isophotes, corresponding to the nE type.
Average values of different ETG- and LSB-types.
The right hand side of Fig. 5 shows the same CMD but additionally the Virgo irregulars are plotted. No obvious separation between BCDs or irregulars can be found, since they cover nearly the same region in the CMD, also shown by their average colours (Table 4; see also Roediger et al. 2011). However, the irregulars tend to have redder colours at fainter magnitudes, with colours comparable to the early-type dwarfs. For the different subtypes of BCDs, the average g − i colours show a tendency for BCDs with regular LSB-components (nE or iE) to be redder than irregular shaped LSB-components (see Table 4).
4.2. Sizes of BCDs
Figure 6 shows the effective radius Reff vs. the absolute magnitude (top panels) and the mean effective r-band surface brightness ⟨μ⟩eff,r (bottom panels) of the LSB-components and the ETGs. The LSB-components are small in size when compared to the total early-type dwarf population, but partly cover the same region as the compact ~half of it, and mostly lie within the 2σ-region of all early types. Towards brighter ⟨μ⟩eff,r the LSB-components become more compact. However, compensating for the systematic errors of the derived parameters would move the fainter LSB-components slightly further towards the early-type dwarf region (black arrow centered on the error bars in the left-hand panels of Fig. 6).
Table 4 summarises the average values of the different subtypes of dwarf galaxies, showing that there are no significant morphological differences between the BCD subtypes. Since nE-BCDs have similarities to the dE(bc)s it could have been expected that their sizes are comparable. However, we find that the nE-BCDs, as well as the other BCD subtypes, are more compact than the dE(bc)s. The difference would be even slightly larger, if the (central) star-forming component of the dE(bc)s was also excluded for the derivation of their photometric parameters. The effect of this exclusion would, however, be much smaller in the dE(bc)s, since their central star-forming regions contribute less significantly to the total light of the galaxies than is the case for the BCDs. Comparing the LSB-components and the dE(bc)s in the Mr − Reff-plane (Fig. 6), one may conclude that dE(bc)s fall into the linear extension of the LSB-components toward larger Reff and brighter magnitudes.
The right hand side of Fig. 6 shows additionally the irregular galaxies. In both diagrams one finds a clear offset of the overall location of LSB-components and irregulars, with a transition region where both types can be found. At a given Reff the LSB-components of BCDs tend to be brighter and have a higher ⟨μ⟩eff,r.
Figure 7 shows the mean surface brightness ⟨μ⟩eff,r as a function of the absolute magnitude Mr for the LSB-components and early-type dwarfs on the left-hand side and additionally for the irregulars and all early types on the right-hand side. The BCDs and irregulars show a different average location in this diagram, but both overlap with the early-type dwarfs. Still, several LSB-components are outside the 2σ-region of the early-type dwarfs.
4.3. Apparent axis ratio of the LSB-component
Fig. 8 Distributions of the axis ratios of Virgo cluster dwarf galaxies. The upper panels show the irregulars and the LSB-components of the BCDs. The early-type dwarfs (labelled “dE”) also include the dE(bc)s, which are additionally plotted in the lower panel. The black vertical lines indicate the mean values. |
The axis ratios of the LSB-components have a mean value and standard deviation of ⟨b/a⟩LSB = 0.62 ± 0.18, compatible with the mean value of BCDs from Sung et al. (1998) (⟨b/a⟩BCD,Sung = 0.67) and very similar to the dE(bc)s (⟨b/a⟩dE(bc) = 0.63 ± 0.20). The mean value of the early-type dwarfs is slightly larger – i.e. they are slightly rounder – with ⟨b/a⟩dE = 0.69 ± 0.19. Similarly, the irregulars have ⟨b/a⟩dE = 0.67 ± 0.18.
Figure 8 shows the distribution of axis ratios of the LSB-components and additionally the distribution of irregulars, early-type dwarfs, and specifically dE(bc)s. While the distribution of the dE(bc)s comes closest to what would be expected from the intrinsic shape of a (thick) disc, i.e. an oblate flattened galaxy, none of the distributions shown differ significantly from each other. In addition to what is displayed in the figure, no noticeable trend with magnitude is found. Thus, no evolutionary connections can be ruled out based solely on the overall shape of the galaxy types.
4.4. Galaxies with uncertain morphological classification
Uncertain BCD candidates and suggested reclassification (Col. 3).
As outlined in Sect. 2, in addition to our primary working sample there are 21 uncertain BCD candidates in the VCC, 19 of which are covered by the SDSS DR5 (Table 5). Using the parameter space (Mr − (g − i)LSB − Reff − ⟨μ⟩eff) of the BCDs described in the previous sections, as well as the visual inspection of images, we are able to assess whether these galaxies should be reclassified as genuine BCDs or not.
VCC 1374, classified as “Im/BCD”, falls within the Reff and ⟨μ⟩eff values of the primary BCDs, but is brighter than nearly all of them. It is also not far from the bright end of the irregular galaxies’ locus. However, its appearance is indeed that of an irregular or Sm-type galaxy with strong star-forming regions, hence appropriately described by “Im/BCD”. The relatively small effective radius of 0.79 kpc is merely a result of its very elongated shape (axis ratio 0.32) implying an edge-on view – the half-light semimajor axis is 1.39 kpc instead.
For the galaxies VCC 0309 and VCC 2037, also classified as “Im/BCD”, we found that their parameters are in the same region where the irregular galaxies are located, which is also consistent with their visual appearance. We thus suggest to reclassify both as “Im”. For the rest of the galaxies classified as “Im/BCD”, a clear distinction between BCDs and irregulars is not possible, because they are located in the overlap zone between these two types. Their class thus remains unchanged.
There are three galaxies initially classified as “Sm / BCD” and “SBm / BCD”, respectively. VCC 1356 is located in the same parameter region as the BCDs, consistent with its visual appearance. We thus suggest to reclassify it to “BCD”. VCC 1725 and VCC 1791 do not show a clear alikeness or separation from the BCDs; their class remains unchanged.
The galaxy VCC 0281 was taken from the early-type sample of Janz & Lisker (2009) and was classified in the VCC as “dS0 or BCD”. It was also described by Lisker et al. (2006b) as a dE(bc). In the (g − i) − Mr-plane, VCC 0281 falls into the same region as several of the brighter BCDs. However, it lies at the edge of the structural parameter range of the BCDs, and well inside the early-type dwarf regime. We thus refrain from a reclassification.
Fig. 9 Image (left) and unsharp mask (right) of galaxy VCC 0213. Both panels have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. The unsharp mask was created using a two pixel Gaussian kernel. |
The classification of VCC 0135, 0655, and 1955 (“Spec/BCD”), as well as VCC 0213 (“dS?/BCD?”), implies the presence of hints for spiral structure. Indeed, VCC 0213 shows spiral structure in the SDSS image, and also for VCC 0135, a weak spiral-like feature is visible. For VCC 0655, an inner ring or spiral-like structure is mimicked by several star-forming knots. No spiral structure can be seen for VCC 1955. In the classification scheme of Sung et al. (2002), they reported that the “Postmerger” BCDs of their sample sometimes showed “spiral structures with a compact off-centred core”. However, to the knowledge of the authors, there is no detailed study about BCDs with spiral structure visible at optical wavelengths. HI observations of BCDs by van Zee et al. (2001) revealed rotation in BCDs, but the optical morphology of the used sample has smooth, symmetrical isophotes.
VCC 0213 is located at the bright and red end of the BCD region. In terms of structure, it has similar parameters as the BCDs. Its red core (Fig. 9) seems to contradict the classification as a BCD. The spiral structure, together with the red core – which appears similar to a bulge – might hint at a spiral galaxy with a greater distance than the one of the GOLDMine data base of 17.0 Mpc. However, the heliocentric velocity according to GOLDMine is vhelio = −165 km s-1, which is consistent with being a Virgo member (cf. Fig. 10) and makes it unlikely that the galaxy is a background object. We thus decide to leave its classification unaltered.
VCC 0135 lies at the edge of the structural locus of the brighter BCDs. However, its (g − i)-colours are too red, and its image looks like that of a low-luminosity lenticular or early-type spiral galaxy. The spectrum shows strong emission lines in the center, where blue features can be seen upon detailed inspection of the “core-region”. Therefore, we refrain from changing its classification.
VCC 0655 seems to be too bright (−18.73 mag) and too large (Reff = 1.49 kpc) compared to the BCDs. However, the image displays several star-forming regions in the central part of the galaxy, which are surrounded by a redder stellar host, typical for BCDs. Furthermore, VCC 0655 has some similarities to the well-known BCD NGC 2537 (Mrk 86, “Bear’s Paw Galaxy”). The sample of BCDs in Papaderos et al. (2008) also shows a wide spread in radius from very compact objects up to objects with . Thus, it cannot be ruled out that VCC 0655 is a bright BCD. We refrain, however, from reclassifying it, since deeper observations are needed to undoubtedly judge whether spiral structure is present.
At first glance, VCC 1955 appears similar to VCC 0655, except that its star-forming region is more centrally concentrated. However, its surface brightness is much lower, and in turn it has an even larger effective radius, too large compared to the other BCDs. Given its overall red colour and regular shape, VCC 1955 would qualify as a large-sized dE(bc), i.e. an early-type dwarf with a blue core. We therefore reclassify this object as “dE(bc)?”, but remark that deeper observations might reveal further structural features.
Finally, as mentioned in Sect. 2, visual inspection of VCC 0429 and VCC 1713 (classified as “?”) suggests a reclassification to “Im / BCD”. Indeed, their parameters fall in the overlap region of BCDs and irregulars. VCC 1411 was reclassified from “pec” to “BCD?” based on visual inspection.
The newly determined morphological types of the galaxies are listed in Table 5. For one galaxy (VCC 1356) we suggest a reclassification to a primary BCD class. For three other objects (VCC 309, 1955, and 2037), a unique type could be assigned, instead of a mixed type. All other galaxies indeed show complex properties that can only be described with a combination of two types.
5. Discussion
From Table 6 it is apparent that various photometric parameters of certain and possible Virgo cluster members agree within 1σ. Therefore, in the following discussion no differentiation based on Virgo membership was applied.
Fig. 10 GOLDMine distance vs. the velocity of the Virgo BCDs. The vertical solid line at vLG = 1900 km s-1 corresponds to the limit from which galaxies belong to the M and W cloud according to Gavazzi et al. (1999). The three BCDs in the box lie in the projected cloud region, but are below or only slightly above this velocity limit. We therefore adopt a distance of 16.5 Mpc for them instead of the GOLDMine distance of 32 Mpc. |
Comparison of the BCDs regarding their membership class.
5.1. Distance dependence of the results
GOLDMine distances of the different Virgo clouds.
The GOLDMine data base provides individual distances to the BCDs, which can be found in Col. 13 in Table 3. The distances are based on their association with various substructures in the projected Virgo cluster area (Gavazzi et al. 1999). These were defined based on the Tully-Fisher relation (Tully & Fisher 1977) for a sample of spiral galaxies and the Fundamental Plane (Djorgovski & Davis 1987) for E/S0 galaxies, revealing separate regions in the parameter space of velocity and distance. Table 7 shows the different distances and the resulting magnitude offsets when adopting the GOLDMine distances instead of d = 16.5 Mpc (m − M = 31.09) mag. The factor by which Reff [kpc] changes due to different distances is given by α(Reff).
Since the allocation of the BCDs to the different substructures of the Virgo cluster only depends on the projected location of the galaxies, the GOLDMine distances can not be taken for granted. For instance, the column depth of the different substructures – and of the main cluster itself – cannot be taken into account.
Fig. 11 Same as Figs. 6 and 7, but instead of a constant distance of 16.5 Mpc, the distances given by the GOLDMine database are used for the LSB-components of the primary sample, except for the three objects (shown with magenta symbols) marked in Fig. 10, which remain at 16.5 Mpc. All other galaxy types also remain at a constant distance of 16.5 Mpc. |
Fig. 12 Parameter space of the LSB-components of our Virgo BCDs and a comparison sample of star-forming dwarf galaxies (SF-dwarfs) in various environments. The black squares correspond to the comparison sample of star-forming dwarf galaxies, while the other galaxy types are colour coded as in the previous figures. |
Figure 10 shows the distance of the Virgo BCDs against the velocity vLG with respect to the Local Group (cf. Fig. 5 of Gavazzi et al. 1999), calculated via vLG = vhelio + 220 km s-1. The heliocentric velocity vhelio was given by GOLDMine. According to Gavazzi et al. (1999) galaxies with vLG > 1900 km s-1 belong to the M or W cloud with a distance of 32 Mpc. As one can see, three BCDs (VCC 0022, 0024, and 0274) have a GOLDMine distance of 32 Mpc, but have velocities either below or roughly at vLG = 1900 km s-1. Therefore, we do not assume a distance of 32 Mpc for these galaxies, but instead a distance of 16.5 Mpc.
Figure 11 shows the (Mr − Reff − ⟨μeff⟩)-plots of the analysed galaxies with the GOLDMine distances. Interestingly, there are only three LSB-components of the BCDs that do not fall into the ETGs’ 2σ-plane, and these BCDs are classified as “iI” and “iI,C”. One of these three BCDs shows a profile flattening of the LSB-component, shifting it more to the 2σ-region of the early types. In the Mr − Reff-diagram, there are two additional BCDs outside the 2σ-region of early-type dwarfs, for which the above velocity criterion was not fulfilled and thus no greater distance was assigned.
The comparison of the late-type dwarfs (BCD(-LSB) and irregular) and the early-type dwarfs reveals another interesting result. The LSB-components of the BCDs and the irregular galaxies form a continuous sequence, which together shares the same location as the early-type dwarf galaxies. This implies that, structurally, the Virgo BCDs could evolve into early-type dwarfs with small radii once star formation ceases in them. Additionally, Virgo irregulars could form the progenitor population of future early-type dwarfs with large radii. Moreover, it should be pointed out that these findings are dependent on the distance of the galaxies. As seen in Fig. 6, there are more BCDs outside the 2σ-region when a constant distance of 16.5 Mpc is adopted. We point out that we have applied the individual GOLDMine distances only to the LSB-components of the primary BCD sample, in order to illustrate the effect on their values. In contrast, this would only have a weak effect on the population of early types: firstly, the vast majority of them are not located in the projected substructure regions that contain galaxies at significantly greater distances, and secondly, those that are located there would be shifted largely along the Mr − Reff-region that they already populate (cf. Table 7).
Extremely metal-poor BCDs from the study of Papaderos et al. (2008).
5.2. Comparison with other BCDs
Table 8 shows the values for a sample of BCDs from Papaderos et al. (2008; hereafter P08) from the SDSS. All given values of P08 were measured in the SDSS g-band and for a first order approximation we assume that Reff is equal in the g and r filters. As one can see, the effective radii of the P08 sample are also very compact with an average of , but with a large scatter.
Since the BCDs of our study are located in the Virgo cluster, one may ask whether they deviate from BCDs in field and group environments. To investigate this question a large sample of 340 emission-line galaxies in various environments (unpublished study by H. T. Meyer, R. Kotulla, P. Papaderos, Y. Izotov, N. G. Guseva and K. J. Fricke) was used. The galaxies were spectroscopically selected, therefore, also non-dwarf galaxies are included in the initial sample. To include only star-forming dwarf galaxies, a g-band magnitude cutoff of Mcutoff = −18.0 mag was applied, resulting in the sample that we use here for comparison. This sample includes BCDs as well as other types of star-forming (dwarf)9 galaxies.
To all galaxies of this comparison sample the same photometric analysis (derivation of profiles and decomposition) was applied. The CMD of the LSB-components of the Virgo BCDs and the comparison sample is shown in Fig. 12. The LSB-components of the Virgo BCDs tend to be redder than those of the comparison galaxies. This could indicate an influence of the denser cluster environment on the Virgo BCDs, making them redder due to the removal of gas. In the structural diagrams of Fig. 12, one can see that the locus of the Virgo BCDs nearly completely overlaps with that of the comparison sample of star-forming dwarf galaxies in the field. The fact that the latter extends to larger effective radii and brighter magnitudes than the Virgo sample may reflect a selection effect. An increase in the sample size of cluster BCDs would therefore be desirable for future investigations.
5.3. Evolutionary connection between different galaxy types
Various observational and theoretical lines of evidence suggest that the evolutionary pathways of late-type dwarfs in a cluster environment have to significantly differ from those in the field. Several mechanisms for starburst ignition in the initial stages of the interaction of a dwarf with the ICM can be envisaged, for example triggered gas collapse by the external ICM pressure (as opposed to gas collapse driven by self-gravity), or strong dissipation and gas cooling in large-scale shocks, followed by collective star formation. For such reasons, it is not even sure that the compact structure of the LSB-component – a typical characteristic of field BCDs – has to be common in systems classified as BCDs in galaxy clusters. For example, external agents (see above) could ignite a BCD-typical starburst in a genuine (that is, relatively diffuse) irregular. Alternatively, an externally triggered destabilisation and inflow of gas could eventually lead to an adiabatic contraction of the LSB-component (Papaderos et al. 1996a), and the opposite may happen after complete removal of the dwarf gas halo, as it sinks deeper into the hostile ICM environment.
Even in the cluster periphery, late-type dwarfs have undergone interaction with the hot ICM and other nearby galaxies for a few hundreds of Myr, a time span on the order of their dynamical timescale. In summary, in a late-type dwarf plunging into the ICM, the starburst that makes it appear as BCD may have a different origin than starbursts in the general population of field BCDs. Consequently, one should be cautious when comparing structural and integral photometric properties of cluster-BCDs with literature data for ordinary BCDs.
As our understanding of the early interaction of late-type dwarfs with the cluster environment is still poor, it would be speculative to generalise conclusions drawn from our Virgo BCD sample to the BCD population as a whole. Nevertheless, the following discussion concentrates on a nearly complete and well-selected sample of cluster BCDs, thereby minimising the potential dilution of any trends that could be caused when combining objects from different environments.
The evolutionary scenario of Davies & Phillipps (1988) describes how a dwarf irregular galaxy can ignite a starburst and become a BCD due to infalling gas from a surrounding reservoir. After several starburst phases and non-starbursting irregular-phases, the galaxy evolves to a passive dwarf galaxy. If this evolution was possible in the Virgo cluster outskirts, we should observe a similar structure of BCDs (the starburst phase) and early-type dwarfs (the end product). While this may serve as a possible evolutionary path from the faint end of the late-type sequence of Kormendy & Bender (2012) to the faint end of the early-type sequence, we note that early-type and late-type galaxies in today’s clusters have spent, on average, a large fraction of their lifetime in different environments (Lisker et al. 2013).
Observational studies of BCDs by Drinkwater & Hardy (1991)10 and Papaderos et al. (1996a) reported that BCDs are very compact objects as compared to the early-type dwarfs. They concluded that only the most compact of the latter may be related to the BCDs. An even more extreme conclusion is reached by Drinkwater et al. (1996), who state that “no blue star-forming progenitors of dE galaxies” exist in the Virgo cluster down to an absolute B-magnitude of MB ≈ −14 mag. At first view, these findings seemed to be supported by our photometric measurements for the galaxies’ entire light: Fig. 4 (top left panel) shows that the majority of BCDs are more compact than even the 2σ-region of early-type dwarfs at a given magnitude.
However, if the gas reservoir that is immediately available to a BCD’s starburst region is limited, thus restricting the starburst phase to a duration of some 100 Myr or less (Thuan 1991; Thornley et al. 2000), a population of non-starbursting counterparts must exist within the same spatial volume. Assuming a velocity of 1000 km s-1, a BCD could only move some 100 kpc during the starburst phase, i.e. some tenths of the Virgo cluster’s virial radius – whichwould mean that the outskirts of the Virgo cluster should hold a substantial number of non-starbursting counterparts to BCDs11. How do these look like?
Fig. 13 Possible descendants of BCDs (left) and of irregulars (right) in the Virgo cluster. The left panels show a bright (top) and a faint (bottom) early-type dwarf that fall into the structural parameter range of the LSB-components of Virgo BCDs, i.e. they are more compact than the average early-type dwarf. The right panels show a bright (top) and faint (bottom) early-type dwarf that fall into the range of the Virgo irregulars, i.e. they are more diffuse than the average early-type dwarf. All images have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. The galaxies and their parameters are, clockwise from top left: VCC 0033 (Mr = −16.85 mag, Reff = 0.73 kpc), VCC 0510 (Mr = −17.00 mag, Reff = 1.64 kpc), VCC 1565 (Mr = −15.49 mag, Reff = 1.40 kpc), and VCC 0068 (Mr = −15.25 mag, Reff = 0.53 kpc). |
To investigate this, we determined the structural properties of the underlying LSB-component of the Virgo BCDs, i.e. without the contribution of the starburst region. The results already show much more overlap with the early-type dwarfs in the parameter space of structure and also colour (see Figs. 5 and 6). When we also take the different distances of the various parts of the Virgo region into account, the LSB-components turn out to fit remarkably well to the early-type dwarf population with small effective radii (Fig. 11). This is in agreement with the conclusions of Micheva et al. (2013) from deep imaging of a sample of luminous BCDs and intermediate-mass galaxies. The left panels of Fig. 13 show examples for such early-type dwarfs that fall within the locus of the LSB-components: these may look like the future “red and dead” descendants of BCDs. In contrast, the right panels show early-type dwarfs that are more diffuse than the LSB-components, and instead overlap with the irregulars.
Fig. 14 The two BCDs (left: VCC 0410, right: VCC 1313) that remain more compact than the early-type dwarfs even after applying the individual distances of GOLDMine. These objects lie outside of the 2σ-region of early-type dwarfs in the ⟨μeff⟩ − Reff diagram of Fig. 11. Both images have a width of 80”, corresponding to 6.3 kpc at a distance of 16.5 Mpc. |
As for the two objects (Fig. 14) that remain more compact than the 2σ-region of early types in the ⟨μeff⟩ − Reff diagram (Fig. 11), their LSB-component colours are also bluer than the 2σ-region of early types. VCC 1313 has hardly any discernible LSB-component within our detection limits. This galaxy seemsto follow the literal meaning of the term “blue compact dwarf” (see Fig. 14, right panels), while most of the other Virgo BCDs rather have the appearance of a (sometimes faint) early-type dwarf or irregular with additional starburst. The non-detection of a red LSB-component down to low surface brightness levels for some “extreme” BCDs was already reported by Drinkwater & Hardy (1991): two BCDs in their study had no surrounding low-surface-brightness envelopes detected to a limit of μ(V) = 27mag/arcsec2. Therefore, they would be classified as i0-BCDs. For all of these, the Next Generation Virgo Cluster Survey (NGVS, Ferrarese et al. 2012) will provide unprecedented insight into their outer structure and older stellar population, if existent.
Fig. 15 Histogram of morphological types in the outer region of Virgo (DM87 > 0.6Mpc). A cutoff of Mr > −16.5mag and ⟨μr⟩eff,r < 23.0mag/arcsec2 was applied to restrict the parameter region to the locus of most of the LSB-components (see Fig. 7), using a distance of 16.5 Mpc for all galaxies. For the red histogram lines, an additional colour cutoff of (g − i) > 0.7 mag was applied. The histogram of the early-type dwarfs (labelled “dE”) also includes the dE(bc)s. Candidate BCDs are included with their respective first-mentioned class (e.g. “Sm/BCD”). The numbers are obtained from the full sample of Virgo cluster galaxies that was used in Weinmann et al. (2011), based on the VCC and adopting a faint magnitude limit of 18 mag in the B-band. For the galaxies that could not be analysed with the SDSS, we use estimations for Mr as described in the appendix of Weinmann et al. (2011), and adopt the average surface brightness and colour of the respective class. |
Are early-type dwarfs the only morphological type that qualifies as non-starbursting counterparts of BCDs? Fig. 15 shows the number distribution of galaxies in the outer region of the Virgo cluster that fall in the same parameter space of magnitude and mean surface brightness as the majority of the LSB-components of the BCDs, displayed separately for blue and red (g − i > 0.7) galaxies. Among blue galaxies, late-type spirals and especially irregular galaxies (Sánchez Almeida et al. 2008) form a population of significant size. However, the majority of the LSB-components of BCDs is already in the red colour regime (red histogram in Fig. 15), where the early-type dwarfs clearly outnumber all other galaxy types.
Vaduvescu et al. (2006) concluded from a near-infrared analysis of Virgo cluster galaxies, decomposed with a different scheme than the one used in this study, that “BCDs and dIs are similar structurally”, which is incompatible with our conclusions (see Fig. 6). We measure irregulars to be systematically more extended than the LSB-components of BCDs at a given magnitude, hence only a small fraction of the irregulars would qualify as non-starbursting BCDs. Vaduvescu et al. fitted structural parameter relations to irregulars, and claimed consistency of the BCDs with these relations. Their Figs. 4 and 5 show large scatter among their BCD sample, with nine of the 18 BCDs having substantially smaller radii than the relation of irregulars. Among the five BCDs that lie significantly above the relation of irregulars, the three largest ones have been classified as ambiguous types between Sm, irregular and BCD in the VCC, and are more similar to Sm/Im by visual inspection of SDSS colour images. This may provide an explanation for the discrepant findings.
Indications that the Virgo BCDs have already experienced some environmental influence are given by the redder LSB-components of the Virgo BCDs as compared to a large BCD sample from various environments (Sect. 5.2), as well as by the fact that the integrated colours of many Virgo BCDs are not dominated by the starburst component (Fig. 4), consistent with their visual appearance but unlike typical field BCDs (see Sect. 1). Due to these differences, we cannot be certain that the commonly assumed short duration of a BCD’s starburst phase also holds for the Virgo BCDs. If their starbursts are indeed relatively short, then many of the early-type dwarfs in the Virgo outskirts with similar magnitudes must be the non-starbursting counterparts to the BCDs. Unless we are witnessing all BCDs in their very last starburst phase, which seems unlikely, this would mean that a number of early-type dwarfs must be able to re-ignite a starburst.
At first glance, this may seem consistent with the recent or even ongoing star formation in the dE(bc)s of the Virgo cluster (Lisker et al. 2006b), as well as with the interpretation of Hallenbeck et al. (2012) that gas detection in some Virgo early-type dwarfs could imply re-accretion of gas. However, De Looze et al. (2013) recently reported from a Herschel far-infrared analysis that BCDs contain on average more dust than early-type dwarfs with residual star formation (“transition-type dwarfs”), and we find that the dE(bc)s are on average less compact and brighter than the LSB-components of the BCDs (see Fig. 6). Future high-sensitivity comparative studies of the stellar and gas distribution and dynamics in those galaxy types are needed for a conclusive interpretation of their possible evolutionary links.
Online material
Appendix A: Error estimation
To obtain a better insight into the errors of the derived parameters of the BCDs, we created a set of artificial BCDs. To mimick the real Virgo BCDs, we combined the light distribution from two components, representing the LSB and starburst component. Firstly, the LSB-components were divided into a bright version (Mr, LSB = −17.0mag), with half-light semi-major axes of ahl = 7.6arcsec and ahl = 15.1arcsec, and a faint version (Mr,LSB = −15.0mag), with ahl = 6.3arcsec and ahl = 12.6arcsec. Secondly, two different Sérsic indices n were used, with n = 1.0 (pure exponential) and n = 1.3. We chose three different axis ratios (b/a = 0.5,0.7,0.9) and two position angles (pa = 60, 90) for the LSB component; the starburst component was chosen to be circular.
To mimick the additional outer tail of the LSB-component of some of the real Virgo-BCDs (see Sect. 3.5) at larger radii, we optionally added an exponential outer tail to the exponential LSB-component, having the same axis ratio and position angle. The tail has twice the half-light semi-major axis of the LSB-component, and can be weak (Δm = 1.2 mag fainter in total brightness) or strong (Δm = 0.75mag). Note that we did not add a tail to the LSB component with n = 1.3, since our intention was to test whether our approach to assume an exponential shape for the LSB component’s analysis, and to account for a possible outer tail, would be able to approximate the n = 1.3 case reasonably well.
For the starburst component, two different effective radii (2.0 and 4.0arcsec) were used to account for the various extents of the starbursts. We created one version in which the starburst is co-aligned with the LSB component, and another version in which it is offset by 0.5ahl, since a fraction of the real BCDs show a clear offset between starburst and LSB component (see Table 1).
With these different parameters we end up with 384 artificial BCDs. We first created noise-free images without the starburst component, to determine the true total (two-Petrosian) magnitude Mr,LSB of the LSB+tail component, as well as the true value of ahl, using the input axis ratios and position angles. We then created images with realistic noise and seeing characteristics, which were analysed in the same manner as the Virgo BCDs (see Sect. 3). The magnitudes and radii derived with LAZY for the LSB+tail component were then compared to the true values, thereby yielding estimates for systematic errors of magnitude and relative radius, as well as for their statistical errors (standard deviation) σMr and σReff. Note that here and in the remainder of this section, we omit the subscript “LSB” from all quantities, for better readability.
Mean values of the derived parameters for the LSB components of the artificial BCDs for different Sérsic indices n.
Table A.1 summarises the results, subdivided by LSB component brightness and profile shape. Based on this approach, the final systematic and statistical errors that we adopt for our BCD analysis are provided in Table A.2, subdivided into bright
and faint LSB components. We refrain here from differentiating between the different Sérsic indices n and between the cases with/without a tail component, as the resulting values were very similar.
The mean effective surface brightness ⟨μ⟩eff can be derived directly from Mr and Reff: (A.3)Its statistical error was calculated assuming independent errors on magnitude and radius: The systematic surface brightness error was calculated by (A.6)The average differences between the input magnitudes and radii and the ones measured by LAZY are almost negligible for the bright galaxies (0.03 mag / 1%) and moderate for the faint galaxies (0.15 mag/11%). However, both add up to a significant systematic effect on the surface brightness of 0.41 mag/arcsec2 for the faint galaxies, while being negligible (0.01 mag/arcsec2) for the bright galaxies. The statistical errors are reasonably small for all quantities, with 0.10 mag for Mr, 6% for Reff, and 0.16 mag/arcsec2 for ⟨μ⟩eff.
We point out that, due to our parameter setup, the fraction of artificial BCDs with inner flattening (see Sect. 3.4) was much larger than for the real Virgo BCDs, where only three BCDs were fitted with a inner flattening. We therefore calculated the above quantities only from those BCDs without inner flattening, to avoid an overestimation of the errors.
Note: The LSB-component of BCDs should not be mixed up with the galaxy type commonly called “Low Surface Brightness galaxies”. To the contrary, a high central surface brightness (μ0(B) < 22 mag/arcsec2) has been found to be a characteristic property of the LSB-host of BCDs (Papaderos et al. 1996a; Gil de Paz & Madore 2005) implying a higher central mass density than in typical irregulars, and the more so, in dSphs and genuine (μ0 ≥ 23.5 mag/arcsec2) LSB galaxies.
Using the online Image List Tool of the SDSS, http://skyserver.sdss3.org/dr8/en/tools/chart/list.asp
Lisker et al. (2007) used so-called asymmetry centering for the early-type dwarfs instead.
This should not be confused with the colour−magnitude diagrams used in the studies of resolved stellar populations. We refrain from using the common term “colour−magnitude relation” (CMR) for Fig. 5, as it may not be clear a priori whether an actual relation exists.
It should be pointed out that the findings of Drinkwater & Hardy (1991) were based on a particular sample of the most compact BCDs in Virgo.
This does not necessarily apply to BCDs that are located at greater distances than the main cluster, which, however, is the minority (Fig. 10).
Acknowledgments
The authors would like to thank the anonymous referee for helpful suggestions. H.T.M., T.L., and J.J. were supported within the framework of the Excellence Initiative by the German Research Foundation (DFG) through the Heidelberg Graduate School of Fundamental Physics (grant number GSC 129/1). H.T.M. was supported by the DFG through grant LI 1801/2-1. J.J. acknowledges support by the Gottlieb Daimler and Karl Benz Foundation. P.P. is supported by Ciencia 2008 Contract, funded by FCT/MCTES (Portugal) and POPH/FSE (EC). He also acknowledges support by the Fundação para a Ciência e a Tecnologia (FCT) under project FCOMP-01-0124-FEDER-029170 (Reference FCT PTDC/FIS-AST/3214/2012), funded by FCT-MEC (PIDDAC) and FEDER (COMPETE). Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the U.S. Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, 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 Max-Planck-Institute for Astronomy (MPIA), the Max-Planck-Institute for Astrophysics (MPA), New Mexico State University, Ohio State University, University of Pittsburgh, University of Portsmouth, Princeton University, the United States Naval Observatory, and the University of Washington.
References
- Adelman-McCarthy, J. K., Agüeros, M. A., Allam, S. S., et al. 2007, ApJS, 172, 634 [NASA ADS] [CrossRef] [Google Scholar]
- Amorín, R., Alfonso, J., Aguerri, J. A. L., Muñoz-Tuñón, C., & Cairós, L. M. 2009, A&A, 501, 75 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bekki, K. 2008, MNRAS, 388, L10 [Google Scholar]
- Bertin, E., & Arnouts, S. 1996, A&AS, 117, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Binggeli, B., & Cameron, L. M. 1991, A&A, 252, 27 [NASA ADS] [Google Scholar]
- Binggeli, B., Sandage, A., & Tammann, G. A. 1985, AJ, 90, 1681 [NASA ADS] [CrossRef] [Google Scholar]
- Binggeli, B., Tarenghi, M., & Sandage, A. 1990, A&A, 228, 42 [NASA ADS] [Google Scholar]
- Binggeli, B., Popescu, C. C., & Tammann, G. A. 1993, A&AS, 98, 275 [NASA ADS] [Google Scholar]
- Blanton, M. R., Dalcanton, J., Eisenstein, D., et al. 2001, AJ, 121, 2358 [NASA ADS] [CrossRef] [Google Scholar]
- Boselli, A., & Gavazzi, G. 2006, PASP, 118, 517 [NASA ADS] [CrossRef] [Google Scholar]
- Cairós, L. M., Caon, N., Vílchez, J. M., González-Pérez, J. N., & Muñoz-Tuñón, C. 2001a, ApJS, 136, 393 [NASA ADS] [CrossRef] [Google Scholar]
- Cairós, L. M., Vílchez, J. M., González Pérez, J. N., Iglesias-Páramo, J., & Caon, N. 2001b, ApJS, 133, 321 [NASA ADS] [CrossRef] [Google Scholar]
- Cairós, L. M., Caon, N., Papaderos, P., et al. 2003, ApJ, 593, 312 [NASA ADS] [CrossRef] [Google Scholar]
- Davies, J. I., & Phillipps, S. 1988, MNRAS, 233, 553 [NASA ADS] [CrossRef] [Google Scholar]
- De Looze, I., Baes, M., Boselli, A., et al. 2013, MNRAS, 436, 1057 [NASA ADS] [CrossRef] [Google Scholar]
- Djorgovski, S., & Davis, M. 1987, ApJ, 313, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Dressler, A. 1980, ApJ, 236, 351 [NASA ADS] [CrossRef] [Google Scholar]
- Drinkwater, M., & Hardy, E. 1991, AJ, 101, 94 [NASA ADS] [CrossRef] [Google Scholar]
- Drinkwater, M. J., Currie, M. J., Young, C. K., Hardy, E., & Yearsley, J. M. 1996, MNRAS, 279, 595 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, B. G., & Elmegreen, D. M. 2010, ApJ, 722, 1895 [NASA ADS] [CrossRef] [Google Scholar]
- Elmegreen, B. G., Zhang, H.-X., & Hunter, D. A. 2012, ApJ, 747, 105 [NASA ADS] [CrossRef] [Google Scholar]
- Ferrara, A., & Tolstoy, E. 2000, MNRAS, 313, 291 [NASA ADS] [CrossRef] [Google Scholar]
- Ferrarese, L., Côté, P., Cuillandre, J.-C., et al. 2012, ApJS, 200, 4 [NASA ADS] [CrossRef] [Google Scholar]
- Gavazzi, G., Boselli, A., Scodeggio, M., Pierini, D., & Belsole, E. 1999, MNRAS, 304, 595 [NASA ADS] [CrossRef] [Google Scholar]
- Gavazzi, G., Boselli, A., Mayer, L., et al. 2001, ApJ, 563, L23 [NASA ADS] [CrossRef] [Google Scholar]
- Gavazzi, G., Boselli, A., Donati, A., Franzetti, P., & Scodeggio, M. 2003, A&A, 400, 451 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gavazzi, G., Fumagalli, M., Cucciati, O., & Boselli, A. 2010, A&A, 517, A73 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gil de Paz, A., & Madore, B. F. 2005, ApJS, 156, 345 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Graham, A. W., & Driver, S. P. 2005, Publ. Astron. Soc. Australia, 22, 118 [Google Scholar]
- Gu, Q., Zhao, Y., Shi, L., Peng, Z., & Luo, X. 2006, AJ, 131, 806 [NASA ADS] [CrossRef] [Google Scholar]
- Gunn, J. E., & Gott, III, J. R. 1972, ApJ, 176, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Guseva, N. G., Papaderos, P., Izotov, Y. I., et al. 2003, A&A, 407, 105 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Hallenbeck, G., Papastergis, E., Huang, S., et al. 2012, AJ, 144, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Hensler, G. 2012, Morphological Mutations of Dwarf Galaxies, eds. P. Papaderos, S. Recchi, & G. Hensler, Astrophys. Space Sci. Proc., 75 [Google Scholar]
- Izotov, Y. I., Chaffee, F. H., Foltz, C. B., et al. 1999, ApJ, 527, [Google Scholar]
- Janz, J., & Lisker, T. 2008, ApJ, 689, L25 [NASA ADS] [CrossRef] [Google Scholar]
- Janz, J., & Lisker, T. 2009, ApJ, 696, L102 757 [NASA ADS] [CrossRef] [Google Scholar]
- Janz, J., Laurikainen, E., Lisker, T., et al. 2012, ApJ, 745, L24 [NASA ADS] [CrossRef] [Google Scholar]
- Janz, J., Laurikainen, E., Lisker, T., et al. 2013, ApJS, accepted [arXiv:1308.6496] [Google Scholar]
- Jedrzejewski, R. I. 1987, MNRAS, 226, 747 [Google Scholar]
- Jerjen, H. 2012, Dwarf Elliptical Galaxies: United and Divided, eds. P. Papaderos, S. Recchi, & G. Hensler, 133 [Google Scholar]
- Kniazev, A. Y., Pustilnik, S. A., Grebel, E. K., Lee, H., & Pramskij, A. G. 2004, ApJS, 153, 429 [NASA ADS] [CrossRef] [Google Scholar]
- Kormendy, J. 1985, ApJ, 295, 73 [NASA ADS] [CrossRef] [Google Scholar]
- Kormendy, J., & Bender, R. 2012, ApJS, 198, 2 [NASA ADS] [CrossRef] [Google Scholar]
- Krüger, H., Fritze-v. Alvensleben, U., & Loose, H.-H. 1995, A&A, 303, 41 [NASA ADS] [Google Scholar]
- Kunth, D., & Östlin, G. 2000, A&ARv, 10, 1 [NASA ADS] [CrossRef] [Google Scholar]
- Lelli, F., Verheijen, M., Fraternali, F., & Sancisi, R. 2012, A&A, 537, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Lisker, T., Grebel, E. K., & Binggeli, B. 2006a, AJ, 132, 497 [NASA ADS] [CrossRef] [Google Scholar]
- Lisker, T., Glatt, K., Westera, P., & Grebel, E. K. 2006b, AJ, 132, 2432 [NASA ADS] [CrossRef] [Google Scholar]
- Lisker, T., Grebel, E. K., Binggeli, B., & Glatt, K. 2007, ApJ, 660, 1186 [NASA ADS] [CrossRef] [Google Scholar]
- Lisker, T., Grebel, E. K., & Binggeli, B. 2008, AJ, 135, 380 [NASA ADS] [CrossRef] [Google Scholar]
- Lisker, T., Weinmann, S. M., Janz, J., & Meyer, H. T. 2013, MNRAS, 432, 1162 [NASA ADS] [CrossRef] [Google Scholar]
- Loose, H.-H., & Thuan, T. X. 1986a, ApJ, 309, 59 [NASA ADS] [CrossRef] [Google Scholar]
- Loose, H.-H., & Thuan, T. X. 1986b, in Star-Forming Dwarf Galaxies and Related Objects, eds. D. Kunth, T. X. Thuan, & J. Tran Thanh van, 73 [Google Scholar]
- Marquart, T., Fathi, K., Östlin, G., et al. 2007, A&A, 474, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Mas-Hesse, J. M., & Kunth, D. 1999, A&A, 349, 765 [NASA ADS] [Google Scholar]
- Mastropietro, C., Moore, B., Mayer, L., et al. 2005, MNRAS, 364, 607 [Google Scholar]
- Mayer, L., Governato, F., Colpi, M., et al. 2001, ApJ, 547, L123 [NASA ADS] [CrossRef] [Google Scholar]
- McQuinn, K. B. W., Skillman, E. D., Cannon, J. M., et al. 2010, ApJ, 724, 49 [NASA ADS] [CrossRef] [Google Scholar]
- Mei, S., Blakeslee, J. P., Côté, P., et al. 2007, ApJ, 655, 144 [NASA ADS] [CrossRef] [Google Scholar]
- Micheva, G., Östlin, G., Bergvall, N., et al. 2013, MNRAS, 431, 102 [NASA ADS] [CrossRef] [Google Scholar]
- Moore, B., Katz, N., Lake, G., Dressler, A., & Oemler, A. 1996, Nature, 379, 613 [NASA ADS] [CrossRef] [Google Scholar]
- Noeske, K. G. 1999, Diploma thesis, Göttingen [Google Scholar]
- Noeske, K. G., Guseva, N. G., Fricke, K. J., et al. 2000, A&A, 361, 33 [NASA ADS] [Google Scholar]
- Noeske, K. G., Papaderos, P., Cairós, L. M., & Fricke, K. J. 2003, A&A, 410, 481 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Noeske, K. G., Koo, D. C., Phillips, A. C., et al. 2006, ApJ, 640, L143 [NASA ADS] [CrossRef] [Google Scholar]
- Östlin, G., Amram, P., Masegosa, J., Bergvall, N., & Boulesteix, J. 1999, A&AS, 137, 419 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papaderos, P., & Östlin, G. 2012, A&A, 537, A126 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papaderos, P., Loose, H.-H., Fricke, K. J., & Thuan, T. X. 1996a, A&A, 314, 59 [NASA ADS] [Google Scholar]
- Papaderos, P., Loose, H.-H., Thuan, T. X., & Fricke, K. J. 1996b, A&AS, 120, 207 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papaderos, P., Izotov, Y. I., Fricke, K. J., Thuan, T. X., & Guseva, N. G. 1998, A&A, 338, 43 [NASA ADS] [Google Scholar]
- Papaderos, P., Fricke, K. J., Thuan, T. X., Izotov, Y. I., & Nicklas, H. 1999, A&A, 352, L57 [NASA ADS] [Google Scholar]
- Papaderos, P., Noeske, K. G., Cairós, L. M., Vílchez, J. M., & Fricke, K. J. 2001, in Dwarf galaxies and their environment, eds. K. S. de Boer, R.-J. Dettmar, & U. Klein, 283 [Google Scholar]
- Papaderos, P., Izotov, Y. I., Thuan, T. X., et al. 2002, A&A, 393, 461 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papaderos, P., Guseva, N. G., Izotov, Y. I., et al. 2006, A&A, 457, 45 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Papaderos, P., Guseva, N. G., Izotov, Y. I., & Fricke, K. J. 2008, A&A, 491, 113 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Peng, C. Y., Ho, L. C., Impey, C. D., & Rix, H.-W. 2010, AJ, 139, 2097 [NASA ADS] [CrossRef] [Google Scholar]
- Petrosian, V. 1976, ApJ, 209, L1 [NASA ADS] [CrossRef] [Google Scholar]
- Recchi, S., & Hensler, G. 2006, A&A, 445, L39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roediger, J. C., Courteau, S., McDonald, M., & MacArthur, L. A. 2011, MNRAS, 416, 1983 [NASA ADS] [CrossRef] [Google Scholar]
- Salzer, J. J., & Norten, S. A. 1999, The Low Surface Brightness Universe, IAU Col., 171, 170 [Google Scholar]
- Salzer, J. J., MacAlpine, G. M., & Boroson, T. A. 1989, ApJS, 70, 479 [NASA ADS] [CrossRef] [Google Scholar]
- Sánchez Almeida, J., Muñoz-Tuñón, C., Amorín, R., et al. 2008, ApJ, 685, 194 [NASA ADS] [CrossRef] [Google Scholar]
- Sandage, A., & Binggeli, B. 1984, AJ, 89, 919 [NASA ADS] [CrossRef] [Google Scholar]
- Schlegel, D. J., Finkbeiner, D. P., & Davis, M. 1998, ApJ, 500, 525 [NASA ADS] [CrossRef] [Google Scholar]
- Searle, L., & Sargent, W. L. W. 1972, ApJ, 173, 25 [NASA ADS] [CrossRef] [Google Scholar]
- Sérsic, J. L. 1968, Atlas de galaxias australes (Cordoba, Argentina: Observatorio Astronomico) [Google Scholar]
- Silich, S. A., & Tenorio-Tagle, G. 1998, MNRAS, 299, 249 [NASA ADS] [CrossRef] [Google Scholar]
- Staveley-Smith, L., Davies, R. D., & Kinman, T. D. 1992, MNRAS, 258, 334 [NASA ADS] [CrossRef] [Google Scholar]
- Sung, E.-C., Chun, M.-S., Freeman, K. C., & Chaboyer, B. 2002, in The Dynamics, Structure & History of Galaxies: A Workshop in Honour of Professor Ken Freeman, eds. G. S. Da Costa, E. M. Sadler, & H. Jerjen, ASP Conf. Ser., 273, 341 [Google Scholar]
- Sung, E.-C., Han, C., Ryden, B. S., Chun, M.-S., & Kim, H.-I. 1998, ApJ, 499, 140 [NASA ADS] [CrossRef] [Google Scholar]
- Tajiri, Y. Y., & Kamaya, H. 2002, A&A, 389, 367 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Telles, E., Melnick, J., & Terlevich, R. 1997, MNRAS, 288, 78 [NASA ADS] [CrossRef] [Google Scholar]
- Thornley, M. D., Schreiber, N. M. F., Lutz, D., et al. 2000, ApJ, 539, 641 [NASA ADS] [CrossRef] [Google Scholar]
- Thuan, T. X. 1985, ApJ, 299, 881 [NASA ADS] [CrossRef] [Google Scholar]
- Thuan, T. X. 1991, in Observations and Models of Blue Compact Dwarf Galaxies, eds. C. Leitherer, N. Walborn, T. Heckman, & C. Norman (Cambridge: Cambridge University Press), 183 [Google Scholar]
- Thuan, T. X., & Martin, G. E. 1981, ApJ, 247, 823 [NASA ADS] [CrossRef] [Google Scholar]
- Tully, R. B., & Fisher, J. R. 1977, A&A, 54, 661 [NASA ADS] [Google Scholar]
- Vaduvescu, O., Richer, M. G., & McCall, M. L. 2006, AJ, 131, 1318 [NASA ADS] [CrossRef] [Google Scholar]
- van Zee, L., Westpfahl, D., Haynes, M. P., & Salzer, J. J. 1998, AJ, 115, 1000 [NASA ADS] [CrossRef] [Google Scholar]
- van Zee, L., Salzer, J. J., & Skillman, E. D. 2001, AJ, 122, 121 [NASA ADS] [CrossRef] [Google Scholar]
- Vigroux, L., Souviron, J., & Vader, J. P. 1984, A&A, 139, L9 [NASA ADS] [Google Scholar]
- Vilchez, J. M. 1995, AJ, 110, 1090 [NASA ADS] [CrossRef] [Google Scholar]
- Weinmann, S. M., Lisker, T., Guo, Q., Meyer, H. T., & Janz, J. 2011, MNRAS, 416, 1197 [NASA ADS] [CrossRef] [Google Scholar]
- Yasuda, N., Fukugita, M., Narayanan, V. K., et al. 2001, AJ, 122, 1104 [NASA ADS] [CrossRef] [Google Scholar]
- Zhao, Y., Gao, Y., & Gu, Q. 2013, ApJ, 764, 44 [NASA ADS] [CrossRef] [Google Scholar]
All Tables
Different subtypes of BCDs according to the morphology of the starburst- and LSB-component based on Loose & Thuan (1986b).
Mean values of the derived parameters for the LSB components of the artificial BCDs for different Sérsic indices n.
All Figures
Fig. 1 Distribution of Virgo cluster galaxies. Red symbols refer to objects with BCD as primary class, black symbols to uncertain BCD candidates, small black dots to early-type galaxies, green asterisks to dE(bc)s, and the black cross to M87. Open symbols denote galaxies that could not be analysed (see text for details). |
|
In the text |
Fig. 2 Colour images from SDSS (upper panels), spectra (middle panels), as well as surface brightness and colour profiles (lower panels) of VCC 1744 (left) and VCC 0130 (right). In the case of VCC 1744 a pure exponential law was assumed to fit the LSB-component, while for VCC 0130 an inner profile flattening was assumed (black dashed lines). The horizontal bar denotes the interval in R⋆over which the LSB-component was fitted. The images in this figure and throughout the paper have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. They were obtained from the online Image List Tool of the SDSS (http://skyserver.sdss3.org/dr8/en/tools/chart/list.asp) and were slightly enhanced at low intensity, using the same correction curve for all. |
|
In the text |
Fig. 3 Surface brightness profile of a BCD (VCC 0001). Left panel: the line corresponds to an exponential fit to the LSB-component and the red data points show the outer tail in the r-band. Right panel: same, but for the g, r and i-filters. Obvious are the outer tails of the profiles at low surface brightness levels. The horizontal bars indicate the radius at which the profile starts deviating from a pure exponential fit. |
|
In the text |
Fig. 4 The Mr − (g − i) − Reff − ⟨μ⟩eff-plane of BCDs and ETGs. No decomposition into starburst and LSB-component was applied. Red squares represent the primary BCD sample, blue squares denote the candidate sample. Black dots represent the sample of early-type galaxies taken from Janz & Lisker (2008, 2009), for which colours were derived from the light within 2Reff. Green asterisks are dE(bc)s from Lisker et al. (2006b). For each (vertical) interval in magnitude or surface brightness, the (horizontal) ± 2σ-region of the ETGs is indicated by the black bars. VCC 0428 and 1313 have colours ((g − i)tot = −0.049 and −0.263mag, respectively), lying outside the colour range of the plot and are therefore only indicated by an arrow. |
|
In the text |
Fig. 5 Left: colour−magnitude diagram of the LSB-components of our primary BCD sample, compared to the Virgo early-type dwarf galaxies (labelled “dE” but including dS0s), for which colours were derived within 2Reff. Red squares correspond to values of the LSB-components, where black crosses show LSB-components with inner profile flattening. The notation “LSB(corr)” indicates that the parameters of the LSB-components are corrected for the outer tail (see Sect. 3.5). The black vectors display the change of the parameters when an inner profile flattening is applied. Green asterisks correspond to dE(bc)s and blue squares to galaxies with a uncertain morphological classification (see Sect. 4.4). Typical errors are shown on the left-hand side. Right: same diagram, but additionally plotted are the irregulars (cyan triangles) and the entire population of early-type galaxies (black dots). In both panels we indicate the 2σ-region of the entire ETG population with vertical bars. |
|
In the text |
Fig. 6 Effective radius Reff vs. absolute magnitude (top panels) and mean effective surface brightness (bottom panels). Symbols are the same as in Fig. 5. The black error bars in the left-hand panels illustrate the statistical errors, while the arrows indicate the systematic error, as derived from analysing artificial BCD images (see the Appendix). The error bars are placed on the ordinate so that they correspond to the average surface brightness and the magnitude of the respective set of artificial BCDs. The arrows point in the direction towards which the measurements would have to be shifted to correct for the systematic error. Note that the upper arrow is almost invisible due to the very small error. |
|
In the text |
Fig. 7 Absolute magnitude Mr vs. mean effective surface brightness ⟨μ⟩eff,r in r-band. Symbols are the same as in Fig. 5. The black error bars in the left panel illustrate the statistical errors, while the arrows indicate the systematic error, as derived from analysing artificial BCD images (see the Appendix). The error bars are placed on the abscissa so that they correspond to the magnitude of the respective set of artificial BCDs. The arrows point in the direction towards which the measurements would have to be shifted to correct for the systematic error. Note that the left arrow is almost invisible due to its very low value. |
|
In the text |
Fig. 8 Distributions of the axis ratios of Virgo cluster dwarf galaxies. The upper panels show the irregulars and the LSB-components of the BCDs. The early-type dwarfs (labelled “dE”) also include the dE(bc)s, which are additionally plotted in the lower panel. The black vertical lines indicate the mean values. |
|
In the text |
Fig. 9 Image (left) and unsharp mask (right) of galaxy VCC 0213. Both panels have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. The unsharp mask was created using a two pixel Gaussian kernel. |
|
In the text |
Fig. 10 GOLDMine distance vs. the velocity of the Virgo BCDs. The vertical solid line at vLG = 1900 km s-1 corresponds to the limit from which galaxies belong to the M and W cloud according to Gavazzi et al. (1999). The three BCDs in the box lie in the projected cloud region, but are below or only slightly above this velocity limit. We therefore adopt a distance of 16.5 Mpc for them instead of the GOLDMine distance of 32 Mpc. |
|
In the text |
Fig. 11 Same as Figs. 6 and 7, but instead of a constant distance of 16.5 Mpc, the distances given by the GOLDMine database are used for the LSB-components of the primary sample, except for the three objects (shown with magenta symbols) marked in Fig. 10, which remain at 16.5 Mpc. All other galaxy types also remain at a constant distance of 16.5 Mpc. |
|
In the text |
Fig. 12 Parameter space of the LSB-components of our Virgo BCDs and a comparison sample of star-forming dwarf galaxies (SF-dwarfs) in various environments. The black squares correspond to the comparison sample of star-forming dwarf galaxies, while the other galaxy types are colour coded as in the previous figures. |
|
In the text |
Fig. 13 Possible descendants of BCDs (left) and of irregulars (right) in the Virgo cluster. The left panels show a bright (top) and a faint (bottom) early-type dwarf that fall into the structural parameter range of the LSB-components of Virgo BCDs, i.e. they are more compact than the average early-type dwarf. The right panels show a bright (top) and faint (bottom) early-type dwarf that fall into the range of the Virgo irregulars, i.e. they are more diffuse than the average early-type dwarf. All images have a width of 80′′, corresponding to 6.3 kpc at a distance of 16.5 Mpc. The galaxies and their parameters are, clockwise from top left: VCC 0033 (Mr = −16.85 mag, Reff = 0.73 kpc), VCC 0510 (Mr = −17.00 mag, Reff = 1.64 kpc), VCC 1565 (Mr = −15.49 mag, Reff = 1.40 kpc), and VCC 0068 (Mr = −15.25 mag, Reff = 0.53 kpc). |
|
In the text |
Fig. 14 The two BCDs (left: VCC 0410, right: VCC 1313) that remain more compact than the early-type dwarfs even after applying the individual distances of GOLDMine. These objects lie outside of the 2σ-region of early-type dwarfs in the ⟨μeff⟩ − Reff diagram of Fig. 11. Both images have a width of 80”, corresponding to 6.3 kpc at a distance of 16.5 Mpc. |
|
In the text |
Fig. 15 Histogram of morphological types in the outer region of Virgo (DM87 > 0.6Mpc). A cutoff of Mr > −16.5mag and ⟨μr⟩eff,r < 23.0mag/arcsec2 was applied to restrict the parameter region to the locus of most of the LSB-components (see Fig. 7), using a distance of 16.5 Mpc for all galaxies. For the red histogram lines, an additional colour cutoff of (g − i) > 0.7 mag was applied. The histogram of the early-type dwarfs (labelled “dE”) also includes the dE(bc)s. Candidate BCDs are included with their respective first-mentioned class (e.g. “Sm/BCD”). The numbers are obtained from the full sample of Virgo cluster galaxies that was used in Weinmann et al. (2011), based on the VCC and adopting a faint magnitude limit of 18 mag in the B-band. For the galaxies that could not be analysed with the SDSS, we use estimations for Mr as described in the appendix of Weinmann et al. (2011), and adopt the average surface brightness and colour of the respective class. |
|
In the text |
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.