A&A 398, 89-100 (2003)
DOI: 10.1051/0004-6361:20021631
J. Schulz - U. Fritze - v. Alvensleben - K. J. Fricke
Universitäts-Sternwarte, Geismarlandstr. 11, 37083 Göttingen, Germany
Received 18 March 2002 / Accepted 5 November 2002
Abstract
HST has opened the possibility to decompose the surface brightness profiles of
galaxies up to significant redshifts and look-back times into
r1/4-bulge and exponential disk components. This should allow to study the
redshift evolution of bulge and disk luminosity contributions and discriminate
between the different formation scenarios for these galaxy components
currently discussed, i.e. decide if star formation in bulges and disks started
at the same time or was delayed in either of the two components. An
indispensable prerequisite for the comparison of bulge-to-disk ratios of
galaxies at different redshifts is to properly account for cosmological band
shift and evolutionary effects.
We present evolutionary synthesis models for both components and add their
spectra in various proportions to obtain the full range of local galaxies'
B-band bulge-to-total light ratios. Bulge star formation is assumed to occur on
a short timescale of 109 yr, disk star formation proceeds at a constant rate.
We study the evolution of the relative light contributions of both
components backward in time and, for a given cosmological model, as a function
of redshift. This allows us to see how far back into the past the locally
well-established correlation between galaxy morphologies and spectral properties
can hold. To cope with the present uncertainty about the formation epochs of
bulge and disk components we present models for three scenarios: bulges and
disks of equal age, old bulges and delayed disk star formation, and old disks
with subsequent bulge star formation.
We quantitatively show the wavelength dependence of bulge-to-total (=B/T)
light ratios for local galaxies. The different star formation timescales for
bulge and disk components lead to B/T ratios that significantly increase from Uthrough I-bands (by factors 4-6 for weak bulge systems Sc) with the
rate of increase slightly depending on the relative ages of the two components.
The redshift evolution of B/T-ratios in various bands U, B, V, I, H is
calculated
accounting both for cosmological and evolutionary corrections assuming a
standard cosmology (
). In
particular, for the two scenarios with old bulges and old or younger disks, the
redshift evolution of B/T-ratios is dramatic in every band and both for galaxies
ending up at
with low and high B-band B/T light ratios. Our
results clearly show that it does not make any sense to compare B/T ratios
measured in one and the same band for galaxies at different redshifts without
fully accounting for evolutionary and cosmological effects. These,
unfortunately, significantly depend on the relative ages of the two components
and, hence, on the galaxy formation scenario adopted.
We also show that simultaneous decomposition of galaxy profiles in several bands
can give direct information about these relative ages and constrain formation
scenarios for the different galaxy components.
Of the wavelength bands we explore (U, B, V, I, H), the
I- and H-bands show the smoothest redshift evolution and, hence, are best
suited for a first order comparison of galaxies over the redshift
range from z=0 to .
Our robust result that -
irrespective
of the respective ages of the bulge and disk stellar components -
I-band B/T-ratios apparently increase with increasing redshift for all
galaxy types with present B/T > 0.1 implies that the scarcity of
bulge-strong systems at
reported by Marleau & Simard (1998) and Aguerri & Trujillo (2002) for HDF and Hawaiian
Deep Field galaxies is further enhanced.
Key words: galaxies: general - galaxies: evolution - galaxies:
formation - galaxies: fundamental parameters -
galaxies: stellar content
The bulge-to-disk ratio was one of the main qualitative classification criteria that led Hubble (1926) to establish his famous Hubble sequence of morphological galaxy types. When de Vaucouleurs (1958) showed that surface brightness profiles of M31 and other spirals could be decomposed into spheroidal components following an r1/4 law and exponential disk components, a quantitative decomposition became possible. Simien & de Vaucouleurs (1986) present bulge-to-disk light ratios for 98 galaxies and investigate their relation to other galaxy parameters, as e.g. the Hubble type.
Conventionally, Bulge/Disk (B/D) or Bulge/Total (B/T) light ratios for local galaxies are measured in the B-band. Different Hubble types of galaxies are mainly characterised by their respective B/T ratios.
Deep HST observations open up the possibility to determine structural parameters
of galaxies to considerable redshifts and to directly assess the interesting
question how galaxy morphologies have evolved cosmologically (cf.
Marleau & Simard 1998; Simard et al. 2002).
Despite considerable effort since many years both from the observational side and from semianalytical and hydrocosmological galaxy formation models, it is not yet clear, how and when bulges formed. They may either have formed together with the halo stars and globular clusters prior to the disks or else later from the disks by internal or external dynamical instabilities such as bar formation and destruction or an interaction or accretion event.
Among the parameters easiest to access in observational investigations of galaxy morphologies and their evolution at the limits of present observational capabilities are the B/D and B/T light ratios. To compare B/D or B/T light ratios among galaxies at different redshifts in an attempt to study their evolution it is important to know the wavelength dependence of these light ratios. From the observational side, a very first systematic comparison of galaxy morphologies obtained from B- and H-band classifications of different spiral types is presented by Eskridge et al. (2002).
In the local Universe, a good correlation of galaxy morphological properties with their respective spectral and photometric properties is well established.
Spheroidal galaxies and the spheroidal components of spirals, i.e. the bulges,
are well described with star formation laws that declined rapidly from high
values in the distant past to the very low or zero values observed at present.
Evolutionary synthesis models usually use a time evolution of their Star
Formation Rate (SFR) of the form
with short e-folding times of order
Gyr (e.g.
Bruzual & Charlot 1993; Guiderdoni & Rocca-Volmerange
1987; Fritze - v. Alvensleben & Gerhard 1994, and many others) or
even describe those objects by single burst populations
(Vazdekis 1999;
Maraston & Thomas 2000).
Star formation in galaxy disks, on the other hand, is a much more steady process
with the SFR density depending on the surface density of HI according to a
Schmidt law
with n close to 1 (cf.
Kennicutt 1998a, 1998b). In evolutionary
synthesis models, therefore, the global SFRs for spiral galaxies are
conventionally ascribed much longer characteristic timescales ranging from
Gyr for Sa to
for pure disk Sd
galaxies.
In the work presented here, we use a very simple and standard evolutionary
synthesis approach to investigate the wavelength dependence of B/T light ratios
in various wavelength bands in their time evolution. We separately model the two
galaxy components, bulge and disk, and then combine them to yield, after 12 Gyr
of evolution, the presently observed average B-band B/T light ratios typical of
different spiral types. We can then study the differences between conventional
B-band B/T ratios and the corresponding light ratios
in
other wavelength bands
.
We find significant differences in (B/T) ratios between the B-band where (B/T) ratios are conventionally
determined for local galaxies and other wavelength bands. We show
quantitatively how these differences depend on galaxy type.
Studying our models backward in time, we explore the time and redshift evolution of B/T-ratios in various bands which is important to take into account when comparing galaxy morphologies at various redshifts.
Finally, we discuss the implications that our results have for the interpretation of the first observational results on the redshift evolution of B/T-ratios.
We describe the bulge component by a SFR
with
Gyr and the disk component with a constant SFR
const. We use our evolutionary synthesis code (cf. Schulz et al. 2002 for details) to
separately calculate the time evolution of spectra, luminosities in various
bands (U...K), mass-to-light ratios, etc. for both components and then, in a
second step, we combine the two components in mass ratios as to yield, by
today's age, the observed average B-band B/T ratios of the respective galaxy
types (Table 1).
Type | (B/T)|B |
E | ![]() |
S0 | 0.57 |
Sa | 0.41 |
Sb | 0.24 |
Sc | 0.094 |
Sd | 0.022 |
We note that both components are modeled in the conventional way, i.e. with
stellar evolutionary tracks/isochrones for solar metallicity. Observational
evidence points to average subsolar metallicities for both bulge and disk stars
(cf. McWilliam & Rich 1994;
Sadler et al. 1996;
Ramírez et al. 2000;
Rocha-Pinto & Maciel 1996, 1998). At
about half-solar metallicity, a value typical
for bulge stars and stars in the solar neighborhood, both components bulge and
disk will become slightly bluer (0.05) in their optical color
evolution.
Currently, there are basically two scenarios for the formation of bulges and
disks. The first is the early collapse scenario with bulge Star Formation ( SF) occurring on a short timescale 1 Gyr from pre-galactic gas early in
the Universe, possibly related to halo star and globular cluster formation.
Disks form later, perhaps after a time delay of a few Gyr from dissipative
infall of left-over halo gas. Elliptical galaxies and the bulges of spirals obey
the same relations e.g. for surface brightness vs. luminosity
(Kormendy 1985), magnesium index vs. central velocity dispersion
(e.g. Bender et al. 1993;
Jablonka et al. 1996), suggesting a
common origin for all spheroidal stellar systems.
Large bulges like those of the Milky Way and M31 seem to host old stellar populations with narrow age dispersions (Peletier 2001; Abraham et al. 1999) in support of this early collapse scenario.
In the second scenario, disks form first and bulges are created secularly from disk material. The onset of bulge formation may have various causes. Internal instabilities in the disk as well as external perturbations may create bars that efficiently funnel stars and gas into the inner parts, providing material for intense SF (e.g. Martin & Friedli 1997; Aguerri 1999). A high central mass concentration can, in turn, dissolve the bar and leave a bulge or, in a further step, increase the mass of a pre-existing bulge (Norman et al. 1996). These processes may happen at any time (and even several times) in the evolution of a spiral galaxy. In this scenario, part of the bulge stars may only be spatially regrouped and have an extended age distribution like that of the disk stars, while others form in the bulge formation event. Clearly, accretion events where a larger disk galaxy swallows a smaller companion can also lead to bulge formation. In this case, again, either a stellar object (dE, dSph) may be accreted and make up the bulge or a gas-rich dwarf galaxy that may fuel a central starburst producing the stars in the bulge of the later remnant. While a mere regrouping of pre-existing disk stars would give the bulge stars a broad age distribution, the formation of a significant population of new stars in the bulge-forming event will occur on short timescales (108-109 yr) and leave a bulge with a stellar population younger than that of the disk but with a narrow age distribution.
Some recent observational results seem to favor the view that while large bulges were formed according to the early collapse scenario, i.e. on a short timescale in the early Universe with the disks around them starting SF later, small bulges may form later than the disks during their secular evolution and by internal instabilities or external accretion (Balcells & Peletier 1994; Prugniel et al. 2001; see also Wyse 1999 for a recent review).
For the very abundant numbers of S0 galaxies in local galaxy clusters, yet a different formation scenario is discussed. Cluster S0 galaxies are found to be most plausibly transformed into what they are from the spiral-rich field galaxy population falling into the cluster potential. Spectroscopic analyses by Jones et al. (2000), our own comparison of their photometric properties with evolutionary synthesis models including SF truncation with or without prior starbursts (Bicker et al. 2002) and dynamical modelling (Moore et al. 1998; Quilis et al. 2000) together draft a scenario where rapid interactions with the cluster potential and environment largely destroy the stellar disks of infalling spirals and the dense and hot intracluster medium sweeps away their HI supply, ultimately leaving remnants with S0 morphologies and SF completely truncated. This, however, is not the subject of the present paper where we want to explore the effects of bulge versus disk SF along the Hubble sequence.
To cope with the uncertainty about the evolutionary sequence of bulge and disk
formation, we investigate 3 different scenarios. In the first one, both
components are assumed to start SF at the same time, i.e. 12 Gyr ago. In
the second one, only the bulge starts forming stars that early, at a redshift
-5 in standard cosmologies (we use
)
and SF in the disk is assumed to only start with some
delay about 7 Gyr ago, i.e. at
.
In the third scenario, the disk
is assumed to start SF first and bulge formation is assumed to take place in a
secondary SF event about 7 Gyr ago, like a central starburst as e.g. following a
disk instability, an external perturbation of the disk, or an accretion event.
We just chose one redshift
for the onset of delayed bulge or disk SF and postpone a detailed study of this parameter to a future paper that directly compares models with observational B/T-decompositions.
Note that our models do not include any dynamical effects, they are simple 1-zone models without spatial resolution. They only go beyond classical global galaxy evolutionary synthesis models by the inclusion of two phases of SF for the disk and the bulge which are superposed at appropriate times.
In Fig. 1 we present the spectral evolution of a galaxy with present-day (B/T)|B = 0.6, classified as S0, in the three different scenarios (Cols. 1 to 3) at three different ages ranging from 6.5 Gyr (1st row) through 9 and 12 Gyr (2nd and 3rd row, respectively). Each panel presents the integrated spectrum of the galaxy from the UV through NIR (solid line) and the decomposition into bulge (dashed line) and disk components (dotted line).
The absolute luminosities of the galaxies at 12 Gyr are the same in all three scenarios. The earlier the bulge forms, the higher its early SF rate has to be to lead to a fixed (B/T)|B = 0.6 at 12 Gyr, the later disk SF starts, the higher its constant rate has to be.
![]() |
Figure 1: Spectra of a galaxy with present-day (B/T)|B = 0.6 at ages of 6.5, 9, and 12 Gyr, respectively. The first column shows a model with equal age bulge and disk components, the second column the model with an old bulge and a younger disk, and the third column the model with old disk and a younger bulge. |
Open with DEXTER |
The normalization to B/T = 0.6 in the B-band at 12 Gyr is clearly seen
in terms of nearly comparable fluxes of the bulge and disk components at
Å.
For the equal age bulge and disk scenario presented in the 1st column we see that at early times 6.5 Gyr, the bulge with its strong early SF
strongly dominates the spectrum at optical and NIR wavelengths. Bulge SF at this
age, however, is essentially finished already and disk SF is largely responsible
for the UV flux. As both components age, the UV contribution of the disk
component forming stars at a constant rate is essentially unchanged since the
number of high mass stars has reached an equilibrium. Towards longer
wavelengths, the flux from the disk slightly increases with age as low mass
stars accumulate. The flux from the bulge component decreases as bulge stars go
on dying out from higher to lower masses.
In the scenario with the old bulge and a disk starting SF only around 5 Gyr (2nd column), the disk contribution at optical and NIR wavelengths at an age of 6.5 Gyr is still fairly small with the bulge dominating very strongly, because of his many small mass stars. Only in the UV again, where the old bulge has not much flux any more at this age, the disk component, although still very small in mass, clearly dominates because of their newly born massive stars. To end up with the same (B/T)|B = 0.6 as in the first scenario, the disk now has to have a higher constant SFR and is seen to catch up in optical and NIR fluxes at ages 9 and 12 Gyr like in the scenario shown in Col. 1.
In the 3rd Col. showing the scenario of an old disk and delayed
bulge SF, the situation is very different. At 6.5 Gyr we see a very young bulge
shortly after its violent SF phase that largely dominates the flux from UV
through NIR. Note the high overall luminosity of a galaxy with such a young
bulge component. Around 5000 Å e.g. the total flux in this scenario is
higher by a factor 3 as compared to both other scenarios. At 9 Gyr, the total
flux at 5000 Å is only higher than in the other scenarios by a factor 1.3, bulge SF has ceased and the UV has become dominated by the disk component
which also gains relative importance at longer wavelengths, a trend that
continues through the age of 12 Gyr.
![]() |
Figure 2: V-band luminosity (upper panels) and (V-I) color evolution (lower panels) for model galaxies with present-day (B/T)|B = 0.6. The left column again shows the model with equal ages of bulge and disk, the middle column shows the model with an old bulge and a younger disk, and the right column the model with an old disk and a younger bulge. |
Open with DEXTER |
In this section we present the photometric evolution of galaxies with
present-day
(B/T)|B = 0.6 in the three different scenarios. In
Fig. 2 we show as an example the V-band luminosity and (V-I) color
evolution for the total system (solid lines) and the decomposition into bulge
and disk components (dotted and dashed lines, respectively). Evolution is shown
from the onset of SF through the present time, corresponding to 12 Gyr,
and beyond. Hence, we zoom into the V-band wavelength region of
Fig. 1 to study it in some more detail. Of course, this can be
done for all other wavelengths from UV through NIR, but we will present only one
example here.
We see in all the top panels of Fig. 2 the strong V-band luminosity
evolution associated with the burst-like bulge SF (
)
producing a stellar population that after strong fading over
12 or 7 Gyr, for the old and younger bulge scenarios, respectively, still contributes 60% of the total light in B today. The
burst-like bulge SF leads to a strong increase in V-luminosity shortly after the
onset of bulge SF and a strong decline starting
1 Gyr later. The
continuous disk SF mode, on the other hand, goes along with a weak but steady
rise in disk V-luminosity from the onset of disk SF to the present.
For the total V-luminosity evolution it hardly makes any difference if disk SF
starts as early as bulge SF or 5 Gyr later. In both the old bulge/younger disk
and equal age cases, the V-luminosity is much higher (about a factor 4) in early
evolutionary stages (3 Gyr) than today. Only in the old disk/younger
bulge case the system starts out at significantly lower V-band luminosity than
today until the onset of the bulge SF sharply raises its V-luminosity by a factor >10. About
2 Gyr later, it starts fading again.
Looking at the (V-I) color evolution (lower panels) we see, in the case
of equal age bulge and disk components, a significant reddening by 0.5 mag
during the first 4 Gyr and a constant
for all times
thereafter. The bulge component further reddens at t > 4 Gyr, but this
is compensated by the increasing importance of the bluer disk.
For the old bulge and younger disk case, the total color evolution, of course,
only shows the strong reddening of the bulge component to
until, at
Gyr, disk SF sets in and acts to bring the (V-I) back to 0.85 within
2 Gyr from whereon it remains
constant.
In the old disk and younger bulge case, as expected, the system reddens much
less and only reaches
after 5 Gyr of evolution. The
onset of the strong burst-like bulge SF then induces a strong and sudden blueing
to
which then turns into a steady reddening to reach
about 3 Gyr later.
Differences in today's total (V-I) color for a galaxy with present-day (B/T)|B = 0.6 are small between the different scenarios. At 12 Gyr, (V-I) = 0.87 in the case of equal age components, and 0.85 and 0.83 for older bulge and older disk scenarios, respectively.
![]() |
Figure 3: Present-day B/T-ratios in different filters for model galaxies with (B/T)|B = 0.1, 0.3, and 0.6, respectively. Top panel: equal age components, middle panel: old bulge and younger disk, bottom panel: old disk and younger bulge. |
Open with DEXTER |
In this section, we investigate the wavelength dependence of B/T ratios for
galaxies of various types, i.e. with various present-day B-band (B/T)|Bratios of 0.1 (
Sc), 0.3 (
Sab), and 0.6
(
S0). In Fig. 3, we present the wavelength
dependence of
values for these three galaxy types in
all three scenarios. All galaxies here are depicted at their present age of 12 Gyr. Note that we use HST filter bands for UBVI throughout.
For all galaxies, there is a clear trend for B/T ratios to increase with
increasing wavelength of observation. There is little difference between the
scenarios with equal age bulge and disk components and with old bulges and
younger disks. E.g., an Sc-type galaxy with B-band
(B/T)|B = 0.1 shows
(B/T)|U = 0.03,
(B/T)|V = 0.13, and
(B/T)|I = 0.17. Its
B/T ratio hence increases by a factor 6 from U through I in both scenarios
with an old bulge.
Type | (B/T)|U | (B/T)|B | (B/T)|V | (B/T)|I |
equal age bulge and disk | ||||
E | 0.72 | 0.9 | 0.93 | 0.94 |
S0 | 0.30 | 0.6 | 0.68 | 0.72 |
Sa | 0.11 | 0.3 | 0.38 | 0.42 |
Sc | 0.03 | 0.1 | 0.14 | 0.16 |
old disk | ||||
E | 0.76 | 0.9 | 0.92 | 0.93 |
S0 | 0.35 | 0.6 | 0.66 | 0.69 |
Sa | 0.13 | 0.3 | 0.35 | 0.39 |
Sc | 0.04 | 0.1 | 0.12 | 0.14 |
old bulge | ||||
E | 0.70 | 0.9 | 0.93 | 0.94 |
S0 | 0.28 | 0.6 | 0.69 | 0.74 |
Sa | 0.10 | 0.3 | 0.39 | 0.45 |
Sc | 0.03 | 0.1 | 0.14 | 0.17 |
Table 2 shows that the trend for B/T to increase with increasing wavelength is stronger for the small bulge galaxies than for those with strong bulge components.
For the old disk and younger bulge scenario, however, the trend of increasing
B/T with increasing wavelength is weaker. From
(B/T)|U = 0.04 through
(B/T)|I = 0.15 it only increases by a factor 4.
All of this is due to the different SF histories and timescales of the disk and bulge components. Disk SF being essentially constant in time from its onset to the present epoch soon produces a stellar population in which the fraction of high mass stars dominating the wavelength region shortward of V keeps constant and the low mass stars dominating longward of V slowly accumulate. Bulge SF being strongly peaked around the time of bulge formation implies that after a very blue high luminosity period bulge stars become collectively old and tend to dominate the longer wavelength regions. Their effect, of course, is the stronger, the higher the final required B/T ratio of the galaxy.
We conclude that as a result of the different SF timescales for bulge and disk components the B/T light ratios of galaxies are predicted to be significantly dependent on the wavelength at which the decomposition is performed. Slightly dependent on the galaxy formation scenario, they systematically increase from U through I with the strongest effect, a factor of 4-6 difference from U to I, seen for the weak bulge galaxy types.
In principle, bulge - disk decomposition of local galaxies at different wavelengths may give information about the relative ages of the bulge and disk stars and, hence, help constrain bulge and disk formation scenarios.
Eskridge et al. (2002) recently presented a first systematic comparison of galaxy classifications in B and H. They decompose their sample of 250 bright spirals of various types from S0 through Sm from the Ohio State University Bright Spiral Galaxy Survey in B- and H-bands independently. They find that, while for extremely bulge-strong (S0) as well as pure disk systems (Sd...Sm) the classifications from B and H do not significantly diverge, spirals from Sa through Scd appear one T-type earlier in H than in B, on average, albeit with large scatter. We see from Fig. 3 and Table 2 that in I our Sab and Sc models models have B/T-ratios very similar to the B-band B/T of Sa and Sbc models, respectively, not much dependent on the scenario adopted and well consistent with Eskridge et al.'s results. As in their observations, our S0 model with
(B/T)|B=0.6 does not change much when observed in I:
.
The large observational scatter precludes us from definitely discriminating between the three different scenarios for bulge and disk formation, although the observations seem to better agree with old bulge models than with the one where the disk is older than the bulge. Their results may, as well, indicate that there is in fact a range of ages for real disks (and bulges) rather than just the two extremes that we study. A detailed quantitative analysis of these data with our models, however, will be the subject of a forthcoming paper.
Comparing galaxy morphologies at different redshifts is a difficult task. Two effects have to be considered, the band-shift effect, related to the cosmological correction, resulting from the redshifting and dimming of the spectrum, and the evolutionary effect due to changes in the SF rate.
As to the first effect, it is immediately clear from a comparison of images of local galaxies at different wavelengths, e.g. in the UV, optical, and NIR, that galaxy morphologies appear significantly different at different wavelengths. E.g., star-forming regions or knots stand out in the UV. To describe quantitatively how galaxy morphologies change due to this band shift effect when going towards higher redshift, however, is very difficult.
To quantify evolutionary effects on the morphological appearance of an individual galaxy as a function of redshift would require to know its individual SF history and, in particular, the SF histories of its components bulge and disk. While the SF histories of individual galaxies may to some extent be different from galaxy to galaxy, in particular in their short term fluctuations, there seem to be reasonable average global SF histories, smoothed out over short term fluctuations, that are characteristic of galaxies of various spectral types. This is a very basic fact on which all evolutionary synthesis models rely. And, in the local Universe, at least, a good correlation is observed between spectral and morphological galaxy types. How far this correlation holds back into the past is an open question and depends on the formation scenario for the bulge and disk components of galaxies.
Careful analysis of observations of galaxy morphologies out to significant
redshifts, that are becoming possible now with HST and large ground-based
telescopes, are a promising way to investigate these questions. Marleau & Simard (1998)
present B/T ratios from an r1/4-bulge and exponential disk
decomposition in the I-band (F814W) of more than 500 galaxies brighter than
I = 26 mag in the HDF North over the redshift range
.
They find the number of bulge dominated galaxies to decrease with
increasing redshift. The apparent lack of pure disk galaxies at
in their data may be due to a selection effect driving the low surface
brightness disk galaxies below the detection/classification limit at those
redshifts. The lack of high surface brightness strong bulge components at redshifts
,
however,
cannot be explained by a selection bias. A detailed analysis of the relative
contributions of bulge and disk components, derived from I-band images over this
large redshift interval requires consideration of both band-shift and
evolutionary effects.
As simplified as they are, our bulge and disk models offer a first tool for
comparison with results from bulge - disk decompositions of distant
galaxies. In this Section, we present a theoretical analysis of the redshift
evolution of B/T ratios as obtained from our models for the three different
galaxy formation scenarios we consider. As for the wavelength dependence of
present-day B/T ratios discussed above, the redshift evolution of
ratios, as well, is a result of the different SF timescales
that apply for spheroidal and disk galaxy components, respectively. The results
we present fully account for both the cosmological and evolutionary corrections
as explained above with the cosmological model
.
We stress that UBVI note observer-frame (not
rest-frame) passbands throughout.
Let us look first at the scenario where both bulge and disk components were
forming early in the Universe at .
Here the situation is dramatic for the progenitors of both types of galaxies,
those ending up today with small and large bulge contributions. In Fig. 4, we present the redshift evolution of B/T ratios at
U, B, V, I, H out to z > 4 for two types of galaxies ending
up today at z = 0 with low and high (B/T)
Sc) and 0.6
,
respectively.
At present, no H-band decompositions into bulge and disk light contributions of distant galaxies have been published. As HST NICMOS observations, however, will soon provide a wealth of spatially resolved data (Balcells 2002, priv. comm.), we include our model H-band predictions in Figs. 4-6 and discuss them in Sect. 6.4.
For the late-type galaxies we first recover at z = 0 the difference of a
factor 2 between its B- and I-band B/T ratios. This difference increases
to a factor >20 out to
due to both a decrease in (B/T)|B and an increase in (B/T)|I. At
,
corresponding to ages
2 Gyr in our cosmology, the active bulge SF phase
acts to rapidly raise the B/T ratios in all bands to a uniform value of 0.5 from
whereon they synchronously continue to increase towards
.
![]() |
Figure 4: Redshift evolution of B/T-ratios in various bands for two different galaxy types with present-day (B/T)|B = 0.1 (top) and (B/T)|B = 0.6 (bottom), respectively. Equal age bulge and disk components. |
Open with DEXTER |
The progenitors of present-day low B/T galaxies, e.g. with
show systematically higher (B/T)|I ratios
at higher redshift, ranging from (B/T)
0.2-0.3 at
to
(B/T)
-0.85 at
-4.5. Hence, including evolutionary and cosmological effects, we find the I-band B/T-ratios of present-day late-type galaxies to rise strongly with redshift.
For bulge-dominated galaxies with (B/T)|B = 0.6 at z = 0, the
situation is even more extreme in this scenario. Out to ,
their
B-band B/T ratios decrease strongly down to (B/T)
at
as a result of the Balmer discontinuity being redshifted towards longer passbands. Their I-band B/T ratios are not affected and slowly increase.
For strongly bulge-dominated galaxies uniformly classified in I over the huge
redshift range from
,
in this scenario, no
significant biases are expected except that bulge contributions are slightly
overestimated at z > 1.
In this scenario, galaxies observed at
should show vastly
different B/T ratios in different bands, varying between
for Sc progenitors and between
for S0 progenitors.
Galaxies at z > 2 are predicted to show the same B/T ratios in all bands
U - I and only high values B/T > 0.5 are expected to be observed at
z > 2 since even low mass bulges in their early evolutionary phases (2 Gyr) strongly dominate the light at all wavelengths. Galaxies observed
at high redshift
2 with B/T ratios of 0.5 and greater may well be the
progenitors of local Sc galaxies (
(B/T)|B = 0.1). They simply are
observed in or shortly after their active bulge SF phase.
In Fig. 5, we present the redshift evolution of B/T ratios in
various bands U, B, V, I, H for two different galaxy types with present-day B-band
(B/T)|B (z=0) = 0.1 (
)
and 0.6 (
), respectively, for the scenario with
bulge formation in the early Universe and disk formation delayed until
.
In this scenario, the evolution of low B/T systems is very weak and smooth until
.
At
,
the disk did not form stars,
hence only pure bulges are expected. With increasing redshift, the difference
between B- and I-band B/T ratios increases steadily from a factor
2 at
z = 0 to a factor > 10 at
.
For high B/T systems, we again observe with increasing redshift a decrease of the U- and B-band B/T-ratios like in the previous scenario and a weak increase of the V- and I-band B/T-values.
With I-band classification over the redshift range
,
a
galaxy ending up with
(B/T)|B = 0.1 at
would be observed to have
at
and
around
.
Strong differences are seen between B/T ratios in different bands for late galaxy
types. Within a given band, however, the redshift evolution is moderate to
.
A galaxy ending up with a strong bulge
(B/T)|B = 0.6 at
would be observed in the I-band to even have a somewhat stronger one at
.
Hence, in this scenario the progenitors of local weak bulge systems, when
classified in I at redshifts
look like having slightly stronger bulge components and this effect tremendously increases, of course, as
one approaches the redshift where disk SF slowly started.
In this scenario (Fig. 6) at redshifts higher than
,
bulge SF did not
occur and only pure disks are expected. At z = 0, we recover the
results from Sect. 6.1 showing the strong wavelength dependence of
the B/T ratio. Next, we note the tremendous differences in the redshift
evolution between the U-band B/T-ratios and those in longer wavelength
bands. For strongly bulge-dominated systems, the B/T-ratios in B, V, and Ievolve quite homogeneously in contrast to the previous scenario. They increase
slowly from
(B/T)|B = 0.6 at z = 0 to
at
in all bands except for U.
For galaxies with weak bulge contributions, e.g.
(B/T)|B = 0.1 at z = 0, the B/T-ratios increase considerably in all bands B, V, I to
,
from whereon their slope has a turning point and is beginning
to decrease to
.
Around
,
the curves for
the B/T in different bands cross each other, i.e. at
the B/Tvalues at shorter wavelengths become higher than those at longer ones in
contrast to the situation at lower redshift.
![]() |
Figure 5: Redshift evolution of B/T-ratios in various bands for two different galaxy types with present-day (B/T)|B = 0.1 and 0.6, respectively. Old bulge and younger disk. |
Open with DEXTER |
![]() |
Figure 6: Redshift evolution of B/T-ratios in various bands for two different galaxy types with present-day (B/T)|B = 0.1 and 0.6, respectively. Old disk and younger bulge. |
Open with DEXTER |
With I-band bulge-disk decomposition, the B/T ratio of a galaxy at
would be overestimated as compared to local I-band B/T ratios by typically
a factor
2 for low B/T systems and less so for higher B/T galaxies, except for the U-band.
In the I-band, the B/T ratio of a local Sc-type galaxy with
(B/T)|B (z=0) = 0.1 increases from
(B/T)|I=0.14 at z=0 with increasing redshift to
(B/T)|I = 0.2 at
and
(B/T)|I = 0.5 at
.
This is due to evolutionary and cosmological effects. At
,
I-band observations sample rest-frame B where the bulge component at
an age of the disk of
5 Gyr is actively forming stars in this scenario.
Hence, the intrinsically relatively low mass bulge components of later spiral types are much more prominent in early evolutionary phases closer to their active SF phase and cosmological corrections further enhance this evolutionary effect.
Systems with
at
-0.9 are in fact
the progenitors of local late-type Sc galaxies with
(B/T)|B = 0.1 at z = 0.
Comparing in Figs. 4-6 the redshift evolution of Johnson's (Bessell & Brett 1988) H-band B/T-ratios with those for the I-band discussed above, we see that both are fairly similar with the H-band evolution being slightly more continuous than that in I. The H-band hence appears slightly better suited for a comparison among galaxies over a range of redshifts than I. Note, however, that already at z = 0 strong bulge systems with
(B/T)|B = 0.6 show
and weak bulge systems (
Sc) with
(B/T)|B = 0.1 show
.
While for systems that are bulge-strong at z = 0 the bulge components of their progenitors at higher redshift are slightly overestimated this effect is much stronger for the progenitors of local weak-bulge systems. H-band decomposition of the progenitors of local Sc galaxies are predicted by our models to yield B/T-light ratios of
(B/T)|H = 0.4 around
in case of equal age bulge and disk or for old disk and younger bulge components and as high as
(B/T)|H = 0.7 for the old bulge and younger disk scenario. For the equally old bulge and disk scenario the H-band B/T-ratio of local weak-bulge (
Sc) progenitor galaxies smoothly increases further at redshifts beyond 0.8 and levels off around
(B/T)|H = 0.8 at
.
A quantitative comparison of bulge and disk decompositions over a range of redshifts - unless restricted to I- or H-band investigations in a first approximation - requires evolutionary as well as cosmological corrections that can be given by our models.
For unresolved galaxies evolutionary and cosmological corrections are conventionally given in terms of magnitude differences, which are normalized to the average
luminosity of locally observed galaxies of the respective spectral types. Bulge and disk subcomponents, however,
do not have such an absolute magnitude scale. We therefore chose to give their
evolutionary and cosmological corrections in terms of luminosity ratios rather
than magnitude differences - like in the case for single burst stellar populations (cf. Schulz et al. 2002) - and call them
and
,
respectively, not to be confused with the
and
conventionally given for galaxies. With the luminosities
for the bulge and disk components chosen as to match
and
the
and
for bulge and disk, respectively, are defined as
z |
![]() |
![]() |
![]() |
![]() |
![]() |
... |
0.03 | 1.0828 | 0.8461 | 1.0511 | 0.9286 | 1.0404 | |
0.05 | 1.1416 | 0.7082 | 1.0826 | 0.8462 | 1.0638 | |
0.08 | 1.2474 | 0.5911 | 1.1300 | 0.7597 | 1.0962 | |
![]() |
z |
![]() |
![]() |
![]() |
![]() |
![]() |
... |
0.03 | 1.1321 | 0.8562 | 1.0909 | 0.9346 | 1.0794 | |
0.05 | 1.2118 | 0.7265 | 1.1494 | 0.8599 | 1.1321 | |
0.08 | 1.3444 | 0.6151 | 1.2568 | 0.7816 | 1.2323 | |
![]() |
z |
![]() |
![]() |
![]() |
![]() |
![]() |
... |
0.03 | 0.9973 | 0.9608 | 0.9901 | 0.9629 | 0.9850 | |
0.05 | 0.9965 | 0.9272 | 0.9856 | 0.9169 | 0.9780 | |
0.08 | 0.9954 | 0.9003 | 0.9791 | 0.8704 | 0.9671 | |
![]() |
z |
![]() |
![]() |
![]() |
![]() |
![]() |
... |
0.03 | 0.9945 | 0.9649 | 0.9824 | 0.9667 | 0.9741 | |
0.05 | 0.9925 | 0.9350 | 0.9741 | 0.9247 | 0.9615 | |
0.08 | 0.9897 | 0.9114 | 0.9621 | 0.8827 | 0.9425 | |
![]() |
In Table 3 we give the evolutionary and cosmological corrections
and
as a function of redshift z for broad band filters
,
separately for our old and younger bulge and old and younger disk models in a cosmology with H0=65,
.
z | t [Gyr] | BDM |
0.0 | 11.92 | 31.15 |
0.5 | 6.957 | 42.267 |
1.0 | 4.535 | 44.130 |
1.5 | 3.115 | 45.305 |
2.0 | 2.189 | 46.179 |
2.5 | 1.542 | 46.877 |
3.0 | 1.065 | 47.461 |
3.5 | 0.702 | 47.961 |
4.0 | 0.417 | 48.400 |
4.5 | 0.188 | 48.790 |
5.0 | 0.0 | 49.141 |
The longstanding controversy about the origin of the different galaxy components bulge and disk, in particular about their respective formation epochs and mechanisms, may soon be solved. High resolution imaging of faint sources, e.g. with HST, opens up the possibility of surface brightness decomposition into r1/4- bulge and exponential disk components for galaxies up to considerable redshifts z > 1. Applied to reasonably large samples of galaxies with observational biases carefully taken into account, this technique will allow to directly study the redshift evolution of both components to significant look-back times. To derive the redshift evolution of bulge and disk luminosity contributions and study how far back in time the locally well established correlation between galaxy morphological and spectrophotometric properties will hold, however, careful consideration of both cosmological band-shift and evolutionary effects are of crucial importance.
We present a first and simplified evolutionary synthesis approach to analyze these two effects. Bulge star formation - like that for any kind of spheroidal system - is conventionally assumed to occur on short timescales (108-109 yr) while the star formation rate in the disk component is safely assumed to be roughly constant in time.
We separately model the spectral evolution from UV through NIR of both
components and superpose them in various proportions as to match by an age of 12 Gyr the average bulge-to-total light ratios as conventionally measured in the B-band for different galaxy types (S0, Sa, ..., Sd). To cope with the full
range of currently discussed formations scenarios for both galaxy components, we
present models in which both bulge and disk star formation is assumed to start
around the same time in the early Universe (at )
and others
in which either of the two is somehow delayed until
.
We can then study in each of the three scenarios the wavelength dependence of present-day bulge-to-total light ratios that results from the different star formation timescales (and possibly ages) of the two components. While it is intuitively clear that galaxy morphologies depend on the wavelength of observation, we can show quantitatively the importance of this effect for the decomposition of surface brightness profiles. It is strongest for weak bulge systems and only slightly dependent on the adopted scenario. For an Sc-type galaxy with B-band B/T = 0.1 e.g., the B/T ratio increases by a factor 4 (old disk and younger bulge) to 6 (old bulge and old or younger disk) from U- through I-band. A very first comparison shows recent observational results for a sample of galaxies decomposed both in B- and H-bands by Eskridge et al. (2002) to be well consistent with our model predictions for I-band decomposition. The large scatter in the observed (B/T)|H vs. (B/T)|B relation, however, precludes us from discriminating between our three different scenarios.
In a second step, we study for a standard cosmological model (
)
the redshift evolution of bulge-to-total
light ratios and find it to be significantly different for the three different
scenarios. The progenitors of local galaxies of a given type, i.e., of a given
B-band B/T light ratio at z = 0, are shown to feature strongly different
B/T ratios when observed over a range in redshift in the same band. The redshift
evolution of B/T ratios in a given band depends in a complex way on galaxy type,
i.e. present-day B/T ratio, on the wavelength band in which the profile
decomposition is performed, and on the formation epochs of the two components.
In particular, we find that in the scenario with old bulge and equally old disk
the progenitors of local weak bulge Sc-type galaxies are expected to show
significantly stronger bulge components at all wavelengths at .
This is also significant in the scenario with old disk and younger bulge at
.
In case of a younger bulge its strong and relatively recent star formation
makes even low present-day B/T systems appear bulge-dominated at .
The I- and H-bands generally show the smoothest redshift evolution in terms of
B/T light ratios and hence are the best wavelengths to compare B/T-ratios over some range in redshifts. Nevertheless, it is important to take into account that the progenitors of any type of local galaxy with
are expected to show stronger bulge contributions at higher redshifts.
So, even in the I- and H-bands, and much more so in other wavelength bands, B/T ratios
of galaxies at different redshifts can only be compared if evolutionary and
cosmological corrections are properly accounted for.
Irrespective of the respective ages of the bulge and disk stellar
components, our models indicate that both for low and high present-day
B/T-galaxies, their I-band B/T-ratios rise with increasing
redshift. I.e., I-band B/T-ratios of galaxies at z > 0 should
overestimate their bulge light contributions as compared to local
galaxies. Marleau & Simard (1998), however, report a scarcity of bulge-strong systems at
from their I-band decomposition of HDF-N galaxy profiles. In view of our results the scarcity of bulge-strong systems at high redshift should even be more dramatic than what they find. The same holds true for the results recently reported by Aguerri & Trujillo (2002) from an I-band B/T-decomposition of galaxies from two Hawaiian Deep Fields. The
decrease in the number density of E/S0 galaxies to z = 0.8 they find should also be seen as a lower limit to their real decrease in view of our results. We recall that while surface brightness dimming may cause a selection bias against disks at high redshift this is not expected for the high surface brightness bulge components.
We show that surface brightness profile decomposition in more than one band and for a sample of galaxies at different redshifts (see Simard et al. 2002 for a first data set) should allow to discriminate between the different scenarios for bulge and disk formation. It will also be useful to eliminate the effects of dust reddening and absorption shown to be very important for the proper bulge-disk decomposition of local galaxies by Balcells & Peletier (1994).
We give electronic tables for the luminosity evolution of bulge and disk subcomponents as well as for their respective evolutionary and cosmological corrections for a standard cosmological model in all scenarios considered. They allow for a detailed quantitative comparison of galaxy B/T ratios in various bands over a wide redshift range.
Clearly, more observations, refined models, and, in particular, profile decompositions at more than one wavelength are required to definitively settle this important issue which has far-reaching consequences for our understanding of the formation and cosmological evolution of galaxies.
We believe that as simplified as our first set of models is at the present stage, they are a useful and indispensable tool for the interpretation of galaxy morphologies at z > 0 and their redshift evolution. Detailed comparison with existing HST data (e.g. Marleau & Simard 1998; Aguerri & Trujillo 2002; Simard et al. 2002) should allow to seriously constrain galaxy formation scenarios.
We caution that all our conclusions result from the different star formation timescales (and possibly different redshifts at which star formation sets in) in the bulge and disk components. They are therefore not able to describe a situation - if it exists - in which a bulge is formed without significant additional star formation by a mere regrouping of disk stars due to dynamical instabilities.
Acknowledgements
We thank our referee, Marc Balcells, for his valuable suggestions that helped to significantly improve the presentation of our results.