The Complete Initial Mass Function Down to the SubSolar Regime in the Large Magellanic Cloud with Hubble Space Telescope ACS Observations^{1}^{1}affiliation: Based on observations made with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 526555. ^{†}^{†}affiliationmark: ^{2}^{2}affiliation: Research supported by the German Research Foundation (Deutsche Forschungsgemeinschaft) and the German Aerospace Center (Deutsche Zentrum für Luft und Raumfahrt).
Abstract
In this photometric study of the stellar association LH 95 in the Large Magellanic Cloud (LMC) we focus on the premain Sequence (PMS) population in order to construct, for the first time, the subsolar initial mass function (IMF) in the LMC. The basis for this investigation consists of the deepest photometry ever performed in the LMC with the Advanced Camera for Surveys (ACS) onboard the Hubble Space Telescope (HST). We improve our catalog of point sources obtained with ACS (Paper I). We carry out a Monte Carlo technique to subtract the contribution of the general field of LMC. We isolate the central region in the observed area of the association, as the part characterized by the highest concentration of PMS stars. We analyze the reddening distribution of the system in order to obtain the mean value of extinction and we study the mass function of its fieldsubtracted population. For this purpose, we introduce a new set of evolutionary models, derived from the calculations on the evolution of PMS stars by Siess and collaborators. We use these models with our observations of LH 95 to derive the IMF of the system. This mass function is reliably constructed for stars with masses down to M, the lowest mass ever observed within reasonable completeness in the Magellanic Clouds. Consequently, its construction offers an outstanding improvement in our understanding of the lowmass star formation in the LMC. The system IMF of LH 95 shows a definite change in its slope for masses M, where it becomes more shallow. In general, the shape of this IMF agrees very well with a multiple powerlaw, as the typical Galactic IMF, down to the subsolar regime. The change in the slope (“the knee”) of our IMF at M also agrees with the average Galactic IMF. As far as the slope of this system IMF is concerned, it is found to be somewhat more shallow than the corresponding “classical” Galactic IMF in the subsolar regime, probably due to unresolved binarity, while for stars with it becomes slightly steeper. We do not find significant differences in the shape of the overall IMF of LH 95 from that of each of the three individual subclusters of the association. This clearly suggests that the IMF of LH 95 is not subject to local variability.
Subject headings:
stars: formation — stars: premainsequence — open clusters and associations: individual(LH95) — Magellanic Clouds — initial mass function1. Introduction
The Magellanic Clouds, the closest undisrupted dwarf galaxies to our own galaxy, host numerous young stellar associations. The study of the lowmass populations of these systems provide important improvements in understanding extragalactic star formation. The investigation of the star formation processes, and the Initial Mass Function (IMF), together with the time dependency of the formation events, are key points in the characterization of stellar populations in any kind of concentrations, from small clusters to entire galaxies.
Regarding the IMF, several studies have been carried out in the past for the Galaxy (e.g. Salpeter, 1955; Scalo, 1986; Kroupa, 2001; Chabrier, 2003), describing the initial numbers of massive stars ( M), stars of intermediate mass ( M M), lowmass stars ( M M) and brown dwarfs ( M). The findings of such investigations support the hypothesis of the universality of the IMF, in the sense that the average IMF measured in a diversity of galactic systems does not show noticeable variations (Kroupa, 2002). This statement turns out to be rather certain for the intermediate and highmass stars, while for lowmass stars there could be a dependence of the IMF slope on metallicity, in the sense that metalrich environments tend to produce more lowmass stars than metalpoor systems (Piotto & Zoccali, 1999; Reylé & Robin, 2001), although several IMF measurements are needed to confirm such a systematic effect. In lowmass studies there are several uncertainties present. For example the study of the IMF in metalpoor galactic globular clusters (Piotto & Zoccali, 1999) is significantly limited by the corrections for dynamical evolution effects. Specifically, phenomena such as mass segregation and unknown binarity fraction change dramatically the presently measured mass distribution with respect to the one at the time of the formation of these systems (Marks et al., 2008).
Under these circumstances, ideal environments for the search of a metallicity dependence of the IMF are the Magellanic Clouds. Their stellar populations are characterized by metallicities of to times lower than in the disk of the Milky Way (Luck et al., 1998), but have a similar star formation rate (SFR) (Westerlund, 1997). Furthermore, the low metal abundances of the Magellanic Clouds suggest that their environments as well as in other galaxies of the Local Group are closer to the ones at the early ages of the universe, when the peak of star formation occurred (, Pei et al., 1999). Moreover, the lower dusttogas ratio in the Magellanic Clouds (Koornneef, 1982; Bouchet et al., 1985), in comparison with the Milky Way, assures a lower extinction, reducing the bias introduced by differential reddening, a well known limitation in the study of young stellar populations.
R.A.  DECL.  

#  (J2000.0)  (J2000.0)  (mag)  (mag)  (mag)  (mag) 
1  05 37 06.04  66 21 37.15  14.839  0.001  15.066  0.001 
2  05 37 19.26  66 21 07.88  15.956  0.001  15.001  0.001 
3  05 36 59.26  66 21 20.66  15.406  0.001  15.570  0.001 
4  05 37 15.12  66 21 44.39  15.503  0.001  15.494  0.001 
5  05 37 05.62  66 21 35.39  15.372  0.001  15.473  0.001 
6  05 36 49.60  66 23 26.16  16.930  0.003  14.873  0.001 
7  05 37 05.53  66 21 59.65  15.632  0.001  15.740  0.001 
8  05 37 04.53  66 22 01.02  15.670  0.001  15.783  0.002 
9  05 37 14.29  66 22 52.00  15.951  0.001  15.990  0.001 
10  05 36 56.68  66 21 10.40  16.828  0.001  15.102  0.001 
11  05 36 53.49  66 21 55.55  17.183  0.002  15.595  0.001 
12  05 36 59.67  66 22 02.10  16.813  0.001  14.982  0.001 
13  05 37 01.86  66 22 35.65  16.752  0.002  15.663  0.001 
14  05 36 57.29  66 21 49.03  16.267  0.001  16.440  0.001 
15  05 36 53.06  66 22 11.28  16.610  0.001  16.736  0.002 
Note. – Magnitudes are given in the Vega system. Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. The spatial resolution of the ACS/WFC is 0.05 .
Most of the star formation occurs in dense cores of Giant Molecular Clouds, where one or more large stellar concentrations, known as OB associations (Ambartsumian, 1947), are formed. However, most of the star formation studies in the Milky Way are carried out in smaller fields, like the young Taurus, Lupus and Chamaeleon starforming regions. Although it has been shown that galactic OB associations host populations of faint, lowmass PMS stars (Preibisch et al., 2002; Sherry et al., 2004; Briceño et al., 2007) their analysis is strongly limited by the contamination of background and foreground evolved populations, requiring detailed measurements of proper motions, timeconsuming spectroscopy and multiepoch photometry to define their membership. On the other hand, the study of OB associations in the LMC is less subject to contamination by field populations, due to the limited distance spread (Caldwell & Coulson, 1986; Cole, 1998). However, in order to identify the faint red PMS stellar component of LMC associations, photometry with instruments of high sensitivity and angular resolution is required, due to the larger distance from us.
Recent findings from the Hubble Space Telescope (HST) confirm the presence of PMS stars in OB associations of the Magellanic Clouds. Gouliermis et al. (2006a) studied this case in the LMC association LH 52, using WFPC2 observations in the  and equivalent bands, and they discovered a population of objects in the colormagnitude diagram (CMD) consistent with TTauri stars. Subsequently, the Advanced Camera for Surveys (ACS) has been used to perform similar studies in other associations of the Magellanic Clouds, such as NGC 346 (Nota et al., 2006; Gouliermis et al., 2006b; Hennekemper et al., 2008; Sabbi et al., 2008) and NGC 602 (Schmalzl et al., 2008) in the SMC. However, these studies could not characterize the subsolar mass function, because of the insufficient detection of such stars.
In the first part of our study ACS observations enabled us to discover a rich population of PMS stars in the starforming region LH 95 in the LMC (Gouliermis et al., 2007, from hereafter Paper I). In this paper we focus on the IMF of these stars, which is reliable for masses down to about M, considering that our photometry allowed to perform precise completeness corrections for stars with masses down to this limit. In § 2 we describe our photometry and the completeness of our data. In § 3 we apply the decontamination of the observed stellar sample in the area of the association LH 95 from the contribution of the general LMC field for the identification of the stellar members of the system. We also investigate the spatial distribution of the identified premain sequence population, in order to define the limits of the main body of the system, and we apply extinction measurements. In § 4 we present our calculations for new PMS evolutionary models in the observational plane. We make use of these models in combination with our photometry for the construction of the system IMF of LH 95 in § 5. We discuss our findings in § 6 and final conclusions are given in § 7.
2. Photometry
In this study we analyze the photometry presented in Paper I, derived from deep observations with the WideField Channel (WFC) of ACS onboard HST within the GO Program 10566 (PI: D. Gouliermis). Two pointings were observed, one centered on the association LH 95 itself, which we refer to as the “system”, and another one about 2 to the west on an empty area, typical of the local LMC field. We refer to the latter as the “field”. For each pointing observations in two photometric bands were available: F555W and F814W, roughly equivalent to the Johnson V and Cousins I bands respectively. The observations and their photometry are described in detail in Paper I. We found about 16,000 stars in the system, about 2,600 of which lie in the PMS region of the CMD, and about 17,000 in the field.
In the present work we enhance the photometric catalog of stars derived in Paper I by searching for eventually missing objects, especially in the lowmass PMS regime. Photometry was obtained using the ACS module of the package DOLPHOT^{1}^{1}1DOLPHOT is an adaptation of the photometry package HSTphot (Dolphin, 2000). The software and its documentation can be obtained from http://purcell.as.arizona.edu/dolphot/. (Ver. 1.0). The selection of stars in our original photometry was based on the quality parameters estimated for each detected source by the package. However, objects which are actual stars may have been rejected as spurious detections in crowded regions, especially in the presence of a nonuniform background due to nebular emission, as well as in the neighborhood of very bright sources. Therefore, we considered the catalog with all the rejected sources and we performed a visual inspection of each of them on the original ACS FITS images constructed with Multidrizzle. With this process we recovered about 900 additional stars. The updated photometric catalog of the association LH 95 includes in total 17,245 stars. A sample of this catalog is shown in Table 1. The corresponding CMD is shown in Fig. 1. Typical uncertainties of our photometry as a function of the magnitude for both filters are given in Figure 2 (left), for both the “system” and the “field” ACS pointings.
2.1. Completeness
The completeness of the data is evaluated by artificial star experiments with the use of lists of almost 400,000 artificial stars created with the utility acsfakelist of DOLPHOT for each of the observed areas of the system and the field. This utility enables to add stars of given magnitudes to all the frames simultaneously, searching for a detection in none, one, or both photometric bands. In this way it is possible to derive automatically a 2dimensional completeness function for both color () and magnitude or . Considering that at any given faint magnitude in one band the color term covers a range of up to magnitudes, as shown in the CMD of Figure 1, allowing the magnitude in the other band and therefore the completeness to vary significantly, this twodimensional approach is necessary for an accurate treatment of the completeness of our data. Therefore, we define with our artificial star technique the 2D completeness as a function of () and , by computing the ratio of the detected sources over the total added stars in a grid of equally spaced bins magnitudes wide in () and and magnitudes in .
Such completeness maps have been constructed for both the system and the field, enabling, via interpolation, an accurate determination of the actual completeness at any given point of the CMD. The right panel of Figure 2 shows the completeness function in both bands for both the system and the field (top and middle respectively). The completeness is found to vary with location, and therefore different regions in the observed fieldofview are characterized by different completeness functions according to their crowding. In the bottom panel of the figure we show the completeness as a function of magnitude for the central region of the association, as it is defined in § 3.2. It is remarkably lower than that in the entire fieldofview, due to higher crowding of faint stars and overdensity of the bright ones.
To avoid further contamination by spurious objects in our catalog, which could have been included through additional stars, we trimmed the data set used in this work excluding sources with mag or mag. At these magnitude limits stars are faint enough to be well below any reasonable completeness limit. The completeness in the “system” CMD at V= mag and I= mag is less than 1%. Consequently, their removal does not affect the quality and completeness of our data.
3. The PreMain Sequence Population of LH 95
3.1. Field subtraction
The observed stellar population in the area of the association LH 95 is naturally contaminated by the general field population of LMC. Therefore, in order to study the stellar members of the association alone, and in particular its PMS stars, the field contamination should be removed. In our case, having at our disposal photometry in only two bands, the possibility to distinguish if a given star is a true member or not is limited, but a statistical approach can provide accurate results.
We use a Monte Carlo technique, which considers that a star belongs to the cluster or the field with a probability, which depends on the local density in the same part of each CMD (of system and field). In order to maximize the statistical accuracy of the result, two iterations are applied. First, the field subtraction is applied to the entire observed fieldofview centered on the system, enabling to study the spatial distribution of the LH 95 stars and to isolate the boundaries of the association. After that, the study is limited on this subregion of the frame, and another field subtraction is carried out on the central region providing an accurate estimate of the membership of the sources included.
Specifically, for the first iteration we consider each star in the catalog of the system and we select an elliptical region in the CMD of the system centered on the position of this star, and the same region in the CMD of the field. Each such region is being chosen with semiaxes mag and mag. This particular size has been chosen as to be wide enough to include a considerable number of stars, and still small enough to achieve a considerably high resolution in the resulting CMD. We compute the probability of each star in the area of the system to belong to the field as:
(1) 
where and are the numbers of stars included in the elliptical region of the CMD of the system and the field respectively. and are the completeness factors, for system and field, measured for exact position of the considered star in the CMD, as explained in § 2.1. is a normalization factor which takes into account differences in the field stellar density between the two regions. From repeating measurements in different parts of the area of the field and that of the system we found that on average the field includes systematically times more main sequence stars than LH 95. Therefore, considering that both areas cover equal surfaces, we set . According to Eq. 1, when the star is considered to be a field star and it is removed; means that there are no field stars in the field CMD in the elliptical neighborhood of the star and therefore it is considered a member of the association. When the star is randomly kept or rejected as a member star with a probability equal to .
The result of the first fieldsubtraction iteration, applied to the entire LH 95 field of view is shown in Fig. 3, where stars flagged to belong to the association are plotted by black dots. As shown in Paper 1, it is evident that the PMS population belongs entirely to LH 95, while all lower main sequence stars turn out to belong to the LMC field. The upper main sequence (UMS) shows an overabundance of stars in the area of the system, demonstrating the presence of the corresponding population in the association itself.
3.2. Topography of the PMS Stars
In this section we study the spatial distribution of the PMS stars in the area of LH 95, as an indication of the places where star formation occurs. For this purpose, we construct surface density maps of the whole area by counting the stars identified as system members after applying the Monte Carlo removal of the field contamination as described in § 3.1. Since we are interested in the PMS cluster properties, we considered only the PMS population, as selected in the region of the CMD shown in Figure 3.
The star counts are performed by dividing the whole ACS fieldofview into a grid of elements, each corresponding to a size of pixels (pc), and by counting the (fieldsubtracted) PMS stars in each of them. The selection of the specific gridelement size, which corresponds to the “resolving efficiency” of the detection of stellar concentrations, was chosen, after several experiments, as the most appropriate for revealing the smallest concentrations with a physical meaning.
In order to remove from this map eventual noise, we applied a smoothing on the density map with a kernel of size pc pixels. The majority of the PMS stars is found to be located within a compact region , with a density well above the threshold of , where is the standard deviation of the background surface density. The derived twodimensional density map is shown in Figure 4 (left), where the density limit, corresponding to the first (lower) isopleth, is chosen to define the statistically significant concentration of the PMS population of the association LH 95.
RA  DEC  Size  

(J2000.0)  (J2000.0)  (pc)  (pc)  (10 M)  (M pc)  (100 Myr)  (Myr)  (Myr)  
LH 95  5 37 06.50  66 22 03.50  33.5  6.5  2.37  0.47  0.9  6.4  433.0 
A  5 37 15.25  66 21 40.87  8.7  1.7  0.33  1.28  2.4  2.2  26.9 
B  5 36 59.78  66 21 36.88  4.7  0.8  0.19  2.70  5.1  0.9  6.7 
C  5 37 03.96  66 22 09.73  12.9  3.0  0.98  1.11  2.1  3.1  95.5 
From here on, we refer to the region confined by the density isopleth of Figure 4 as the central region of LH 95, and we focus our subsequent analysis on the stellar population included within the boundaries of this region, as the most representative of LH 95. In Figure 4 (left) can be seen that within the central region there are smaller substructures characterized by a higher projected density of PMS stars. We select the three most prominent ones, isolating circular areas within the central region as shown in Figure 4 (right). We refer to these substructures as subcluster A, B and C respectively. The positions and sizes of the selected circular regions around the subclusters, shown in Figure 4 (right), are given in Table 2 with those for the whole central region. We measure the total mass, , included in each stellar concentration, assuming that all systems follow a mass function similar to the standard Galactic field IMF (Kroupa, 2002), and extrapolating their stellar content down to 0.08 M. We then estimate the stellar density, , and the disruption time, , of each cluster following the method by Gouliermis et al. (2002, their section 6). The latter is given as (Spitzer, 1958):
(2) 
Moreover, the dynamical status of a stellar system is defined by two additional timescales, the crossing and the twobody relaxation time (Kroupa, 2008), which are given as:
(3) 
respectively. The threedimensional velocity dispersion of the stars in the cluster, , is given as
(4) 
where is the star formation efficiency (SFE) and the halfmass radius of the cluster. In order to make a rough estimation of the aforementioned timescales for the subclusters and the whole central region of LH 95, we apply the formulas of Eq. 3.
The SFE in several nearby Galactic gasembedded clusters has been found to range typically from 10% to 30% (Lada & Lada, 2003) or 20% to 40% (Kroupa, 2008), value which increases with time while the gas is removed. However, the average age of the system is 4 Myr (see §5), slightly greater than the typical time necessary to remove most of the gas (Lada & Lada, 2003). Therefore, and taking into account the very low optical extinction (see §3.2), LH 95 should be considered as an emerging cluster, meaning that it should be at the process of separation from the parental cloud. Furthermore, according to Wilking & Lada (1983), a high SFE is required for a bound young cluster to emerge from its parental cloud. As a consequence, we consider a value of as more adequate for our system. We derive, thus, values of of the order of 2 km s for our objects and we provide the additional estimated structural parameters for each system also in Table 2. We find that the crossing time for the whole region is greater than the age of the system, and this is quite consistent with the observed subclustering of the system.
From a comparison between the characteristics of the subclusters and the central region, as they are given in Table 2, one can see that there is a significant fraction of mass outside the immediate regions of the subclusters. This can also be seen in the maps of Figure 4, where a prominent population of PMS stars is easily distinguished away from the subclusters, spread in the area among them within the central region. We estimate that the fraction of distributed stellar mass over the total corresponds to 40%, with the remaining 60% being clustered. This fraction of distributed PMS stars in LH 95 is consistent with statistical analyses of several Galactic nearby starforming regions that predict 60% (see review by Allen et al., 2007). However, it is not clear if this distributed population is the result of the star formation or a merging process. Detailed simulations do predict that subclusters may merge to form a larger one (e.g. Fellhauer et al., 2006), but our data do not allow us to verify if this takes place in LH 95. In any case, the stellar density and disruption time of all considered systems, including the whole central region, show that all are rather compact and none of them is under disruption. Subclusters A and B, located at the east and west of the central region, are particularly compact and roughly spherical shaped, while subcluster C, close to the southern limit of the central region, presents a more elongated and extended distribution of PMS stars. This subcluster corresponds to the visibly prominent main part of the association (see Figure 1 in Paper I).
3.3. LH 95 Central region
For the detailed study of the PMS population in LH 95, while minimizing the contamination by other stars, we focus on the central region of the system. Specifically, we consider the entire population included in the central region and we apply, for a second time, the Monte Carlo technique for field subtraction, taking into account the ratio between the surface of the observed LMC field and that of the central region.
It should be noted that this second application of the field subtraction technique to the central region is not equivalent of simply isolating the stars, which are found to be located within the boundaries of the region after the first iteration of the Monte Carlo subtraction.
Naturally, the central region has a higher density of system stars (specifically PMS stars), and therefore a lower relative contamination by field stars. Consequently, the actual fraction of stars which are marked to be field stars is lower, reducing the probability of (statistically probable) wrong membership assignments. The number of field stars included in the elliptical region of the CMD is computed based on the photometry in the entire field and scaled according to the area occupied by the LH 95 central region, which we found to be of the total ACS fieldofview.
The result is shown in Figure 5, where we plot the CMDs of the whole population included within the boundaries of the central region (left panel), the one of the whole area of the field, sampled randomly and downscaled according to the area coverage of the central region (central panel), and the corresponding CMD of the central region after the field contribution has been removed (right panel). From this figure it is evident that the observed stellar population of the low main sequence (LMS) can be entirely considered as field population. On the other hand the association has a prominent upper main sequence (UMS) component, and all premain sequence stars, easily distinguished by redder colors, are solely present in LH 95 alone. Furthermore, the second application of our field subtraction technique, limited to the central region, allows us to reduce considerably the contamination by objects not located in the PMS and UMS parts of the CMD, which are still present as residuals from the first application of the method on the entire fieldofview (Figure 3).
In the CMDs of Figure 5 can be noted that the field population has a tighter sequence than the PMS population, which may indicate an agespread among the PMS stars of LH 95 (see also Figure 9). Indeed, the locations of lowmass PMS stars in the CMD of young stellar systems of our Galaxy often show a widening, which may be evidence for an agespread in the system (e.g. Palla & Stahler, 2000). However, several characteristics of these PMS stars (being T Tauri stars), such as variability and circumstellar extinction, can cause considerable deviations of their positions in the CMD (Sherry et al., 2004), which may be misinterpreted as an agespread. As a consequence, detailed simulations of the characteristics of PMS stars are required to quantify the effect of these characteristics on their positions in the CMD and to conclude on any true agespread (Hennekemper et al., 2008) with the use of populationsynthesis techniques.
In the following sections we focus our study on the isolated, field subtracted population in the LH 95 central region, shown in Figure 5 (right).
3.4. Interstellar Extinction
With photometry available only in two bands it is not possible in general to measure the visual extinction independently for every star of the system. It is, however, possible to perform a statistical approximation obtaining the average value of reddening as well as its distribution. For this purpose we consider the UMS stars ( mag) in the central part of LH 95, after performing field subtraction, and we compare the position observed in the CMD to that expected according to the Padova grid of evolutionary models in the ACS photometric system (Girardi et al., 2002). Considering that LH 95 is a very young association ( 10 Myr) we assume that all UMS stars have a age equal to the youngest available isochrone in the Padova models, for a metallicity of (typical value for the LMC; Kontizas et al. 1993) and distance modulus mag, compatible with several distance measurements for the LMC, available in the literature (e.g. Alcock et al. (2004)).
For every star we compute the intersection between the reddening vector and the isochrone, and we obtain, thus, the reddening value in terms of . Since the ACS photometric system is slightly different from the standard JohnsonCousins system, we compute the exact reddening parameters for our ACS bands using the typical galactic extinction law of Cardelli et al. (1989), parameterized by a value of . We consider a set of template spectra taken from the NextGen catalog (Hauschildt et al., 1999), and we apply the extinction curve with an arbitrary which we choose to be equal to 1, and we measure the extinction and by means of integration of the original and reddened spectra within the ACS filter profiles. In this way we compute and .
The derived reddening distribution, shown in Figure 6 (middle), appears to have a tail extended to higher values (filled histogram). This tail appears due to a number of UMS objects, which show a high color excesses. Taking into account all objects, including the ones with high color excess, the average reddening has a value mag, corresponding to an optical extinction of . However, these red objects are being suggested to be Herbig Ae/Be (HAeBe) stars (Gouliermis et al., 2002). Such stars of intermediate mass, aged between 1 and 10 Myr, occupy the faint end of the UMS with mag. Being in their premain sequence phase, they have an intrinsic color higher than that of MS stars in the same luminosity range. Consequently, if indeed these are not MS but HAeBe stars, which is possible due to the youthfulness of the system, then the derived reddening would be biased towards higher values.
We estimate the mean reddening without considering these stars (darker histogram in Figure 6  middle) and we derived a lower mean value of mag, corresponding to an optical extinction of mag. In Figure 6 (right) both values are applied to the 4 Myr isochrone model (the youngest available in the grid of models by the Padova Group). From this figure it is evident that the difference between the two values is not dramatic, producing more or less the same shift of the model in the CMD. As a consequence, and in order to avoid any biases due to the presence of young stellar objects in LH 95, we consider an optical reddening of mag as the most representative for the region. We repeated the construction of the reddening distributions shown in Figure 6 (middle) and the determination of the mean reddening values using different isochrones, spanning an age range between 4 and 7 Myr and assuming different distances varying by 0.2 mag around the assumed distance modulus. We found that the mean reddening remains practically unchanged within 0.02 magnitudes, a difference small enough to confirm the accuracy of our estimation.
4. New set of observational PMS models
In the last decades several theoretical models for the premain sequence stellar evolution have been computed. Some of the most popular, for PMS masses above the hydrogen burning limit, are the models of Siess et al. (2000), Palla & Stahler (1999), D’Antona & Mazzitelli (1994, 1997) and Swenson et al. (1994). These models are generally expressed in terms of physical quantities, such as the effective temperature and the total bolometric luminosity, describing the evolution of stars in the HR diagram. However, one of the primary aims of stellar evolution theory is the explanation of the observed photometric data of stars in order to extract their masses and ages from their magnitudes and colors, and therefore a conversion between physical and observable quantities of the models is required. Such a conversion is specifically described by Siess et al. (2000), who include the transformation between and to colors and magnitudes in the Cousins system and infrared bands for their models. These authors use simple relations versus color, as well as bolometric corrections from either Siess et al. (1997) or Kenyon & Hartmann (1995). These conversion tables, derived from observations of stellar clusters, are valid for solar metallicity dwarf stars, but the discrepancies can be significant when one deals with populations with different stellar parameters. In particular, the age of a star is related to the stellar radius and therefore to the surface gravity, which introduces differences in the spectral behavior, and therefore a population of different ages or metallicities may require different conversion relations. This issue can be important for e.g. cold Mtype stars, for which the broad molecular absorption bands dominate the optical spectra, affecting integrated colors and color corrections.
A more thorough method to analyze this issue is to make use of synthetic atmosphere models, performing photometry directly on synthetic spectra. Girardi et al. (2002) applied this method for their evolutionary models for evolved populations, which were converted in several photometric systems using a grid of atmosphere models described by three parameters (, , ). This method was also applied for stars of subsolar mass and brown dwarfs by Baraffe et al. (1998, 2001). For low and intermediatemass stars, however, there are no similar conversions of evolutionary models into the observable plane for various photometric systems available. As a consequence, and in order to have an accurate transformation of the PMS evolutionary models of Siess et al. (2000) into the observational plane, we apply a conversion of these models following an approach similar to that of Girardi et al. (2002). We construct, thus, observational models (both tracks and isochrones) for four assumed metallicities (, , , ) and several photometric systems (specifically Johnson, Cousins, HST WFPC2, HST ACS, HST STIS, HST NICMOS, Galex, Strömgren, SDSS, Tycho and 2MASS).
4.1. The Method
We apply our conversion method first by collecting a grid of atmospheric models in the , , space. These models are shown in Fig. 7. We selected spectra from NextGen (Hauschildt et al., 1999) for low and mediumtemperature stars. The selected grid covers temperatures 2600 K 8200 K in steps of 200 K and surface gravities in cgs units, spaced in 0.5 dex. For higher temperatures we used the grid of models by Kurucz (1979, 1993), with gravities between and . This grid covers entirely the range of stellar parameters in the evolutionary models of Siess et al. (2000). The considered metallicities cover the values , representative of stellar populations in the Magellanic Clouds, as well as in many regions of the Galaxy. We resampled the NextGen models, the resolution of which is much higher than required for our purposes of broadband synthetic photometry, to the wavelength values of the Kurucz models, and we normalized them to the same units (erg cm s Å).
The absolute magnitude in a photometric band of a star with a spectral energy distribution , and stellar radius is given by
(5) 
where is a reference spectrum that gives a known apparent magnitude ; in the Vegamag photometric system, which uses the flux of Lyr as reference, and the zeropoints are close to zero^{2}^{2}2Although the historical Vegamag photometric system was supposed to have zeropoints defined so that to impose the apparent magnitude of Vega equal to zero, improvements in the modeling of the spectrum of Vega, in its measure and in the definition of standard throughputs, have introduced small corrections..
We assume an extinction , since we are interested in constructing a general observable grid of evolutionary models for unreddened objects. Eq. (5) can be written then as:
(6) 
where
(7) 
The latter term can be directly calculated having the synthetic spectrum of any considered star, a calibrated Vega spectrum, the band profile of the photometric filter considered and the eventual zeropoint . For the computation of we used the Gensynphot code included in the Chorizos package (MaízApellániz, 2004). Given a set of stellar parameters (, , and ) for a sample of points (stars) on an evolutionary track or an isochrone, the corresponding spectra are computed by interpolation on the Cartesian grid of synthetic models of stellar atmospheres. Then, the integral at the numerator of Eq. (7) is computed. The final normalization is performed using a recent reference spectrum of Vega (Bohlin, 2007) with the corrected zeropoints as reported by MaízApellániz (2007).
Following this method we converted all the original evolutionary models of Siess et al. (2000) into absolute magnitudes, for the metallicities of and corresponding to respectively. Considering that the original Siess et al. models were computed only for and 0.02, we extrapolated the theoretical isochrones computed for to lower metallicities with the use of the atmosphere models computed for such metallicities. Naturally, a more detailed treatment of lowmetallicity tracks and isochrones requires both evolutionary and atmospheric models to be available at the same considered metallicity. However, for the present analysis in the LMC, for which , the uncertainties are quite low given the small difference between the metallicity of the theoretical isochrones (0.01) and that of the atmosphere models (0.008).
4.2. Comparison with previous conversions
With our construction of the observational plane for PMS evolutionary models, we verified that a single temperaturecolor relation is not sufficient for stars with different physical parameters, especially surface gravity. In Figure 8 (left) we present a comparison between the , colorcolor diagrams for different isochrones, as they were derived from the method of Siess et al. (2000) (red dotted lines) and from our method with the use of atmosphere models (blue continuous lines). In this figure it is shown that the isochrones of Siess et al. (2000) are distributed along one single sequence, due to the use of monolithic temperaturecolor relations, which are independent of the surface gravity. However, our observational isochrones show a spread for high color terms (low ), because differences in ages imply differences in the surface gravities, which produce variations in the integrated colors according to the atmosphere models we used. The major differences arise for Mtype stars, in which the optical spectra are dominated by strong oxide absorption bands. Figure 8 (right panel) shows these differences on the CMD. As shown in this figure, the differences between our method and that of Siess et al. (2000) are mostly evident at the red end of the models, where our computations predict redder colors, and for the youngest models, where the stellar surface gravities are lower than those for main sequence stars. In Figure 8 the example of is shown.
5. The Initial Mass Function
The mass function (MF) is defined as the number distribution of stars as a function of mass. While a general observational approach measures the so called Present Day Mass Function (PDMF), a cornerstone for understanding how stars form is the Initial Mass Function (IMF), which is the mass distribution according to which stars are born. For a premain sequence population younger than the time required by the most massive stars to disperse it, or to evolve into postMS evolutionary phases, the observed mass function coincides with the IMF.
Generally, the IMF is parameterized as follows:
(8) 
namely, approximating it with a powerlaw (e.g. Salpeter, 1955) or with a series of powerlaws, with exponents changing in different mass ranges (Scalo, 1986; Kroupa, 2002). In this section we describe the “average” IMF derived for the stellar populations of the central region of LH 95, as it has been defined in § 3.2.
The , CMD of Figure 5 (right panel) is shown again in Figure 9, with our observable PMS evolutionary tracks (left panel) and three indicative PMS isochrones (right panel), computed in § 4, overlayed. In the right panel of the figure the ZAMS from our grid of models is also plotted along with the youngest MS isochrone from the Padova grid of evolutionary models (Girardi et al., 2002) corresponding to . Both isochrones fit to each other very well demonstrating that the bright MS stars of LH 95 have ages of 4 Myr. As far as the PMS stars of LH 95 concerns, their positions in the CMD corresponds to an agespread that covers ages between 1.5 and 10 Myr according to our PMS models. A simple statistics on the distribution of the ages of these PMS stars, based on their CMD positions, results to a normal distribution of ages which peaks at about 4 Myr. However, if this agespread is true or not can be the subject of debate, as we discuss in § 3.3. Therefore, based on the PMS population of LH 95 we define an indicative age for LH 95 of . Corresponding PMS isochrones for ages of 1.5, 4 and 10 Myr from our grid of observable PMS models are overlayed on the CMD in the right panel of Figure 9. A comparison of this age, which also represents the individual subclusters of LH 95, with the corresponding dynamical crossing time, estimated in § 3.2, shows that all subclusters, as well as the whole central region have an age comparable to their crossing times, indicating that these systems are under ‘mixing’ process if not already ‘mixed’ and therefore we probably do not observe them close to their initial morphology.
5.1. Construction of the IMF
Based on our new set of PMS models computed above for the ACS photometric bands and the average LMC metallicity (), we assigned a value of mass to all the stars identified in the region of LH 95, by comparing our observed data with the evolutionary tracks of these models. For this comparison we assumed a distance modulus of mag, and an average reddening of mag, as it is derived in § 3.4. According to Figure 9, our photometry reveals PMS stars with masses as low as 0.2 M. Our sample, though, as we discuss later, is complete for masses down to 0.43 M. For the PMS stars we assign a mass measurement to each observed star according to their positions in the CMD by interpolating between the PMS evolutionary tracks shown in Figure 9.
The position of the brightest main sequence stars is not covered by our PMS evolutionary models, limited for masses below 7M. Therefore, we derive the masses of these stars directly from the Padova 4 Myr isochrone – already used for the determination of the average reddening in § 3.4 – shifted in the direction of the reddening vector to reach individually each one of these stars. The good agreement between this Padova isochrone and the corrected ZAMS from the Siess models is evident from Figure 9. While we deredden the most massive MS stars according to their individual position in the CMD, their mass estimation is based solely on photometry. A more accurate estimation of the mass for massive stars certainly requires spectroscopy, because for such hot stars the broadband colors are degenerate (Massey et al., 1995), and naturally this effect may influence the mass estimation for the brightest MS stars in our sample and the corresponding IMF. We discuss the influence of the use of photometry alone to the slope of the derived IMF in §5.2.
With the masses of both the MS and the PMS stars in our sample available the construction of their IMF is in principle straightforward, by binning the stars according to their masses and fitting the derived distribution. However, there are considerations, which affect the construction and interpretation of the IMF, and which should be taken into account. We discuss these issues below.

An important issue to be considered is the decontamination process of the stellar sample in the region of the association from the field populations. The Monte Carlo method for the field subtraction, as applied in § 3.3, is based on probabilities and therefore it is a stochastic process. As a consequence, the number of stellar members of the system, which are counted in an arbitrary mass bin is not determined in an unambiguous manner. Moreover, the statistical significance of this number is not determined by the fraction of stars considered as system members over the total number of stars, but it depends also on the probability that all the observed stars in the bin are system members. This dependence affects specifically the MS populations with 19 /mag 22, which as it is shown in Figure 5 are well mixed with field populations, implying a higher uncertainty in the determination of the fraction of system members. In the faint part of the CMD, system and field populations are in general well distinguished, and therefore the uncertainties in the numbers of PMS stars that are true systemmembers are insignificant. However, there is an important possibility that a fraction of the binary sequence, evident in the CMD for the LMC field population (see Figures 5 and 9), contaminates the selected sample of PMS starmembers of the system. We consider this bias and report on its effect in § 5.2

The fitting process of the IMF. Typical linear regression methods, which are used for obtaining a functional form of an observed IMF, such as the method, consider the presence of measurement errors and are based on the assumption that the uncertainty associated to each point follows a Gaussian distribution. However, in this case the uncertainty in the number of counts within a mass bin is the overall effect of both the Poissonian error that naturally comes from the counting process, and the uncertainty that arises from the field subtraction, as discussed above. Although the first can be well approximated with a Gaussian for large , the latter depends on the positions of stars and tracks in the CMDs. Moreover, in the high mass regime, where the number statistics is poor, the Poissonian error cannot be treated as Gaussian, requiring the asymmetry of the distribution to be taken into account.

The completenesscorrection process. In the lowmass regime corrections for the incompleteness of the stellar sample are required to estimate the actual number of stars. However, given the width of the lowmass evolutionary tracks in the CMD shown in Figure 9 and the 2dimensional variation of the completeness itself in the CMD (in both magnitudes and colors; see § 2.1), this correction is not unique for all the stars counted in a considered mass bin.
We construct the system IMF of the association LH 95, with a special care in addressing the above issues, as follows. We first determine a grid of logarithmic mass bins, for which we choose variable width in , that increases from dex at low masses to dex for the higher masses. The counting of the stars is made in variablesized bins because such bins yield very small biases, which are only weakly dependent on the number of stars, in contrast to uniformly binned data (MaízApellániz & Úbeda, 2005). For each bin a number of stars is derived. In order to achieve a statistically correct sampling of points for the constructed IMF, we simulate for each bin the probability distribution expected for every with a sample of 1000 points randomly selected from a Poisson distribution with mean .
We then consider the position in the CMD of all the stars in each bin, and we compute for each star its completeness (the ratio between the number of sources detected and the total number of sources) as a function of both magnitude and color according to the completeness measurements described in § 2.1, for the central region of LH 95. Let us consider the star counted in the bin. The star is located in a specific position of the CMD, where completeness is , meaning that the total number number of stars corresponding to is simply . Consequently, for the total number of stars counted in the bin, , the completeness correction is:
(9) 
which is a quantity .
We multiply our Poissonian sample associated to each bin by , and we divide by the width of the bin in order to express numbers in units of number of stars per solar mass^{3}^{3}3It is worthnoting that the completeness correction should be applied after simulating the Poisson distribution, and not before. For example, if one star of completeness is counted into the bin, the true number of sources associated to the sample is , with an uncertainty of and not .. In order to take into account the additional error introduced by the stochasticity of the Monte Carlo fieldsubtraction (§ 3.3), we iterate the technique of field removal times, each time repeating the procedure explained in the previous paragraphs, namely, counting the number of stars in each bin, sampling a Poisson random population associated to it, and computing and applying the completeness correction.
The outcome of this approach is that for each massbin, instead of a single pointerror combination, we have a population of 200,000 random IMF points, derived according to the true statistical distribution of the expected number of stars in the bin. Clearly, the mean and the standard deviation of each of these distributions represent the IMF data points and error bars. The main advantage of this technique is that we override the nonGaussian behavior of the massbin error bars. Specifically, now for every massbin we do not have one point associated with an error that cannot be treated by the regression theory, but a numerous set of points with no uncertainties associated to them, and consequently with arbitrary equal weights for all massbins. Naturally, this allows a more accurate fitting of the mass function.
In Figure 10 we visualize the concept of this approach. For every massbin the set of random points, representing the distribution of the expected is drawn with tiny colored dots. In this plot an artificial spread along the axis (a random Gaussian shift with of the bin width) is introduced to better highlight the distributions along the axis. The black dots with error bars show the average and standard deviation of each sample of points. The derived numbers are reliable for 0.35 M, whereas for the lower masses (first five mass bins from the left) the completeness corrections are not sufficient. This limitation is evident from the CMD of Figure 9, where it can be seen that for the evolutionary tracks for very low masses ( 0.35 M), the completeness is extremely low, so that no stars can be detected, limiting the ability to recover the true number distributions for such masses.
5.2. IMF model fitting
A commonly used procedure to derive the functional form of an observed IMF is to fit a predefined function to the individual numbers of stars originally counted per massbin taking into account the corresponding counting errors. However, in our case, we have simulated the statistical distribution of per massbin taking into account the additional errors due to incompleteness and field subtraction. Within our statistical approach we perform 200 repetitions of the Monte Carlo field subtraction method and for each of them we map the Poissonian statistics 1000 times. The outcome is a nonanalytical probability distribution for each massbin that represents the true expected statistical uncertainties due to the issues discussed in §5.1. This process provided a large sample of 200,000 points per massbin, and naturally the bestfitting function can be derived more accurately with the use of these data, rather than the original 14 bins. The only assumption that we make is that the functional form of the IMF of LH 95 can be correctly represented by powerlaws, according to the parameterization given in Eq. (8), and therefore we search for the best multiple powerlaw which may reproduce our data. As seen in Figure 10 the IMF seems to be more shallow for the low masses ( M) than for intermediate and high masses, implying that it could be reproduced by a twophase powerlaw. We investigate, however, whether a twophase powerlaw is sufficient or a threephase powerlaw, such as that proposed by Scalo (1998), is required to represent our IMF.
In order to find the best fit of two and threephase powerlaws to our IMF we applied a LevenbergMarquard nonlinear least square minimization technique (Levenberg, 1944; Marquardt, 1963). For this fit we consider only the massbins of the IMF for which the completeness correction allowed us to reconstruct the actual number of stars, meaning for . The slopes of the powerlaws, as well as the position of the break points along the abscissa are the free parameters of our fit algorithm. All 200,000 points of each massbin are considered of equal weight. However, since a twophase powerlaw is a subcase of the threephase powerlaw, the sumofsquares of the latter is always lower. Therefore, the results of the fit algorithm for each case are normalized in a rigorous manner by performing the socalled statistical “F test”.
This test compares the relative increase in the residual sum of squares by reducing the complexity of the model with the relative increase of degrees of freedom, namely, the differences between the number of data points and the free parameters. Factor is, thus, computed as
(10) 
where and are the residual sum of squares for the two and threephase powerlaw best fit models, and are the corresponding number of free parameters, and the total number of points, which, in our case, is the effective number of considered mass bins (14) and not the total number of variables used in the fit (200,000 14), given that the additional multiplicity has been introduced to map the non Gaussian errors. Typically, if the simpler model is valid (which means that there is no need to add parameters), then , whereas if the complex model represents the distribution better, then . We obtain . The probability, , that the null hypothesis, which is “the improvement found adding additional parameters is solely due to chance” (or “the simpler model is good enough”) should be rejected is given by the value of the cumulative Snedecor’s distribution with degrees of freedom. In our case it is , while, generally, a value of is required to state with sufficient significance, that the complex model is better. Under these circumstances, we conclude that our system IMF is best approximated by a twophase powerlaw.
We compute the uncertainty associated to this result, by utilizing a sampling technique. For each massbin we pick randomly one out of the 200,000 points from the corresponding distribution of in that bin, and we perform a twophase powerlaw fit as described above, obtaining the two measurements for the slope, and the point (, ) where the change of the slope occurs. After repeating this process for a substantial amount of times, we obtain a statistically significant sample for the values of the two slopes and the point of change, and consequently the relative confidence level of our fit. Figure 11 shows the result of our fit to the system IMF of LH 95. The best twophase powerlaw fit is shown as a solid line, while the uncertainty of this fit is demonstrated by the grey shaded areas for confidence of 95% (dark grey) and 99% (light grey) respectively. In Table 3 we give the system IMF slopes and their uncertainties with the corresponding mass ranges. The 1 error in the determination of the break point is 0.07 dex.
Another bias that should be considered in the construction of the system IMF of LH 95 is the possible contamination of our selected sample of PMS stars by a binary sequence of the LMC field MS population, discussed in § 5.1. In order to quantify the effect of this bias we select this part of the CMD and remove the corresponding stars. Specifically, we remove all stars bluer than the borders defined by the points (F555WF814W, F814W) (0.9, 22.9), (1.65, 24.3) and (2.2, 26.9) on the CMD. The number of these stars changes according to the stochastic Monte Carlo field subtraction; on average they count to 82 members, less than 4% of the selected PMS population. Then, we reapply the analysis described in §§ 5.1 and 5.2 for the construction and modelfitting of the system IMF. We find that this new IMF does not differ significantly from the one shown in Figure 11, and its shape is not different from that described with the slopes of Table 3. Specifically, this IMF has 99%confidence slopes for , and for , indistinguishable from the slopes of Table 3. Our analysis on the subclusters of LH 95 showed that the contamination of our sample by a MS binary sequence does not alter also the shape of their individual system IMFs, with slopes as well indistinguishable from those given in Table 4.
The lack of spectroscopy and the use of broadband photometry alone for the mass estimation of massive MS stars is found to lead to steeper IMF slopes (Massey et al., 1995), and this may affect the overall intermediate and highmass slope of the IMF of LH 95. However, considering that in their analysis Massey et al. (1995) show that this effect is significant for stars more massive than 15 M, the only massbin in our IMF that can be biased is the most massive one, comprising stars with /M . In addition, the IMF derived from photometry alone for highmass stars in starforming regions of the Magellanic Clouds is found to be also quite shallow and not steep, and within the observed IMF slope variations in such regions (see discussion in §6.2). Under these circumstances, the corresponding uncertainty in our derived system IMF slope due to the use of photometry alone should not be expected to be larger than those given in Table 3.
mass range  

1.05  
2.05 
Note. – A Salpeter (1955) IMF would have a slope of .
6. Discussion
The IMF constructed in the previous section is the most complete extragalactic system IMF down to the subsolar regime derived so far. Specifically, our correction for incompleteness in the derivation of actually provides a reliable IMF down to , as discussed in §5.1. Probably, the most interesting result concerning this IMF is the change in its slope, which we verified statistically for stars in the the subsolar mass range. The IMF slope, , decreases for these stars by one unit. An analogous behavior is generally found in related studies in the Galaxy. Comprehensive studies which summarize the average properties of the galactic IMF are those of Scalo (1986), revisited in Scalo (1998) and Kroupa (2001, 2002). The latter investigations are rather complete, since they cover stellar mass distributions in a wide mass range, from brown dwarfs with to the most massive stars. According to Kroupa (2002), the IMF slope (Eq. 8) changes from in the substellar mass range, to for masses between 0.08 M and 0.5 M, for 0.5 M 1 M, and for stars of higher masses. This IMF is generally characterized as the Galactic average, in the sense that it is reasonably valid for different regions of the Galaxy. A comparison between the system IMF we construct here for LH 95 and that of Scalo (1998), as well as the universal IMF of Kroupa (2001) and the average Galactic IMF of Kroupa (2002), is shown in Figure 12. We scaled all other IMFs to match our mass distribution at . In Figure 12 an overall agreement between the various IMFs is quite evident.
The mass limit of the breaking point in our system IMF ( 1 M) coincides with the “standard Galactic field IMF”, which has been reported to have a change of slope at (Kroupa, 2002), in agreement also with the knee of the IMF derived by Scalo (1998). In addition, there is a second knee in the Kroupa IMF, which occurs at the lower masslimit of 0.5 M. Since we cannot constrain the system IMF below we do not identify this second knee. The slope of our IMF for stars with 1 M falls within the range of slopes found by both Scalo and Kroupa. This is in line with the recent finding that the system IMF is essentially equal to the stellar IMF for intermediate and massive stars (Weidner et al., 2008). Our system IMF, however, is somewhat steeper, as an indication of a bottomheavy IMF, but the use of photometry alone for the estimation of the masses for the brightest stars ( 15 M) may bias the slope to steeper values (see §5.2). As far as the slope of our IMF in the subsolar regime is concerned, it is found slightly more shallow than that of these authors. Specifically, our IMF is by  more shallow than the stellar IMF (as seen from the theoretical plots of Kroupa, 2001), making the resulting stellar IMF for LH 95 indistinguishable to the Galacticfield stellar IMF (e.g. Kroupa, 2002). It should be noted that both Scalo and Kroupa IMFs extend down to masses as low as the hydrogen burning limit ( 0.08 M), well below our detection limit. In general, one can conclude that the system IMF we derive for LH 95 is in agreement with the average galactic IMF in the entire mass range we studied, from 0.43 M to 20 M.
The highmass IMF of LH 95 was constructed earlier from direct CCD imaging in with the 1m telescope at Siding Spring Observatory by Gouliermis et al. (2002). These authors provide the fieldsubtracted completenesscorrected IMF of the system for stars with masses down to M, and they find that this IMF can be fitted with a single powerlaw with slope for stars with masses up to about 17 M. This result indicates a more shallow slope than what we found for intermediate and high mass stars (). However, in our analysis we included masses as low as , extending the mass range investigated in Gouliermis et al. (2002). Indeed, if we exclude the mass bins from 1 to 2 M () from our analysis we derive a more shallow value for for the intermediate and highmass stars. The effect of a steeper IMF for the mass range between 1 and 2 M is evident in Figs. 10 and 11. After excluding stars of M from the considered stellar sample and performing a single powerlaw fit on the resulting , we find an IMF slope of , comparable to that of Gouliermis et al. (2002).
Concerning our derived IMF for the central region of LH 95 it is important to assess the accuracy of our findings, in relation with the possible sources of uncertainties. One of the most important error sources normally comes from the incomplete stellar samples. However, in our case we consider the completeness as it has been derived from artificial star tests in the central part of LH 95, where the crowding is higher (Figure 2b), and we performed rigorous corrections for incompleteness taking into account variations of completeness for all stellar positions in the CMD. As far as the verylow mass stars is concerned, we do not consider those with 0.35 M, where our completeness correction is limited.
6.1. Stellar Binarity
An important source of uncertainty in the construction of the IMF is also the presence of unresolved binaries. Unresolved binarity, can modify the shape of the observed mass function so that it appears more shallow than what it actually is. For this reason, we stress that the IMF reported in Table 3 must be considered as the system IMF of our region, and not a stellar IMF as which one refers when the effect due to binarity is removed. Clearly it is not possible to correct this effect without the knowledge of the actual binary fraction of the studied population, but one can estimate the overall effect on the final mass distribution. Sagar & Richtler (1991) applied Monte Carlo simulations on artificial stellar populations, in order to investigate the change in the mass function slope due to binarity for stars in the mass range . These authors found that the result depends on both the binary fraction and the original slope of the mass function. Specifically, for an original slope and a binary fraction Sagar & Richtler (1991) found a slope change , while for and , . On the other hand, Kroupa et al. (1991) in their investigation on the lowmass IMF in the solar neighborhood, showed that binaries can make a major effect in mass function determinations below 1 M. Consequently, the expected changes in the slope due to unresolved binarity may introduce uncertainties in the derived slopes given in Table 3, but not large enough to justify the intrinsic change of slope that we observe in our IMF for LH 95 at 1 M. In any case, as discussed by Kroupa (2002), most of the studied mass functions in the literature are not corrected for binarity, mostly due to the high uncertainty in the binary fraction. The Galactic IMF derived by (Kroupa, 2001, 2002), however, has been corrected for unresolved multiple systems. A more recent study of the binary effect on the massivestar IMF by Weidner et al. (2008) shows that a powerlaw IMF is not significantly affected even by large numbers of unresolved binaries or higher order multiples.
In order to analyze the effect of unresolved binarity in the measured slope of the system IMF in the subsolar regime, we apply a MonteCarlo technique adopted from the work of Sagar & Richtler (1991) and Kroupa et al. (1991, 1993). For this method we assume a predefined set of functions for the IMF, and in order to be consistent with those we find for the system IMF, as well as with the Galactic IMF of Kroupa (2002), we consider threephase power laws with break points at 0.3 M and 1 M. We limit our mass samples in the stellar regime with 0.08 M. We let the slopes , , within the three selected mass ranges to change respectively in the range for , for and for . For each of these mass distributions, we generate a sample of 50,000 stars positioned in the CMD along our 4 Myr PMS isochrone (corresponding to the age of the system). We assume binary fraction, , between 0 and 1, and we randomly pair couples of stars from the simulated distributions. We then compute the new position in the CMD of these binary systems and assign new masses as if they were single stars. The outcoming mass distributions are fitted with a 3phase power law imposing the same position of the break points.
We isolate the models for which the derived system IMF in the range is compatible with the one measured for LH 95 within the confidence interval of the latter. Naturally there are multiple solutions, given the lack of knowledge about the actual IMF of LH 95, its shape at lower masses, and the binarity fraction appropriate for this system. However, we can constrain the uncertainties due to unresolved binarity in our measured slope, , for the subsolar regime . The IMF slope within this mass range is mainly affected by the value of . If we assume to be that of the Galacticfield IMF, the slope in this mass range changes by for . This change decreases to for , and for .
A change in the original IMF for masses larger than Solar does not affect our results by more than , whereas changing the slope at very low masses () affects significantly the slope , diminishing it when becomes much steeper or shallower than the value of Kroupa (2002) IMF. Even in this case, the effect is always less than . Therefore, we conclude that for a reasonable binary fraction of , unresolved binarity can bias our IMF slope, , in the subsolar regime on average by no more than 0.3 units. As a consequence unresolved binarity does not affect significantly the measured IMF for LH 95 within the estimated errors, and thus the corrected stellar IMF for LH 95 remains indistinguishable from the Galactic IMF.
6.2. Metallicity effects
Recent findings suggest that variations in the metallicity of a starforming cloud might cause variations in the IMF slope for stars with masses M, in the sense that Galactic regions with higher appear to produce more lowmass stars (Kroupa, 2002). This is further supported by the measured IMF in Galactic open clusters (Barrado y Navascués et al., 2001), globular clusters (Piotto & Zoccali, 1999) and old and relatively metalpoor thick disk stars (Reylé & Robin, 2001). Kroupa (2002) proposes a systematic metallicity dependence of the IMF exponent of the form , where . If we apply this dependence to the LMC, for which , we find that the slope of the average IMF in the LMC for masses M should be flattened to the value . The flattening of the IMF of the LMC general field for stars with M has been discussed by Gouliermis et al. (2006a), who found only indications of such a change in the IMF slope, because of observational limitations. Our data, on the other hand, being deeper, show a constant IMF slope in the lowmass regime down to the detection limit of M with no indication of flattening for M. Consequently, we cannot verify any metallicity dependence of the subsolar IMF in the LMC, at least based on the formula of Kroupa (2002). Considering that this result is based on a single stellar system, more studies of lowmass populations in the Magellanic Clouds may provide a better constraint on any systematic metallicitydependence of the lowmass IMF.
As far as the highmass regime is concerned, while the IMF of LH 95 seems to be steeper than that derived from groundbased observations by Gouliermis et al. (2002), our analysis (§6) shows that if we consider only stars more massive than M the IMF of LH 95 has a Salpeterlike slope, comparable to that of Gouliermis et al. (2002). While, as discussed in §5.2, the use of photometry alone leads to steeper IMFs, there are several photometric studies of associations in the LMC, which derive a highmass IMF that is not steeper, but more shallow than the typical Salpeter IMF (e.g. Hill et al., 1994; Oey, 1996; Dolphin & Hunter, 1998). Taking into account all available studies for the intermediate and highmass stars with 7 M, both spectroscopic and photometric (e.g. Massey et al., 1989a, b; Hill et al., 1994; Massey et al., 1995; Hill et al., 1995; Massey & Hunter, 1998) one can see that the corresponding IMF slopes in starforming regions of the LMC are found clustered around the value , similar to the Galactic average IMF. The system IMF we derive here for LH 95 for stars within the same mass limits does not differ significantly from this value. Extensive studies have also shown that the highmass IMF is independent of metallicity in the Galaxy and the Magellanic Clouds (see review by Massey, 2003), and our results also suggest a weak metallicity dependence of the highmass IMF, providing further support to its universality, in the sense that the average IMF of a galaxy seems to be independent from environmental factors, and that any individual stellar system is not expected to share the exact IMF. Indeed, the standard galactic IMF, as presented by, e.g., Scalo (1986) and Kroupa (2002), is the average over a set of systems, each of them showing slightly different mass distributions.
6.3. Variability of the IMF within LH 95
The clustering behavior of the PMS stars in LH 95 is discussed in § 3.2, where three high concentrations of such stars within the main body of the system are identified as the subclusters of the association LH 95 (given the names cluster A, B and C respectively). In order to study how the IMF may vary from one subcluster to the other, and how the IMF of each individual concentration of PMS compares to the “overall” IMF of the system, we constructed the IMF for each of the three subclusters.
We applied the same method used in §5.1 and 5.2, isolating each time the regions of the subclusters, as shown in Figure 4b, iterating the Monte Carlo field subtraction, computing the distribution of number of stars in mass bins and correcting for incompleteness. For the model fitting, we used the same LevenbergMarquard technique, with the difference that, in this case, we constrained the position of the break point between the two powerlaws to be the same as for the general IMF, namely, at . This choice is reasonable considering that we are mainly interested in eventual changes in the derived slopes of the IMF.
The constructed IMF for each subcluster is shown in Figure 13. The overall IMF for the whole system is also overlayed, after it is normalized in order to match the subclusters’ best fit on the break point. In each case the IMF for the subsolar masses is the same as the general IMF, as discussed in §6 for the entire system of LH 95. Moreover, the derived twophase models that fit the IMF of each subcluster and the corresponding uncertainties, given in Table 4, are compatible to that of the overall IMF given in Table 3.
subcluster A  
subcluster B  
subcluster C 
Note. – slopes are in terms of according to Equation 8, uncertainties are confidence intervals.
The IMF constructed for each subcluster in LH 95 suggests that the mass distribution of the PMS stars in this association does not vary much with location. Moreover, the spatial distribution of the intermediate and highmass stars within the whole of the system does not show any indication that the system is mass segregated. On the other hand if the individual subclusters are considered, the more massive stars do occupy the central parts of each of them, providing clear evidence that mass segregation in LH 95 probably occurs in rather small clustering scales.
According to recent results by Weidner & Kroupa (2006) on the integrated galaxial stellar IMF, the overall IMF of a composite system should be steeper than the one it would have if it was not composite, in the intermediate and highmass regime. Taking into account the fact that LH 95 is resolved into at least three individual subclusters (§3.2), it may be considered as a scaleddown case of the scenario of these authors, and therefore as a composite system. It should be noted that also this may explain the somewhat steep slope of the IMF we derive for LH 95.
7. Conclusions
We present a comprehensive investigation of the stellar populations of the starforming association LH 95, located on the northwestern edge of the supergiant shell LMC 4 in the Large Magellanic Cloud, using of deep HST/ACS photometry in the filters and (proxies for standard and respectively). Our findings can be summarized as follows.
We construct a new set of observational PMS models (both isochrones and evolutionary tracks) by starting from the original computations by Siess et al. (2000) and performing synthetic photometry on artificial stellar spectra. The latter were derived from interpolation of atmosphere models from the NextGen (Hauschildt et al., 1999) and Kurucz (Kurucz, 1993) grids. This approach enabled several improvements considering observed data for PMS stars. Specifically:

We computed absolute magnitudes for each point of the original theoretical models in 88 different photometric bands, including the ACS/WFC system, which is useful for our observations. The availability of observational PMS evolutionary models in any given photometric system allows us to avoid the application of any color corrections to the observed data, or to approximate observations in the ACS/WFC system with any standard photometric system.

Our computations, being similar to that of Girardi et al. (2002), do not rely on a single dependence on effective temperature for the computation of colors and magnitudes, but they also take into account the dependencies on the stellar surface gravities (i.e. in luminosity), and therefore they provide a more accurate treatment of the young stages of PMS evolution.

These observational models are extrapolated to lower metallicities than those of Siess et al. (2000) by computing our synthetic photometry for the average LMC and SMC metallicities of and , respectively.
We enhanced the catalog of stars detected in the observed field around the association LH 95 with ACS (Gouliermis et al., 2007) by identifying new point sources, which were previously rejected by the photometry mostly in crowded regions. Therefore, about additional stars are included in our photometric catalog (Table 1), which now includes in total 17,245 stars. Our deep ACS photometry allowed us to analyze the properties of the abundant PMS population found in LH 95, which shows an average age of 4 Myr. We performed star counts on our stellar sample of the area of the association and we selected the central part of the observed field, where a statistically significant concentration of PMS stars is found, as the most representative region of LH 95. The 3 isopleth was selected to define the boundaries of the system. The PMS stars within the association are found to be mostly concentrated in three distinctive subclusters, named cluster A, B and C respective, providing evidence of the existence of stellar subgroups within a single association.
We estimate the total stellar masses and corresponding dynamical timescales for the system itself as well as the three subclusters, and we find that the system cannot be considered dynamically evolved, and it does not experience any disruption. Moreover, the fact that LH 95 is a bound cluster with its gas partially removed implies that the star formation efficiency is probably higher than the typical values for young clusters.
We estimated the stellar contribution of the general background field of the galaxy in the area of the system from our ACS observations of a nearby control field. We then performed field subtraction from the area of the system with the use of Monte Carlo simulations.
With our evolutionary models we estimated the interstellar extinction in LH 95 and we computed the masses of all the members of the system. Our stellar sample reaches masses as low as 0.2 M, and we are able to correct the incompleteness accurately down to 0.43 M. For the first time, this allows the construction of the extragalactic IMF down to the subsolar mass regime.
We derive the system IMF by counting stars within bins of variable width, which increase with mass, and correcting for incompleteness based on artificial star experiments for the stars included in each mass range. We derive the actual statistical distribution of the counting errors associated to each mass bin via a Monte Carlo technique, in order to take into account uncertainties introduced by the contamination from the field population. The statistical F test for model selection allowed us to determine a twophase powerlaw as the most appropriate model to fit our IMF. We detect a change in the slope of the IMF (the “knee” of the IMF) for masses lower than 1 M. A LevenbergMarquard algorithm for least square fitting was applied to determine the powerlaw exponents and the point where the slope changes. The derived system IMF slope is found to be in the subsolar regime, and at intermediate and high masses, with a of confidence.
This system IMF, once corrected for unresolved binarity, is statistically compatible, within the studied mass range, with the average Galactic IMF, concerning its slope and the mass limit where it changes (1 M). As a consequence, we cannot consider the metallicity difference between the LMC () and the Galaxy important enough for a significant change in the IMF slope. Certainly similar data on such stellar systems will show if this is a normal behavior of the IMF in this galaxy. We find no significant differences in the shape of the overall IMF of LH 95 from that of each of the three individual PMS subclusters of the association. This clearly suggests that the IMF of LH 95 is not subject to local variability.
References
 Alcock et al. (2004) Alcock, C., et al. 2004, AJ, 127, 334
 Allen et al. (2007) Allen, L., et al. 2007, Protostars and Planets V, 361
 Ambartsumian (1947) Ambartsumian, V. A. 1947, In Stellar Evolution and Astrophysics, Armenian Acad. of Sci.
 Baraffe et al. (1998) Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403
 Baraffe et al. (2001) Baraffe, I., Chabrier, G., Allard, F., Hauschildt, P. 2001, From Darkness to Light: Origin and Evolution of Young Stellar Clusters, 243, 571
 Barrado y Navascués et al. (2001) Barrado y Navascués, D., Stauffer, J. R., Bouvier, J., & Martín, E. L. 2001, ApJ, 546, 1006
 Bohlin (2007) Bohlin, R. C. 2007, The Future of Photometric, Spectrophotometric and Polarimetric Standardization, 364, 315
 Bouchet et al. (1985) Bouchet, P., Lequeux, J., Maurice, E., Prevot, L., & PrevotBurnichon, M. L. 1985, A&A, 149, 330
 Briceño et al. (2007) Briceño, C., Hartmann, L., Hernández, J., Calvet, N., Vivas, A. K., Furesz, G., & Szentgyorgyi, A. 2007, ApJ, 661, 1119
 Cardelli et al. (1989) Cardelli, J. A., Clayton, G. C., Mathis, J. S. 1989, ApJ, 345, 245
 Caldwell & Coulson (1986) Caldwell, J. A. R., & Coulson, I. M. 1986, MNRAS, 218, 223
 Cole (1998) Cole, A. A. 1998, ApJ, 500, L137
 Chabrier (2003) Chabrier, G. 2003, PASP, 115, 763
 D’Antona & Mazzitelli (1994) D’Antona, F., Mazzitelli, I. 1994, ApJS, 90, 467
 D’Antona & Mazzitelli (1997) D’Antona, F., Mazzitelli, I. 1997, Memorie della Societa Astronomica Italiana, 68, 807
 Dolphin & Hunter (1998) Dolphin, A. E., & Hunter, D. A. 1998, AJ, 116, 1275
 Dolphin (2000) Dolphin, A. E. 2000, PASP, 112, 1383
 Fellhauer et al. (2006) Fellhauer, M., Kroupa, P., & Evans, N. W. 2006, MNRAS, 372, 338
 Girardi et al. (2002) Girardi, L., Bertelli, G., Bressan, A., Chiosi, C., Groenewegen, M. A. T., Marigo, P., Salasnich, B., Weiss, A. 2002, A&A, 391, 195
 Gouliermis et al. (2002) Gouliermis, D., Keller, S. C., de Boer, K. S., Kontizas, M., & Kontizas, E. 2002, A&A, 381, 862
 Gouliermis et al. (2006a) Gouliermis, D., Brandner, W., & Henning, Th. 2006a, ApJ, 636, L133
 Gouliermis et al. (2006b) Gouliermis, D. A., Dolphin, A. E., Brandner, W., & Henning, Th. 2006b, ApJS, 166, 549
 Gouliermis et al. (2007) Gouliermis, D. A., Henning, Th., Brandner, W., Dolphin, A. E., Rosa, M., Brandl, B. 2007, ApJ, 665, L27 (Paper I)
 Hauschildt et al. (1999) Hauschildt, P. H., Allard, F., Ferguson, J., Baron, E., & Alexander, D. R. 1999, ApJ, 525, 871
 Hennekemper et al. (2008) Hennekemper, E., Gouliermis, D. A., Henning, Th., Brandner, W., & Dolphin, A. E. 2008, ApJ, 672, 914
 Hill et al. (1994) Hill, R. J., Madore, B. F., & Freedman, W. L. 1994, ApJ, 429, 204
 Hill et al. (1995) Hill, R. S., Cheng, K.P., Bohlin, R. C., O’Connell, R. W., Roberts, M. S., Smith, A. M., & Stecher, T. P. 1995, ApJ, 446, 622
 Kenyon & Hartmann (1995) Kenyon, S. J., Hartmann, L. 1995, ApJS, 101, 117
 Kontizas et al. (1993) Kontizas, M., Kontizas, E., & Michalitsianos, A. G. 1993, A&A, 269, 107
 Koornneef (1982) Koornneef, J. 1982, A&A, 107, 247
 Kroupa et al. (1991) Kroupa, P., Gilmore, G., & Tout, C. A. 1991, MNRAS, 251, 293
 Kroupa et al. (1993) Kroupa, P., Tout, C. A., & Gilmore, G. 1993, MNRAS, 262, 545
 Kroupa (2001) Kroupa, P. 2001, MNRAS, 322, 231
 Kroupa (2002) Kroupa, P. 2002, Science, 295, 82
 Kroupa (2008) Kroupa, P. 2008, in The Cambridge Nbody Lectures, S. Aarseth, C. Tout, R. Mardling (eds), Lecture Notes in Physics Series, Springer Verlag (arXiv:0803.1833)
 Kurucz (1979) Kurucz, R. 1979, ApJS, 40, 1
 Kurucz (1993) Kurucz, R. 1993, ATLAS9 Stellar Atmosphere Programs and 2 km/s grid. Kurucz CDROM No. 13. Cambridge, Mass.: Smithsonian Astrophysical Observatory, 1993., 13
 Lada & Lada (2003) Lada, C. J., & Lada, E. A. 2003, ARA&A, 41, 57
 Levenberg (1944) Levenberg, K. 1944, The Quarterly of Applied Mathematics 2, 164
 Luck et al. (1998) Luck, R. E., Moffett, T. J., Barnes, T. G., III, & Gieren, W. P. 1998, AJ, 115, 605
 MaízApellániz (2004) MaízApellániz, J. 2004, PASP, 116, 859
 MaízApellániz & Úbeda (2005) Maíz Apellániz, J., & Úbeda, L. 2005, ApJ, 629, 873
 MaízApellániz (2007) MaízApellániz, J. 2007, The Future of Photometric, Spectrophotometric and Polarimetric Standardization, 364, 227
 Marks et al. (2008) Marks, M., Kroupa, P., & Baumgardt, H. 2008, MNRAS, 386, 2047
 Marquardt (1963) Marquardt, D. 1963, SIAM Journal on Applied Mathematics 11, 431
 Massey et al. (1989a) Massey, P., Silkey, M., Garmany, C. D., & DegioiaEastwood, K. 1989a, AJ, 97, 107
 Massey et al. (1989b) Massey, P., Parker, J. W., & Garmany, C. D. 1989b, AJ, 98, 1305
 Massey et al. (1995) Massey, P., Lang, C. C., DegioiaEastwood, K., & Garmany, C. D. 1995, ApJ, 438, 188
 Massey & Hunter (1998) Massey, P., & Hunter, D. A. 1998, ApJ, 493, 180
 Massey (2003) Massey, P. 2003, ARA&A, 41, 15
 Nota et al. (2006) Nota, A., et al. 2006, ApJ, 640, L29
 Oey (1996) Oey, M. S. 1996, ApJS, 104, 71
 Palla & Stahler (1999) Palla, F., Stahler, S. W. 1999, ApJ, 525, 772
 Palla & Stahler (2000) Palla, F., & Stahler, S. W. 2000, ApJ, 540, 255
 Pei et al. (1999) Pei, Y. C., Fall, S. M., & Hauser, M. G. 1999, ApJ, 522, 604
 Piotto & Zoccali (1999) Piotto, G., & Zoccali, M. 1999, A&A, 345, 485
 Preibisch et al. (2002) Preibisch, T., Brown, A. G. A., Bridges, T., Guenther, E., & Zinnecker, H. 2002, AJ, 124, 404
 Reylé & Robin (2001) Reylé, C., & Robin, A. C. 2001, A&A, 373, 886
 Sabbi et al. (2008) Sabbi, E., et al. 2008, AJ, 135, 173
 Sagar & Richtler (1991) Sagar, R., & Richtler, T. 1991, A&A, 250, 324
 Salpeter (1955) Salpeter, E. E. 1955, ApJ, 121, 161
 Scalo (1986) Scalo, J. M. 1986, Fundamentals of Cosmic Physics, 11, 1
 Scalo (1998) Scalo, J. 1998, The Stellar Initial Mass Function (38th Herstmonceux Conference), 142, 201
 Schmalzl et al. (2008) Schmalzl, M., Gouliermis, D. A., Dolphin, A. E., & Henning, Th. 2008, ApJ, 681, 290
 Sherry et al. (2004) Sherry, W. H., Walter, F. M., & Wolk, S. J. 2004, AJ, 128, 2316
 Siess et al. (1997) Siess, L., Forestini, M., Dougados, C. 1997, A&A, 324, 556
 Siess et al. (2000) Siess, L., Dufour, E., Forestini, M. 2000, A&A, 358, 593
 Spitzer (1958) Spitzer, L. Jr. 1958, ApJ, 127, 17
 Swenson et al. (1994) Swenson, F. J., Faulkner, J., Rogers, F. J., & Iglesias, C. A. 1994, ApJ, 425, 286
 Weidner & Kroupa (2006) Weidner, C., & Kroupa, P. 2006, MNRAS, 365, 1333
 Weidner et al. (2008) Weidner, C., Kroupa, P., & Maschberger, T. 2008, accepted for publication in MNRAS (arXiv:0811.3730)
 Westerlund (1997) Westerlund, B. E. 1997, The Magellanic Clouds (Cambridge: Cambridge Univ. Press)
 Wilking & Lada (1983) Wilking, B. A., & Lada, C. J. 1983, ApJ, 274, 698