A&A 373, 190-198 (2001)
DOI: 10.1051/0004-6361:20010585
R. Behrend - A. Maeder
Geneva Observatory,
1290 Sauverny, Switzerland
Received 3 April 2001 / Accepted 12 April 2001
Abstract
We perform calculations of pre-main sequence evolution
of stars from 1 to
with growing accretion rates
.
The values of
are taken equal to a constant fraction
of the rates of
the mass outflows observed by Churchwell (1998) and Henning
(2000). The evolution of the various stellar parameters is
given, as well as the evolution of the disc luminosity;
electronic tables are provided as a supplement to the articles. Typically, the
duration of the accretion phase of massive stars is
and there is less than
difference in the time necessary to form a 8 or
star.
If in a young cluster all the proto-stellar cores start to
accrete at the same time, we then have a relation M(t) between the masses of the new stars and the time t of their
appearance. Since we also know the distribution of stellar masses at the end of star
formation (IMF), we can derive the star formation history
N(t). Interestingly enough, the current IMF implies two
peaks of star formation: low mass stars form first and high mass
star form later.
Key words: stars: formation - stars: evolution - stars: pre-main sequence - stars: statistics - accretion, accretion discs - Hertzsprung-Russel (HR) and C-M diagrams
Two scenarios have been proposed to explain the formation of massive stars.
Star coalescence (Bonnell et al. 1998; Stahler et al. 1998)
and the accretion scenario (Beech & Mitalas 1994; Bernasconi & Maeder 1996).
The merging of stars requires very high stellar densities
(10
)
to be efficient (Bonnell et al.
1998; Henning 2001).
The imaging of new-born stars with the HST (Burrows et al. 1996), the IRAM interferometer (Dutrey 1998), the VLA (Wilner & Ho & Kastner & Rodríguez 2000) and the VLT (Brandner 2000, among others) clearly shows discs (or remnants of discs), favouring the accretion scenario for low and intermediate mass stars. The jets and/or outflows that are sometime observed are closely related to accretion discs; they can contribute to solve the angular momentum problem (Tomisaka 2000).
Accretion rates of the order of 10-5 to
explain very well low mass stars (Palla 1998). But such values fail to describe the formation
of massive stars, because with such rates the massive stars would leave the ZAMS
(i.e. would have burned a significant part of their hydrogen) before being fully formed,
as pointed out by Nakano (1988) and Norberg & Maeder (2000).
Nakano (2000) studied the implications of a huge constant accretion rate of
10
for an observed
protostar.
The reality could well lie between these two models. Can an accretion rate that depends on the stellar mass (or luminosity) describe the star formation from low to very high masses? This is the question of the growing accretion rates.
A nice observational relation between the outflow mass rates and the
stellar bolometric luminosities in ultra compact
regions was established by
Churchwell
(1998) over the ranges 1 to
and 10-6 to
.
Independent observations by Henning et al.
(2000) confirm its validity.
Churchwell's relation can bring some new constraints for the scenarios of the formation of massive stars. A crucial question is to know if the accretion scenario is compatible or not with this relation. The relation can bring something new to the description of the pre-main sequence evolution of stars and for the explanation of the rapid birth of massive stars. The relation can also be a used as a strong check for the computer simulations of discs around stars and the jets/outflows, and offers the possibility to fix uncertain parameters of these models (e.g. viscosity).
The principal properties of stars with accretion discs on the birthline are discussed. Tracks for stars leaving the birthline are calculated. Links between the IMF and the processes that end the stellar mass evolution are also presented.
![]() |
Figure 1:
The
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
The matter which escapes from the disc (in the central region) is divided in two major flows:
one part falls on the star while the other part is ejected in the form of jets and outflows.
Let us denote by
the mass rate transferred inwards from the disc, by
the mass rate accreted by the star, and by
the mass rate in the jets/outflow.
The relation, found by Churchwell (1998) and confirmed by Henning et al.
(2000), between the bolometric luminosity L of
stars and the outflow mass rate
is single valued
(at least in the error bars), monotonously increasing and valid over a very wide
range of luminosity: from 1 to
.
A polynomial fit of the observations by Churchwell (1998) gives
Due to the great number of changing physical parameters which can have an effect on f, f is probably non-constant during the evolution of accreting stars. The last three values of f cited above suggest that f could be a decreasing function of the stellar mass, but this is far from being firmly established. For this preliminary study, we simplify this huge astrophysical problem, retaining only a ratio f that is constant during the whole accretion period.
The "best'' f is presently unknown. A way to assign a more or less realistic value to f is as follows.
The birthline is the path in the HR diagram that continuously accreting stars follow.
Stars on the birthline are difficult to observe,
because they are hidden by the accretion disc and a dense cocoon of interstellar matter. The star becomes
visible once a sufficient amount of the surrounding matter has been dispersed.
It follows that the young stars are observed generally between the birthline and the zero
age main-sequence (ZAMS). Before hydrogen ignition, gravitational
contraction and deuterium burning are the main energy sources for the star; the
localization of the birthline depends on the deuterium accretion rate. We will choose as the value for which the
computed birthline fits "at the best'' the upper envelope of the observations of young stars in
the pre-main sequence part of the HR diagram (as in Norberg & Maeder 2000, denoted hereafter N&M).
The evolution of the pre-main sequence models with growing accretion rate
is computed with the Geneva code (Meynet & Maeder 2000), under
the assumption of a non-rotating star, at Z=0.02 and at deuterium concentration
.
Winds are not taken into account.
The computations begin with a fully convective
protostar which is taken
to be
old, in agreement with models of low mass star
formation (Palla 1998). The mass of the protostar grows at a rate
.
We compare the models with accretion rates parameterized by the bolometric luminosity of
the star to the model of mass parameterized accretion rate found by N&M:
where units are
and
.
The M- and L-parameterized accretion rates are different in the sense that the
first one depends only on one parameter (the mass) of the star, and the second one depends
on many involved quantities (mass, profiles of chemical composition and internal energy,
etc.) and thus is sensitive to the history of the star (e.g. opacity via chemical abundances,
diffusive processes). For this reason, the memory of the initial luminosity of the star in
the M-parameterization is rapidly lost. This is not the case for the L-parameterized
accretion; if the initial star is over-/under-luminous, the luminosity of the accreted
deuterium is also above/below the normal value and the birthline is also
over-/under-luminous. This sensitivity to the initial conditions for the L-parameterized
accretion remains small with respect to the possible variations of the birthlines with f,
for stars initially on the Hayashi band.
Processes that depend on their past can "easily'' oscillate or even be chaotic, for physical
reasons (Hairer et al. 1987) or for numerical reasons (Bergé et al. 1988).
Apparently the calculations do not reveal strong signatures of instabilities, except maybe
some oscillations near the surface of accreting stars above
.
![]() |
Figure 2:
A birthline in
![]() ![]() ![]() ![]() |
Open with DEXTER |
The first step consists of the estimate of a plausible .
The second step is the prediction of physical parameters for the stars and discs using
the
found. Connections with the IMF is postponed to the next section.
Before hydrogen ignition, the main energy sources for the star are gravitational
contraction and deuterium burning. The localization of the birthline is then governed
essentially by the deuterium accretion rate on the star. The higher ,
the higher the
birthline in the HR diagram.
![]() |
Figure 1 shows birthlines for different values of ,
and the
observations compiled in N&M.
While
seems to be a little too low,
give reasonable
envelopes to the observations.
Two statistical biases in the estimate of the envelope are possible, with consequences
for the determination of
.
1) The stars which constrain the birthline are usually all on one side of the birthline (a
nonsignificant exception is nevertheless visible in Fig. 8).
Because the density of stars near the birthline is small, the probability that one misses a
crucial observation is relatively large. The effect of this bias is an underestimated
luminosity of the envelope of the observation and thus also of the birthline. Unfortunately, this occurs in
the upper pre-main sequence part of the HR diagram where the sensibility of the computed
birthline to
is the greatest. The deficiency of stars in this region is due to the
very short time spent by a star in this region, as will be shown in Fig. 8.
2) The stars leave the birthline faster than they reach the ZAMS. In other words, the median of the
positions of the stars for which the error bars constrain the birthline (by their envelope)
is shifted towards the ZAMS; this effect vanishes as error bars go to zero.
With these two facts in mind, our preference is for the highest value:
,
i.e. f=1/3.
The ratio
supported by Fig. 1 is in good agreement
with theoretical estimates of
,
as mentioned above.
Properties of stars on the birthline and the relative luminosity of the disc to that of the star,
for
,
are summarized in Table 1, while Fig. 2
locates different phases of the star evolution in the HR diagram.
![]() |
Figure 3:
Location of the convective and radiative zones. The relative radius of the
boundaries between the convective and radiative zones are drawn for the birthline with
![]() |
Open with DEXTER |
The -burning begins at the center of the star. Nascent stars of mass below
(
)
are fully convective (Fig. 3),
and the deuterium accreted is carried inside the
star where it is burnt. As the temperature rises, the deuterium
burns farther from the center. At
,
a radiative core is created,
and the accreted deuterium is then burnt in a shell at the periphery of this core;
the star inflates (second hump in Fig. 4) and the luminosity grows.
If the relation of the L-parameterized accretion is right, and in particular its single valued
character, the shell
-burning could play a major role in the beginning of star formation,
raising the stellar luminosity and inducing via Churchwell's relation, a real boost
in the accretion rate, shortening the time needed to reach high mass stars. The convective
-burning
shell becomes thin at the surface of the star around
(
)
(slightly depending on the mass accretion rate); the star becomes fully radiative.
![]() |
Figure 4: M - R relation for stars evolving on the birthlines. |
Open with DEXTER |
The luminosity of the -burning is governed by the accretion rate;
the star inflates to balance a excess of power, or contracts during a shortage of deuterium.
The R(M) relation (Fig. 4) is thus very sensitive to
during the pre-main
sequence evolution of the star and especially during the fully radiative phase.
![]() |
Figure 5: Mass evolution for stars on the birthlines. |
Open with DEXTER |
The consequences of the luminosity's growth during the shell -burning are clearly visible
in the
-
diagram (Fig. 5) in the form of an
elbow near
(
)
with a rather abrupt change of slope in
the curves, for the L-parameterized accretion.
Below the bend, the much smoother curve of the M-parameterized accretion of N&M is a little bit higher
than the one for the L-parameterized accretion, and lower afterwards. Rather than physical, the origin
of this difference resides in the choice of the parameterization for
.
Because both laws were adjusted on the same observations, they globally look alike and make similar
predictions, but they can locally differ.
![]() |
Figure 6:
![]() ![]() |
Open with DEXTER |
![]() |
Figure 7:
Relative luminosity of disc to that of the star as a function of stellar bolometric luminosity for
stars evolving on the birthlines. For the model of N&M,
![]() |
Open with DEXTER |
![]() |
Figure 8:
Tracks in
![]() ![]() ![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
One half of the gravitational energy gained during the rotational
infall is transformed into kinetic energy; the other half is used to heat up the disc and is
eventually radiated. Chemical energies (sublimation of dust and grains, ionization) are small
with regards to gravitational
energy and are neglected in the computations. The disc luminosity,
,
is thus of the order of
,
as given by
Hartmann (1998a, Sect. 5.4).
It appears in Fig. 7 that during the accretion phase,
the luminosity of the disc is always greater than the star luminosity.
On the birthline, the star encounters three phases of nearly constant luminosity. In these
phases, the radius of the star shrinks while the mass continues to grow; according to the
previous formula, the luminosity of the disc grows. That explains the three vertical steps
in the
-
diagram. The step
around
corresponds to the arrival of the star on the ZAMS and is
the only feature really sensitive to
.
The probability of having good observations of
stars in this region of the diagram is very small and we conclude that the relative
luminosity of the disc to that of the star cannot help us to fix
.
The smaller variability of the ratio of disc to star luminosities for the
L-parameterized accretion than for the M-parameterized accretion has the same origin as
in the
-
diagram.
After the accreting star has joined the ZAMS, at
,
it evolves in a very similar manner to the one
calculated in N&M. They used
and the linear fit in Churchwell's
data implies
using the L(M) relation on the ZAMS.
An accretion rate of
is obtained for a
star.
Though it seems very high, this rate is still in the permitted domain of accretion with
spherical symmetry (Wolfire & Cassinelli 1987); if the accretion is made via a disc,
even higher rates are permitted due to the less important cross-section of infalling material
to radiations from the star and from the shock, and to the enhanced cooling
(Nakano 1988; Henning 2001).
Once the disc is no longer fed by its environment, for example by the disappearance of the parent cloud by
the outflows of the star or jets/ouflows and the radiation of neighbouring hot luminous stars (Nakano et al. 1995),
the disc get thinner by accretion and evaporation. The transition from active to passive discs
is probably a smooth function of time. As its mass decreases, the residual disc plays a
smaller and smaller role in the evolution of the star.
Both the mass of discs and the typical duration of the transition are still very uncertain;
in the computations, we neglect their effects, using an instantaneous stopping of the
accretion at time
.
Figure 8 shows the resulting theoretical tracks for some values of
(in units of
).
The following considerations can partly justify the assumption that from a certain stage,
the disc has no more significant influence on the central star.
For low and intermediate mass stars, the disc is usually lighter
than the central star (Hartmann 1998a, Sects. 6.3 and 9.7; 1998b; Bernasconi 1997).
Observations of the residual mass transfer rate of the disc around T Tauri stars is
described by Hartmann (1998b); they can be approximately represented by
(units are
and
), where
is a relative scale for the ages.
Integrating
from
the
corresponding to a disc mass transfer rate
at
to
infinity, one finds a crude estimate of the mass of the disc at the moment of the
stopping of the accretion:
(units are
and
).
is
probably out of the validity domain for
,
but
remains less than
of the mass of the star,
for
.
Bernasconi (1997) studied the residual mass accretion of dying discs,
but his models contain many arbitrary parameters; for intermediate mass stars, the difference between
the tracks with and without residual accretion are nonetheless much smaller than the error bars of the
present day observations.
The time
corresponds approximately to the time at which the cocoon
of interstellar matter vanishes, thus the star becomes visible in the optical bands only after this
moment; the location of the first optical visibility of stars above
is very close to the ZAMS.
For the stars studied in this paper (mass on the ZAMS >
),
there are four regions in the pre-main sequence part of the HR diagram, as shown in
Fig. 8. 1) For
below
(i.e.
),
the stars are fully convective; they leave the birthline with a vertical Hayashi-like
displacement in the HR
diagram. No more accreting stars between
and
(i.e.
)
have a radiative core surrounded by a
convective
-burning shell; if the star already has a radiative core before
the end of accretion, the track closely follows the birthline (with an increasing luminosity) for a short period of time after
;
if a radiative core is absent, one is swiftly created after the
end of accretion and the star leaves the birthline with a decreasing luminosity.
3) For
above
,
the stars are fully radiative, except
on 4) a small band above the ZAMS where stars already have a convective
-burning core.
Electronic tables for the tracks are available via the anonymous ftp server of the Geneva observatory; the file is located at ftp://obsftp.unige.ch/pub/evol/prems/AA2001BM.
The initial mass function (IMF) represents the density of probability that a star ends its accretion with a given mass. The IMF results from many physical phenomena (fragmentation of the parent cloud, vanishing of the reservoir of interstellar matter, accretion rate and its efficiency, binarity, etc.). But under the assumption that the accretion rate does not depend on the local conditions of the cloud and disc, the IMF can be seen as the result of the mechanisms that end the feeding of the disc with interstellar matter (in the case of a small ratio of disc to star masses).
The IMF is usually noted .
is the relative number
of new-born stars in the mass interval
centered on m; the IMF is
normalized so that
.
The assumptions made for the models of this paper imply that the star mass evolution
on the birthline is a single valued and monotonously increasing function M(t).
Thus, from the observed IMF and using the M(t) relation (as in Fig. 4),
one can compute the density of probability
of the age
of stars
when the accretion is abruptly stopped:
![]() |
Figure 9:
Density of probability of
![]() ![]() ![]() ![]() ![]() |
Open with DEXTER |
For the L-parameterized accretion, the break in the mass evolution graph located
near
is great enough to generate a sharp hump (an order of magnitude) in
the density of probability of
,
.
For the M-parameterized accretion of N&M which gives a monotonously
decreasing
,
a bump is also
distinguishable, but it is smoothed by the nature of the law used. In
comparison to the L-parameterized accretion rate, the M-parameterized accretion rate
is indeed
a little higher for masses in the range 1.6 to
and lower for the other stars
(Fig. 5).
The presence of the peak for the L-parameterized accretion does not depend on the ratio
of accreted and ejected masses, if this quantity
does not evolve with the mass and luminosity of the star.
If theses peaks are confirmed, it would be most interesting to investigate the various
physical processes at their origin.
One possible way to explain theses curves could be as follows. For stars below
,
the
fractal parent cloud is probably dissipated mostly by neighbouring stars. For stars above
,
the
star and disc luminosities could grow high enough to scatter the parent cloud.
Numerical experiments have been done with mass accretion rates proportional to the
luminosity parameterized outflow rate
Acknowledgements
The authors express their thanks to Dr. G. Meynet and to P. Norberg for their help and fruitful discussions during the use of the Geneva stellar evolution code.