abstract: Simulating images representative of neurodegenerative diseases is important for predicting patient outcomes and for validation of computational models of disease progression. This capability is valuable for secondary prevention clinical trials where outcomes and screening criteria involve neuroimaging. Traditional computational methods are limited by imposing a parametric model for atrophy and are extremely resource-demanding. Recent advances in deep learning have yielded data-driven models for longitudinal studies (e.g., face ageing) that are capable of generating synthetic images in real-time. Similar solutions can be used to model trajectories of atrophy in the brain, although new challenges need to be addressed to ensure accurate disease progression modelling. Here we propose Degenerative Adversarial NeuroImage Net (DaniNet) — a new deep learning approach that learns to emulate the effect of neurodegeneration on MRI by simulating atrophy as a function of ages, and disease progression. DaniNet uses an underlying set of Support Vector Regressors (SVRs) trained to capture the patterns of regional intensity changes that accompany disease progression. DaniNet produces whole output images, consisting of 2D-MRI slices that are constrained to match regional predictions from the SVRs. DaniNet is also able to maintain the unique brain morphology of individuals. Adversarial training ensures realistic brain images and smooth temporal progression. We train our model using 9652 T1-weighted (longitudinal) MRI extracted from the Alzheimer's Disease Neuroimaging Initiative (ADNI) dataset. We perform quantitative and qualitative evaluations on a separate test set of 1283 images (also from ADNI) demonstrating the ability of DaniNet to produce accurate and convincing synthetic images that emulate disease progression.
author: Daniele Ravi; Daniel C. Alexander; Neil P. Oxtoby; for the Alzheimer's Disease Neuroimaging Initiative\*
title: Degenerative Adversarial NeuroImage Nets:
  Generating Images that Mimic Disease Progression

# Introduction

Neurodegenerative diseases are strongly age-related, and the incidence of these diseases is expected to rise in the next few years due to increasing life expectancy. Creating tools that assist our limited understanding of neurodegeneration can help to answer some of the open questions in this field, e.g., What are the causes? How can we improve subtype classification?

Disease progression modelling  analyzes clinical and image-based biomarkers to map out longitudinal changes during chronic diseases. However, modelling temporal neurodegeneration on full resolution MRI is still a major challenge. To tackle this problem, a few simulators have been proposed in the literature . The recent SimulAtrophy  uses a computational model based on fluid mechanics. This approach combines two deformation fields: one from a biophysical model and the other obtained by non-rigid registration of two real images. For a given baseline MRI, the combined deformation field is used to impose the desired level of atrophy and generate the simulated image. SimulAtrophy is extremely resource-demanding and is not scalable to high-resolution images. Exploiting the success of deep learning, Bowles et al.  proposed a framework based on Generative Adversarial Networks (GANs) to model and manipulate MRI directly through the technique of image arithmetic. Although the system is able to introduce or remove atrophy patterns from regions in the brain, it presents a few important limitations: disease progression is modelled linearly and morphological changes are the same across all patients. A more advanced deep learning framework developed to predict face ageing was proposed in . Our idea is to use a similar framework to generate realistic images that also preserve biological/anatomical constraints associated with disease progression.

In this paper we present a new framework that overcomes the limitations of existing approaches. Firstly, similar to  we employ adversarial training to ensure a high level of realism in the synthetic MRI. Secondly, we avoid imposing predefined atrophy patterns, prefering to implement novel biological constraints that model neurodegeneration. These constraints control the MRI intensity in each brain region through a set of localised SVRs learnt directly from the data. Finally, we design DaniNet to handle non-imaging characteristics, which we employ here to condition upon age and diagnosis. Experiments evaluate the performance of the various components of DaniNet and demonstrate the ability to produce accurate and realistic synthetic images that emulate disease progression.

# Methods

In this section we present the main blocks which compose DaniNet. A colour-coded work-flow of our system is depicted in Fig. .

The first pre-processing block (shown in blue in Fig. ) extracts a normalized slice $x$ from the MRI $I$. The second component (shown in grey) is a Conditional Deep Autoencoder (CDA) composed of two deep neural networks: an encoder $E$ that embeds $x$ in a latent space $Z$, and a generator $G$ that projects the vectors produced by $E$ back to the original manifold. Before this projection, the latent vector $z$ is conditioned with two variables: i) $d$ — a numerical representation \[0-3\] of diagnosis (i.e. cognitively normal, subjective memory concern, early/late mild cognitive impairment, Alzheimer's disease); and ii) $a \in [1,A]$ — an index describing age, binned into $A$=10 groups. This age discretisation is important for computing deformation loss (block in yellow) to learn morphological changes along the progression. The third component (shown in green), consists of two discriminator networks: i) $D_z$ that drives $E$ to produce $z$ with a uniform prior and smooth temporal progression; and ii) $D_{b}$ that drives $G$ to produce realistic *brain* neuroimages. Finally, in the blocks coded in orange, we present the proposed biological constraints used to model disease progression. In the next sub-sections, we will describe each of these components in detail.

## Pre-processing

The pre-processing block removes irrelevant variation on the data and ensures that intensity values in each voxel only decrease. In an ideal scenario, intensity of T1-weighted MRI decreases with age since tissue density (having high intensity) will reduce, while water content (having low intensity) will increase. However, due to scanner variability and other sources of noise, this may not be the case in practice. To handle these variabilities, each input image is normalised using the following pre-processing steps: i) linear co-registration to a 1mm isotropic MNI template using FLIRT-FSL, ii) skull-stripping using BET-FSL, and iii) intensity normalisation to zero mean and unit standard deviation. Images where pre-processing failed were not included in the training set.

## Conditional Deep Autoencoder

A CDA is an extension of the Deep Autoencoder (DA) with the ability to integrate conditional variables inside the generator $G$. The advantage of training a single end-to-end CDA over training separate DAs is that it avoids overfitting individual groups when longitudinal data are missing. For the architecture of $E$ and $G$ we follow . The output of the encoder $E$ is a feature vector $z \in \mathbb{R}^s$ that preserves the brain morphology of the current participant. The output of $G$ are synthetic images defined as $g_a=G(E(x),a,d)$. The aim of $G$ is to learn the mapping between the linear transition of the latent vector conditioned on $a$ and $d$ and the non-linear transition in the original manifold.

## Adversarial Training

Again following , DaniNet includes two discriminator networks that are trained adversarially with the CDA. The first discriminator $D_z$ guides $E$ to generate $z$ with a uniform distribution $\mathbb{U}$ to ensure temporal smoothness, as demonstrated in . Specifically, $D_z$ is trained to distinguish the vector generated by the encoder $E$ from samples extracted from $\mathbb{U}$. On the other side, $E$ is trained in a zero-sum game with the purpose to fool $D_z$. The objective function used for this adversarial training is

$$\label{dz_loss}
\min\limits_{E}\max\limits_{D_z}\mathbb{E}_{z^*}\big[ \log D_{z}(z^*) \big]+\mathbb{E}_{x}\big[1- \log D_{z}(E(x)) \big],$$ where $\mathbb{E}$ is the expectation operator, $z^*$ is a vector sampled from $\mathbb{U}$, $D_{z}$ estimates the probability that a vector comes from $\mathbb{U}$, and $E(x)$ is the latent vector obtained from $x$.

The second discriminator $D_{b}$ guides the generator $G$ to produce realistic brain images. This discriminator is trained to distinguish the synthetic images $g_a$ from the real images in the training set via the following objective function

$$\label{dbrain_loss}
\min\limits_{G}\max\limits_{D_{b}}\mathbb{E}_{x}\big[\log D_{b}\big(x)\big) \big] +
\mathbb{E}_{x}\big[1-D_{b}\big(G(E(x),a,d)\big)\big],$$ where $D_{b}$ estimates the probability that a slice is extracted from a real MRI.

## Biological Constraints

To capture the patterns of image intensity changes that accompany disease progression, our framework uses two separate loss functions during training. These loss functions impose biological constraints that mimic neurodegeneration by ensuring decreased intensity (tissue density) that is consistent with disease progression.

The first loss function operates at the voxel level. Given a synthetic output $g_a$ in age group $a$, we impose that corresponding voxels of each preceding image $g_i$ with $i<a$ have higher-intensity values. Likewise, later images $g_j$ with $j>a$ have lower-intensity values. We express this constraint as $$\label{L_vox}
L_{\textrm{vox}}=\frac{1}{MN(A-1)}\bigg[\displaystyle\sum_{p=1}^{a-1}\textrm{sgn}(g_a-g_p)+
\displaystyle\sum_{p=a+1}^{A}\textrm{sgn}(g_p-g_a)\bigg],$$ where $M$ and $N$ are the number of rows and columns in each slice, and sgn is the sign function.

The second loss function operates at a regional level to improve spatial consistency across neighbouring voxels. We pre-train a set of SVRs to learn how to predict the rate of intensity progression in fixed, overlapping, data-driven regions mask obtained from hierarchical clustering of local intensity of the training images. On the considered axial slice the clustering algorithm detects 128 regions with an average size of 233 voxels. We train one SVR per region, where the input features are the age at baseline, the age at the follow-up, and diagnosis $d$. We restrict the SVR training to monotonically decreasing data by removing time-points where regional intensity increases. The loss function used to constrain the image sequence $g_i$ (with $i \in [1,A]$) through the regional SVR predictions is

$$\begin{gathered}
\label{L_reg}
L_{\textrm{reg}}=\frac{1}{R(A-1)}\displaystyle\sum_{i=1}^{R}\bigg[\displaystyle\sum_{p=1}^{a-1}\bigg(\textrm{SVR}_i(p,a,d)-\frac{\sum[g_a\cdot r_i]+\epsilon}{\sum[g_p\cdot r_i]+\epsilon}\bigg)+
\\
\displaystyle\sum_{p=a+1}^{A}\bigg(\textrm{SVR}_i(a,p,d)-\frac{\sum[g_p\cdot r_i]+\epsilon}{\sum[g_a\cdot r_i]+\epsilon}\bigg)\bigg],
\end{gathered}$$ where $R$ is the number of regions, $r_i$ is the i-th region-mask, $SVR_i(p,a,d)$ is the prediction of intensity rate progression for the region $i$, and $\epsilon=0.1$ avoids division by 0.

The pre-trained SVRs supply a set of constraints on the generator's outputs for ensuring accurate emulation of neurodegeneration. The process used to compute the rate of regional intensity progression is summarized in Fig. .

## Deformation Loss

Our final loss function ensures consistency with the individual subject over time by minimising the difference between individual input image sequences and the corresponding output sequences. For each input image $x$, the difference is minimised between $x$ and a weighted average of two outputs from the nearest age bins:

$$\label{def}
L_{\textrm{def}}=L_{2}\bigg(x , \big[g_a\alpha+g_{a+1}(1-\alpha)\big]\bigg),$$ where $\alpha$ reflects the distance between input age and the group $a$ and $a+1$.

In summary, the loss functions $L_{\textrm{reg}}$ and $L_{\textrm{vox}}$ ensure monotonic intensity change to mimic neurodegeneration, and $L_{\textrm{def}}$ together with $D_{b}$ ensure realistic brain morphology.

## Training Details, Parameters and Evaluation

Data used in the preparation of this article were obtained from the ADNI database (adni.loni.usc.edu). The ADNI was launched in 2003 as a public-private partnership, led by Principal Investigator Michael W. Weiner, MD. The primary goal of ADNI has been to test whether serial magnetic resonance imaging (MRI), positron emission tomography (PET), other biological markers, and clinical and neuropsychological assessment can be combined to measure the progression of mild cognitive impairment (MCI) and early Alzheimer's disease (AD).

Training data (MRI-slices: 9852; participants: 876) and testing (MRI-slices: 1283; participants: 179) are pre-processed T1-weighted MRI from the ADNI dataset. This includes data from multiple sites, scanners, and pre-processing pipelines. We included them all to improve robustness to such variability, and hence increase generalizability. The participants were aged between 63 and 87 years old, with the following distribution of clinical disease stage: 28% cognitively normal, 4% subjective memory concern, 54% mild cognitive impairment and 14% Alzheimer's disease. Each participant has on average 4.7 MRI spanning 3 years. We select participants in the test set having at least one follow-up visit two years after baseline, to allow sufficient time for observable neurodegeneration to occur.

The architectures of each network $E$, $G$, $D_{b}$, $D_Z$ are based on the implementation proposed in . The size of the latent space $Z$ is fixed to 200. The network's parameters of $E$ and $G$ are trained to minimize a total loss defined as $$\label{Tot}
L_{\textrm{tot}} = E_z * 0.0003 + G_b * 0.0003 + L_{\textrm{vox}} * 0.0003 +L_{\textrm{reg}} * 0.0003 + L_{\textrm{def}} * 0.2$$ where $E_z$ is the cross entropy obtained by the discriminator $D_z$ on the generated latent vectors, and $G_b$ is the cross entropy obtained by the disctiminator $D_b$ on synthetic images $g_a$. To train our system we use the stochastic gradient descent solver, ADAM ($\alpha$ = 0.0002, $\beta$`<!-- -->`{=html}1 = 0.5) on an NVIDIA GTX TITAN-X GPU card. The training procedure converged after 3000 epochs of random mini-batch of 100 slices with a size of 128$\times$`<!-- -->`{=html}128. We also align each test individual to the model by baseline age and diagnosis, then personalise DaniNet with a transfer learning step ($T$ component, see results section ) that involves an additional 200 training iterations on the single baseline image. Personalisation is essential to tune the model for the specific morphology of the individual's brain, ready to simulate realistic, personalised followup MRI. To note, only a single MRI from each test subject's first visit is used to personalise the trained model.

In our evaluation, input $x$ are from baseline visits. Follow-up data is the ground truth $y$ against which we evaluate generated images $g_i$ using complementary quantitative and qualitative analysis. The quantitative analysis involves the use of the Structural Similarity Index matrix (SSIM). Since this matrix is not able to quantify the actual realism of the obtained images we also perform a complementary qualitative analysis via a survey of clinicians and medical imaging experts that evaluate the synthetic images based on their perception.

```latex
\begin{table*}[t]\scriptsize
\caption{Performance of DaniNet on different configurations.}
\label{table:Configuration}
\begin{center}
\resizebox{0.90\textwidth}{!}{%
\begin{tabular}{c|>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}>{\centering\arraybackslash}p{1cm}}
 \multirow{ 4}{*}{\rotatebox{90}{SSIM}}& ~\cite{zhang2017age} & $C$ & $T$ & $C$-$T$ &$P$ &$P$-$C$ &$P$-$T$ & $P$-$C$-$T$\\
\cline{2-9}
&0.43 $\pm$ 0.07&0.44 $\pm$ 0.06&0.59 $\pm$ 0.14&0.60 $\pm$ 0.14&0.44 $\pm$ 0.08&0.48 $\pm$ 0.08&0.61 $\pm$ 0.15&\textbf{0.62 $\pm$ 0.16} \\

\end{tabular}
}
\end{center}
\end{table*}
```

# Results

Our first experiment is designed to show that DaniNet improves similarity between generated image and ground truth, with respect to the baseline approach developed for face ageing . Specifically, we evaluate the contribution obtained by three components of DaniNet: i) $P$ obtained exploiting $L_{\textrm{vox}}$ and $L_{\textrm{reg}}$, ii) $C$ obtained conditioning the system with the clinical diagnosis, and iii) $T$ obtained when the learned model is transferred before testing on new data.

For this purpose, the SSIM results obtained using different compositions of $P$, $C$, and $T$ are reported in Table . Clearly, the best results are obtained when the framework includes all three components. The baseline model produces inferior SSIM. Transfer learning $T$ is the feature that provides the highest contribution. A moderate contribution is provided by $P$ and $C$. These improvements were assessed statistically using a paired t-test and all p-values were less than 0.0001. Visual results for the configurations $C$-$T$, $P$-$C$-$T$ and baseline  are shown in Fig. . The error maps, shown in the last row of the figure, confirm our findings. A full simulation with the optimal configuration $P$-$C$-$T$ is shown in Fig.  where neurodegeneration is apparent in the progression (left to right), most obviously in ventricular expansion and cortical thinning.

Our second experiment tests the qualitative performance of DaniNet based on human visual perception. In a survey, we asked 25 medical imaging experts (neurologists and computer scientists) to evaluate on average 36 sets of three images each, randomly selected from the test set, generated using two different configurations of DaniNet ($C$-$T$ and $P$-$C$-$T$) and the baseline . These were shown to the user in a random order. The input $x$ and the related follow-up $y$ were also displayed on the screen as references for the participants. For each set of images, the user was asked to select the closest synthetic image to $y$. Results from the survey confirm that DaniNet is a considerable improvement of the baseline approach. The medical imaging experts selected the configuration $P$-$C$-$T$ $26\pm5$ (mean $\pm$ std) times ($72\pm14\%$), the configuration $C$-$T$ $7\pm4$ times ($19\pm11\%$), and the baseline $1\pm1$ times ($3\pm3\%$). For only $2\pm3$ outputs ($6\pm8\%$) the users were not happy with any of the generated synthetic images.

# Conclusion and Future Work

We have proposed and evaluated (quantitatively and qualitatively) a novel deep-learning framework that is able to learn how to emulate the effect of neurodegenerative disease progression on structural MRI. The framework produces personalised, realistic output images through a combination of biological constraints, transfer learning, and conditioning upon both fixed and variable non-imaging characteristics. To the best of our knowledge, we are the first to propose a simulator that imitates realistic neurodegeneration by imposing biological constraints. In future work we will extend the framework to simulate entire 3D-MRI with the aid of low-memory techniques . Additionally, focal brain pathologies, such as glioma and white matter hyperintensities, are not currently modelled by our framework; future work will extend to consider such effects. Finally, although we have demonstrated our framework for modelling atrophy in MRI, we believe that this solution can be used with different image modalities (e.g. PET, CT, etc.) and to model disease progression on other organs (e.g lung, retina, etc.).

# Acknowledgement

The authors would like to thank NVIDIA Corporation for the donation of the Titan Xp GPU used for this research.

Data collection and sharing for this project was funded by the ADNI (National Institutes of Health Grant U01 AG024904) and DOD ADNI (Department of Defense award number W81XWH-12-2-0012). ADNI is funded by the National Institute on Aging, the National Institute of Biomedical Imaging and Bioengineering, and through generous contributions from the following: AbbVie, Alzheimer's Association; Alzheimer's Drug Discovery Foundation; Araclon Biotech; BioClinica, Inc.; Biogen; Bristol-Myers Squibb Company; CereSpir, Inc.; Cogstate; Eisai Inc.; Elan Pharmaceuticals, Inc.; Eli Lilly and Company; EuroImmun; F. Hoffmann-La Roche Ltd and its affiliated company Genentech, Inc.; Fujirebio; GE Healthcare; IXICO Ltd.; Janssen Alzheimer Immunotherapy Research & Development, LLC.; Johnson & Johnson Pharmaceutical Research & Development LLC.; Lumosity; Lundbeck; Merck & Co., Inc.; Meso Scale Diagnostics, LLC.; NeuroRx Research; Neurotrack Technologies; Novartis Pharmaceuticals Corporation; Pfizer Inc.; Piramal Imaging; Servier; Takeda Pharmaceutical Company; and Transition Therapeutics. The Canadian Institutes of Health Research is providing funds to support ADNI clinical sites in Canada. Private sector contributions are facilitated by the Foundation for the National Institutes of Health (www.fnih.org). The grantee organization is the Northern California Institute for Research and Education, and the study is coordinated by the Alzheimer's Therapeutic Research Institute at the University of Southern California. ADNI data are disseminated by the Laboratory for Neuro Imaging at the University of Southern California.

This project has received funding from the European Union's Horizon 2020 research and innovation programme under grant agreement No. 666992.

EPSRC grant EP/M020533/1 supports DCA's work on this topic. The NIHR UCLH Biomedical Research Centre also supports this work.


abstract: The early evolution of star clusters in the Small Magellanic Cloud (SMC) has been the subject of significant recent controversy, particularly regarding the importance and length of the earliest, largely mass-independent disruption phase (referred to as "infant mortality"). Here, we take a fresh approach to the problem, using an independent, homogeneous data set of $UBVR$ imaging observations, from which we obtain the SMC's cluster age and mass distributions in a self-consistent manner. We conclude that the (optically selected) SMC star cluster population has undergone at most $\sim 30$ per cent ($1\sigma$) infant mortality between the age range from about $(3-10)$ Myr, to that of approximately $(40-160)$ Myr. We rule out a 90 per cent cluster mortality rate per decade of age (for the full age range up to $10^9$ yr) at a $> 6 \sigma$ level. We independently affirm this scenario based on the age distribution of the SMC cluster sample.
author: Richard de Grijs$^{1,2}$[^1] and Simon P. Goodwin$^1$
$^1$ Department of Physics & Astronomy, The University of Sheffield, Hicks Building, Hounsfield Road, Sheffield S3 7RH
$^2$ National Astronomical Observatories, Chinese Academy of Sciences, 20A Datun Road, Chaoyang District, Beijing 100012, China
date: Received date; accepted date
title: Star cluster "infant mortality" in the Small Magellanic Cloud (Redivivus)

# Introduction

The early evolution of the star cluster population in the Small Magellanic Cloud (SMC) has been the subject of considerable recent attention and vigorous debate (e.g., Rafelski & Zaritsky 2005; Chandar, Fall & Whitmore 2006; Chiosi et al. 2006; Gieles, Lamers & Portegies Zwart 2007). The key issue of contention is whether the SMC's star cluster system has been subject to the significant early cluster disruption processes observed in "normal", interacting and starburst galaxies commonly referred to as "infant mortality" and "infant weight loss". Chandar et al. (2006) argue that the SMC has been losing up to 90 per cent of its star clusters per decade of age, at least for ages from $\sim 10^7$ up to $\sim 10^9$ yr, whereas Gieles et al. (2007) conclude that there is no such evidence for a rapid decline in the cluster population, and that the decreasing number of clusters with increasing age is simply caused by fading of their stellar populations. They contend that the difference between their results was due to Chandar et al. (2006) assuming that they were dealing with a mass-limited sample, whereas it is actually magnitude-limited. In fact, this is not entirely correct; Chandar et al. (2006) analyse the full magnitude-limited sample and conclude that it is approximately surface-brightness limited. They then compare the cluster age distribution of the full sample (expressed in units of ${\rm d}N_{\rm cl} / {\rm d}t$, i.e., the number of clusters per unit time period) to that of a subsample for masses $\ge 10^3$ M$_\odot$ (which they do not analyse in the same manner), and suggest both to be similar, although the latter is much flatter[^2], hence giving rise to the discrepancy between their results and those of Gieles et al. (2007). Both studies are based on the same data set, the Magellanic Clouds Photometric Survey (MCPS; Zaritsky, Harris & Thompson 1997).

The main contribution of this paper to the ongoing debate is two-fold: (i) We revisit the SMC's early star cluster evolution using an alternative approach; and (ii) we use independently obtained ages and masses based on an independent data set, i.e., the $UBVR$ photometric survey of the Magellanic Clouds by Massey (2002), originally analysed by Hunter et al. (2003). We conclude that there is indeed only marginal evidence for infant mortality in the SMC star cluster sample, supporting the careful analysis of Gieles et al. (2007). In Sect. we first briefly introduce the concept of cluster infant mortality. We discuss our observational data and the basic analysis leading to the age and mass estimates in Sect. . In Sect. we justify our choice of age ranges to construct cluster mass functions (CMFs). Finally, in Sect. , we present our case for the absence of significant cluster infant mortality.

# Cluster infant mortality

Observations of increasing numbers of interacting and starburst galaxies, including the Antennae system, M51 and NGC 3310, show a significantly larger number of young ($\lesssim 10-30$ Myr) star clusters than expected from a simple extrapolation of the cluster numbers at older ages, taking into account the observational completeness limits and the effects of sample binning, and under the additional, simplifying assumption that the star cluster formation rate (CFR) has been roughly constant over the host galaxy's history (e.g., de Grijs et al. 2003b; Whitmore 2004; Bastian et al. 2005; Fall, Chandar & Whitmore 2005; Mengel et al. 2005; Chandar et al. 2006; see also Whitmore, Chandar & Fall 2007 for a presentation of earlier results, and de Grijs & Parmentier 2007 for a review). This significant overdensity remains, even in view of the presence of a recent burst of star cluster formation in many of these galaxies.

These observations have prompted a flurry of activity in the area of star cluster disruption processes. This has led to suggestions that cluster systems appear to be affected by a disruption mechanism that acts on very short time-scales ($\lesssim 10-30$ Myr) and which may be mass-independent – at least for masses in excess of $\sim 10^4$ M$_\odot$ (e.g., Fall et al. 2005; Bastian et al. 2005; Fall 2006). This fast disruption mechanism, which is thought to effectively remove around 50 (Goodwin & Bastian 2006; although their sample is very probably biased), 70 (Bastian et al. 2005; Mengel et al. 2005) or even 90 per cent (Lada & Lada 1991; Whitmore 2004; Whitmore et al. 2007) of the youngest clusters from a given cluster population, is thought to be the rapid removal of the intracluster gas on a time-scale of $\sim 5$ Myr, the signatures of which have been seen in several clusters (Bastian & Goodwin 2006). The observational effect resulting from this rapid gas removal has been coined cluster "infant mortality" (Lada & Lada 2003); it was originally reported in the context of the number of very young embedded clusters, compared to their older, largely gas-free counterparts in the Milky Way.

The general consensus emerging from recent studies is that rapid gas removal from young star clusters is likely to leave the clusters super-virial and hence lead to the rapid disruption of many clusters (see, e.g., Goodwin 1997a,b; Boily & Kroupa 2003a,b; Goodwin & Bastian 2006; see also de Grijs & Parmentier 2007 for a review). This leaves surviving clusters more susceptible to destruction (Vesperini & Zepf 2003; Bastian et al. 2005; Fall et al. 2005).

As described by Goodwin & Bastian (2006, and references therein) the effect of gas removal is to rapidly decrease the potential well in which the stars reside. The cluster will expand in an attempt to return to virial equilibrium. If the virial ratio of the stars *after* gas expulsion is $\ga 3$ the cluster will be unable to return to an equilibirum and will be destroyed[^3]. Clusters with a higher effective star-formation efficiency than around $30$ per cent will survive, but may undergo significant "infant weightloss" (loosing in excess of 50 per cent of their initial mass in some cases). The signature of infant weightloss has been observed in several young clusters (Bastian & Goodwin 2006). The time-scale over which a cluster will be destroyed, or attain a new (lower-mass) equilibrium configuration is $10 - 40$ Myr (depending on the effective star-formation efficiency and the cluster mass).

The $10 - 40$ Myr time-scale of gas removal-induced infant mortality and infant weightloss has important consequences for the analysis and interpretation of the data in this paper. Clusters undergoing expansion will have decreasing surface brightnesses, thus reducing their chances of being detected as they grow older. However, some clusters will recollapse after $10 - 40$ Myr, which may bring them back into the sample. In addition, the speed at which clusters are lost from the sample would be expected to depend on their (initial) mass. Lower-mass clusters which are initially only just above the detection limit will drop out of the sample very quickly, whilst larger clusters may remain in the sample (albeit with a lower surface brightness) for longer. It is also almost impossible (without extensive observations of the surface brightness profiles and/or dynamical state of the cluster) to determine which clusters that are present in this age range will survive gas expulsion, and which are headed for destruction. Thus, the interpretation of the numbers and mass function of clusters in the age range $10 - 40$ Myr is fraught with problems (in addition to these "physical" problems, age determinations for clusters in this age range also cause problems; see below).

# A homogeneous photometric database

The basis for our detailed re-analysis of the SMC star cluster system is provided by the $UBVR$ broad-band spectral energy distributions (SEDs) of Hunter et al. (2003), based on Massey's (2002) CCD survey of the Magellanic Clouds.

In a series of recent papers, we developed a sophisticated tool for star cluster analysis based on broad-band SEDs, AnalySED, which we tested extensively both internally (de Grijs et al. 2003a,b; Anders et al. 2004) and externally (de Grijs et al. 2005), using both theoretical and observed young to intermediate-age ($\lesssim 3 \times
10^9$ yr) star cluster SEDs, and the galev "simple" stellar population (SSP) models (Kurth et al. 1999; Schulz et al. 2002). The accuracy has been further increased for younger ages by the inclusion of an extensive set of nebular emission lines, as well as gaseous continuum emission (Anders & Fritze-v. Alvensleben 2003). We concluded that the *relative* ages and masses within a given cluster system can be determined to a very high accuracy, depending on the specific combination of passbands used (Anders et al. 2004). Even when comparing the results of different groups using the same data set, we can retrieve any prominent features in the cluster age and mass distributions to within $\Delta \langle \log( {\rm Age / yr} )
\rangle \le 0.35$ and $\Delta \langle \log( M_{\rm cl} / {\rm M}_\odot
) \rangle \le 0.14$, respectively (de Grijs et al. 2005), which confirms that we understand the uncertainties associated with the use of our AnalySED tool to a very high degree.

In de Grijs & Anders (2006) we presented newly and homogeneously redetermined age and mass estimates for the entire Large Magellanic Cloud (LMC) star cluster sample covered by the Massey (2002) data. Based on the comparison of our results in de Grijs & Anders (2006) with those published previously in a range of independent studies (mostly based on spectroscopic or isochrone analyses), and additionally on a detailed assessment of the age-metallicity and age-extinction degeneracies, we concluded that our broad-band SED fits yield reliable ages, with statistical *absolute* uncertainties within $\Delta\log( \mbox{Age/yr}) \simeq 0.4$ overall. Here, we extend this to the SMC cluster sample, using the same age-dating technique as described above.

Our cluster age and mass determinations assume an average metallicity of $Z = 0.008$ (where Z$_\odot = 0.020$), and a mean foreground extinction $E(B-V) = 0.08$ mag. We will justify both of these choices below. In Fig. we show the distribution of our SMC cluster sample in the log(age) vs. log(mass) plane; the adopted 50 per cent completeness limit is overplotted. We have also indicated the regions in this plane from which we have drawn statistically complete subsamples, which we will discuss in detail in Sect. .

The determination of the 50 per cent completeness limit of the SMC cluster data is in essence based on a close inspection of the cluster photometry contained in Hunter et al.'s (2003) fig. 11. These authors selected their sample from the catalogue of Pietrzyński et al. (1998), matched to the observational field of view of the Massey (2002) data. Therefore, our completeness is that of this catalogue; Hunter and her team did not quantify the completeness levels themselves (D. Hunter, priv. comm.), although they discuss an observed fading limit. However, for our analysis it is important to understand the sample incompleteness affecting our observations. As such, we adopted the conservative approach that the present-day SMC cluster luminosity function (CLF; see Hunter et al.'s fig. 11) is best represented by a power-law function in luminosity. Based on this assumption, we used the same observational data as used by Hunter et al. (2003) to determine the 50 per cent completeness limit at $M_V
\sim -4.5 \pm 0.2$ mag (based on a power-law fit to the clusters brighter than $M_V = -5$ mag; varying this lower limit by a few tenths of a magnitude does not result in a significant change), i.e., at the same level as Hunter et al.'s (2003) observed fading limit. We note that if the underlying CLF is *not* a single power law over the entire observed luminosity range, the limit we adopt following this approach is in fact a lower limit. In the latter case the observations will likely be more complete than estimated here. Since the adopted 50 per cent completeness limit describes the lower envelope of the distribution of our SMC cluster sample very well, we are confident that our approach is reasonable. In addition, we point out that a variation in the magnitude limit of 0.2–0.4 mag will shift the mass limit by at most 0.1–0.2 dex, which clearly is still within our range of uncertainties. The magnitude of the shift expected when going from the 50 to the 90 per cent completeness limit in the SMC disc is of the same order, $\sim 0.5$ mag. Finally, we point out that it is most likely that the completeness of our SMC cluster sample is in fact determined by the $U$-band observations. From a direct comparison of the $U$-band and the $V$-band data, we derive a 50 per cent completeness in the $U$ band at $M_U = -5.0 \pm 0.3$ mag.

Chiosi et al. (2006) recently analysed the star-formation history in the SMC in detail using an independently selected star cluster sample. Where possible, they derive the extinction towards individual clusters based on colour-magnitude diagram analysis, and otherwise assume a mean extinction $E(B-V) = 0.08$ mag, following Tumlinson et al. (2002) and Hunter et al. (2003; see also Rafelski & Zaritsky 2005). We have adopted the same average extinction value to our SMC cluster sample, using the Calzetti attenuation law (Calzetti 1997, 2001; Calzetti et al. 2000; Leitherer et al. 2002) with $R_V =
4.05$.

Rafelski & Zaritsky (2005) obtained SMC cluster ages of a small cluster sample on the basis of three sets of models, for metallicities of $Z=0.001, 0.004$ and 0.008. They concluded that some of the lowest-metallicity models could be rejected and adopted $Z=0.008$ as an appropriate mean metallicity for their SMC cluster sample. Chiosi et al. (2006) also adopted this metallicity, but for their younger sample clusters. For the older ($\gtrsim 1-2$ Gyr) clusters, they assumed $Z=0.004$, as was also done by Hunter et al. (2003). However, as shown in Fig. , the large majority of our sample clusters (and in particular the subpopulations we will analyse in more detail below) are younger than $\sim 1$ Gyr, so that we adopt $Z=0.008$ as the mean metallicity for our SMC cluster sample.

# The cluster mass function

In de Grijs & Anders (2006) we presented the cumulative CMFs of the LMC star clusters younger than certain age limits. We concluded that while the older cluster (sub)samples are characterised by CMF slopes consistent with the $\alpha \simeq 2$ slopes generally observed in young star cluster systems – where $\alpha$ is defined as $N(M_{\rm
cl}) \propto M_{\rm cl}^{-\alpha}$ – the youngest mass and luminosity-limited LMC cluster subsets show shallower slopes (at least below masses of a few $\times 10^3$ M$_\odot$). We noted that we could not disentangle the unbound from the bound clusters at the youngest ages. This is what we set out to do here for the SMC cluster system.

In order to achieve this goal, we present the CMFs for subsets of our SMC cluster sample in Fig. , where the cluster subsamples were selected based on their age distributions. The age limits used to generate the different panels in Fig. are shown as the vertical dash-dotted lines in Fig. ; the various subsets are also cross-linked between the figures.

A closer look at Fig. reveals that, because of the variation in the observational detection limit as a function of age, the lower-mass limits of our cluster subsamples differ. Thus, the CMFs presented in Fig. are statistically complete above different mass limits, as indicated by the horizontal dash-dotted lines in Fig. and the vertical dotted lines in Fig. .

In all panels of Fig. , we have overplotted CMFs with the canonical slope of $\alpha = 2$ (corresponding to a slope of $-1$ in units of ${\rm d} \log(M_{\rm cl}/{\rm M}_\odot) / {\rm d}
\log (N_{\rm cl})$, used in these panels). We have only shifted and scaled these lines vertically, as justified below.

*We emphasise that for the star cluster infant mortality study performed here, we need to choose the age ranges of our cluster subsamples carefully*, for both physical reasons and also because of the discrete nature of the model isochrones. Regarding the latter, it is well known that broad-band SED fitting results in artefacts in the cluster age distribution. This is predominantly caused by specific features in the SSP models, such as the onset and presence of red giant branch or asymptotic giant branch (AGB) stars at, respectively, $\sim 10$ and $\sim 100$ Myr (e.g., Bastian et al. 2005). Alternatively, both the age-metallicity and the age-extinction degeneracies will affect the resulting cluster age distributions, thus also leading to artefacts in the data (e.g., de Grijs et al. 2003b; Anders et al. 2004). We have attempted to avoid placing our age range boundaries around ages (and, where possible, have taken account of the uncertainties in age in doing so) where the effects of such artefacts might seriously impede the interpretation of the results. For instance, one can see a clear artefact in the cluster age distribution (which we will refer to as a "chimney") at $\log( {\rm Age / yr} )
\simeq 7.2 (\simeq 16$ Myr); the average uncertainties for these ages are of order a few Myr, so that we decided to limit our youngest cluster subsample to clusters younger than 10 Myr. If, instead, we had adopted an age limit at $\log( {\rm Age / yr} ) = 7.17 (15$ Myr), we would have had marginally better statistics, but our analysis would be affected by the unknown effects of the age uncertainties associated with this chimney (see Goodwin et al., in prep., for a detailed discussion of the issues involved).

The rationale for adopting as our youngest subsample all clusters with ages $\le 10$ Myr is that at these young ages, the vast majority of the star clusters present will still be detectable, even in the presence of early gas expulsion (e.g., Goodwin & Bastian 2006) – as long as they are optically conspicuous. The CMF of this subsample is shown in Fig. a.

Fig. c includes our sample clusters with ages in excess of 40 Myr, up to 160 Myr. While the upper age limit ensures the full inclusion of the clusters affected by the onset of the AGB stage, its exact value is rather unimportant for our analysis, and it was mainly determined by the need to have reasonable statistics in this (and the upper) age range in Fig. d. The lower age limit of this subsample is crucial, however. As shown by Goodwin & Bastian (2006), most dissolving clusters will have dispersed by an age of $\sim 30$ Myr, while the surviving clusters will have returned to an equilibrium state by $\sim 40$ Myr, when some of the early expansion will have been reversed, depending on the effective star-formation efficiency. This latter age is therefore a good lower boundary to assess the surviving star cluster population.

We explicitly exclude any star clusters aged between 10 and 40 Myr from our analysis. In this age range, which is shown in Fig. b for completeness, it is likely that dissolving star clusters that will not survive beyond about 30–40 Myr might still be detectable and therefore possibly contaminate our sample. In addition, this is the age range in which early gas expulsion causes rapid cluster expansion, before settling back into equilibrium at smaller radii; because of the expanded nature of at least part of the cluster sample, we might not be able to detect some of the lower-luminosity (and hence lower-mass) clusters that would again show up beyond an age of $\sim 40$ Myr. At the same time, the effects of "infant weightloss" (Weidner et al. 2007) will further confuse the analysis in this age range (see Section for details).

# Is there evidence for cluster infant mortality in the SMC?

## Young and intermediate-age clusters

In Fig. a, we have included the best-fitting CMF slope (dash-dotted line), in addition to the canonical $\alpha = 2$ CMF slope (dashed line). Both slopes are the same, within the uncertainties. This also shows that the SMC's CMF at the youngest ages is consistent with an $\alpha = 2$ power law down to cluster masses of $\sim 125$ M$_\odot$, within the (Poissonian) uncertainties.

In the simplest case, in which the cluster formation rate has remained roughly constant throughout the SMC's evolution (see, e.g., Boutloukos & Lamers 2003, their fig. 10; see also Gieles et al. 2007), the number of clusters would simply scale with the age range covered. In Figs. b, c and d, we show the canonical $\alpha=2$ CMF scaled from the best-fitting locus in Fig. a by the difference in (linear) age range between the panels. The main uncertainties introduced by this method are (i) fluctuations caused by small-number statistics in the youngest age range (the effects of which will be enhanced when scaling the young-age CMF to a greater age range), and (ii) the exact length of the youngest age range (especially considering the necessarily short extent of our youngest age bin). While our galev SSP models start at an age of 4 Myr, the actual ages of a small subset of our sample clusters might be somewhat younger. This introduces an artificial concentration of clusters at our youngest model age, as can be seen in Fig. . It is, unfortunately, not straightforward to remedy this situation based on broad-band imaging observations alone. However, we note that it may take up to 3–5 Myr for an embedded cluster to clear a cavity in its natal gas cloud for its stars to become visible at optical wavelengths (see, e.g., Plante & Sauvage 2002 for a review of embedded young massive star cluster observations). Therefore, we adopt the conservative working assumption that our youngest age range runs from 3–10 Myr. The scaling from the youngest age bin to that covering \[40,160\] Myr (Fig. c, i.e., the most important age range for our CMF comparison in the context of our infant mortality analysis) is therefore a factor of $\sim 17$, or $\Delta \log (N_{\rm cl}) =
1.234$.

Despite the caveat regarding the absence of embedded star clusters in our youngest subsample, we argue that this has a negligible effect on the CMF presented in Fig. a, because of their very small number. Recent, homogeneous observations of the SMC using the *Spitzer Space Telescope* in a number of mid-infrared passbands (Bolatto et al. 2007) have shown that the vast majority of the embedded sources are low-mass ($\ll 100$ M$_\odot$) young stellar objects rather than more massive clusters and associations (we point out that for the comparison done here, we are mostly interested in clusters with masses greater than $10^3$ M$_\odot$). The possible exceptions to this rule are few, and include the four youngest $\sim
3$ Myr-old SMC clusters, NGC 299, NGC 346, NGC 376, and NGC 602 (e.g., Sabbi et al. 2007).

The scaled canonical CMF in Fig. c is an almost perfect fit to the observational CMF. Although the best-fitting CMF slope is ${\rm d} \log(M_{\rm cl}/{\rm M}_\odot) / {\rm d} \log
(N_{\rm cl}) = -0.82 \pm 0.18$, this compares to ${\rm d} \log(M_{\rm
cl}/{\rm M}_\odot) / {\rm d} \log (N_{\rm cl}) = -1.01 \pm 0.20$ if we ignore the lowest-mass clusters at $\log(M_{\rm cl}/{\rm M}_\odot) \le
3.2$, where there may be residual incompleteness effects (see the selection area in Fig. compared to the age-dependent detection limit).

This very good match between the observed CMF for the age range from 40–160 Myr (Fig. c) and the scaled CMF from Fig. a implies that *the SMC cluster system has not been affected by any significant amount of cluster infant mortality for cluster masses greater than a few* $\times 10^3$ M$_\odot$. Based on a detailed assessment of the uncertainties in both the CMFs and the age range covered by our youngest subsample, we can limit the extent of infant mortality between the youngest and the intermediate age range to a maximum of $\lesssim 30$ per cent ($1\sigma$). We rule out a $\sim 90$ per cent (infant) mortality rate per decade of age at a $>6 \sigma$ level. This result is in excellent agreement with that of Gieles et al. (2007); it is, however, in direct contradiction to the claim of Chandar et al. (2006) that the SMC cluster system has undergone sustained destruction at very high rates (up to 90 per cent per decade in logarithmic age) for the full age range up to $\sim 1$ Gyr, although we note that Chandar et al. (2006) do not include the youngest SMC clusters in their analysis.

As an important caveat, we remind the reader that the main underlying assumption leading to this result is the notion that the SMC's CFR has been approximately constant over the time-scale of $\sim 1$ Gyr. If this were seriously in error, in order for this to give rise to the result reported here, the SMC's average CFR must have been significantly enhanced in the $40 - 160$ Myr-old age range, by up to an order of magnitude, compared to that at present. There is no clear evidence, in either the current data set or the earlier work by Boutloukos & Lamers (2003; see also Gieles et al. 2007), nor in the age distribution of the field stars (Chandar et al. 2006; Chiosi & Vallenari 2007), to suggest that this is the case. In fact, if anything, we might expect an enhanced CFR around the time of the last close encounter between the SMC and the LMC, some $200 - 500$ Myr ago (see, e.g., Heller & Rohlfs 1994; Gardiner & Noguchi 1996; see also Chiosi et al. 2006, but see Chiosi & Vallenari 2007 for an alternative interpretation), i.e., significantly longer ago than the age range probed by our intermediate-age clusters in Fig. b.

For completeness, we also show the best-fit power-law CMF, as well as the scaled canonical CMF, in Fig. b for the age range between 10 and 40 Myr. Although *we strongly caution against placing too much significance on the analysis of the CMF in this age range*, for the reasons outlined in Section , it is interesting to note that the scaled canonical CMF does in fact seem to describe the extremities of this CMF reasonably well. However, at intermediate masses (a few $\times 10^3 - 10^4$ M$_\odot$) the observed CMF exceeds the theoretical prediction for a constant cluster formation rate. Although the effects of cluster expansion and infant weightloss most likely contribute to confusing the emerging picture, the main cause of this discrepancy is owing to the artificial chimney at $\log( {\rm Age / yr} ) \simeq 7.2$. The cluster numbers in this age range are dominated by this artefact, so it is important that we understand in which sense this affects our results. Because of the discreteness of the isochrones in our SSP models around this age, and the tendency for the broad-band fitting routine to iterate to a local minimum $\chi^2$ solution, most (but not all) of the clusters in this chimney should have been assigned somewhat greater ages. Because their ages have been underestimated (by up to 0.1–0.2 dex in logarithmic age), the associated masses have also been underestimated, by a similar amount. Unfortunately, until more detailed SSP models become available, there is no easy way out. However, a qualitative assessment suggests that if we were able to correct for this chimney, the derived masses of at least a fraction of the clusters affected would be greater, and thus that the apparent excess in Fig. b would be redistributed towards greater masses. The result would be a smoother CMF, more akin to the scaled canonical $\alpha=2$ CMF of Fig. a.

The alternative interpretation, i.e., that the star cluster formation rate in the SMC has undergone a significant increase in the age range between 10 and 40 Myr appears to be effectively ruled out by the complementary analysis of Chiosi et al. (2006). In fact, these authors find evidence suggesting the contrary, i.e., that the SMC cluster population has seen periods of enhancement during the first $\sim 15$ Myr, and at around 90 Myr – this implies that in the age range of interest here, the cluster formation rate found by Chiosi et al. (2006) was in fact *reduced* with respect to the most recent cluster forming episode (roughly equivalent to the cluster sample shown in Fig. a). Equivalently, neither Rafelski & Zaritsky (2005), nor either Chandar et al. (2006) or Gieles et al. (2007) find an episode of increased cluster formation at a few $\times 10^7$ yr, despite the fundamentally different conclusions drawn by the authors of the latter two studies.

## The oldest sample clusters

Finally, in Fig. d we show the combined SMC CMF for clusters from 160 Myr up to 1.0 Gyr, as well as the scaled canonical CMF. The latter matches the highest-mass ($\log(M_{\rm cl}/{\rm
M}_\odot) \ga 4.8$) part of the observed CMF, but significantly overpredicts the number of lower-mass ($\log(M_{\rm cl}/{\rm M}_\odot)
\lesssim 4.8$) clusters. This flattening of the CMF with respect to our intermediate age range (Fig. c) evidences the increased importance of mass-dependent cluster disruption, as described in detail by, e.g., Lamers et al. (2005; see specifically their fig. 4). However, we note that the calculations of Lamers et al. (2005) are based on coeval stellar populations, while our oldest age bin contains a mixture of differently aged clusters. Nevertheless, it is apparent from fig. 4 of Lamers et al. (2005) that the flattening of the CMF increases as a function of age, starting at early ages. In a mixed intermediate-age population, the individual roughly coeval subpopulations giving rise to the overall CMF therefore lead to a flattened CMF with respect to the initial cluster mass distribution (which we have shown in the SMC to be equivalent to the canonical $\alpha=2$ power-law CMF; see Fig. a). A comparison of our observations with the Lamers et al. (2005) calculations is therefore relevant to first order. A similar flattening of the CMF with increasing age is also predicted as owing to the effects of the underlying galactic tidal field (e.g., Vesperini & Heggie 1997; Baumgardt & Makino 2003). We will discuss these older clusters in more detail in a follow-up paper (Goodwin et al., in prep.), in which we will discuss the entire evolutionary sequence of the star cluster systems of both Magellanic Clouds, and where we aim to understand the *physics* driving the early evolution of these star clusters systems.

## Additional supporting evidence for little early disruption

In Fig. we show the SMC cluster age distribution expressed in units of ${\rm d}N_{\rm cl} / {\rm d}t$, i.e., the number of cluster per unit time-scale (for which we adopt $10^6$ yr). To first order, our age distribution is similar to that based on the Rafelski & Zaritsky (2005) data used by both Chandar et al. (2006) and Gieles et al. (2007). It is also roughly similar to the age distribution derived independently by Chiosi et al. (2006).

Gieles et al. (2007) analysed their cluster age distribution very carefully, and found little evidence for infant mortality in the SMC cluster system. They showed that the decline in ${\rm d}N_{\rm cl} /
{\rm d}t$ in their sample could be attributed entirely to evolutionary fading of the cluster population, irrespective of which SSP models are used. They derived that for a magnitude-limited sample, as also discussed in this paper, the decline in ${\rm d}N_{\rm cl} / {\rm d}t$ as a function of age is graphically described by a slope of $-0.72$, if the cluster ages are based on the galev SSP models. In Fig. , we show the expected effects of evolutionary fading of the cluster population as the dashed arrow. It is immediately clear that the decline in the age distribution up to $\log( \mbox{Age yr}^{-1} ) \simeq 7.8$ can indeed be entirely attributed to evolutionary fading. This is supported by the mass-limited subsamples shown in Fig. : for all mass ranges, they show an essentially flat age distribution up to $\sim
10^8$ yr. This is in support of the results of both Gieles et al. (2007), and Chandar et al. (2006; their fig. 1), although the latter authors favoured a different interpretation.

The age distribution at $\log( \mbox{Age yr}^{-1} ) \simeq 8.2 \pm
0.3$ (covering the age range from about 80 to 320 Myr) falls below the fading line, however. For a constant cluster formation rate over this entire period, and if we normalise the age distribution at our youngest data point, we would need the SMC cluster sample to have suffered from $\sim 20-50$ per cent disruption in order to match the observations. We can firmly rule out a constant $\sim 90$ per cent disruption rate per decade in age, up to an age of 1 Gyr. The expected effects of evolutionary fading combined with a 90 per cent disruption rate are shown as the dash-dotted arrow in Fig. . The arrow clearly does not fit the observed age distribution, if we require it to pass through our youngest data point. We note, however, that the slope of this latter arrow is very similar to that of the age distribution of the full SMC sample for ages in excess of a few $\times 10^8$ yr, when secular disruption is likely to take over. We will discuss these older age ranges in detail in Goodwin et al. (in prep.).

In summary, we set out to shed light on the controversy surrounding the early evolution and disruption of star clusters in the SMC. We embarked on a fresh approach to the problem, using an independent, homogeneous data set of $UBVR$ imaging observations, from which we obtained the cluster age distribution in a self-consistent manner. We conclude that the optically selected SMC star cluster population has undergone at most $\sim 30$ per cent ($1\sigma$) "infant mortality". Using the age distribution of the SMC cluster sample in units of the number of clusters observed per unit time-scale, we independently confirm this scenario. Gieles et al. (2007) reached a similar conclusion.

# acknowledgments

We thank Peter Anders for performing the SED cluster fits using his AnalySED tool. We are also grateful to the anonymous referee for a careful reading of this paper and for suggesting useful additional lines of investigation. RdG acknowledges the discussions at the ESO workshop on "12 Questions in Star and Massive Star Cluster Formation" organised by Markus Kissler-Patig and Tom Wilson. RdG and SPG also acknowledge partial financial support from the Royal Society in the form of a UK-China International Joint Project, as the revised version of this paper was prepared during a collaborative visit to China. We are grateful to an anonymous independent referee for suggesting a number of improvements to the paper that strengthened the main results. This research has made use of NASA's Astrophysics Data System Abstract Service.

[^1]: E-mail: email@example.com

[^2]: Although Chandar et al. (2006) suggest that their sample is roughly mass limited, they also note that the mass-limited subsample, constrained to clusters with masses $\log( M_{\rm cl} / {\rm M}_\odot ) \ge 3.5$, shows a flatter age distribution.

[^3]: This is usually given in terms of a star-formation efficiency of $\sim 30$ per cent. However, as noted by Goodwin & Bastian (2006), it is the virial ratio of the stars after gas expulsion that is the crucial parameter, and this can only be related to the star-formation efficiency in a simple way if the stars and gas were initially in virial equilibrium with one another. For this reason, Goodwin & Bastian (2006) use the term "effective star formation efficiency".


Astronomy Letters, Vol. 32, No. 5, pp. 302-307, 2006

**Interacting Galaxies in Sparsely Populated Groups**
O. V. Melnyk [^1]
*Astronomical Observatory, Kiev National University, 3 Observatorna str., Kiev, 04053 Ukraine
email@example.com*

**Abstract**
We calculated the median values of the following parameters for 87 groups of galaxies with three to eight components: the mean rms velocity of the galaxies in the group, $S_{v}$ = 166 km/s; the harmonic mean radius, $R_{h}$ = 29 kpc; and the mass-to-light ratio, $M_{vir}/L = 33 M_{\odot}/L_{\odot}$. The $M_{vir}/L$ ratio depends on the population of the system, while $S_{v}$ does not depend on $R_{h}$. To ascertain the relationship between the activity of galaxies in groups and their morphological composition and the effect of the environment on the evolutionary processes in groups, we consider the fraction of galaxies with UV excess in the sample of interacting galaxies in groups (6%), single peculiar galaxies (8%), and isolated galaxies (4%) and their morphological composition. We also show that the number of active galaxies decreases with increasing population of the group of galaxies, while the frequency of occurrence of early-type ($E/S0$) galaxies increases. **Key words : galaxies, groups of galaxies, interacting galaxies.**

# Introduction

Vorontsov-Velyaminov, Arp, and other authors began an active discussion of the nature of interacting galaxies in the 1950s. The definition of interacting galaxies as objects embedded in a common luminous shell, interpenetrating, and with distortions of their shapes (asymmetry, tails, and bridges) belongs to Vorontsov-Velyaminov (1958). Photographic, spectroscopic, and photometric observations of interacting galaxies were carried out in the 1970s at the Special Astrophysical Observatory of the Russian Academy of Sciences and the Crimean and Byurakan Astrophysical Observatories. Interacting galaxies that form groups or are components of compact groups of galaxies (see, e.g., Karachentseva and Karachentsev 2000, Karachentseva et al. 1987, Hickson et al. 1989) clusters of galaxies, are encountered; there are also isolated peculiar objects among interacting galaxies (Karachentseva 1973). Compact Hickson groups (Hickson et al. 1989) that contain a significant fraction of interacting galaxies enjoy the greatest popularity from both observational and theoretical points of view.

Studying the relationship between the activity of galaxies in groups and their morphological composition is important for understanding the evolutionary processes in close groups. Comparison of the kinematic parameters and mass-to-light ratios for groups containing different numbers of galaxies with those for isolated galaxies provides no less important information. To this end, we compare the occurrence of active galaxies with UV excess in groups with different numbers of members containing interacting galaxies from the catalog by Vorontsov-Velyaminov et al. (2001) and in the sample of isolated galaxies by Karachentseva (1973). In this paper, we also consider observational, kinematic, and virial properties of these systems.

The catalog of interacting galaxies (Vorontsov-Velyaminov 2001) is a combination and continuation of the two previous parts of the catalog published by Vorontsov-Velyaminov in 1959 and 1977. The catalog contains 2014 interacting galaxies and is homogeneous to $15^{m}$ and $\delta > -45^{\circ}$. However, it also includes faint galaxies (to $20^{m}$); the radial velocities of the objects reach 27 000 km/s. Among the 2014 interacting galaxies of the catalog, 446 (22%) and 73 (3%) were classified as nests N and chains Ch, respectively; pairs P are most numerous, 1230 (61%); and the remaining 14% are peculiar galaxies or M51-type interacting galaxies with close companions.

# The Sample

Systems classified by the type of interaction as nests N or chains Ch with three or more members (2/3 and 1/3 of the systems in the sample contain nests and chains of galaxies, respectively) were selected from the catalog of Vorontsov-Velyaminov et al. (2001). Close galaxies, but without any traces of distortions of their shapes were not considered to be interacting galaxies according to the criterion of Vorontsov-Velyaminov (1958) and, accordingly, were not included in the catalog. Therefore, to identify a "bound" group of galaxies, we corrected the number of components in the system by visually checking each system on DSS images and galaxy parameters in the LEDA and NED databases. This selection was necessary, since, according to Zasov and Arkhipova (2000), more than a third of the nests of galaxies identified by Vorontsov-Velyaminov are actually single galaxies with several sites of star formation or pairs. The triplet VV143 (HCG18), which is actually a giant irregular galaxy (Plana et al. 2000), is a dramatic example. Thus, if an interacting system belonged to a more densely populated Hickson group, then the components were included in the group; conversely, if a galaxy with a radial velocity differing significantly from the radial velocities of other components was detected in the system, then it was excluded from the group (as, e.g., in the Stephan Quintet VV288). The sample contains a total of 30 groups identified with compact Hickson groups. The final sample contains 335 galaxies in 87 groups consisting of three to eight components with known observational parameters (radial velocities, sizes, apparent magnitudes, and morphological types) for all components.

# The morphological composition and activity of interacting galaxies

We compared the parameters of interacting galaxies in groups with two samples of galaxies. One of these samples consists of single peculiar galaxies that were classified in the catalog by Vorontsov-Velyaminov (2001) as En (enigmatic, peculiar), K (comet-like), H (with vast HII regions), and R (with rings); nests N or chains Ch, which may be barred spirals or irregular galaxies with large bright knots (Zasov and Arkhipova 2000); and M-galaxies with a satellite in the spiral arms. There are a total of 533 such galaxies in this sample; the radial velocities and apparent magnitudes are known for 471 of them. For comparison, we also selected 1022 from the 1051 galaxies of the catalog of isolated galaxies by Karachentseva (1973) (CIG), which have the status of isolated galaxies in the NED database; the radial velocities, apparent magnitudes, and morphological types are known for 878 of them.

Figures 1 and 2 show the distribution of galaxies in morphological type for the sample of interacting galaxies in the selected groups and the comparison samples and their differential luminosity functions.

The fraction of early-type $E/S0$ (the LEDA classification -5..0) galaxies among the interacting galaxies in groups is 47%. The fractions of early-type galaxies among the single peculiar and isolated CIG galaxies are 7% and 19%, respectively. The sample of single peculiar galaxies has more late-type and irregular galaxies and more low-luminosity galaxies.

The morphological type is luminosity-independent in all three samples, i.e., galaxies of various morphological types are represented in the entire range of luminosities.

Figure 3 shows the distribution of galaxies in ($U - B$) and ($B - V$) color indices (according to LEDA). The number of galaxies with known color indices in each sample is $N_{(U - B)}$ = 30 and $N_{(B - V)}$ = 41 for interacting galaxies in groups, $N_{(U - B)}$ = 77 and $N_{(B - V)}$ = 109 for single peculiar galaxies, and $N_{(U - B)}$ = 66 and $N_{(B - V)}$ = 92 for isolated CIG galaxies. The availability of a color index for a galaxy does not depend on its luminosity and morphological type.

Because of the poor statistics, we cannot analyze in detail the differences in color index, for example, for the same morphological type depending on the galaxy membership in a certain sample. However, it follows from Fig. 3 that the interacting galaxies in groups are redder than the isolated galaxies and considerably redder than the peculiar galaxies. Thus, the $(U - B)$ and $(B - V)$ color distributions for the two samples from the catalog by Vorontsov-Velyaminov (2001) confirm the performed morphological classification: an excess of early-type objects, i.e., galaxies of an older population among the interacting galaxies in groups, and an excess of irregular, blue objects among the single peculiar galaxies.

We then compared the fractions of Markarian and Kazarian galaxies with UV excess (according to LEDA and NED) in the samples under consideration.We found that active galaxies account for 6%, 8%, and 4% of the interacting galaxies in groups, single peculiar galaxies, and the CIG sample, respectively. The fractions of active galaxies among groups with different numbers of members are 9%, 7%, and 3% for triplets, quartets, and groups with five to eight galaxies, respectively. Active galaxies with UV excess are encountered rarely among the early-type objects in all samples, 2% for interacting galaxies in groups and $\sim$ 1% among single-galaxy samples. Therefore, most of the active objects in all samples are Sa..Irr galaxies. Galaxies with emission-line spectra, Seyferts, and LINERs are the dominant types of activity in all samples; galaxies with vast HII regions, active star formation, and absorption-line spectra are also encountered. We found no dominance of these activity types when comparing the isolated and interacting galaxies. The activity of a galaxy can be produced by both internal and external processes: the fall of a companion to the galaxy and tidal phenomena (see Byrd and Valtonen (2001) and references therein). The high percentage of active galaxies among the single peculiar galaxies and their occurrence among the isolated CIG galaxies could be the result of both factors.

According to Karachentsev and Karachentseva (1974) and Byrd and Valtonen (2001), there is an excess of galaxy pairs in which both components are Markarian galaxies. We considered the fractions of active Markarian and Kazarian galaxies in the galaxy triplets of two samples: interacting triplets (ITs) and the combined sample of Northern and Southern triplets (NSTs) (Karachentseva et al. 1987; Karachentseva and Karachetnsev 2000); the rms velocity of the member galaxies is $S_{v} < 300$ km/s for all triplets. Note that the IT groups with $S_{v} > 300$ km/s contain no active galaxies with UV excess.

The columns in Table 1 list the expected (E) (for a random distribution) and observed (O) numbers of triplets, $N_{tr}$, without galaxies with UV excess, with one UV galaxy, with two UV galaxies, and with three UV galaxies for the IT and NST samples.
Table 1. Number of triplets with UV galaxies.

| Sample | $N_{tr}$ | $N_{tr}$ without UV |  | $N_{tr}$ with 1 UV |  | $N_{tr}$ with 2 UV |  | $N_{tr}$ with 3 UV |  |
|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|
| 2-10 | total | E | O | E | O | E | O | E | O |
| IT | 42 | 30 $\pm$ 5 | 32 | 11 $\pm$ 3 | 7 | 1 $\pm$ 1 | 2 | 0 | 1 |
| NST | 86 | 79 $\pm$ 9 | 78 | 10 $\pm$ 3 | 6 | 0 | 1 | 0 | 1 |

It can be seen from Table 1 that the expected number of triplets for a random distribution without active galaxies and with two UV galaxies is smaller than their observed number in both samples, although the difference is within the error limits. The expected number of triplets with one UV galaxy is greater than their observed number, while the expected number of triplets with three UV galaxies is smaller than their observed number. Three galaxies with UV excess have two triplets: VV672 and KTG82 (the VV2002 pair). Thus, among the triplets, as in double galaxies (Karachetnsev and Karachentseva 1974; Byrd and Valtonent 2001), there is an excess of groups where all components are active galaxies with UV excess.

The occurrence of active galaxies in interacting groups is also related to their population: the more components in the group, the fewer galaxies with UV excess it contains (the fraction of active galaxies is only 3% in groups of five to eight members). These data may be indicative of different evolutionary stages of these systems. Coziol et al. (2004) suggested an evolutionary theory of compact groups based on the morphological type and spectral indices of galaxy activity determined for 27 compact Hickson groups of galaxies. The evolved groups, i.e., those in which a merger has already occurred, contain elliptical galaxies that exhibit no activity, while at the first merger stage, there are late-type galaxies with active star formation in the groups.

# Kinematic and virial parameters

We calculated the kinematic and virial parameters for 87 groups in the same way as was done by Karachentseva and Karachentsev (2000) and Vavilova et al. (2005) for groups with different numbers of galaxies. Table 2 lists the median parameters: $n_{m}$ is the number of galaxies in the group; $N_{gr}$ is the number of groups; $N(E/S0)$ is the ratio of the number of early-type galaxies to the total number of galaxies in the group: $N(E/S0)$ = $n(E/S0)/n_{m}$; $\langle V_{LG} \rangle$ are the mean radial velocities of the groups corrected for the solar motion as prescribed by Karachentsev and Makarov (1996), in km/s; $S_{v}$ are the rms velocities of the galaxies, in km/s; $R_{h}$ is the harmonic mean radius of the group, in kpc; $\tau$ is the dimensionless crossing time of the system, in units of the Hubble time 1/$H_{0}$; $M_{vir}$ is the virial mass of the group, in $10^{12}$ $ÌM_{\odot}$; $L$ is the luminosity, in $10^{10}$ $L_{\odot}$; and $M_{vir}/L$ is the virial mass-to-light ratio, in $M_{\odot}/L_{\odot}$.

Table 2. Median kinematic and virial parameters of groups.

|  |  |  |  |  |  |  |  |  |  |
|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|:--:|
| $n_{m}$ | $N_{gr}$ | $N(E/S0)$ | $\langle V_{LG} \rangle$, | $S_{v}$, | $R_{h}$, | $\tau$, | $M_{vir}$, | $L$, | $M_{vir}/L$, |
|  |  |  | km/s | km/s | kpc | 1/$H_{0}$ | $10^{12}$ $M_{\odot}$ | $10^{10}$ $L_{\odot}$ | $M_{\odot}/L_{\odot}$ |
| 3 | 51 | 0.33 | 7291 | 104 | 26 | 0.04 | 0.99 | 6.16 | 19 |
| 4 | 19 | 0.50 | 7265 | 237 | 38 | 0.02 | 5.33 | 8.52 | 53 |
| 5-6 | 10 | 0.78 | 8628 | 385 | 33 | 0.02 | 11.8 | 9.21 | 99 |
| 7-8 | 7 | 0.55 | 9004 | 384 | 60 | 0.02 | 24.6 | 13.3 | 151 |
| 3-8 | 87 | 0.43 | 7673 | 166 | 29 | 0.03 | 2.29 | 7.53 | 33 |

It follows from Table 2 that the velocity dispersion of the group and its mass, luminosity, and virial mass-to-light ratio increase with increasing group population. The mass-to-light ratio $M_{vir}/L$ depends strongly on the rms velocity and does not depend on the harmonic mean radius. Figure 4 shows that the luminosity of a group increases only slightly with $M_{vir}/L$.

Triplets are found in the entire $M_{vir}/L$ range; the spread in $M_{vir}/L$ is greatest for them. Almost all of the triplets located in Fig. 4 to the right in the region of groups consisting of four to eight components have rms velocities higher than 300 km/s. In addition, such high $S_{v}$ are typical of variously populated groups with $M_{vir}/L >$ 40 $M_{\odot}/L_{\odot}$. The similarity of $M_{vir}/L$ and $S_{v}$ for these triplets to the corresponding parameters for more densely populated groups may stem from the fact that these triplets are located in clusters and/or have more distant companions, for example, as in Hickson groups where dwarf galaxies were found at distances of several group radii (Hunsberger et al. 1998).

According to Coziol et al. (2004), the evolution level of a group increases with system mass and with galaxy velocity dispersion in the group and elliptical galaxies dominate in evolved groups. Tovmassian et al. (2004) also showed that the velocity dispersion in a group increases with number of $E/S0$ galaxies, with the frequency of occurrence of early-type galaxies being an increasing function of the group population and, accordingly, its mass. We found a similar dependence for the sample of interacting galaxies in groups, $\log S_{v} = 0.50 N(E/S0) + 1.96$, with the correlation coefficient R = 0.44 and the standard deviation SD = 0.36. However, the rms velocity of the group galaxies depends on the mean radial velocity, as can be seen from the linear regression equation: $\log S_{v} = 0.83 \log\langle V_{LG} \rangle - 1.03$, R = 0.43, and SD = 0.36. Note that these two dependences have equal correlation coefficients. Thus, the morphological classification of galaxies may be subject to distance selection, since the spiral arms in distant galaxies are more difficult to distinguish.

Numerical simulations suggest (see, e.g., Zheng et al. (1993) and references therein) that the product of a galaxy merger is an elliptical galaxy. Observations of elliptical galaxies in compact Hickson groups showed that their properties differ from those of elliptical field galaxies (de Oliveira and Hickson 1994). The frequency of occurrence of $E/S0$ galaxies as a function of the population in interacting groups is 0.43, which is much higher than that in isolated CIG galaxies, 0.19. This result is acceptable for the merger hypothesis, especially if it is considered that these galaxies could have already undergone mergers in the past (Zheng et al. 1993; Hunsberger et al. 1998). This may be evidenced by the observational data on binuclear galaxies in compact groups, which could be the result of previous mergers (Amram et al. 2004; Bentoni and Buson 2000), and the presence of double nuclei in Markarian galaxies and their tendency to be found in pairs (Keel and van Soest 1992).

# Conclusion

Based on the catalog by Vorontsov-Velyaminov (2001), we compiled a sample of 87 interacting galaxies with known radial velocities and with three to eight components. The median parameters for these galaxies are: the rms velocity $S_{v}$ = 166 km/s, the harmonic mean radius $R_{h}$ = 29 kpc, the system crossing time $\tau$ = 0.03 1/$H_{0}$, the virial mass $M_{vir}$ = 2.29 $\times  10^{12}$ $M_{\odot}$, the luminosity $L$ = 7.53 $\times 10^{10}$ $L_{\odot}$, and the mass-to-light ratio $M_{vir}/L$ = 33 $M_{\odot}/L_{\odot}$.

The interacting galaxies in groups and the single peculiar galaxies differ in morphological composition, as confirmed by their color indices: bluer galaxies are typical of the sample of single galaxies, while redder galaxies are typical of the galaxies in groups. The samples of interacting and peculiar galaxies contain about twice as many galaxies with UV than CIG galaxies.

The fraction of active galaxies in interacting groups correlates with the group population: the more components in a group, the fewer galaxies with UV excess it contains (9% in triplets, 7% in quartets, and only 3% active galaxies in groups of five to eight members). The frequency of occurrence of $E/S0$ galaxies in interacting groups is 0.33, 0.50, and 0.67 for triplets, quartets, and groups with five to eight members, respectively. This is much higher than that in isolated CIG objects (0.19) and in peculiar galaxies (0.07), where the contribution of active galaxies is 4% and 8%, respectively. Thus, the contribution of early-type galaxies increases and the number of active galaxies decreases with increasing group population. The velocity dispersion of the group, its mass, luminosity, and virial mass-to-light ratio also increase with increasing group population. An excess of groups where all components are active galaxies with UV excess are observed in triplets. These results may suggest that interacting galaxies in differently populated groups are characterized by a certain morphological composition and are at different evolutionary stages.

**ACKNOWLEDGMENTS.** I am grateful to V.E. Karachentseva and I.B.Vavilova for helpful discussions of the results and valuable remarks that contributed to the improvement of this work. I used the LEDA (http://leda.univ-lyon1.fr), NED (http://nedwww.ipac.caltech.edu), and DSS images (http://archive.eso.org/dss/dss).

[^1]: O. V. Mel'nik