text:
# On the recent star formation history of the Milky Way disk

## 1 Introduction

Throughout the Milky Way’s history, a non-negligible fraction of star formation has apparently occurred in starburst-like events (see Majewski, 1993 or Freeman & Bland-Hawthorn, 2002, for a review). On the other hand, observations suggest that present-day star formation in the disk of our Galaxy takes place in stellar groupings rather than in isolation. These stellar aggregates form from molecular clouds, and in the disk of the Milky Way they appear in two types: bound and unbound. Unbound, short lived stellar groupings are called associations; bound, long lived stellar groups are known as open clusters. Open clusters can also be formed out of the remains of rich stellar associations (Kroupa et al., 2001). As most of the field stars appear to have been formed in the so-called clustered mode (i.e., in clusters or associations), not in the dispersed mode (i.e., in isolation), it is a natural choice to consider these stellar clumpings as the de facto units of star formation in the disk of our Galaxy (Clarke et al., 2000). The idea of open clusters being fundamental units of star formation in the Galactic disk is, however, controversial (Meyer et al., 2000). It has been argued that bound open clusters cannot contribute significantly to the field star population of the Galactic disk because they are rare and long lived (e.g. Roberts, 1957). In contrast, most young embedded clusters are thought to evolve into unbound stellar associations, which comprise the majority of stars that populate the Milky Way disk (Lada & Lada, 1991). Nevertheless, if star clusters are the elementary units of star formation they can, in principle, be used to derive the star formation history, recent and old. Unfortunately, stellar associations evolve and dissolve in a time-scale of $``$ 50 Myr (Brown, 2002), therefore they cannot be used to study the star formation history of the Galactic disk. On the other hand, open clusters are comparatively long lived objects that may serve as excellent probes into the structure and evolution of the Galactic disk.

The problem of deriving the star formation history of the Milky Way has been considered by a number of authors using different techniques (Twarog, 1980; Scalo, 1987; Barry, 1988; Gomez et al., 1990; Marsakov et al., 1990; Noh & Scalo, 1990; Meusinger, 1991; Soderblom et al., 1991; Micela et al., 1993; Díaz-Pinto et al., 1994; Rocha-Pinto & Maciel, 1997; Chereul et al., 1998; Isern et al., 1999; Lachaume et al., 1999; Rocha-Pinto et al., 2000a; Hernandez et al., 2000a; Bertelli & Nasi, 2001; Gizis et al., 2002; Just, 2002, 2003). The majority of these studies suggest that the disk of the Milky Way Galaxy has not experienced a smooth and constant star formation history but a bursty one with several episodes of enhanced star formation. Although most of these papers are restricted to the study of the star formation history in the solar neighbourhood, the star sample birth sites are in fact distributed over a larger range of distances because of orbital diffusion, and so they can provide an estimate of the global star formation rate and their conclusions can be extrapolated to the entire Milky Way disk.

In this paper we revisit the topic of the recent star formation history of the Milky Way disk by using data from the Open Cluster Database (Mermilliod, 2003), the New Catalogue of Optically Visible Open Clusters and Candidates (Dias et al., 2003a, b) as well as Hipparcos data (ESA, 1997) to construct the open cluster age distribution. The age distribution diagram is then interpreted using numerical results from an extensive library of $`N`$-body calculations carried out during the last ten years. This method permits the reconstruction of the star formation history of the Milky Way disk over the last 2 Gyr with a time resolution of 0.05 Gyr using very few a priori assumptions. Our main objective is to understand the recent star formation history of the entire Galactic disk although only the solar circle ($`\pm `$ 1.75 kpc) can be studied due to the incompleteness of the available open cluster sample (see Section 2).

This paper is organized as follows: in Section 2, we present the properties of the Galactic open cluster sample. The raw histogram of open cluster number vs. age is also included in this Section as well as a discussion of the potential limitations of our approach. In Section 3, we summarize some relevant results on the evolution and dissolution of realistic $`N`$-body star cluster models. This Section also includes a comparison between published open cluster ages and number of stars with the results of the numerical models. In Section 4, we present our method as well as results for the Galactic disk. A closer view to the last 0.2 Gyr and the derived open cluster initial mass function are considered in Section 5. The star formation history is presented and discussed in Section 6; corrections are also discussed here. A detailed comparison with results from other authors is carried out in Section 7. An alternative interpretation of our results in the context of dynamical merger histories is included in Section 8. Open questions and conclusions are summarized in Section 9.

## 2 An open cluster sample

For many years, the Lyngå catalogue (1985, 1987) has been the classical reference source for open cluster data. An exhaustive analysis of the properties of the Milky Way open cluster system using this catalogue was carried out by Janes et al. (1988). Unfortunately, the Lyngå catalogue is not the most updated source of data in this field. The Open Cluster Database (WEBDA) created and maintained by J.-C. Mermilliod (1988, 1992a, 1992b, 1993, 1995, 1996, 2003) includes all the data already covered in Lyngå’s catalogue and many more. The latest update of the Open Cluster Database (October 2003) includes 1731 open clusters with ages for 616 objects (36%). In this database we have found 581 clusters with age $``$ 2 Gyr. Slightly less complete is the New Catalogue of Optically Visible Open Clusters and Candidates (NCOVOCC). The first version of this catalogue was published in Dias et al. (2002). The October 2003 version (Dias et al., 2003b), includes 1637 objects, 603 (37%) with published ages. This new catalogue updates the previous catalogues of Lyngå (1987) and of Mermilliod (1995) included in the WEBDA database. There are however non-negligible differences between the two databases as we can see from Figs. 1 and 2. NCOVOCC includes 569 clusters of age $``$ 2 Gyr. On the other hand, Salaris et al. (2004) have recently re-determined ages for a sample of 71 of the oldest open clusters (all of them older than 0.67 Gyr) using a morphological age indicator, metallicities, and Hipparcos parallaxes. In general, the new ages are older than those from WEBDA and NCOVOCC. In order to study the impact of different age determination techniques on our results we have modified the NCOVOCC data with the new age determinations from Salaris et al. (2004). The new subsample consists of 559 open clusters younger than 2 Gyr. The age-distance diagram for this new sample appears in Fig. 3.

In this paper we use the three subsamples presented above to construct an age distribution diagram. As explained below, we choose 2 Gyr in order to obtain statistically relevant results. Our samples contain objects from dense, very young clusters to old and almost dissolved moving groups. The age distribution diagram for the three open cluster samples considered is shown in Fig. 4. Although a detailed analysis of the different features present in this histogram is postponed to Section 4, the figure suggests that several epochs of enhanced star formation may have taken place during the last 2 Gyr.

The age-distance diagrams for the three samples appear in Figs. 1, 2, and 3. These figures indicate that the open cluster data used in this paper provide an incomplete sample; unfortunately, it is difficult to estimate the degree of incompleteness of this sample due to the irregular distribution of Galactic open clusters. The sample appears to be rather incomplete for open clusters more distant than $``$ 3.5 kpc, therefore our conclusions can only be rigorously applied to the solar circle: an annulus of Galactocentric distance $`R_{}\pm 1.75`$ kpc, where $`R_{}`$ is the solar Galactocentric distance ($``$ 8.5 kpc, e.g., Zombeck, 1990). On the other hand, the Hipparcos Catalogue contains numerous selection biases which are difficult to characterize. It is, however, reasonable to assume that the sample is fairly complete (166 clusters) for objects in the solar neighbourhood (distances to the Sun of 1 kpc or less) although selection effects are likely to be very important for clusters located in the outskirts of the Galactic disk as well as in the region behind the Galactic center and for small clusters with Galactocentric distance $`<8.5`$ kpc as they are projected against the dense stellar fields of the southern Milky Way. On the other hand, intermediate age and old open cluster samples are incomplete as clusters with higher $`|z|`$ are easier to detect than clusters with lower $`|z|`$. In this paper, we are interested in studying only the recent star formation history of the Milky Way ($``$ 2 Gyr), therefore any effects arising from lack of completeness of the sample are likely to be less important than those from age errors. Besides, open cluster samples are biased toward objects where the total luminosity is dominated by ultra-luminous O and B stars (young clusters) or red (or blue) giants (old clusters) as these systems are easier to identify. However, if O and B stars were the only stars the cluster would not be stable since the massive stars are destined to become supernovae, returning their mass to the interstellar medium. Hence, the existence of low mass stars is critical for the survival of open clusters. Old, dynamically depleted open clusters are very difficult to identify although they may exist in large numbers (de la Fuente Marcos, 1998). Faint, population depleted open clusters are very difficult to identify against the field stars. Therefore we are working with an incomplete, volume-limited sample that is also flux-limited.

On the other hand, and as for any other physical parameter, open cluster age determinations are affected by errors. In a large sample, these errors are very likely to be non-homogeneous as different methods have been used by different authors to calculate the ages. The disparity in open cluster ages (when considering different authors, see WEBDA for multiple examples) partly reflects the critical dependence of cluster isochrone fitting on the adopted reddening (often large for open clusters), even small reddening errors can create significant errors in the derived age and metallicity. However, determining effective temperatures and metallicities for cluster turn off stars directly through echelle spectroscopy is free from systematic errors and subject only to uncertainties in the model atmospheres. Some ages for clusters in the samples considered have been determined using the first technique, but some others have been found using the second one or even other indirect methods. It is relatively difficult, specially in papers older than about 5 years, to find published estimations of the errors associated with open cluster age determinations. In order to evaluate the impact of age errors in our analysis we have compiled errors for 52 objects within the studied age range from published sources: NGC 6134 and NGC 3680 (Bruntt et al., 1999); NGC 2141 (Carraro et al., 2001); NGC 2158 (Grocholski & Sarajedini, 2002); NGC 2112 (Carraro et al., 2002); Berkeley 104, Berkeley 60, King 15, NGC 381, Berkeley 64, King 6, NGC 1348, Berkeley 23, NGC 2259, NGC 2304 (Ann et al., 2002); NGC 1663 (Baume et al., 2003); NGC 3990 (Prisinzano et al., 2004); Pismis 19 (Carraro & Munari, 2004); IC 166, NGC 752, King 5, NGC 1254, NGC 1278, NGC 1817, NGC 2158, NGC 2194, NGC 2192, NGC 2236, NGC 2266, Berkeley 30, NGC 2324, NGC 2354, NGC 2355, NGC 2360, Haffner 6, Melotte 71, Pismis 2, NGC 2660, NGC 2849, NGC 3680, NGC 4815, NGC 5822, IC 4651, IC 4756, Berkeley 42, NGC 6802, NGC 6827, NGC 7044, NGC 2477, NGC 7789, NGC 2204, Hyades, and Praesepe (Salaris et al., 2004). Fig. 5 shows the uncertainty as a function of the age for this open cluster sample. The average age error is about 22% and the plot suggests that younger ages are affected by larger errors (for a recent discussion on the problems associated with age determinations in young open clusters see, e.g., Piskunov et al., 2004). If this error sample can be considered as representative of the error range for the entire sample, the youngest cluster (age $`<`$ 0.5 Gyr) errors are very likely in the range 50-150 Myr with for older cluster errors in the range 150-250 Myr.

In order to test the reliability of the features found in the open cluster age distribution we have generated synthetic age distribution diagrams from the original data. Each simulation was composed of the following steps:

1) The original sample from Mermilliod (2003) is the input

2) A randomly distributed error within a given percentage of the input value is added/subtracted to each value in the original sample.

3) The simulated open cluster sample is binned at 0.050 Gyr intervals. In this way the real age of the open clusters in the sample is shifted randomly according to a given average error (15, 25, and 35 %).

4) The process is repeated 175,000 times per error value and we calculate the average number of objects per age bin. The final, averaged, sample is plotted in Fig. 6.

We have carried out three sets with 175,000 simulations to study this effect. These simulations appear to suggest that for age errors of about 15% or lower, most of the features in the age distribution remain relatively unaltered. For errors in age determination larger than about 25% only the youngest features are recovered.

## 3 The life span of open cluster models

Fig. 4 is nothing more than the actual age distribution, a raw histogram of open cluster numbers vs. age, of relatively young, nearby open clusters. In order to interpret this distribution additional tools are needed. Experimental Stellar Dynamics may help us in this analysis by matching cluster dynamics and cluster data to provide a better view of the star formation history of the Milky Way disk. This section summarizes some relevant results on the evolution and dissolution of realistic $`N`$-body star cluster models. Some technical details on the modeling process are also discussed as well as the initial mass functions used in the calculations. After an open cluster is formed and its massive stars have expelled the gas it will evolve steadily, losing stars on a relatively long time-scale. Distant encounters among cluster stars together with the Galactic tidal field cause the slow and gradual evaporation of the cluster. Close encounters and stellar evolution can speed up the process, specially for small, sparsely populated clusters. Primordial binaries are also very important (although less than stellar evolution) as they dominate completely the dynamical evolution of poor clusters. Violent, very fast cluster disruption can be induced during encounters with giant molecular clouds (Terlevich, 1983, 1985, 1987). The disintegration time-scale of open clusters depends on both population and Galactocentric distance; but, for a given Galactocentric distance, the initial population ($`N_o`$) as well as the initial mass function are the main physical parameters to consider in order to analyze the long term evolution of the cluster.

The numerical results discussed here come from an extensive set of simulations (de la Fuente Marcos, 1995, 1996a, 1996b, 1997a, 1997b; de la Fuente Marcos & de la Fuente Marcos, 2000, 2002) calculated with the standard $`N`$-body code NBODY5 (Aarseth, 1985, 1994, 1999, 2003) for clusters located in the solar neighbourhood. These calculations include the effects of stellar evolution, the Galactic tidal field, primordial binaries, and realistic initial mass functions (hereafter IMF). NBODY5 consists of a fourth-order predictor-corrector integration scheme with individual time steps. It utilizes an Ahmad-Cohen (1973) neighbour scheme to facilitate calculation of the gravitational forces, and handles close encounters via two-, three-, four-, and chain regularization techniques (Kustaanheimo & Stiefel, 1965; Aarseth & Zare, 1974; Mikkola, 1985; Mikkola & Aarseth, 1993). Five different IMFs were used to generate stellar masses. Two standard power law models (Salpeter, 1995; Taff, 1974) and three modern IMFs due to the Miller & Scalo (Miller & Scalo, 1979; Eggleton et al., 1989), Kroupa (Kroupa & Tout, 1992; Kroupa et al., 1990, 1991, 1993) and Scalo (Scalo, 1986). The IMFs used in our calculations differ mainly in the number of low-mass and high-mass stars generated. The Salpeter IMF is given reasonably well by the approximation

$$\xi (M)0.03M^{2.35},$$
(1)
where $`\xi (M)\mathrm{d}M`$ is the number of stars in the mass interval $`M`$ to $`M+\mathrm{d}M`$. The Taff IMF is given by

$$\xi (M)0.03M^\alpha ,$$
(2)
with $`\alpha =2.50`$ for clusters with $`N100`$ and $`\alpha =2.65`$ for clusters with $`N>100`$. For generating masses with the Miller & Scalo IMF we use the function (Eggleton et al., 1989)

$$M(\mathrm{X})=\frac{0.19\mathrm{X}}{(1\mathrm{X})^{0.75}+0.032(1\mathrm{X})^{0.25}},$$
(3)
where X is uniformly distributed in the interval $`[0,1]`$. For X $``$ 0.7–0.999, this gives a power-law expression with $`\alpha =7/3`$ and for X $`<<`$ 1 it gives an approximately constant value. For X $`>`$ 0.999, the IMF slope is increased to 5. The Kroupa IMF is given by (Kroupa et al., 1993)

$$\xi (M)=\{\begin{array}{cc}0.035M^{1.3}\hfill & \text{if 0.08 }\text{ M }<\text{ 0.5},\hfill \\ 0.019M^{2.2}\hfill & \text{if 0.5 }\text{ M }<\text{ 1.0},\hfill \\ 0.019M^{2.7}\hfill & \text{if 1.0 }\text{ M }<\text{ }\mathrm{}.\hfill \end{array}$$
(4)
We use the mass-generating function

$$M(\mathrm{X})=0.08+\frac{\gamma _1\mathrm{X}^{\gamma _2}+\gamma _3\mathrm{X}^{\gamma _4}}{(1X)^{0.58}},$$
(5)
as a convenient representation of these relations, where X is uniformly distributed in the interval $`[0,1]`$ and the $`\gamma `$-parameters are given by $`\gamma _1=0.19`$, $`\gamma _2=1.55`$, $`\gamma _3=0.050`$ and $`\gamma _4=0.6`$. No masses less than 0.08 $`M_{}`$ are generated by this formula. It is based on Scalo’s (1986) initial mass function. For the Scalo IMF we use the mass-generating function

$$M(\mathrm{X})=\frac{0.3\mathrm{X}}{(1X)^{0.55}},$$
(6)
where X is uniformly distributed in the interval $`[0,1]`$. It is similar to Kroupa IMF except for $`M<`$ 0.16; Scalo IMF gives a smaller number of stars with those masses. The Kroupa IMF is the most recent IMF we have used, it generates a greater number of low-mass stars and a smaller number of heavy stars. On the contrary, the Salpeter IMF gives a smaller number of low-mass stars and a greater number of massive stars. Among these, the Taff IMF gives a smaller number of massive stars than the Salpeter one but greater than the Miller & Scalo one. The Scalo IMF is intermediate between Miller & Scalo and Kroupa IMFs. For low-mass stars, the Miller & Scalo IMF gives a greater number than the Salpeter and Taff ones but smaller than Scalo and Kroupa IMFs. In addition, three different star densities (uniform, $``$ 1/$`r^2`$, Plummer) have been used to generate the initial positions of the stars in the cluster model. Since the initial membership of a star cluster is the main parameter with regard to cluster life-time for a given Galactocentric distance, stellar evolution (through the IMF) is the only process able to induce a differential behaviour as a function of $`N`$.

Averaging cluster life-times for the entire library of models as a function of $`N`$ to fit a power law to the data we obtain $`\tau =0.011N^{0.68}`$ (in Gyr) with a correlation coefficient of $`r=0.995`$ (Fig, 7). This slope is very similar to the slopes found by Baumgardt (2001) and Baumgardt & Makino (2003) although the details of our simulations are rather different. The above evolutionary track is only strictly valid for cluster models with $`N>`$ 100-150 members. Unfortunately, for smaller clusters the characteristic life-time is not very well defined as the energy released in a single supernova event may be larger that the binding energy of the entire cluster (de la Fuente Marcos, 1993). For clusters with $`N<50`$, life-times are in the range 5-40 Myr and for $`N`$ = 75, $`\tau <`$ 80. In general, the fluctuation range for the life-time of small-$`N`$ cluster models is very large and depends strongly on the initial mass function, the spatial distribution, and the binary fraction. For richer clusters, simulations indicate that, in the solar neighbourhood and neglecting interactions with giant molecular clouds and the effect of the cluster gas on the early stages of the evolution of the cluster, an open cluster has to include about 200-400 stars (at least) in order to survive for about 0.5 Gyr, 400-700 to last 0.7 Gyr, 700-1000 to be detectable after 0.9 Gyr and 1000-2000 to survive for about 1.3 Gyr. If a cluster is still observable after 2 Gyr, its initial population was at least 3000 stars. In terms of cluster masses and for an average stellar mass of 0.4 $`M_{}`$, the oldest open clusters considered in this paper (2 Gyr) were born with masses in excess of 1.2 $`\times 10^3M_{}`$ and very likely around $`2\times 10^3M_{}`$.

The astrophysical interpretation of the solid line in Fig. 7 is rather clear: it is the open cluster death line. Open clusters well below the line are young and they evolve towards the line as they lose stars into the field. As it evolves, an average open cluster spends most of its life-time close to the death line moving parallel to it, in the population vs. age diagram. Open clusters become part of the stellar field when they cross the line. It is, however, very difficult to find open cluster memberships and masses among the published literature. Star counts are in many cases unreliable and proper motions are not always available to confirm membership. There is however an observational parameter that may help in open cluster membership determinations: the integrated magnitude of the cluster.

The integrated absolute visual magnitude of a star cluster can be written as:

$$M_v(\mathrm{cluster})M_v(\mathrm{sun})=2.5\mathrm{log}\frac{L(\mathrm{cluster})}{L(\mathrm{sun})},$$
(7)
where the absolute visual magnitude of the Sun is 4.82 (e.g. Zombeck, 1990) and the luminosity of the cluster could be, in principle, estimated using the mass-luminosity relationship $`LM^\beta `$, with $`\beta 3.8`$ (Popper, 1980) for stars in the main sequence. Unfortunately, the main contribution to the intrinsic brightness of a relatively young cluster is not coming from the stars on the main sequence but from red (or blue) subgiants or giants. Taking this fact into account the luminosity of the cluster can be written as:

$$L(\mathrm{cluster})=N_{NMS}<L_{NMS}>+N_{MS}<L_{MS}>,$$
(8)
where $`N_{MS}`$ and $`N_{NMS}`$ are the number of stars in the main sequence and outside the main sequence, respectively. Their relative contribution to the integrated luminosity of the cluster is obtained multiplying by the corresponding average luminosity. This average luminosity can be obtained from the average mass and the mass-luminosity relationship for the contribution from the main sequence. For subgiants or giants it is also possible to obtain an average luminosity using data for standard stars (e.g. Zombeck, 1990). In our calculations we will consider an average luminosity for stars outside the main sequence of 565 $`L_{}`$. Combining Eqs. 7, 8, $`N=N_{MS}+N_{NMS}`$, and considering the average stellar mass $`M`$ we obtain:

$$N=\frac{1}{M^{3.8}+\frac{N_{NMS}}{N}<L_{NMS}(L_{})>}10^{\frac{(M_v(\mathrm{cluster})M_v(\mathrm{sun}))}{2.5}}.$$
(9)
The fraction of stars outside the main sequence can be estimated assuming an IMF and taking into account that an estimation of the life-time in the main sequence for a given star of mass $`M`$ and luminosity $`L`$ is $`M/L`$. The relationship given by Eq. 9 can be calibrated using a few clusters with well established memberships (e.g. Hyades, Pleiades). For the age range studied in this work the fraction of stars outside the main sequence is in the range 0.5-10%. The final value of $`N`$ is affected by our choice for this fraction but the main contribution to the uncertainty in $`N`$ is coming from the observational determination of the integrated magnitude. For Galactic open clusters this uncertainty could be $`\pm 0.5`$ mag (Battinelli et al., 1994; Lata et al., 2002). This translates into 20% errors in $`N`$. We have applied this approach to estimate the current population of a small sample of open clusters (Hyades, Pleiades, IC 4756, NGC 2204, NGC 2506, NGC 7789, NGC 1245, NGC 381, NGC 2324, and NGC 2192) with age, error in age, and integrated absolute visual magnitude available from published sources (WEBDA; Battinelli et al., 1994; Gray, 1965; Pandey et al., 1989; Lata et al., 2002; Sagar et al., 1983; Spassova & Beav, 1985). The results are also plotted in Fig. 7, the dashed line represents the power law fit to the open cluster data (grey circles), $`\tau =9.2N^{0.70}`$ with a correlation coefficient of $`r=0.941`$. Only one open cluster younger than 500 Myr has been included because all of them are well below the limiting line. The power law fit includes all the plotted points. If the youngest cluster is not considered then the correlation coefficient is essentially 1. In any case, Fig. 7 appears to confirm that observational data and results from $`N`$-body simulations are fully consistent within the error limits.

## 4 Method and results

In this section we provide a detailed description of our open cluster age distribution method. As stated in Section 2, our results are mainly sensitive to the errors in age determination but also to the degree of completeness of the open cluster sample. These two effects have been analyzed in Section 2 and are discussed in detail below. On the other hand, any errors in theoretical stellar models also propagate into the results obtained because the published open cluster ages always make reference to theoretical stellar models, sometimes through direct isochrone comparison, others through the use of morphological features found in the cluster color-magnitude diagram.

### 4.1 Method and assumptions

The open cluster catalogues offer high-quality data for a relatively large number of objects in the Milky Way disk, both in the solar neighbourhood and beyond. In particular, the age data can be used to construct an age distribution for the Milky Way open clusters, Fig. 4. Assuming that the sample considered is representative of the entire Galactic disk, the star formation rate can be derived from its age distribution, as the number of open clusters in each age bin is, by hypothesis, correlated with the number of open clusters initially born at that time. Once an age distribution is available, it is in principle easy to recover the star formation history that gave rise to the observed age distribution by using results from realistic $`N`$-body simulations. This method permits the reconstruction of the global star formation history of the Galactic disk with a time resolution of 0.050 Gyr over the last 2 Gyr. Our method does not assume any a priori structure or condition on the star formation rate and it basically consists of two steps: (i) Construct a representative sample of open clusters. In principle, the optimal approach will be to construct a volume-limited open cluster sample in the solar vicinity (see Section 6). The cluster sample spans a very large range of ages. (ii) Construct the age distribution diagram for the sample. (iii) Infer the star formation history from the diagram. Observations of bursts of star formation in other galaxies indicate that intense star formation is always associated with production of large star clusters. The life span of larger clusters is longer; therefore, an usually high number of open clusters at a given time interval can be interpreted as the result of an event of enhanced star formation at that given age. This is probably the main a priori assumption maintained throughout the paper although it appears to be strongly supported by observational evidence. If star cluster masses are sampled from an open cluster initial mass function, larger numbers translate into increased probability of formation of large star clusters. Luminous, and therefore rich, young star clusters are found whenever there is vigorous star formation, whether it be in galaxy mergers or starburst galaxies. Galaxies with very active star formation have proportionally more of their stars in clusters than in the field, with some of them devoting as much as 15-20% of their luminosity to clusters (Larsen & Richtler, 2000). Presently available observational data strongly suggests that the cluster forming frequency is highest during violent bursts of star formation (van den Bergh, 2000).

### 4.2 Results

The resulting star formation history comes directly from the age distribution (Figs. 4, 8), in an approach which assumes that the most frequent ages of the open clusters indicate the epochs when the star formation was more intense if we take into account that star clusters are the elementary units of the star formation process. We have not included any volume correction because the number of clusters more distant than about 6 kpc is negligible for the two samples considered. The volume effect as well as evolutionary corrections are considered in Section 6). We refer basically to the clumps of clusters at about 0.35, 0.70, 1.13, 1.50, and 1.93 Gyr. These clumps will be identified as burst 1, 2, 3, 4, and 5, respectively. The actual existence of burst 1, 2 and 3 is not significantly (see the detailed discussion below) affected by uncertainties in the cluster ages as errors are smaller than 0.3 Gyr for all clusters (younger than 1.2 Gyr) in the samples considered. However, for clusters older than 1.2 Gyr age errors are larger, therefore the statistical significance of bursts 4 and 5 is not, in principle, as well established as for 1, 2, and 3.

On the other hand, if we assume that the Galactic disk is not currently undergoing a burst of star formation and that the age distribution for stars younger than 0.2 Gyr is representative of the monotonic (quiescent) star formation rate in the Milky Way disk (i.e. it has not been polluted by an enhanced episode of star formation), then the fraction of open clusters that survives for up to about 0.25 Gyr is about 92%. This corresponds to the characteristic life-time of clusters with $`N150`$. In other words, observations indicate that only 8% of all the open clusters are born with stellar populations larger than about 200 stars or masses $`>80M_{}`$. This feature can be interpreted as implying that star formation takes place preferentially in clusters of that size. This result, that most of the stars appear to have been formed in small star clusters, has already been pointed out by Kroupa (1995a, b, c) and de la Fuente Marcos (1997b). However, a single 3,000 stars open cluster can produce as many field stars as 60 small (50 stars) clusters and, as pointed out above, the number of large clusters violently destroyed during cluster-giant molecular cloud interactions, even though difficult to estimate, can be very significant. On the other hand, during an episode of enhanced star formation even much larger clusters (up to 40,000 stars) can be formed. Therefore it may also be possible that most of the field stars have been formed in large clusters, however we do not observe many large clusters today and the previous possibility remains rather speculative. In any case, the cumulative age distribution for the full open cluster sample shows a departure from the predictions of constant formation rate and exponentially declining dissolution rate.

### 4.3 Quiescent vs. burst of star formation

Assuming that the Galactic disk has never been polluted by any extragalactic stellar populations, our previous analysis suggests that superimposed on a relatively small level of quiescent star formation, mainly in small-$`N`$ star clusters, the star formation rate has experienced at least 5 episodes of enhanced star formation lasting about 0.2 Gyr each, with production of larger clusters. This cyclic behavior appears to show a period of 0.4$`\pm `$0.1 Gyr. Although the observed cyclic behavior can be easily explained as triggered by density waves (Hernandez et al., 2000a; Martos et al., 2004), the analysis in this section is also exploring the possibility of periodic, tidally triggered star formation bursts as an alternative (or possibly concurrent) scenario to explain the enhancements observed in the open cluster age distribution. There is evidence for a recent pericentric passage of the Sagittarius Dwarf Spheroidal Galaxy $``$ 50 Myr ago (Johnston et al., 1999). If this close interaction was able to tidally trigger star formation in the disk of the Milky Way, its effects should currently be observable and it could be considered as burst 0. There are in fact observational signatures of unusually rich young open clusters (or even protoglobulars) in the disk of the Milky Way: e.g. Cyg OB2 (Knödlseder, 2000; Knödlseder et al., 2002; Comerón et al., 2002; Hanson, 2003). We will however neglect this possibility in our subsequent discussion.

#### 4.3.1 The 0.35 Gyr burst

The most recent star formation bursts are also the most likely bursts to have occurred, since they took place in the very recent past, and so are less affected by the age errors. Unfortunately, our analysis in Section 2 appears to suggest that age errors are larger for younger clusters. The youngest (0.25-0.45 Gyr) burst found in this work appears around the same age in all the three samples, lasting 0.2 Gyr that is the typical duration of bursts found in previous work. This burst is temporally coincident with a perigalactic passage of the Small Magellanic Cloud (hereafter SMC) with the Milky Way as predicted by Lin et al. (1995). Harris and Zaritsky (2004) have studied the star formation history of the SMC and have found four major bursts of star formation: 0.06, 0.4, 2.5, and 8.4 Gyr. Models by Zaritsky and Harris (2004) indicate that the 0.4 Gyr peak in the SMC stellar age function coincides with one of the epochs of closest approach between SMC and the Milky Way. During the interaction, the elevation in the star formation rate was significant ($`>`$ factor of 2) over the quiescent rate of star formation in the SMC.

#### 4.3.2 The 0.70 Gyr burst

This burst appears at 0.6-0.8 Gyr, lasting also 0.2 Gyr, in WEBDA and NCOVOCC but appears to be wider, 0.6-1.0 Gyr, for the sample including Salaris et al. (2004), shifting the peak to 0.8 Gyr. Johnston et al. (1999) have found that the Sagittarius Dwarf Galaxy had a pericentric passage $``$ 0.7 Gyr ago. This close interaction is temporally coincident with burst 2.

#### 4.3.3 The 1.13 Gyr burst

From WEBDA, this burst started 1.3 Gyr ago, lasting about 350 Myr. However, NCOVOCC gives the same starting age but lasting 0.3 Gyr. The sample corrected with Salaris et al. (2004) data differs also only in the duration, 0.2 Gyr. This burst is temporally coincident with a perigalactic passage of the SMC with the Milky Way (Lin et al., 1995).

#### 4.3.4 The 1.50 Gyr burst

The three samples indicate that the burst started 1.6 Gyr ago, lasting 0.2 Gyr. This burst appears to be coincident with a close encounter between the Magellanic Clouds and the Galaxy (Murai & Fujimoto, 1980; Gardiner et al., 1994; Lin et al., 1995). On the other hand, the Sagittarius Dwarf Galaxy had another pericentric passage around the same time (Johnston, 1998; Johnston et al., 1999).

#### 4.3.5 The 1.93 Gyr burst

This is the oldest star formation burst identified in this work and it is also the most affected by the age errors. WEBDA and NCOVOCC indicate that this burst started 2 Gyr ago, lasting 0.15 Gyr. The third sample suggests a burst lasting 0.1 Gyr. This is the only burst that does not appear to coincide with any predicted or modelled interaction with a neighbour galaxy.

In spite of the time coincidences, interpreting the observed periodicity in the age distribution of young, nearby open clusters as evidence in support of a tidally triggered star formation history is somewhat speculative. On the other hand, a close interaction with the Magellanic Clouds cannot be compared dynamically to that of the Sagittarius or the Canis Majoris dwarf galaxies.

Using our method without any corrections it is, however, difficult to estimate the evolution of the absolute star formation rate with the age of the disk. Nevertheless, it is possible to attempt a rough estimate: the number of clusters in bursts 1, 2, 3, 4, and 5 from the WEBDA sample is 73, 31, 28, 10, and 5 respectively, for the Dias et al. sample is 73, 30, 28, 12, and 4 and for Salaris et al. is 66, 49, 12, 12, and 4. Data seem to suggest that for WEBDA and Dias et al. burst 2 was likely the strongest. If under quiescent star formation only 2% of open clusters survive for more than 2 Gyr, during burst 2 the star formation increased by at least a factor 4. Salaris et al. data also indicate that burst 2 was the most important, followed by burst 4. During burst 2 the number of open clusters formed increased by a factor 8 at least.

### 4.4 Limitations

As pointed out before, the star formation history inferred for the Milky Way disk in the previous section is affected by selection effects and cluster age errors, therefore some features in the derived age distribution could be caused by the incompleteness of the sample. As a result of these limitations, the actual amplitude of the bursts is likely higher. Working with a flux-limited sample makes it very difficult to estimate the evolution of the absolute star formation rate with the age of the disk. Nevertheless, in this section we try to assess the accuracy of our results.

The age errors affect considerably the duration of the star formation events, since they tend to scatter the ages of the star clusters originally born in a burst. We can expect that this error could smear out peaks and fill in gaps in the age distribution. In Section 2, we have tested our results using synthetic age distribution diagrams from the original data (see Fig. 6). Simulations appear to indicate that for age errors of about 15% or lower bursts 2, 3, 4, and 5 remain relatively unaltered both in terms of age and duration, however burst 1 seems to split into two bursts, one at 0.2 Gyr and another one at 0.5 Gyr. For errors in age determination larger than about 25% only bursts at 0.2, 0.5, and 0.7 are recovered. However, our error analysis in Section 2 appears to indicate that the older the cluster, the lower the error. Therefore it could be possible that errors for clusters older than 1.5 Gyr are in the range 15-25% and the five bursts can be recovered.

Our derivation of the star formation history of the Milky Way disk uses a sample of dynamically bound, disk star clusters or open clusters. This sample is assumed as not polluted by unbound star clusters or associations. However, associations dissolve on a time-scale of $``$ 50 Myr so even if some degree of contamination remains, the contribution to the final conclusions is negligible. It is only significant for the first bin in our age distribution diagram. On the other hand, it is known that a typical cluster in the outer disk will survive about twice as long as one in the inner disk (Janes et al., 1988). The outer disk group is an homogeneous group whose members are dissolving on a time-scale on the order of 3-4 Gyr. This contribution is however negligible as we restrict our analysis to clusters younger than 2 Gyr. In our previous analysis we have assumed that this age distribution is depopulated from extragalactic objects, so no open clusters have been incorporated into the galactic disk from accreted galaxies. An alternative, plausible but rather speculative scenario with excess clusters coming from accreted dwarf galaxies is considered in Section 8.

Limitations on the applicability of our method to a wider time range appear in connection with the small number of known open clusters older than 2 Gyr. With a resolution of 0.05 Gyr, our results are not statistically significant for an age older than 2 Gyr as the number of objects per age bin is too low. The need to provide reliable results limits the age range over which we can derive the star formation history to 0-2 Gyr, with the time resolution considered throughout the paper. However, our method can be applied to a wider time range if the time resolution is lowered (de la Fuente Marcos & de la Fuente Marcos, 2004).

## 5 The last 0.2 Gyr

Contrary to what is sometimes assumed in the early stages of the star formation history of the Milky Way disk, which were ruled by accretion events and early mergers of satellite dwarf galaxies and star clusters that contributed to the formation of a substantial fraction of the old and metal-deficient stars in the Galactic disk, during the last 0.2 Gyr the star formation rate (see Fig. 9) has been remarkably stable with apparently no events of enhanced star formation.

Fig. 9 appears to indicate that the open cluster sample younger than 0.2 Gyr can be considered representative in the analysis of quiescent star formation in the Milky Way disk. The characteristic time-scale to assemble and give birth to an open cluster seems to be $``$ 10 Myr with (at least) about 80 new open clusters being formed every 10 Myr or 8 clusters/Myr within about 3 kpc from the Sun. About 80% of newly formed clusters are destroyed during the first 20 Myr of life.

For a time resolution of 5 Myr (Fig. 9, small panel) a very interesting feature is visible: clusters younger than 5 Myr are scarce. In principle, this feature may be interpreted as a minimum in the star formation rate, however young clusters are still embedded in their parent molecular cloud and this makes them difficult to detect. Dust obscuration appears to be the natural explanation for the origin of this feature. On the other hand, there is a sharp decrease in the number of open clusters observed for ages older than 15-20 Myr. There is an initial very rapid rate of cluster dissolution, gradually declining thereafter. Again the most simple interpretation is to assume that an enhanced star formation period started $``$ 20 Myr ago. However, this feature can be better explained if we consider stellar evolution in open clusters. For a metallicity $`Z`$ = 0.02, stellar evolution models (e.g. Schaller et al., 1992; Bressan et al., 1993) indicate that stars with masses above 10 $`M_{}`$ and older than about 20 Myr have already exploded as supernovae. For a small cluster ($`N<`$ 75), the energy released during one supernova event can be equivalent to the binding energy of the entire cluster. Therefore, early supernovae appear to be responsible for the violent destruction of about 80% of young open clusters.

In the following section, we use this information to generate a power-law open cluster initial mass function with slope -2.7.

### 5.1 Open Cluster Initial Mass Function

The frequency distribution of stellar masses at birth, the so-called stellar initial mass function (hereafter SIMF), is a fundamental parameter to study the stellar mass spectrum. Salpeter (1955) used the observed luminosity function for the solar neighbourhood and theoretical evolutionary times to derive a SIMF which may be approximated by a power-law:

$$n(m)m^\alpha ,$$
(10)
where $`n(m)`$ is the number of stars per unit mass interval. The original value of $`\alpha `$ found by Salpeter is 2.35, for masses between 0.4 and 10.0 $`M_{}`$. The determination of the SIMF is a cornerstone in Astrophysics, since the SIMF determines the evolution, surface brightness, chemical enrichment, and baryonic content of all galaxies and it is also a fundamental link between stellar and galactic evolution. For recent reviews on the SIMF see Kroupa (2002) and Chabrier (2003).

As many of the stars in the disk of the Milky Way appear to be the result of a clustered mode of star formation, the open cluster initial mass function (hereafter OCIMF) is also a fundamentally important distribution function to study cluster formation, chemical evolution of galaxies, and star formation in general. The fact that the SIMF and the OCIMF can be intimately related was first pointed out by Reddish (1978). However, deriving the OCIMF is very difficult as cluster mass determinations are affected by several observational and theoretical biases: obscuration by dust, flux limits, incomplete sampling of cluster haloes, lack of information on the evolutionary state of the cluster, and uncertainties in stellar models and cluster dissolution rates. Although the OCIMF is not well constrained it can be written as a power law of the type

$$N(M)d(M)M^\alpha dM$$
(11)
for $`M_{min}<M<M_{max}`$, $`N(M)`$ is the number of open clusters per unit mass interval. The cumulative distribution will be

$$f(M)M^{\alpha +1}.$$
(12)
If the process by which open clusters are formed is also able to produce individual stars, it appears reasonable that both SIMF and OCIMF should show similar slopes within the observational errors. Reddish (1978) found a slope equal to -2.2 from a sample of 72 open clusters. Vázquez and Feinstein (1989) estimated the OCIMF for the Milky Way disk using a sample of 130 clusters from Lyngå’s (1987) catalogue. In order to find the distribution, they used the values of the masses quoted in the catalogue to obtain a slope $`\alpha =2.74\pm 0.09`$ for all the clusters in their sample. They also found that the slope is different for clusters with Galactocentric distances smaller than the Sun’s and clusters with larger Galactocentric distances, although this may be an artifact due to incompleteness of the data. The topic of the distribution of the masses of star clusters has been re-visited for various star-forming galaxies (e.g. Bik et al., 2003; Elmegreen et al., 2000; Whitmore et al., 1999). Kroupa and Boily (2002) have recently approached the subject from a theoretical point of view focusing mainly on large clusters. They found that the entire Galactic population II stellar spheroid can be generated if star formation proceeded via embedded clusters distributed like a power-law mass function with slope in the range (1.9, 3.6).

In this section we will infer the OCIMF for the quiescent Milky Way disk without considering the values of the cluster masses provided by the catalogues but taking into account the life-time scale obtained from simulations in Section 3 and the open cluster sample described in Section 2 to derive the OCIMF for the Milky Way disk. This method to construct an OCIMF has been outlined in Kroupa and Boily (2002). In Section 4.2 it has been evidenced that the disk of the Milky Way may have experienced several episodes of enhanced star formation over the last 2 Gyr, this implies that the OCIMF cannot be derived from the total sample of clusters, because the starburst will bias the OCIMF towards massive clusters. However, if we restrict our analysis to clusters with ages less than about 200 Myr, these clusters can be used to derive the OCIMF. This choice implies that the star formation rate has remained almost constant during the last 200 Myr; this assumption appears to be reasonably supported by our previous analysis. For an average stellar mass of 0.4 $`M_{}`$, we obtain that $``$ 65% of clusters have masses $`20M_{}`$, about 17% with $`20M_{}<M<40M_{}`$, 10% with $`40M_{}M<60M_{}`$, 6% with $`60M_{}<M<75M_{}`$, and $``$ 2% have masses $`>75M_{}`$, with a likely average of 150 $`M_{}`$. Observations and simulations seem to suggest that $`M_{min}2M_{}`$ ($`N_{min}=5`$) and $`M_{max}16,000M_{}`$ ($`N_{max}=40,000`$). These results are consistent with an open cluster population being sampled from a power-law OCIMF with $`\alpha 2.7`$ with a correlation coefficient $`r=0.99`$ (see Fig. 10). The error associated with this determination can be estimated from the standard deviation obtained for the life-time of simulated clusters and is $`20\%`$. In this analysis we are not taking into consideration that a number of clusters (those younger than $`60`$ Myr) may be still in the evolutionary stage driven by gas expulsion. This simplification is likely to have a minor impact on the overall results from this section because our results indicate that a significant fraction of open clusters formed during quiescence are in fact small and therefore short-lived.

## 6 From the age distribution to the star formation history

In the previous sections we have assumed that the open cluster sample under study is representative of the entire Galactic disk, therefore the evolution of the star formation rate can be inferred from its age distribution, since the number of open clusters in each age bin has to be correlated with the number of objects initially at that time as a result of dissolution processes. The most reliable transformation of the open cluster age distribution into history of the star formation rate comprises two intermediate corrections, namely the volume and evolutionary corrections. As stated in Section 4, a scale height correction can be neglected because the studied age range is $`<2`$ Gyr. The applied corrections are explained in the following sections.

### 6.1 Evolutionary correction

As pointed out in Section 3, a correction due to the dynamical evolution of the cluster (cluster disintegration) is needed because our sample includes clusters with different initial populations ($`N_o`$). The more massive clusters have a life expectancy higher than the short-lived, small-$`N`$ clusters, thus the latter would be missing in the older age bins. It is however possible to correct for this effect using the OCIMF obtained in the previous section as well as the results from Section 3.

The corrections are given by the following formalism. The number of clusters born a time $`t`$ ago ($`t`$ = 0 is present time) is the current number observed at age $`t`$ divided by the fraction of surviving clusters expected for an age $`t`$, $`f`$. The fraction of surviving clusters changes with $`t`$, therefore the life-time scale obtained in Section 3 is used to estimate the correct fraction of surviving clusters to be applied. In detail: for the second bin we have corrected by the same number in bin 1; bins 3-4 have been corrected using $`f`$ = 0.18; $`f`$ = 0.08 has been used in bins 5-9; bins 10-16 have been corrected by $`f`$ = 0.02. After 0.8 Gyr, we have assumed that no cluster formed during quiescent star formation epochs may have survived; therefore, the correction is made in a slightly different way: we add the reference level at $`t`$ = 0 to the standard correction with $`f`$ = 0.02. The age histogram from WEBDA and the corrected distribution appear in Fig. 11. The star formation history outlined in Section 4 is somewhat altered:

#### 6.1.1 The 0.35 Gyr burst

This burst appears at 0.20-0.35 Gyr and it seems to be not very prominent.

#### 6.1.2 The 0.70 Gyr burst

This burst appears at 0.6-0.8 Gyr, lasting also 0.2 Gyr. It appears to be much stronger than burst 1.

#### 6.1.3 The 1.13 Gyr burst

It started 1.3 Gyr ago and lasted 0.35 Gyr. It was the strongest. The unusually high level of enhanced star formation could be the result of the combined action of density waves and a strong tidal interaction with the SMC.

#### 6.1.4 The 1.50 Gyr burst

This burst appears at 1.4-1.6 Gyr, lasting also 0.2 Gyr. It appears to be similar in strength to burst 2.

#### 6.1.5 The 1.93 Gyr burst

This burst started 2 Gyr ago, lasting 0.15 Gyr. Its activity was lower than for bursts 2 and 4 but much higher than 1.

### 6.2 Volume correction: the solar neighbourhood

Since our original open cluster sample is not volume-limited, there could be a bias in the relative number of objects in each age bin: open clusters with different initial populations and metallicities have different integrated luminosities, thus the volume of space sampled varies from cluster to cluster. In principle, our method will produce optimal results for volume-limited open cluster samples. In this section we consider a cluster sample in the solar neighbourhood. The solar neighbourhood is defined as a volume centered on the Sun that is much smaller than the overall size of the Milky Way galaxy and yet large enough to contain a statistically useful sample of stars (see, e.g., Binney & Tremaine, 1987). The appropriate size of the volume depends on which stars or objects are going to be investigated: for white dwarfs, which are both common and faint, it may consist of a sphere of radius 10 pc centred on the Sun, while for the bright but rare O and B stars, the solar neighbourhood may be considered to extend as far as 1 kpc from the Sun. In our analysis we will consider the latter in order to retain a statistically significant sample of open clusters. There are 161 open clusters in our sample within 1 kpc of the Sun (and younger than 2 Gyr).

If we restrict our analysis to clusters closer than 1 kpc, the age distribution plotted in Fig. 12 is obtained. Burst 1 is clearly identified at an age of 0.35 Gyr with a duration of about 0.2 Gyr. Burst 2 is also present for an age of about 0.7 Gyr with a duration of about 0.2 Gyr. Burst 3 appears at 1.15 Gyr and is 0.1 Gyr wide. Unfortunately, it is difficult to recover burst 4 and 5 as the number of clusters older than 1.2 Gyr within 1 kpc from the Sun is very low (3 objects). In fact, our results appear to suggest that the solar neighbourhood experienced quiescent star formation at 1.2-2 Gyr. From Fig. 12, the star formation rate history in the solar neighbourhood seems to be slightly different, with not a well defined period and narrower maxima. It is, however, easy to reconcile these apparent differences. The age distribution features interpreted as bursts 1, 2, and 3 are clearly identified in Fig. 12. Burst4 and 5 may also exist but there are too few objects to confirm the hypothesis. On the other hand, uncertainties in the open cluster ages are smaller for closer clusters, therefore maxima are narrower and show better contrast against steady, quiescent star formation.

Fig. 13 shows the solar neighbourhood sample with evolutionary correction as described in the previous section. In this case the second burst is the strongest, followed by the third, the fifth, and the first. The fourth one is now the least important.

## 7 Comparison with other authors

The problem of deducing the star formation history of our Galaxy has been primarily studied through chemical evolution models although other methods have also been used: Hertzsprung-Russell diagram inversion, stellar evolutionary tracks (Twarog, 1980; Meusinger, 1991; Bertelli & Nasi, 2001), chromospheric activity as measured by Ca II and K emission (Barry, 1988; Soderblom et al., 1991; Rocha-Pinto et al. 2000a, b), stellar kinematics (Gomez et al., 1990; Marsakov et al. 1990; Chereul et al., 1998; Just, 2002, 2003), the main-sequence luminosity function (Scalo, 1987), the white dwarf luminosity function (Noh & Scalo, 1990; Diaz-Pinto et al., 1994; Isern et al., 1999), combining the metallicity distribution and age-metallicity relation of G dwarfs (Rocha-Pinto & Maciel, 1997) and the distribution of coronal emission as measured by X-ray luminosities (Micela et al., 1993; Lachaume et al., 1999). Most of these studies have inferred an star formation history that is non-monotonic with time.

Although the majority of studies in this field use samples of stars in the solar neighbourhood with no stars more distant than about 100 pc being considered, it does not mean that the star formation history derived can only be applied to stars born in the solar neighbourhood. Nearby stars older than about 0.2 Gyr. come from birth sites which span a large range in Galactocentric distances. Wielen (1977) showed that the orbital diffusion coefficient deduced from the observed increase of velocity dispersion with age implies that such stars have suffered an rms azimuthal drift of about 2 kpc for an age of 0.2 Gyr. Considerable, but smaller, drift should occur also in the radial direction. Wielen et al. (1996), on the basis of the Sun’s metallicity and the radial metallicity gradient in the Galactic disk, estimated that the Sun has migrated outward by 1.9$`\pm `$0.9 kpc in the past 4.5 Gyr. Sellwood and Binney (2002) have discussed the dynamics of radial migration due to transient spiral arms, and the relationship between radial migration and disk heating. They estimate that old stars formed in the solar neighbourhood should now be scattered nearly uniformly throughout an annulus of 8.5$`\pm `$4.0 kpc. In this sense the star formation history inferred for nearby stars is a measure of the global star formation history, at least at the Sun’s Galactocentric radius. It is, however, not clear how the diffusion process affects entire star clusters not standalone stars and it is likely that star clusters experience little or no (azimuthal, radial) drift over the life-time.

### 7.1 Hernandez et al. (2000a)

Hernandez et al. (2000a) have derived the star formation history of the solar neighbourhood over the last 3 Gyr using data from the Hipparcos catalogue to construct color-magnitude diagrams. They treated these diagrams using advanced Bayesian analysis techniques (Hernandez et al., 1999; Hernandez et al., 2000) to deduce the star formation rate history of this region. They recovered the star formation history with an unprecedented time resolution of 0.05 Gyr. This high resolution makes it easy to compare our present results with those in their papers (Hernandez et al. 2000a, b, 2001). Their results indicate that the local star formation rate has an oscillatory component of period $`0.5`$ Gyr superimposed on a small level of constant star formation activity. Their inferred star formation history appears to be compatible with the observed distribution of stars. This cyclic behaviour is interpreted as the result of repeated encounters with the Galactic arm density pattern. In their work, they find that the last episode of enhanced star formation happened about 0.5 Gyr ago. Another burst is found at about 1.3 Gyr and yet another one at 1.9 Gyr. Their Fig. 4 (Hernandez et al., 2000a) however suggests a double peak for the youngest burst and also for the one at 1.3 Gyr. In fact, their results appear to indicate one burst at about 0.4 Gyr (our burst 1) and another one at 0.7 Gyr (our burst 2). The 1.3 Gyr burst may be the convolution of burst 3 and 4 found in our analysis. Their third burst coincides in time with our fifth one. Therefore, our derivation of the recent star formation history both in the solar neighbourhood and the Galactic disk can be considered as fully consistent with the results found by Hernandez et al. (2000a) in spite of the fact that a completely different technique has been used in its derivation.

### 7.2 Other authors

Rocha-Pinto et al. (2000a) have derived a star formation history from the chromospheric activity-age distribution of a large sample comprising 552 stars with Hipparcos parallaxes in the solar neighbourhood (stars more distant than 80 pc were omitted), and have found intermittency in the star formation rate over 14 Gyr. With a time resolution of 0.4 Gyr, their history of the star formation rate plotted in Fig. 2 (Rocha-Pinto et al., 2000a) and Table 1 (Rocha-Pinto et al., 2000b) indicates that the solar neighbourhood experienced a burst of star formation 0-1 Gyr ago (burst A). An epoch of low star formation took place 1-2 Gyr ago and a new, very strong burst was experienced 2-5 Gyr ago (burst B). An additional burst was found 7-9.5 Gyr ago (burst C). Although their time resolution makes a direct comparison impossible, their burst A is compatible with a convolution of our bursts 1 and 2, and possibly 3. There are, however, strong discrepancies in the 1-2 Gyr range.

Noh and Scalo (1990) using the white dwarf luminosity function found peaks of star formation at 0.3 Gyr and 1.8 Gyr, that are compatible with our bursts 1 and 5. Barry (1988) using a technique similar to the one in Rocha-Pinto et al. (2000a), has also found three bursts using a volume-limited sample of 115 F-G stars, the most recent peaked in the last 0.5 Gyr (burst A). An epoch of low star formation took place 1-4 Gyr ago and a new burst was experienced 4-6 Gyr ago (burst B). The oldest burst was found 7-11 Gyr ago (burst C). His burst A appears to be a convolution of our bursts 1 and 2. Using the main-sequence luminosity function, Scalo (1987) found signatures of starbursts 5-6 and 0.3 Gyr ago. His 0.3 Gyr burst coincides with our burst 1. Twarog (1980) found that during the last 4 Gyr the star formation rate has remained more or less constant but there was a sharp increase from 4 to 8 Gyr ago. Although the time resolution in these studies is rather low (0.2-1.0 Gyr), their results seem to indicate that bursts 1, 2, and 5 have also been identified by other authors.

## 8 Age distribution maxima: enhanced star formation vs. accretion origin

Starburst events are often seen in interacting galaxies, and it seems reasonable to infer a causal connection between galaxy interactions (mergers and close encounters) and starburst activity. It is now becoming increasingly clear that average galaxies engulf and devour smaller satellite galaxies as part of a more general process of hierarchical merging. Our own galaxy, the Milky Way, may have been incorporating dwarf galaxies into its disk since it was formed. The Galaxy is encircled by satellite galaxies that appear confined to one of two great streams across the sky (Lynden-Bell & Lynden-Bell, 1995). The most well-known of these are the Magellanic Clouds and the associated HI Magellanic stream. All of these are expected to merge with the Milky Way in the distant future, largely due to the dynamical friction from the extended halo. The first strong evidence in favour of a hierarchical merging scenario for the Milky Way was presented by Ibata et al. (1994, 1995, 1997) when they discovered the Sagittarius dwarf galaxy, a low mass dwarf spheroidal galaxy about 25 kpc from the Sun in an advanced state of disruption that is being absorbed by the disk of the Milky Way. In fact, our results appear to indicate that the Sagittarius Dwarf Galaxy may have had some role in inducing enhanced star formation in the Milky Way disk. For about 10 years, the Sagittarius Dwarf Galaxy remained as the only strong evidence that the Milky Way is incorporating satellite galaxies, but another object has been added recently. The Canis Major Galaxy (Martin et al., 2004; Bellazzini et al., 2004) appears to be a close relative of the Sagittarius Dwarf Galaxy in terms of structural parameters. This new finding seems to confirm that part of the Galactic disk could be extragalactic in origin. On the other hand, recent numerical models of disk galaxy formation (Abadi et al., 2003a, b) indicate that disrupted galaxy satellites might have contributed a significant fraction of the old stars in the disk of the Milky Way. Navarro et al. (2004) have confirmed that the Arcturus Group, a dynamically-coherent group of 10-12 Gyr old, metal-deficient stars located in the solar neighbourhood around the star Arcturus, are very likely of extragalactic origin. These authors suggest that the Group is part of the tidal debris of a 10<sup>8</sup> $`M_{}`$ disrupted satellite accreted by the Milky Way about 8 Gyr ago.

In this scenario, the peaks in the open cluster age distribution can be interpreted as signatures of merger events. This is however a non-straightforward interpretation as the variations found in the star formation rate appear to be periodic in nature, and mergers are expected to occur randomly in time. However, the actual sequence of events driving enhanced star formation episodes could be even more complex, with both tidal interactions and accretion events contributing to the global star formation history, as the unusual strength of burst 3 suggests.

For the main open cluster sample used in this work (Mermilliod, 2003), we find 97 clusters in the age range 0.2-0.6 Gyr and 80 objects in the age range 0.6-2.0 Gyr. There is a clear excess of clusters older than 0.6 Gyr. Strong star bursts are able to explain this feature but it could also be possible that some of these clusters may have been incorporated by the Galactic disk during accretion events in the recent past. In any case, the accretion scenario can only be important for coplanar interactions with the disk of the Milky Way and large accreted star clusters, otherwise the star clusters surviving the merger would end up either destroyed or in the Galactic halo.

### 8.1 The Canis Major galaxy role

The chain of dynamical events that provoked the formation of the recently discovered Canis Major Dwarf Galaxy remnant (Martin et al., 2003; Bellazzini et al., 2003) has been proposed as the process that built up the Monoceros Ring (Newberg et al., 2002; Yanny et al., 2003; Ibata et al., 2003; Crane et al., 2003; Rocha-Pinto et al., 2003; Majewski et al., 2004). An unusually high spatial density of open clusters has been found in the field of the Canis Major galaxy. The list of open clusters includes, at least: NGC 2204, 2243, and 2477, Tombaugh 1 and 2, Berkeley 20 and 33, Melotte 66, and AM-2. Bellazzini et al. (2003) have also found that the position, distance, and stellar population of the old open clusters AM-2 and Tombaugh 2 strongly suggest that they are part of the Canis Major galaxy. Both are older than 4 Gyr and therefore they have not been included in our age distributions for clusters younger than 2 Gyr. Nevertheless, if the physical association between these open clusters and the merged dwarf galaxy is confirmed, it will be a clear example of intergalactic stellar pollution, with foreign star clusters being absorbed by the Galactic disk. If true, it may have a major impact on the determination of the star formation history of the Milky Way, with a fraction of the observed stars and open clusters being actually outsiders coming from other star formation histories.

### 8.2 Searching for dwarf galaxy remnants

In this paper we have shown that open clusters can be used as tracers to investigate the star formation history of the Milky Way disk. If a number of dwarf galaxies with their cohorts of star clusters have been captured and disrupted in the past, open clusters absorbed from these galaxies can be used to find the trail left by the parent galaxy. Searching for groups of coeval old open clusters with similar distances located in the same region of the sky can help to identify the bulk of the accreted satellite.

## 9 Discussion and conclusions

Does the Galaxy form stars continuously, or in bursts separated by epochs of relative quiescence? If star formation occurs in bursts, what processes mediate the bursts? These are just two of the many questions that we have tried to answer in this paper. We have applied a new method to carry out an objective reconstruction of the star formation history of the Milky Way disk over the last 2 Gyr. A sample composed of 581 open clusters with known ages and distances was used in the derivation of this star formation history. Peaks in the age distribution diagram of the cluster sample were interpreted as signatures of star formation bursts. Our results indicate that the analyzed star formation rate presents two components: periodic episodes of enhanced star formation superimposed on a quiescent star formation level. A constant star formation history during the last 2 Gyr can therefore be discarded. A uniform rate of cluster formation and an exponentially declining dissolution rate would have produced a single population of clusters that would follow a straight line in a logarithmic plot with a slope given by the characteristic life-time of the population. Evidence for at least five epochs of enhanced star formation during the time interval studied has been found. The recent star formation history derived in this work is consistent with star formation rate histories deduced using a range of other techniques. Interpreting the age distribution diagram suggests that quiescent (non-enhanced) star formation has proceeded in discrete, highly-correlated regions of activity producing clusters which dissipate their parent molecular cloud on a time-scale of a few Myr ($`<`$ 10 Myr). The majority of open clusters formed during the non-enhanced star formation periods are destroyed and their members become field stars on a time-scale $`<`$ 20 Myr. However, there appears to be a long-lasting component as well, since some of the clusters are able to survive for over 0.2 Gyr. During bursts of star formation, the huge number of clusters formed make it possible to produce larger than average objects that could survive for several Gyr. A number of simulations were done in order to estimate the impact of the age errors on the features found in the age distribution diagram. Even for age errors larger than 30% some of these features (the younger ones) are still clearly identifiable. On the other hand, we have examined the possibility that the Galactic bursts are coeval with features in the star formation history of the Magellanic Clouds and the Sagittarius Dwarf Galaxy, as well as with close encounters between them and the Milky Way. Although the degree of uncertainty is large, there are several coincidences that suggest tidal interactions can play some role as inductive forces of bursts of star formation in the Milky Way disk. Distant galactic encounters may trigger significant star formation but quiescent star formation in the disk of the Milky Way appears to be characterized by a large number of small, short lived open clusters that contribute a small fraction of the total number of stars formed ($`<`$ 150 per cluster) with a small number of large, relatively long lived clusters that may contribute most of the stars in the disk.

There is, however, a cyclic behaviour in the burst sequence that may be better explained by the density wave hypothesis (Lin & Shu, 1964) for the presence of spiral arms in late-type galaxies. A model like the one outlined in Hernandez et al. (2000a) and developed in Martos et al. (2004) can explain easily the 0.4 Gyr periodicity that we detect. Heating by spiral structure can explain the main features of the age-velocity dispersion relation and other parameters of the velocity distribution in the solar neighbourhood, like why the stars in a single velocity-space moving group have a wide range of ages (De Simone et al., 2004). For a pattern speed $`\mathrm{\Omega }_p`$ = 20 km s<sup>-1</sup> kpc<sup>-1</sup> (e.g. Martos et al., 2004) and a Galactocentric distance $`R_{}`$ = 8.5 kpc it implies an orbital period of about 1 Gyr for the Sun. If the enhanced star formation episodes are, in fact, due to the interactions with the spiral arms it means that our Galaxy has two arms. This has been recently suggested by Martos et al. (2004) using a completely different approach. On the other hand, studies on the nature of star formation triggering in large disks indicate that a number of mechanisms may operate concurrently: gravitational instabilities, shocks between colliding clouds of gas, enhanced pressure of the interstellar medium, strong stellar winds, supernova shocks (star formation triggered by nearby bursts of star formation), density waves, shear forces produced by differential rotation, and interactions or mergers with other galaxies. Star formation triggered by previous star formation events is a self-propagating process which may continue over a much longer period of time, hence it is possible that overall a relatively large fraction of star formation is triggered.

On the other hand, if the star formation history inferred from the age distribution of open clusters is consistent with that deduced for individual stars then one can consider this as an argument in favour of most of the observed stars being formed in some kind of star clusters. Since only a subsample of the young clusters are likely to survive, an obvious question is whether most of the field stars in a galaxy are originally formed in clusters. Even in these young star forming regions many of the field stars are from clusters that have already been dissolved, hence the true percentage of stars that were originally in clusters is even higher, and might conceivably be $``$ 100%. It is tempting to suggest that perhaps the majority of stars are formed in groups and clusters, and that the field stars are simply the remnants of the fainter, less dense clusters which have dissolved. In fact, infrared observations show that most star formation occurs in embedded star clusters within Giant Molecular Clouds. These clusters have masses in the range 50-10<sup>3</sup> $`M_{}`$ (Lada & Lada, 2003). Nevertheless, it is still possible to consider the existence of two channels in the star formation process: the dispersed mode and the clustered mode, by modifying slightly our concepts of dispersed and clustered. The dispersed mode would be represented by the numerous short lived ($`\tau <`$ 30 Myr) small star clusters that are part of the quiescent star formation rate. These are different from classical open clusters and may be called Fast Living Clusters or flusters in opposition to the fully fledged standard open clusters as the Hyades or Pleiades. The clustered mode would be represented by the relatively long lived open clusters. In our analysis of the last 0.2 Gyr of star formation, the number of clusters in the age range 30-100 remains relatively constant as well as the fraction in the age range 100-200. This feature suggests a real gap in the mass distribution of disk star clusters: very few clusters are formed in the mass range 20-30 $`M_{}`$ ($`N`$ = 50-75). To explain the very good matching between the star formation history inferred from individual stars and that from open star clusters we have to assume that although most of open clusters ($`>`$ 90%) are poor, and short lived ($`N`$ 200) most of the stars ($`>`$ 90%) are formed in larger clusters ($`N>`$ 200). This conclusion, however, appears not to be supported by current observations.

The third largest galaxy in the Local Group, Triangulum galaxy (M33, NGC 598) may be undergoing star formation events similar to the ones described in this paper. M33 is a spiral galaxy about half the size of the Milky Way. Chandar et al. (1999) have found 44 young clusters in M33 with ages $``$ 100 Myr and masses in the range $`6\times 10^2`$ to $`2\times 10^4M_{}`$. Currently, M33 appears to be forming many young compact clusters but less massive than an older cluster population formed in a previous starburst with characteristic mass $`10^5M_{}`$.

Although our conclusions are uncertain as they are based on a relatively small (36%) subsample of a larger but still incomplete sample composed of almost 2000 open clusters, they are similar to the ones previously found by a number of other independent studies. This confirms the validity of our approach to infer star formation histories (recent and old), not just for the Milky Way disk but for any other galaxy. On the other hand, if the conclusions obtained in this paper are correct, the Galactic disk will start another episode of enhanced star formation within the next very few Myr if it is not already happening right now.

We thank the referee whose comments greatly improved this paper. We are very grateful to S.J. Aarseth for providing his computer code and for his comments on the manuscript. We thank the Department of Astrophysics of Universidad Complutense of Madrid for providing excellent computing facilities at the Centro de Proceso de Datos in Moncloa. The authors thank S. van den Bergh, L. O. Lodén, J.-C. Mermilliod, G. Carraro, and I. Platais for very helpful and interesting discussions. We also thank M. Martos, X. Hernández and associated group for sharing with us their results on the Galactic spiral pattern and its rotation speed in advance of publication. In preparation of this paper we made use of the Open Cluster Database, operated at the Institut d’Astronomie de l’Université de Lausanne, Switzerland, the New Catalogue of Optically Visible Open Clusters and Candidates, operated at the Universidade de São Paulo, Brazil, the NASA Astrophysics Data System and the ASTRO-PH e-print server.


format: text

text:
# 1 Introduction

## 1 Introduction

Colour Reconnection is the term used for strong interactions between colour singlet parton systems of different origin . The effect can influence the evolution of two nearby parton showers. The kinematics of the hadrons coming from both systems can therefore be perturbed. In the reaction $`\mathrm{e}^+\mathrm{e}^{}\mathrm{W}^+\mathrm{W}^{}\mathrm{q}_1\overline{\mathrm{q}}_2\mathrm{q}_3\overline{\mathrm{q}}_4`$ the colour singlets $`\mathrm{q}_1\overline{\mathrm{q}}_2`$ and $`\mathrm{q}_3\overline{\mathrm{q}}_4`$ formed by the decay products of both W bosons, could rearrange themselves to new colour singlets $`\mathrm{q}_1^{}\overline{\mathrm{q}}_4^{}`$ and $`\overline{\mathrm{q}}_2^{}\mathrm{q}_3^{}`$. Because the flow of the colour quantum numbers, reflecting the particle dynamics at short distances, controls the particle distributions in the final state, one could expect a change in these distributions after the colour rearrangement. In these $`\mathrm{W}^+\mathrm{W}^{}`$ events produced at LEP2, the separation distance between the $`\mathrm{W}^\pm `$ decay vertices is around $`\tau _\mathrm{W}1/\mathrm{\Gamma }_\mathrm{W}0.1`$ fm, while the fragmentation scale of the W bosons is around 1 fm. Hence there is a significant space-time overlap of both W hadronization regions where Colour Reconnection could occur. The resulting kinematic structure of a $`\mathrm{W}^+\mathrm{W}^{}`$ event will be different from the situation without this reconnection. The presence of similar effects was found in hadronic B meson decays $`\mathrm{B}\mathrm{J}/\psi +\mathrm{X}`$ where the colour interference between the two original colour singlets ($`\overline{\mathrm{c}}+\mathrm{s}`$ and c+spectator) suppresses this decay .

The probability to rearrange can be enhanced by gluon exchange between both W decay systems. Within the perturbative parton shower model it was shown that the Colour Reconnection or interference probability is negligible . Colour transmutations between partons from different W bosons can only occur from the interference of at least two emitted gluons. This interference piece is suppressed by the effective number of colours as $`1/\mathrm{N}_\mathrm{c}^2=1/9`$ compared to the $`𝒪(\alpha _\mathrm{s}^2)`$ non-reconnection emissions. Also the effects of a finite W width restrict the energy range of primary gluons generated by the alternative rearranged dipoles, because both W bosons do not necessarily decay at the same time.

Within the framework of the Lund string fragmentation implemented in JETSET the colour fields of both W boson strings can overlap in space-time. When the strings are described in their simplest Lorentz invariant way by a linear confinement potential, the event probability $`𝒫_\mathrm{i}`$ to reconnect is related to the volume, $`𝒪_\mathrm{i}`$, of the string overlap as:

$$𝒫_\mathrm{i}=1\mathrm{e}^{\kappa 𝒪_\mathrm{i}}.$$
(1)
The SK1 model parameter $`\kappa `$, with a fixed value for each event, is unknown and can only be tuned or measured from experimental data. For simplicity only one reconnection per event was allowed.

Within ARIADNE the Dipole Cascade Model is also followed by a string fragmentation according to the Lund model. As a criterion for Colour Reconnection one can therefore use the string length $`\mathrm{\Lambda }`$ defined in the momentum space

$$\mathrm{\Lambda }=\underset{\mathrm{i}=1}{\overset{\mathrm{n}1}{}}\frac{\mathrm{ln}(\mathrm{p}_\mathrm{i}+\mathrm{p}_{\mathrm{i}+1})^2}{\mathrm{m}_0^2}$$
(2)
where the sum covers all $`\mathrm{n}1`$ string pieces for a string connecting $`\mathrm{n}`$ partons and $`\mathrm{m}_0`$ is a typical hadronic mass scale (usually $`\mathrm{\Gamma }_\mathrm{W}`$). Reconnection is allowed when the string length after reconnection is shorter than the original one. In the ARIADNE 2 (AR2) model, reconnection between strings from different W bosons only happens below a fixed energy scale (usually $`\mathrm{\Gamma }_\mathrm{W}`$), while in the ARIADNE 3 (AR3) version reconnection is allowed at all energy scales which includes the perturbative phase. Therefore it is only the AR2 model which is relevant for this study. Within an event multiple inter-string reconnections and self-reconnections inside a single string could occur.

Before the formation of the clusters in the HERWIG fragmentation process, reconnection is allowed with a fixed probability ($``$ 1/9) when it reduces the sum of the squared sizes of the formed clusters. The size of a cluster is defined as the space-time separation of the production vertices of the partons forming the cluster. Multiple reconnections and self-reconnections within a single parton shower are allowed.

## 2 Particle flow measurement

The string effect predicts a higher particle production rate in the region between jets originating from the same $`\mathrm{W}\mathrm{q}\overline{\mathrm{q}}`$ decay (‘intra’ region), compared to the regions between jets from different W decays (‘inter’ region). When Colour Reconnection is present, particles tend to migrate from the intra-W to the inter-W region. Therefore an observable which counts the particles in these different regions could be sensitive to the different Colour Reconnection algorithms. The ratio $`\mathrm{R}`$ of the integrated particle density or particle flow in the intra region to the integrated density in the inter region quantifies this behavior.

The value obtained from the LEP combined data is compared to the different models in Figure 1. In the SK1 interpretation of Colour Reconnection the data prefers a value of $`\kappa `$ = 1.18, and the 68% confidence level lower and upper limits are 0.39 and 2.13 respectively. This corresponds to a reconnection probability of 49% in this model at $`\sqrt{\mathrm{s}}`$ = 189 GeV. The extreme SK1 scenario with 100% reconnection probability disagrees with the data at a level of 5.2 $`\sigma `$. The log-likelihood as a function of $`\kappa `$ is shown in Figure 2. All four experiments have observed a very weak sensitivity to the ARIADNE and HERWIG Colour Reconnection models, which does not coincides with the SK1 prediction of the effect. The expected value of $`\mathrm{R}`$ from those fragmentation models in the hypothesis that no Colour Reconnection is present, differs from the measured LEP combined data value by 3.1 and 3.7 $`\sigma `$ for respectively the ARIADNE and HERWIG model. This may indicate an underestimation of the systematic uncertainties in these particle flow analyses. The dominating systematic uncertainties are the fragmentation modeling in the simulation and the possibility for Bose-Einstein Correlations in the final state.

## 3 The $`𝚫𝐦_𝐖`$ method

It has been shown that the $`\mathrm{m}_\mathrm{W}`$ observable  inferred from hadronic decaying $`\mathrm{W}^+\mathrm{W}^{}`$ events at LEP2 by the method of direct reconstruction, is influenced when changing the value of $`\kappa `$. A second method is therefore based on the observation that two different $`\mathrm{m}_\mathrm{W}`$ estimators have different sensitivity to the parametrized Colour Reconnection effect. Hence the difference between them is an observable with information content about $`\kappa `$.

It is observed that mostly low momentum particles and particles in inter-jet regions are affected by Colour Reconnection and hence influencing the measurant, $`\mathrm{m}_\mathrm{W}`$. An alternative analysis was designed which neglects these particles in the reconstruction of the momenta of the four primary partons. This by applying an optimized cone algorithm to cluster the particles in jets, rather than using inclusive clustering algorithms like DURHAM.

An indirect measurement of the SK1 model parameter $`\kappa `$ is possible from the direct measurement of the difference $`\mathrm{\Delta }\mathrm{m}_\mathrm{W}(\mathrm{i},\mathrm{j})`$ in reconstructed $`\mathrm{m}_\mathrm{W}`$ between the CR-sensitive standard (i) and CR-less-sensitive alternative (j) analysis. The observable $`\mathrm{\Delta }\mathrm{m}_\mathrm{W}(\mathrm{std},\mathrm{R}_{\mathrm{cone}})`$ where $`\mathrm{R}_{\mathrm{cone}}`$ is around 0.5 rad, was found to be most sensitive (4.3 $`\sigma `$ sensitivity to the full SK1 model). The systematic uncertainty is dominated by possible effects like Bose-Einstein Correlations or detector related discrepancies in the energy flow reconstruction. The log-likelihood as a function of $`\kappa `$ obtained by DELPHI is shown in Figure 2.

After rescaling the effect according to the overall reconnection probability difference between the SK1 and the HERWIG model, also similar sensitivities were observed for both models. Only a negligible sensitivity to the ARIADNE implementation was found.

The correlation between the $`\mathrm{\Delta }\mathrm{m}_\mathrm{W}(\mathrm{std},\mathrm{R}_{\mathrm{cone}}=0.5\mathrm{rad})`$ and standard $`\mathrm{m}_\mathrm{W}`$ estimators was found to be around 11%.

## 4 Feasibility study for the $`𝐦_𝐖`$ measurement

The log-likelihood information about $`\kappa `$ from both measurements is combined in Figure 2. Within the SK1 model a significant amount of Colour Reconnection is found.

The $`\mathrm{m}_\mathrm{W}`$ estimators applied by the four LEP Collaborations are equally sensitive to the effect of Colour Reconnection. Therefore one can expect to infer the same amount of likelihood information from all four datasets. In this hypothesis and assuming that the $`\mathrm{\Delta }\mathrm{m}_\mathrm{W}`$ analysis would reveal the same most likely value for $`\kappa `$, the total log-likelihood obtained by LEP2 would be Figure 3. For $`\mathrm{m}_\mathrm{W}`$ measured in the fully hadronic channel the minimum of this log-likelihood corresponds to a shift of 106 MeV/c2 at $`\sqrt{\mathrm{s}}`$ = 200 GeV with a 68% CL of \[77;144\] MeV/c2 (where -2log$``$=1). These numbers should be compared with the statistical uncertainty of 28 MeV/c2 in the absence of systematic uncertainties . Assuming that the SK1 model is the true model of Colour Reconnection, one could calibrate for this shift. The remaining uncertainty due to Colour Reconnection would be about 35 MeV/c2 (neglecting ARIADNE and HERWIG). This would reduce the total uncertainty on $`\mathrm{m}_\mathrm{W}`$ measured in the fully hadronic channel from about 110 to 60 MeV/c2 and increase the relative weight of this channel from 9 to 29% in the combination with the semi-leptonic channel.

## 5 Conclusion

Two basically uncorrelated measurements observe the same amount of Colour Reconnection in the LEP2 data according to the SK1 phenomenological model. Within the direct measurement of $`\mathrm{m}_\mathrm{W}`$ this effect was not corrected for and therefore induced a significant systematic uncertainty. A simple but efficient calibration method was proposed to benefit from the inferred information about the Colour Reconnection effect in the $`\mathrm{m}_\mathrm{W}`$ measurement. The uncertainty on the LEP2 value of the direct measurement of $`\mathrm{m}_\mathrm{W}`$ would reduce from 42 to 39 MeV/c2.