The poorly constrained cluster disruption time-scale in the Large Magellanic Cloud

TL;DR

This study uses simulations and observations to estimate the cluster disruption timescale and formation history in the Large Magellanic Cloud, finding a lower limit of 1 Gyr for disruption and an increasing formation rate over time.

Contribution

It provides new constraints on the cluster disruption timescale and formation history in the LMC using a combination of Monte Carlo simulations and homogeneous observational data.

Findings

Disruption timescale for 10^4 M_sun clusters is at least 1 Gyr.

Cluster formation rate has increased from 0.3 to 20-30 clusters per Myr over 5 Gyr.

The cluster formation history depends on the assumed initial cluster mass function shape.

Abstract

We use Monte-Carlo simulations, combined with homogeneously determined age and mass distributions based on multi-wavelength photometry, to constrain the cluster formation history and the rate of bound cluster disruption in the Large Magellanic Cloud (LMC) cluster system. We evolve synthetic star cluster systems formed with a power-law initial cluster mass function (ICMF) of spectral index $\alpha =-2$ assuming different cluster disruption time-scales. For each of these disruption time-scales we derive the corresponding cluster formation rate (CFR) required to reproduce the observed cluster age distribution. We then compare, in a Poissonian $\chi^2$ sense, model mass distributions and model two-dimensional distributions in log(mass) vs. log(age) space of the detected surviving clusters to the observations. Because of the bright detection limit ($M_V^{\rm lim} \simeq -4.7$ mag) above which the observed cluster sample is complete, one cannot constrain the characteristic disruption time-scale for a $10^4$ M$_\odot$ cluster, $t_4^{\rm dis}$ (where the disruption time-scale depends on cluster mass as $t_{\rm dis} = t_4^{\rm dis} (M_{\rm cl} / 10^4 {\rm M}_\odot)^0.62$), to better than $t_4^{\rm dis} \ge 1$ Gyr. We conclude that the CFR has increased from 0.3 clusters Myr$^{-1}$ 5 Gyr ago, to a present rate of $(20-30)$ clusters Myr$^{-1}$. For older ages the derived CFR depends sensitively on our assumption of the underlying CMF shape. If we assume a universal Gaussian ICMF, then the CFR has increased steadily over a Hubble time from $\sim 1$ cluster Gyr$^{-1}$ 15 Gyr ago to its present value. If the ICMF has always been a power law with a slope close to $\alpha=-2$, the CFR exhibits a minimum some 5 Gyr ago, which we tentatively identify with the well-known age gap in the LMC's cluster age distribution.