This paper discusses a new method to generate self-coherent initial conditions for young substructured stellar cluster. The expansion of a uniform system allows stellar sub-structures (clumps) to grow from fragmentation modes by adiabatic cooling. We treat the system mass elements as stars, chosen according to a Salpeter mass function, and the time-evolution is performed with a collisional N-body integrator. This procedure allows to create a fully-coherent relation between the clumps' spatial distribution and the underlying velocity field. The cooling is driven by the gravitational field, as in a cosmological Hubble-Lema\^itre flow. The fragmented configuration has a `fractal'-like geometry but with a self-grown velocity field and mass profile. We compare the characteristics of the stellar population in clumps with that obtained from hydrodynamical simulations and find a remarkable correspondence between the two in terms of the stellar content and the degree of spatial mass-segregation. In the fragmented configuration, the IMF power index is ~0.3 lower in clumps in comparison to the field stellar population, in agreement with observations in the Milky Way. We follow in time the dynamical evolution of fully fragmented and sub-virial configurations, and find a soft collapse, leading rapidly to equilibrium (timescale of 1 Myr for a ~ 10^4 Msun system). The low-concentration equilibrium implies that the dynamical evolution including massive stars is less likely to induce direct collisions and the formation of exotic objects. Low-mass stars already ejected from merging clumps are depleted in the end-result stellar clusters, which harbour a top-heavy stellar mass function.

We analyse a hydrodynamical simulation of star formation. Sink particles in the simulations which represent stars show episodic growth, which is presumably accretion from a core that can be regularly replenished in response to the fluctuating conditions in the local environment. The accretion rates follow $\dot{m} \propto m^{2/3}$, as expected from accretion in a gas-dominated potential, but with substantial variations over-laid on this. The growth times follow an exponential distribution which is tapered at long times due to the finite length of the simulation. The initial collapse masses have an approximately lognormal distribution with already an onset of a power-law at large masses. The sink particle mass function can be reproduced with a non-linear stochastic process, with fluctuating accretion rates $\propto m^{2/3}$, a distribution of seed masses and a distribution of growth times. All three factors contribute equally to the form of the final sink mass function. We find that the upper power law tail of the IMF is unrelated to Bondi-Hoyle accretion.

Stars form in regions of very inhomogeneous densities and may have chaotic orbital motions. This leads to a time variation of the accretion rate, which will spread the masses over some mass range. We investigate the mass distribution functions that arise from fluctuating accretion rates in non-linear accretion, $\dot{m} \propto m^{\alpha}$. The distribution functions evolve in time and develop a power law tail attached to a lognormal body, like in numerical simulations of star formation. Small fluctuations may be modelled by a Gaussian and develop a power-law tail $\propto m^{-\alpha}$ at the high-mass side for $\alpha > 1$ and at the low-mass side for $\alpha < 1$. Large fluctuations require that their distribution is strictly positive, for example, lognormal. For positive fluctuations the mass distribution function develops the power-law tail always at the high-mass hand side, independent of $\alpha$ larger or smaller than unity. Furthermore, we discuss Bondi-Hoyle accretion in a supersonically turbulent medium, the range of parameters for which non-linear stochastic growth could shape the stellar initial mass function, as well as the effects of a distribution of initial masses and growth times.

In the scope of the star formation process, it is unclear how the environment shapes the initial mass function (IMF). While observations of open clusters propose a universal picture for the IMF from the substellar domain up to a few solar masses, the young association eta Chamaeleontis presents an apparent lack of low mass objects (m<0.1 Msun). Another unusual feature of this cluster is the absence of wide binaries with a separation > 50 AU. We aim to test whether dynamical evolution alone can reproduce the peculiar properties of the association assuming a universal IMF. We use a pure N-body code to simulate the dynamical evolution of the cluster for 10 Myr, and compare the results with observations. A wide range of values for the initial parameters are tested in order to identify the initial state that would most likely lead to observations. In this context we also investigate the influence of the initial binary population on the dynamics and the possibility of having a discontinuous single IMF near the transition to the brown dwarf regime. We consider as an extreme case an IMF with no low mass systems (m<0.1 Msun). The initial configurations cover a wide range of initial density, from 10^2 to 10^8 stars/pc^3, in virialized, hot and cold dynamical state. We do not find any initial state that would evolve from a universal single IMF to fit the observations. Only when starting with a truncated IMF without any very low mass systems and no wide binaries, can we reproduce the cluster core properties with a success rate of 10% at best. Pure dynamical evolution alone cannot explain the observed properties of eta Cha from universal initial conditions. The lack of brown dwarfs and very low mass stars, and the peculiar binary properties (low binary fraction and lack of wide binaries), are probably the result of the star formation process in this association. (abridged)

We propose a functional form for the IMF, the L_3 IMF, which is a natural heavy-tailed approximation to the log-normal distribution. It is composed of a low-mass power law and a high mass power-law which are smoothly joined together. Three parameters are needed to achieve this. The standard IMFs of Kroupa-2001a, Kroupa-2002a and Chabrier-2003b (single stars or systems) are essentially indistinguishable from this form. Compared to other 3-parameter functions of the IMF, the L_3 IMF has the advantage that the cumulative distribution function and many other characteristic quantities have a closed form, the mass generating function, for example, can be written down explicitly.

Context. Brown dwarfs represent a sizable fraction of the stellar content of our Galaxy and populate the transition between the stellar and planetary mass regime. There is however no agreement on the processes responsible for their formation. Aims. We have conducted a large survey of the young, nearby cluster IC 348, to uncover its low-mass brown dwarf population and study the cluster properties in the substellar regime. Methods. Deep optical and near-IR images taken with MegaCam and WIRCam at the Canada-France-Hawaii Telescope (CFHT) were used to select photometric candidate members. A spectroscopic follow-up of a large fraction of the candidates was conducted to assess their youth and membership. Results. We confirmed spectroscopically 16 new members of the IC 348 cluster, including 13 brown dwarfs, contributing significantly to the substellar census of the cluster, where only 30 brown dwarfs were previously known. Five of the new members have a L0 spectral type, the latest-type objects found to date in this cluster. At 3 Myr, evolutionary models estimate these brown dwarfs to have a mass of ~13 Jupiter masses. Combining the new members with previous census of the cluster, we constructed the IMF complete down to 13 Jupiter masses. Conclusions. The IMF of IC 348 is well fitted by a log-normal function, and we do not see evidence for variations of the mass function down to planetary masses when compared to other young clusters.

The current knowledge on the stellar IMF is documented. It appears to become top-heavy when the star-formation rate density surpasses about 0.1Msun/(yr pc^3) on a pc scale and it may become increasingly bottom-heavy with increasing metallicity and in increasingly massive early-type galaxies. It declines quite steeply below about 0.07Msun with brown dwarfs (BDs) and very low mass stars having their own IMF. The most massive star of mass mmax formed in an embedded cluster with stellar mass Mecl correlates strongly with Mecl being a result of gravitation-driven but resource-limited growth and fragmentation induced starvation. There is no convincing evidence whatsoever that massive stars do form in isolation. Various methods of discretising a stellar population are introduced: optimal sampling leads to a mass distribution that perfectly represents the exact form of the desired IMF and the mmax-to-Mecl relation, while random sampling results in statistical variations of the shape of the IMF. The observed mmax-to-Mecl correlation and the small spread of IMF power-law indices together suggest that optimally sampling the IMF may be the more realistic description of star formation than random sampling from a universal IMF with a constant upper mass limit. Composite populations on galaxy scales, which are formed from many pc scale star formation events, need to be described by the integrated galactic IMF. This IGIMF varies systematically from top-light to top-heavy in dependence of galaxy type and star formation rate, with dramatic implications for theories of galaxy formation and evolution.

We apply two different algorithms to search for mass segregation to a recent observational census of the rho Ophiuchi star forming region. Firstly, we apply the Lambda_MSR method, which compares the minimum spanning tree (MST) of a chosen subset of stars to MSTs of random subsets of stars in the cluster, and determine the mass segregation ratio, Lambda_MSR. Secondly, we apply the m-Sigma method, which calculates the local stellar surface density around each star and determines the statistical significance of the average surface density for a chosen mass bin, compared to the average surface density in the whole cluster. Using both methods, we find no indication of mass segregation (normal or inverse) in the spatial distribution of stars and brown dwarfs in rho Ophiuchi. Although rho Ophiuchi suffers from high visual extinction, we show that a significant mass segregation signature would be detectable, albeit slightly diluted, despite dust obscuration of centrally located massive stars.

The fraction of star formation that results in bound star clusters is influenced by the density spectrum in which stars are formed and by the response of the stellar structure to gas expulsion. We analyse hydrodynamical simulations of turbulent fragmentation in star-forming regions to assess the dynamical properties of the resulting population of stars and (sub)clusters. Stellar subclusters are identified using a minimum spanning tree algorithm. When considering only the gravitational potential of the stars and ignoring the gas, we find that the identified subclusters are close to virial equilibrium (the typical virial ratio Q_vir~0.59, where virial equilibrium would be Q_vir~0.5). This virial state is a consequence of the low gas fractions within the subclusters, caused by the accretion of gas onto the stars and the accretion-induced shrinkage of the subclusters. Because the subclusters are gas-poor, up to a length scale of 0.1-0.2 pc at the end of the simulation, they are only weakly affected by gas expulsion. The fraction of subclusters that reaches the high density required to evolve to a gas-poor state increases with the density of the star-forming region. We extend this argument to star cluster scales, and suggest that the absence of gas indicates that the early disruption of star clusters due to gas expulsion (infant mortality) plays a smaller role than anticipated, and is potentially restricted to star-forming regions with low ambient gas densities. We propose that in dense star-forming regions, the tidal shocking of young star clusters by the surrounding gas clouds could be responsible for the early disruption. This `cruel cradle effect' would work in addition to disruption by gas expulsion. We suggest possible methods to quantify the relative contributions of both mechanisms.

By analysing models of the young massive cluster R136 in 30 Doradus, set-up using the herewith introduced and publicly made available code McLuster, we investigate and compare different methods for detecting and quantifying mass segregation and substructure in non-seeing limited N-body data. For this purpose we generate star cluster models with different degrees of mass segregation and fractal substructure and analyse them. We quantify mass segregation by measuring, from the projected 2d model data, the mass function slope in radial annuli, by looking for colour gradients in radial colour profiles, by measuring Allison's Lambda parameter, and by determining the local stellar surface density around each star. We find that these methods for quantifying mass segregation often produce ambiguous results. Most reliable for detecting mass segregation is the mass function slope method, whereas the colour gradient method is the least practical in an R136-like configuration. The other two methods are more sensitive to low degrees of mass segregation but are computationally much more demanding. We also discuss the effect of binaries on these measures. Moreover, we quantify substructure by looking at the projected radial stellar density profile, by comparing projected azimuthal stellar density profiles, and by determining Cartwright & Whitworth's Q parameter. We find that only high degrees of substructure affect the projected radial density profile, whereas the projected azimuthal density profile is very sensitive to substructure. The Q parameter is also sensitive to substructure but its absolute value shows a dependence on the radial density gradient of the cluster and is strongly influenced by binaries. (abridged)

Recent analyses of mass segregation diagnostics in star forming regions invite a comparison with the output of hydrodynamic simulations of star formation. In this work we investigate the state of mass segregation of 'stars' (i.e. sink particles in the simulations) in the case of hydrodynamical simulations which omit feedback. We first discuss methods to quantify mass segregation in substructured regions, either based on the minimum spanning tree (Allison's Lambda), or through analysis of correlations between stellar mass and local stellar surface number densities. We find that the presence of even a single 'outlier' (i.e. a massive object far from other stars) can cause the Allison Lambda method to describe the system as inversely mass segregated, even where in reality the most massive sink particles are overwhelmingly in the centres of the subclusters. We demonstrate that a variant of the Lambda method is less susceptible to this tendency but also argue for an alternative representation of the data in the plane of stellar mass versus local surface number density. The hydrodynamical simulations show global mass segregation from very early times which continues throughout the simulation, being only mildly influenced during sub-cluster merging. We find that up to approx. 2-3% of the "massive" sink particles (m > 2.5 Msun) are in relative isolation because they have formed there, although other sink particles can form later in their vicinity. Ejections of massive sinks from subclusters do not contribute to the number of isolated massive sink particles, as the gravitational softening in the calculation suppresses this process.

The aim of this work is to test to what extent the star cluster population of a galaxy can be utilised to constrain or estimate the star formation history, with the Large Magellanic Cloud as our testbed. We follow two methods to extract information about the star formation rate from star clusters, either using only the most massive clusters (Maschberger & Kroupa 2007) or using the whole cluster population, albeit this is only possible for a shorter age span. We compare these results with the star formation history derived from colour-magnitude diagrams and find good overall agreement for the most recent approximately 1 Gyr. For later ages, and especially during the "cluster age gap", there is a deficiency of star clusters in relation to the star formation rate derived from the colour-magnitude diagram. The star formation rates following from the whole cluster population lie a factor of approximately 10 lower than the star formation rates deduced from the most massive clusters or from the colour-magnitude diagram, suggesting that only approximately 10 per cent of all stars form in long-lived bound star clusters.

We attempt to evaluate whether the integrated galactic IMF (IGIMF) is expected to be steeper than the IMF within individual clusters through direct evaluation of whether there is a systematic dependence of maximum stellar mass on cluster mass. We show that the result is sensitive to observational selection biases and requires an accurate knowledge of cluster ages, particularly in more populous clusters. At face value there is no compelling evidence for non-random selection of stellar masses in low mass clusters but there is arguably some evidence that the maximum stellar mass is anomalously low (compared with the expectations of random mass selection) in clusters containing more than several thousand stars. Whether or not this effect is then imprinted on the IGIMF then depends on the slope of the cluster mass function. We argue that a more economical approach to the problem would instead involve direct analysis of the upper IMF in clusters using statistical tests for truncation of the mass function. When such an approach is applied to data from hydrodynamic simulations we find evidence for truncated mass functions even in the case of simulations without feedback.

We undertake a systematic analysis of the early (< 0.5 Myr) evolution of clustering and the stellar initial mass function in turbulent fragmentation simulations. These large scale simulations for the first time offer the opportunity for a statistical analysis of IMF variations and correlations between stellar properties and cluster richness. The typical evolutionary scenario involves star formation in small-n clusters which then progressively merge; the first stars to form are seeds of massive stars and achieve a headstart in mass acquisition. These massive seeds end up in the cores of clusters and a large fraction of new stars of lower mass is formed in the outer parts of the clusters. The resulting clusters are therefore mass segregated at an age of 0.5 Myr, although the signature of mass segregation is weakened during mergers. We find that the resulting IMF has a smaller exponent (alpha=1.8-2.2) than the Salpeter value (alpha=2.35). The IMFs in subclusters are truncated at masses only somewhat larger than the most massive stars (which depends on the richness of the cluster) and an universal upper mass limit of 150 Msun is ruled out. We also find that the simulations show signs of the IGIMF effect proposed by Weidner & Kroupa, where the frequency of massive stars is suppressed in the integrated IMF compared to the IMF in individual clusters. We identify clusters through the use of a minimum spanning tree algorithm which allows easy comparison between observational survey data and the predictions of turbulent fragmentation models. In particular we present quantitative predictions regarding properties such as cluster morphology, degree of mass segregation, upper slope of the IMF and the relation between cluster richness and maximum stellar mass. [abridged]

Many objects studied in astronomy follow a power law distribution function, for example the masses of stars or star clusters. A still used method by which such data is analysed is to generate a histogram and fit a straight line to it. The parameters obtained in this way can be severely biased, and the properties of the underlying distribution function, such as its shape or a possible upper limit, are difficult to extract. In this work we review techniques available in the literature and present newly developed (effectively) bias-free estimators for the exponent and the upper limit. The software packages are made available as downloads. Furthermore we discuss various graphical representations of the data and powerful goodness-of-fit tests to assess the validity of a power law for describing the distribution of data. As an example, we apply the presented methods to the data set of massive stars in R136 and the young star clusters in the Large Magellanic Cloud. (abridged)

The study of young stellar populations has revealed that most stars are in binary or higher order multiple systems. In this study the influence on the stellar initial mass function (IMF) of large quantities of unresolved multiple massive stars is investigated by taking into account stellar evolution and photometrically determined system masses. The models where initial masses are derived from the luminosity and colour of unresolved multiple systems show that even under extreme circumstances (100% binaries or higher order multiples) the difference between the power-law index of the mass function of all stars and the observed mass function is small (~0.1). Thus, if the observed IMF has the Salpeter index alpha = 2.35 then the true stellar IMF has an index not flatter than alpha = 2.25. Additionally, unresolved multiple systems may hide between 15 and 60% of the underlying true mass of a star cluster. While already a known result, it is important to point out that the presence of a large number of unresolved binaries amongst pre-main-sequence (PMS) stars induces a significant spread in the measured ages of these stars even if there is none. Also, lower-mass stars in a single-age binary-rich cluster appear older than the massive stars by about 0.6 Myr.

We have made a new compilation of observations of maximum stellar mass versus cluster membership number from the literature, which we analyse for consistency with the predictions of a simple random drawing hypothesis for stellar mass selection in clusters. Previously, Weidner and Kroupa have suggested that the maximum stellar mass is lower, in low mass clusters, than would be expected on the basis of random drawing, and have pointed out that this could have important implications for steepening the integrated initial mass function of the Galaxy (the IGIMF) at high masses. Our compilation demonstrates how the observed distribution in the plane of maximum stellar mass versus membership number is affected by the method of target selection; in particular, rather low n clusters with large maximum stellar masses are abundant in observational datasets that specifically seek clusters in the environs of high mass stars. Although we do not consider our compilation to be either complete or unbiased, we discuss the method by which such data should be statistically analysed. Our very provisional conclusion is that the data is not indicating any striking deviation from the expectations of random drawing.

Most if not all stars form in star clusters. Thus the distribution of star clusters preserves the information on the star formation history of a galaxy. Massive clusters form only during episodes of high star formation activity whereas periods of low star formation activity cannot produce them. We present here the method of Maschberger & Kroupa (2007) to derive the star formation history of a galaxy from its star-cluster content.

Star formation happens in a clustered way which is why the star cluster population of a particular galaxy is closely related to the star formation history of this galaxy. From the probabilistic nature of a mass function follows that the mass of the most-massive cluster of a complete population, M_max, has a distribution with the total mass of the population as a parameter. The total mass of the population is connected to the star formation rate (SFR) by the length of a formation epoch. Since due to evolutionary effects only massive star clusters are observable up to high ages it is convenient to use this M_max(SFR) relation for the reconstruction of a star formation history. The age-distribution of the most-massive clusters can therefore be used to constrain the star formation history of a galaxy. The method, including an assessment of the inherent uncertainties, is introduced with this contribution, while following papers will apply this method to a number of galaxies.