Free Access
Issue
A&A
Volume 534, October 2011
Article Number A10
Number of page(s) 8
Section Galactic structure, stellar clusters and populations
DOI https://doi.org/10.1051/0004-6361/201117062
Published online 21 September 2011

© ESO, 2011

1. Introduction

A central question in astrophysics is whether the initial mass function (IMF) is universal or if it is a function of star formation environment. The answer is crucial for our understanding of the buildup of galaxies and provides clues to how stars and brown dwarfs are formed. The stellar part of the IMF appears to be universal based on the comparison of the field IMF (Chabrier 2003) to clusters within 2 kpc (e.g. Meyer et al. 2000). Bastian et al. (2010) reached the same conclusion with the possible exception of the Taurus Dark Cloud (Luhman et al. 2009). For more massive Galactic clusters and for low-metallicity clusters in the Large Magellanic Clouds the IMF is well represented by a power-law to the sensitivity and crowding limit of  ~1 M (e.g. Stolte et al. 2006; Gouliermis et al. 2007; Brandner et al. 2008; Andersen et al. 2009).

The situation is less clear for the brown dwarf part of the IMF. Briceño et al. (2002) and Luhman (2004) suggested a higher ratio of stars to brown dwarfs in Taurus than the Orion nebula cluster (ONC). However, Guieu et al. (2006) utilized a larger area survey of Taurus and found the ratio of stars to brown dwarfs to be similar for the two regions (cf. Luhman et al. 2009). Andersen et al. (2008) combined the measured IMF’s for six nearby star-forming regions and the Pleiades and showed that there is no strong evidence for variations in the derived ratio of low-mass stars to brown dwarfs for the regions and that the combined ratio for the 7 regions showed that the underlying IMF was falling into the brown dwarf regime.

Conversely, a handful of studies for OB associations over large areas of the sky have indicated an IMF that does not show evidence for a falling mass function in the brown dwarf regime (e.g. Lodieu et al. 2007). The reason for the apparent discrepancy is not clear but could be due to either differences in the star formation process or due to dynamical evolution of the clusters. For several of the clusters compiled by Andersen et al. (2008) only the central regions were studied (e.g. a 0.6 pc square box for the ONC and 0.24 pc square box for Mon R2) and the measured IMF in the centre can thus be biased if mass segregation is present. Large field of view surveys of young star forming regions are needed to directly compare with the OB associations.

One nearby region of particular interest is the ONC. So far, the cluster content has only systematically been studied deep into the brown dwarf regime for the central  ~5′ square (Hillenbrand & Carpenter 2000; Muench et al. 2002; Slesnick et al. 2004), corresponding to  ~0.6 pc square for a distance of 414 pc (Menten et al. 2007). Larger area, less deep, studies of the stellar population have indicated at least the higher mass objects in the ONC is mass segregated (Hillenbrand & Hartmann 1998). Thus, the low-mass IMF derived for the central area can be biased towards higher mass objects and may not be representative of the cluster as a whole. However, no studies have combined deep observations in the near-infrared together with a large field of view coverage such that the brown dwarf content can be studied in detail as a function of radius and to determine the global ratio of stars to brown dwarfs. Here we present HST/NICMOS Camera 3 F110W and F160W (similar to the ground-based J and H bands) non-contiguous observations of the ONC area. The observations cover the cluster from the very centre out to a radius of 2.55 pc and were obtained as part of the HST Orion Treasury Project (GO 10246, PI Robberto).

The paper is structured as follows. Section 2 presents the observations and data reduction. The source detection and the photometry is discussed in Sect. 3. In Sect. 4 we derive the ratio of low mass stars to brown dwarfs as a function of radius in the cluster and we discuss the evidence for the lower mass cluster content to be mass segregated. In Sect. 5 we compare the findings for the ONC with the similar results from other nearby star forming regions. Finally, we conclude in Sect. 6.

2. Observations

Regions of the ONC were observed using the HST/NICMOS Camera 3 during the HST Treasury Program covering the ONC (GO10246). The non-contiguous nature of the NICMOS Camera 3 observations stems from the fact they were observed in parallel with the ACS and WFPC2 Treasury Project observations. The field-of-view of the NICMOS Camera 3 is 51\hbox{$\farcs$}2 × 51\hbox{$\farcs$}2 and the pixel scale is 0\hbox{$\farcs$}2. Since historically, the NICMOS Camera 3 is slightly out of focus, there is a slight degradation of the point spread function (PSF) which in turn improves the sampling of the PSF.

A total of 102 regions were covered in both the F110W and F160W bands. The regions were typically build up as a 4 piece mosaic in the F160W band and a 5 piece mosaic in the F110W band utilizing the field offset mechanism (FOM, see below). In one point, a single 51\hbox{$\farcs$}2 × 51\hbox{$\farcs$}2 image was observed in both filters. The integration time of each image was 192 s in F160W and 256 s in F110W, resulting in a total integration time per pixel that varies between 192–762 s for the F160W observations and 256–1280 s for the F110W observations, with the central parts of each sub-mosaic being imaged for the full amount of time. Around 55% of the mosaics are covered for 192 s in the F160W band, and 40% for 380 s. For the F110W filter, around 45% is covered for 256 s, and 40% for 512 s.

The large dithering pattern was possible due to the FOM, which however introduces some vignetting. Fortuitously, a 3 × 3 NICMOS Camera 3 mosaic of the Trapezium was obtained in 1998 (Luhman et al. 2000). We were able to make ratios of one vignetted FOM-dithered exposure to that unvignetted exposure. A smooth function was fit to the vignetted edge in the ratio image and then expanded into a 2D vignetting correction image (one per filter for F110W and F160W) that was applied to all affected exposures. The raw images have further been processed with the appropriate dark and flat frames as described on the STScIs web pages. They were further corrected for bad pixels determined from sigma clipping of the median stacked image.

Offsets were determined between the individual images of each region by identifying stars in common and determining the centroid. The frames were then combined using the drizzle software that takes the distortion of the NICMOS Camera 3 frames into account and combines the frames with sub-pixel accuracy. All the source detection, astrometry, and photometry is being done on the final mosaic for each region.

3. Source detection and photometry

We describe the source detection and photometry process. The photometry is compared to the previous work by Robberto et al. (2010) and is converted into the 2MASS system. Completeness corrections are determined and the fractional coverage of the ONC is described.

3.1. Source detection and photometry

Source detection is complicated by strong nebulosity, especially close to the centre of the ONC, and by the under-sampled PSF of the NICMOS Camera 3. A ring median filtered image was subtracted from each mosaiced image in order to remove the large scale nebulosity. The inner and outer annuli for the ring filter were 7 and 11 pixels, (1\hbox{$\farcs$}4–2\hbox{$\farcs $}2). This makes source detection easier but the photometry was performed on the frames with the nebulosity maintained. The PSF is better sampled in the F160W images than the F110W images due to the larger full width at half maximum at longer wavelengths. Since the F160W data are also deeper than the F110W data, we have used the F160W data for the source detection.

Table 1

Sources detected in the NICMOS camera 3 mosaics.

Point sources were visually identified in each individual frame. The radial profile of each source was inspected in order to exclude spurious detections due to nebulosity. We detect a total of 2054 objects. The positions of the objects detected in the F160W images were then used for the photometry in the F110W band. The photometry was performed with an aperture radius of 2.5 pixels (0\hbox{$\farcs$}5) and a sky annulus between 10 and 15 pixels (2″–3″) for both filters. Zero points were adopted from the HST calibration web pages. Bright isolated stars in different regions were used to determine aperture corrections. The random error on the aperture correction was 0.02 mag for both bands. The adopted zero points, including the aperture corrections, are 22.50 mag for the F110W band observations and 21.66 mag for the F160W band observations.

Objects in common with the dataset by Robberto et al. (2010) were used to calibrate the astrometry of each mosaic. For one of the mosaics, no ISPI star was present within the field-of-view. Coordinates for the stars in this field were obtained from the header information and is expected to be accurate to  ~1″, the average accuracy of the header coordinates for the mosaics that could be calibrated with ISPI. The astrometry is in agreement with the ISPI astrometry to within 0\hbox{$\farcs$}1 for 60% of the sources and within 0\hbox{$\farcs$}2 for almost 90% of the sources. The sources and their photometry are presented in Table 1.

thumbnail Fig. 1

Upper and lower left: a comparison of the NICMOS photometry converted into the 2MASS system with the photometry obtained with ISPI, in the 2MASS system (Robberto et al. 2010). Only objects with photometric errors less than 0.1 mag in both filters are shown. Upper and lower right: the NICMOS photometric errors as a function of magnitude.

3.2. Integrity of the photometry

The NICMOS photometry has been calibrated into the 2MASS system utilizing the ISPI data which was transformed into the 2MASS system (Robberto et al. 2010). We chose relatively well exposed sources, not located close to the edge of the NICMOS mosaics. Least square linear fits were performed to the F160W band magnitudes and the F110W − F160W color. H2MASS=F160W+a0HJ2MASS=F110W+a0J+a1J(F110WF160W),\begin{eqnarray} {H_{\rm 2MASS}} & = & \filter{F160W}+a_{0H}\\ {J_{\rm 2MASS}} & = & \filter{F110W}+a_{0J}+a_{1J}*(\filter{F110W}-\filter{F160W}), \end{eqnarray}where a0H = −0.25 ± 0.01, a0J = −0.26 ± 0.01, and a1J = −0.35 ± 0.01, respectively. Figure 1 shows the ISPI photometry as a function of the NICMOS photometry for both the J and H bands. All sources with a nominal photometric error smaller than 0.1 mag in each band both in this study and the ISPI observations are shown. Further, the photometric error for the NICMOS sources is shown for both filters. The agreement with the ISPI photometry is in general good but there is some scatter. This is to be expected due to the difference in spatial resolution making it easier to separate starlight and nebulosity in our observations but also due to our poorly sampled PSF, especially in the J band. In any case, our photometry is significantly deeper than the previous ISPI study, which reached 3σ limits of J ~ 19.5 mag, and H ~ 18 mag. The faintest sources identified in this study have magnitudes of J ~ 22.0 mag, and H ~ 21.0 mag for the areas with limited nebulosity. The discussion below describes in more detail how the completeness limit depends on the distance from the cluster centre, and hence on the amount of nebulosity.

3.3. Completeness corrections

We have performed completeness experiments to determine the fraction of stars detected as a function of object magnitude. Bright isolated stars from several frames were used to create a model for the PSF. Artificial point sources were then placed randomly within the observed frames, with F160W magnitudes randomly chosen between 15 and 20 mag. In each frame around 80% of the detected number of stars were added as artificial stars. The frames are not crowded and the artificial star experiments are being performed to test the detection limit as a function of nebulosity. The relatively large number of stars added per frame does not influence this limit. A color of F110W − F160W = 2 mag, was adopted for all the artificial stars, corresponding to the colour of the faintest objects within the extinction limited sample defined below (see Fig. 5).

The artificial stars were recovered by hand in the F160W frames in a similar manner to the real star detection. Photometry was then performed in both the F160W and the F110W frames. An artificial star was considered detected if the magnitude error in both filters was less than 0.3 mag, similar to the maximum error for some objects used in the extinction limited sample.

thumbnail Fig. 2

The completeness correction for different annuli as a function of F160W band magnitude. From the upper left is shown the completeness corrections for the inner annulus (0.35–0.7 pc), the middle annulus (0.7–1.1 pc), the outer annulus (1.1–1.5 pc), and in the lower-right a comparison of the three annuli.

Figure 2 shows the completeness corrections for three different annuli as a function of F160W band magnitude. The 50% completeness limit clearly varies across the survey area. The presence of strong nebulosity close to the cluster centre results in a brighter 50% completeness limit, as expected. In particular, within the central 0.35 pc radius, the strong nebulosity and undersampled PSF makes the detection of fainter point sources difficult and the completeness corrections are substantial. On the other hand, farther than 0.35 pc from the centre the completeness is high down to F160W = 18 mag. For comparison, a 1 Myr 20 Mjup object at the distance of the ONC has H band magnitude of 16.4 mag. Thus, even at 0.35 pc we can detect essentially all 1 Myr objects down to 20 Mjup seen through an extinction of AV = 7 mag (see below), adopting the models by Baraffe et al. (1998). At larger radial distances objects below 10 Mjup are detected. Due to the incompleteness within 0.35 pc we focus our analysis on the objects at larger radii. The current survey thus complements the observations by Hillenbrand & Carpenter (2000) and Luhman et al. (2000) that both cover the central 0.6 pc square region. For all radii the completeness limit is higher than the detection limit of the data, which is due to both the nebulosity and to the structure in the nebulosity that can hide in particular faint sources.

3.4. Fraction of the ONC covered

Since the observations were obtained in parallel, the field-of-view covered is not contiguous. The very centre of the cluster is not covered and the majority of the observations are to the south and south-south-west of the cluster centre. Figure 3 shows the location in the ONC of each individual mosaic. For reference, the mosaic observed by Luhman et al. (2000) is inserted and the Trapezium stars marked. Figure 4 shows the fractional area of the ONC covered as a function of distance from the cluster centre.

thumbnail Fig. 3

F160W band 26′ × 33′ mosaic of the observed region within the ONC. For reference, the central 140″ × 140″ region observed by Luhman et al. (2000) is included. The location of the Trapezium stars are marked by small boxes in the mosaic. North is up and East is to the left.

thumbnail Fig. 4

The fraction coverage of the ONC by the NICMOS Camera 3 observations. The fraction covered is shown as a function of radial distance from the centre (θ1 Ori C).

The total area imaged is 162 square arcminutes, corresponding to 2.35 pc2 for a cluster distance of 414 pc. The coverage extends in radius from 0.05 pc out to 2.55 pc from θ1 Ori C. Due to the large region of sky covered with a small filling factor, the coverage is relatively low, 20–50% within 0.5 pc, 10–30% between 0.5 pc and 1.5 pc, and less at larger distances as shown in Fig. 4.

4. Results and analysis

We now present the J − H versus H color–magnitude diagrams for the surveyed region and estimate the contamination due to field stars. The mass of each object is derived from its de-reddened H band magnitude. We create an extinction limited sample to ensure an unbiased selection of cluster members as a function of object mass before determining the ratio of low-mass stars to brown dwarfs for the ONC as a function of radius.

thumbnail Fig. 5

The color–magnitude diagram for the surveyed region. Overplotted are a 1 Myr isochrone (Baraffe et al. 1998) shifted to the distance of the ONC, the same isochrone reddened by AV = 7 mag, and the 50% completeness limit for the data. Masses for a 1 Myr population of 1.0, 0.08, 0.02, and 0.01 M are indicated on the lower left panel. The region used as control field is shown in the lower right corner (2.3–2.6 pc). The typical photometric and colour error are shown in each diagram as a function of object magnitude.

4.1. The color–magnitude diagram

Figure 5 shows the color–magnitude diagram for the surveyed region outside 0.35 pc. The data have been divided into the annuli that will be used below to derive the ratio of low-mass stars to brown dwarfs. Overplotted are the 50% completeness limits, a 1 Myr isochrone from Baraffe et al. (1998) converted into the 2MASS system using the same procedure as in Andersen et al. (2006) with the filter transformations by Carpenter (2001), and the same isochrone reddened by AV = 7 mag.

Most of the objects from 1 M to the detection limit below 0.01 M are located between the 1 Myr isochrone shifted by the distance modulus to the ONC, and the same isochrone reddened by AV = 7 mag. The objects on the blue side of the 1 Myr isochrone are likely foreground stars or a significantly older population.

4.2. Extinction limited samples

One complication in calculating the IMF in young star forming regions is the presence of differential extinction. Objects are reddened by different amounts due to patchy obscuration along the line of sight. As evident from Fig. 5, this is the case for the ONC. Higher mass objects can be detected through larger amounts of extinction than lower mass objects due to their intrinsic brightness. As illustrated in e.g. Andersen et al. (2009) this can potentially introduce a spurious deficit in the derived number of low mass objects. One approach to avoid this bias is to construct an extinction-limited sample, which is obtained by considering only objects within a certain range of extinction while taking the completeness limit into account.

The lower mass limit for the extinction-limited sample is set by a combination of depth of the data, completeness corrections, and maximum extinction allowed. The choice of an extinction limit is a compromise between completeness and the low mass limit (lower extinction) and maximizing the sample size (higher extinction). The extinction for each individual object is determined by shifting the object in the color-magnitude diagram back to the assumed isochrone, in a similar manner as in Andersen et al. (2006). We find that more than 80% of the objects are affected by an extinction of AV = 7 mag or less and use this limit for our extinction limited sample.

We have examined whether the extinction distribution is similar at different distances from the centre of the ONC. Figure 6 shows the cumulative distribution of derived extinction for objects in the magnitude range 10.5 < H < 18 mag, similar to the magnitude range from where the IMF is derived. There is no obvious trend of the extinction distribution as a function of radius. A two-sided Kolmogorov-Smirnov test between the cumulative distributions for the different radii indicates that we cannot rule out that the underlying distribution is the same. The maximum deviations between the extinction distributions for different radii were 0.12 with the significance levels varying between 0.37 and 0.62. Thus, since the distribution of extinction appears to be similar as a function of radius, the same extinction limit is used for all radii.

thumbnail Fig. 6

The derived extinction for sources in the ONC within the magnitude limits 10.5 < H < 18 mag, which corresponds to the range from where the ratio of stars to brown dwarfs is determined.

The mass of each object is then determined by sliding it along the reddening vector in the colour-magnitude diagram back to the isochrone. Since only J and H band observations are available we have no knowledge if an object is associated with a circumstellar disc that would give rise to an infrared excess. With JHK photometry the disk is typically taken into account by de-reddening the objects to the T Tauri locus in the colour–colour diagram (Meyer et al. 1997).

Not taking disks explicitly into account will tend to overestimate the derived mass compared to the actual mass. Thus, a strong radial dependence on the disk fraction could in principle mimic a decreasing ratio of stars to brown dwarfs as a function of radius. Hillenbrand et al. (1998) find a modest decrease of the disk fraction, from some 65% showing strong evidence of a disk in the centre to 45% at a radius of 1.1 pc. The ratio determined from the central 0.3 pc suffers less from the effects of disks. Slesnick et al. (2004) determined the spectral types through spectroscopy and although veiling is still a potential issue, the use of J band spectra reduces this effect since the photosphere for late type objects peak in this wavelength range. However, the bolometric magnitudes still depend on the de-reddening in the colour–magnitude diagram. Using different de-reddening approaches provide difference in the luminosity of Δlog L = 0.05 and is not expected to have an influence on the ratio of stars to brown dwarfs. The colour–magnitude diagrams in Fig. 5 further indicate this is not likely to be the source of the decrease in the ratio of stars to brown dwarfs as a function of radius. We would expect to see the strongest effect is in the 0.35–0.7 pc bin where few objects are located close to the star-brown dwarf division line. For objects at the star-brown dwarf boundary the excess is J − H of 0.1 mag or less (Meyer et al. 1997), which after de-reddening will correspond to maximum shift in H band magnitude of 0.16 mag. At most 4 objects would be redefined from stars to brown dwarfs from 0.35–0.7 pc assuming all objects had disks and all objects were affected by 0.16 mag of excess. Several objects at larger distances could also potentially be re-classified, 1 object from 0.7–1.1 pc, and 3 objects from 1.1–1.5 pc. Neither will change the ratios significantly and according to Hillenbrand et al. (1998) half or less of the objects would show evidence for an infrared excess. We thus conclude it is unlikely the decrease in the ratio of stars to brown dwarfs is due to a lower disk fraction at larger radii.

Table 2

Ratio of stars to brown dwarfs in an extinction limited sample with AV = 7 mag as a function of lower mass limit.

4.3. Field star contamination

Due to the relatively large area covered we expect some field star contamination. The contamination can be separated into two types, namely that of foreground and background stars. Since the ONC is located at only 414 pc, we expect the foreground field star contamination to be relatively minor. In addition, the foreground objects will have relatively blue colors compared to the cluster members (see below).

Background stars are a larger concern. Two effects limit the contamination. First, the ONC is located at (l,b) = (209.07, − 19.49), out of the Galactic plane and away from the inner Galaxy. Further, the ONC and surrounding region is associated with and partly in front of a giant molecular cloud complex acting as a shield for background objects. By comparison with the 13CO line map from Bally et al. (1987) it is evident that all the observed region is located within the shaped filament in Orion A-North, in agreement with strong emission seen in the IRAS 100 μm images. Based on the strength of the CO emission the extinction averaged over a 40″ beam is well above AV = 15 mag, in agreement with the low number of reddened stars seen even far from the centre of the ONC.

We have estimated the field contamination in the following way. Since the cluster surface density at larger radii decreases as r-2 (Hillenbrand & Hartmann 1998), we expect the outer region (2.0–2.6 pc) to have a relatively small number of cluster members. Thus, the outer regions can be used as a conservative over-estimate of the field star contamination. However, as we indicate below there appears to be an increase in the relative number of brown dwarfs at larger radii. Any cluster population present at large radii will therefore result in an over subtraction of the brown dwarf content. The derived ratio of stars to brown dwarfs after correction based on the outer regions will thus be an upper limit.

A Galactic model could in principle be used to predict the field star contamination along the line of sight of the ONC (e.g. Robin et al. 2003). However, in practice this is difficult due to the lack of calibration of the models for faint magnitudes. The model can be used as a test of whether the control field reproduces the expected number of foreground field stars. Using the Galactic model by Robin et al. (2003) we find the predicted number of foreground field objects towards the control area (2.0–2.6 pc) to be 8.9. The measured number of blue objects (J − H < 0.8 mag) within the magnitude range 10 < H < 19 mag is 10, suggesting that the outer annulus traces the field star population well. Based on the agreement with the foreground population and the similarity in the extinction distribution as a function of radius, we adopt the outer annulus, 2.0–2.6 pc as the control region for both foreground and background stars.

4.4. The low-mass initial mass function

From the extinction limited samples we can construct an unbiased estimate of the IMF down to the completeness limit of the data. As illustrated in Fig. 2 the completeness in our inner bin decreases sharply below H = 18 mag but individual objects are still detected down to H = 20 mag. Thus, by restricting the sample to an extinction limit of AV = 7 mag and a lower mass limit of 20 MJup the sample is expected to be complete across the whole mass range for all radii. The completeness tests showed the inner annulus is incomplete for the 10–20 MJup mass range and this part of the mass function is thus an upper limit. In addition, as described below, the field contamination is substantial for the lowest mass range. We thus mainly focus the discussion of the IMF down to 30 MJup but present the corresponding ratios for the samples extending to 10 and 20 MJup.

Following Meyer et al. (2000) and Andersen et al. (2008) we have determined the ratio of low-mass stars (0.08–1 M) to brown dwarfs both for the cluster as a whole, and for 3 separate annuli outside a 0.35 pc radius. The results are summarized in Table 2 and the ratio of stars to brown dwarfs as a function of radius is shown in Fig. 7 where the ratios both with and without field stars subtracted are shown.

thumbnail Fig. 7

The ratio of low-mass stars (0.08–1.0 M) to brown dwarfs (0.03–0.08 M) for the ONC as a function of radius. The plus signs with error bars indicate the field star corrected measurements whereas the squares indicate the measurements before field star subtraction. The inner data point is from Slesnick et al. (2004) and the horizontal line is the average ratio with the 1σ error interval provided by the dashed lines.

The ratio determined at 0.35–0.7 pc is consistent with the results by Slesnick et al. (2004) for the central 0.3 pc radius of the cluster. However, for larger distances there is a drop in the ratio of stars to brown dwarfs. Outside a radius 0.7 pc the ratio drops by almost a factor of two, although the error bars for the individual annuli are relatively large. The total ratio of stars to brown dwarfs for the cluster as a whole has been calculated by correcting the source counts in the different annuli by the fraction of the cluster covered by that annulus (derived from Fig. 4). As shown by Hillenbrand & Hartmann (1998) the surface density is well represented by a King profile, and thus extends as r-2 outside the core radius. Thus, even though the surface density is lower at larger radii, this is compensated by an increase in the area in each annulus. The low ratio far from the centre drives the global ratio to a relatively low value of R03 = 2.4 ± 0.2, lower than determined in the cluster centre.

The drop in the ratio of stars to brown dwarfs as the distance from the centre increases is also seen for the sample with a mass limit of 20 MJup or 10 MJup instead of 30 MJup but due to the stronger contamination of field stars as the sample is extended to lower masses we focus the analysis on the higher mass limit data. Whereas the field contamination is some 10% for the brown dwarf content down to a mass limit of 30 MJup, it increases at lower masses and hence fainter magnitudes. This is a consequence of the underlying luminosity function for the field star population which extends as a power-law to faint magnitudes.

Could the lower ratio of stars to brown dwarfs as a function of distance from the ONC centre be due to age effects? Adopting the new lower distance to the ONC (Menten et al. 2007) and a more extensive sample of members where the spectral type has been determined through spectroscopy (Hillenbrand 1997) or through narrow band imaging of TiO absorption features, Da Rio et al. (2010) suggest that the cluster might be older than the previous estimate of 1 Myr. Recent work has further suggested an age gradient across the ONC (Reggiani et al. 2011), where the centre is dominated by a younger population whereas a wider age distribution is seen at larger radii with older stars having a larger average distance from the centre.

The effect of an older cluster at the same distance is that the pre-main sequence stars and brown dwarfs gets dimmer, an effect that is stronger for lower masses and it could thus replicate a decrease of the ratio of stars to brown dwarfs if the age adopted for the outer parts of the cluster is too low. Reggiani et al. (2011) find an average age for the (stellar) population between 1 pc and 1.5 pc to be log (age = 6.3). We have therefore also calculated the ratios of stars to brown dwarfs in the outer annuli for an age of 2 Myr. We find the ratio of low-mass stars to brown dwarfs only changes slightly from the values determined for a 1 Myr isochrone, R03(2 Myr,d = 1.1 − 1.5 pc)) = 1.1 ± 0.8, and R02(2 Myr,d = 1.1 − 1.5 pc) = 0.8 ± 0.6 and this will thus not change the conclusion that the ratio of stars to brown dwarfs decrease as a function of radius.

5. Discussion

We present the constraints that our observations place on the global IMF for the ONC and compare the results with the field IMF and results from other nearby star forming regions. Evidence for variations in the IMF as a function of environment is discussed.

5.1. Variations in the ratio of stars to brown dwarfs as a function of radius

The decreasing ratio of stars to brown dwarfs as a function of radius suggests that the cluster is mass segregated. A similar conclusion was found by Hillenbrand & Hartmann (1998) for the massive stars in the cluster centre. However, there was no indication for mass segregation for the lower mass population as their study only extended to 0.3 M and thus compared high mass stars to low-mass stars. In contrast, we find the sub-stellar members to be segregated compared to the stars as a whole, which are dominated by the lower mass members.

Hillenbrand & Hartmann (1998) concluded that it is possible for the high-mass segregation to have an dynamical origin. This is not the case for the brown dwarfs observed here with the current cluster surface density profile. The half mass relaxation time was estimated in Hillenbrand & Hartmann (1998) to be 6.5 Myr which is less than the range of ages inferred for the cluster and thus unlikely to cause the variations in the ratio of stars to brown dwarfs through dynamical 2 body evolution (cf. Bonnell et al. 1998).

A decrease of the ratio of stars to brown dwarfs as a function of radius has been observed in some open clusters as well, e.g. IC 2391 (Boudreault & Bailer-Jones 2009) where the brown dwarfs were found to have a flatter distribution than the low-mass stars. Similarly, for IC 348 and central parts of the ONC, Kumar & Schmeja (2007) use the minimum spanning tree to suggest the brown dwarfs are less centrally concentrated than the stars. A flatter distribution for brown dwarfs compared to stars is what would be expected from an ejection scenario where a fraction of the brown dwarfs are due to ejections in multiple systems during the earliest phases of star formation (e.g. Reipurth & Clarke 2001; Bate et al. 2003).

More recent simulation by e.g. Bate (2009); Bonnell et al. (2011) does not predict an overabundance of brown dwarfs further away from the centre when the simulations end. However, the simulations typically end after 1–1.5 free fall times which, depending on the specific simulation corresponds to 0.2–0.5 Myr, significantly younger than the current ONC. Thus, a more extended brown dwarf distribution could be generated by subsequent dynamical interactions in the core leading to a substantial unbound population or a higher overall velocity dispersion.

5.2. Comparison with other star forming regions

A smaller ratio of the number of stars to brown dwarfs for the whole cluster compared to the centre region has several implications. Bate (2009) predicted a value of R03 = 1.25, based on SPH simulations of a spherical collapsing cloud producing a total of 1254 objects, even lower than we determine here for the ONC. It is worth noticing that our value of R03 = 2.4 ± 0.2 is not consistent with a Chabrier IMF, which predicts too few brown dwarfs (R03(Chabrier) = 4.95). The probability of the ratio to be drawn from an underlying population with a Chabrier distribution is less than 1%. The ratio we determine is in better agreement with a flat or rising IMF in linear units into the brown dwarf regime. We stress however, that our IMF is determined from a photometric sample only and that spectroscopic confirmation of the brown dwarf candidates is highly desirable.

Andersen et al. (2008) presented a compilation of well determined ratios of stars to brown dwarfs in nearby young star forming regions and the Pleiades. The measurements did not show evidence for a variable underlying IMF as a function of star formation environment. Results for some OB associations, e.g. Upper Sco, and σ Ori, have indicated that the derived mass function is richer in brown dwarfs compared to the young cluster sample (e.g. Lodieu et al. 2007; Caballero et al. 2007). The derived low-mass slope is only weakly falling or flat into the brown dwarfs regime (a slope of  − 0.4 ± 0.2 in logarithmic units). For both the upper Sco survey and for σ Ori, the mass function was determined over a large field of view, several square degrees and the associations are thus probed to much larger radii than e.g. previous studies of the ONC. However, although for some of the embedded clusters only the central regions have been observed (Mon R2, ONC, NGC 2024), others have much wider regions surveyed and the ratio of stars to brown dwarfs is still relatively high (e.g. Taurus, IC 348).

The current finding of a possible change in the ratio of stars to brown dwarfs at larger radii suggests further studies of other nearby regions to larger radii than previously surveyed are needed. Also, obtaining kinematic data for the individual low-mass members to better characterise any mass segregation would be very valuable. Current near-infrared large field of view capabilities both from the ground and space are well suited for these studies.

6. Conclusions

We have presented deep NICMOS III F110W and F160W non-contiguous observations of the ONC region. The data reach J = 22 mag, and H =  21 mag, deeper than previous large scale high spatial resolution surveys. From the data we derive the ratio of low-mass stars (0.08–1 M) to brown dwarfs as a function of radial distance from the centre. The ratio of stars to brown dwarfs is lower than in the centre, falling from a value of 3.7 in the centre to 1.1 ± 0.8 in the annulus 1.1–1.5 pc. This shows the low-mass content of the cluster is mass segregated to large radii. The age of the cluster is much shorter than the half mass relaxation time suggesting the mass segregation for the low-mass population is primordial. The number of stars to brown dwarfs for the cluster as a whole is found to be N(0.081 M)N(0.030.08 M)=2.4±0.2\hbox{$\frac{N(0.08{-}1~M_\odot)}{N(0.03{-}0.08~M_\odot)}=2.4\,\pm\,0.2$}, which is in better agreement with a flat or rising low-mass IMF than a Chabrier IMF.

References

  1. Andersen, M., Meyer, M. R., Oppenheimer, B., Dougados, C., & Carpenter, J. 2006, AJ, 132, 2296 [NASA ADS] [CrossRef] [Google Scholar]
  2. Andersen, M., Meyer, M. R., Greissl, J., & Aversa, A. 2008, ApJ, 683, L183 [NASA ADS] [CrossRef] [Google Scholar]
  3. Andersen, M., Zinnecker, H., Moneti, A., et al. 2009, ApJ, 707, 1347 [NASA ADS] [CrossRef] [Google Scholar]
  4. Bally, J., Lanber, W. D., Stark, A. A., & Wilson, R. W. 1987, ApJ, 312, L45 [NASA ADS] [CrossRef] [Google Scholar]
  5. Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403 [NASA ADS] [Google Scholar]
  6. Bastian, N., Covey, K. R., & Meyer, M. R. 2010, ARA&A, 48, 339 [NASA ADS] [CrossRef] [Google Scholar]
  7. Bate, M. R. 2009, MNRAS, 392, 590 [NASA ADS] [CrossRef] [Google Scholar]
  8. Bate, M. R., Bonnell, I. A., & Bromm, V. 2003, MNRAS, 339, 577 [NASA ADS] [CrossRef] [Google Scholar]
  9. Bonnell, I. A., Bate, M. R., & Zinnecker, H. 1998, MNRAS, 298, 93 [NASA ADS] [CrossRef] [Google Scholar]
  10. Bonnell, I. A., Smith, R. J., Clark, P. C., & Bate, M. R. 2011, MNRAS, 410, 2339 [NASA ADS] [CrossRef] [Google Scholar]
  11. Boudreault, S., & Bailer-Jones, C. A. L. 2009, ApJ, 706, 1484 [NASA ADS] [CrossRef] [Google Scholar]
  12. Brandner, W., Clark, J. S., Stolte, A., et al. 2008, A&A, 478, 137 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Briceño, C., Luhman, K. L., Hartmann, L., Stauffer, J. R., & Kirkpatrick, J. D. 2002, ApJ, 580, 317 [NASA ADS] [CrossRef] [Google Scholar]
  14. Caballero, J. A., Béjar, V. J. S., Rebolo, R., et al. 2007, A&A, 470, 903 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  15. Carpenter, J. M. 2001, AJ, 121, 2851 [NASA ADS] [CrossRef] [Google Scholar]
  16. Chabrier, G. 2003, PASP, 115, 763 [NASA ADS] [CrossRef] [Google Scholar]
  17. Gouliermis, D. A., Henning, T., Brandner, W., et al. 2007, ApJ, 665, L27 [NASA ADS] [CrossRef] [Google Scholar]
  18. Guieu, S., Dougados, C., Monin, J.-L., Magnier, E., & Martín, E. L. 2006, A&A, 446, 485 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  19. Hillenbrand, L. A. 1997, AJ, 113, 1733 [NASA ADS] [CrossRef] [Google Scholar]
  20. Hillenbrand, L. A., & Carpenter, J. M. 2000, ApJ, 540, 236 [NASA ADS] [CrossRef] [Google Scholar]
  21. Hillenbrand, L. A., & Hartmann, L. W. 1998, ApJ, 492, 540 [NASA ADS] [CrossRef] [Google Scholar]
  22. Hillenbrand, L. A., Strom, S. E., Calvet, N., et al. 1998, AJ, 116, 1816 [NASA ADS] [CrossRef] [Google Scholar]
  23. Kumar, M. S. N., & Schmeja, S. 2007, A&A, 471, L33 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  24. Lodieu, N., Hambly, N. C., Jameson, R. F., et al. 2007, MNRAS, 374, 372 [NASA ADS] [CrossRef] [Google Scholar]
  25. Luhman, K. L. 2004, ApJ, 617, 1216L [NASA ADS] [CrossRef] [Google Scholar]
  26. Luhman, K. L., Rieke, G. H., Young, E. T., et al. 2000, ApJ, 540, 1016 [Google Scholar]
  27. Luhman, K. L., Mamajek, E. E., Allen, P. R., & Cruz, K. L. 2009, ApJ, 703, 399 [NASA ADS] [CrossRef] [Google Scholar]
  28. Menten, K. M., Reid, M. J., Forbrich, J., & Brunthaler, A. 2007, A&A, 474, 515 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  29. Meyer, M. R., Calvet, N., & Hillenbrand, L. A. 1997, AJ, 114, 288 [NASA ADS] [CrossRef] [Google Scholar]
  30. Meyer, M. R., Adams, F. C., Hillenbrand, L. A., Carpenter, J. M., & Larson, R. B. 2000, Protostars and Planets IV, 121 [Google Scholar]
  31. Muench, A. A., Lada, E. A., Lada, C. J., & Alves, J. 2002, ApJ, 573, 366 [NASA ADS] [CrossRef] [Google Scholar]
  32. Reipurth, B., & Clarke, C. 2001, AJ, 122, 432 [NASA ADS] [CrossRef] [Google Scholar]
  33. DaRio, N., Robberto, M., Soderblom, D. R., et al. 2010, ApJ, 722, 1092 [NASA ADS] [CrossRef] [Google Scholar]
  34. Reggiani, M., Robberto, M., Da Rio, N., et al. 2011, A&A, in press [arXiv:1108.1015] [Google Scholar]
  35. Robberto, M., Soderblom, D. R., Scandariato, G., et al. 2010, AJ, 139, 950 [NASA ADS] [CrossRef] [Google Scholar]
  36. Robin, A. C., Reylé, C., Derrière, S., & Picaud, S. 2003, A&A, 409, 523 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  37. Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 [NASA ADS] [CrossRef] [Google Scholar]
  38. Slesnick, C. L., Hillenbrand, L. A., & Carpenter, J. M. 2004, ApJ, 610, 1045 [NASA ADS] [CrossRef] [Google Scholar]
  39. Stolte, A., Brandner, W., Brandl, B., & Zinnecker, H. 2006, AJ, 132, 253 [NASA ADS] [CrossRef] [Google Scholar]

All Tables

Table 1

Sources detected in the NICMOS camera 3 mosaics.

Table 2

Ratio of stars to brown dwarfs in an extinction limited sample with AV = 7 mag as a function of lower mass limit.

All Figures

thumbnail Fig. 1

Upper and lower left: a comparison of the NICMOS photometry converted into the 2MASS system with the photometry obtained with ISPI, in the 2MASS system (Robberto et al. 2010). Only objects with photometric errors less than 0.1 mag in both filters are shown. Upper and lower right: the NICMOS photometric errors as a function of magnitude.

In the text
thumbnail Fig. 2

The completeness correction for different annuli as a function of F160W band magnitude. From the upper left is shown the completeness corrections for the inner annulus (0.35–0.7 pc), the middle annulus (0.7–1.1 pc), the outer annulus (1.1–1.5 pc), and in the lower-right a comparison of the three annuli.

In the text
thumbnail Fig. 3

F160W band 26′ × 33′ mosaic of the observed region within the ONC. For reference, the central 140″ × 140″ region observed by Luhman et al. (2000) is included. The location of the Trapezium stars are marked by small boxes in the mosaic. North is up and East is to the left.

In the text
thumbnail Fig. 4

The fraction coverage of the ONC by the NICMOS Camera 3 observations. The fraction covered is shown as a function of radial distance from the centre (θ1 Ori C).

In the text
thumbnail Fig. 5

The color–magnitude diagram for the surveyed region. Overplotted are a 1 Myr isochrone (Baraffe et al. 1998) shifted to the distance of the ONC, the same isochrone reddened by AV = 7 mag, and the 50% completeness limit for the data. Masses for a 1 Myr population of 1.0, 0.08, 0.02, and 0.01 M are indicated on the lower left panel. The region used as control field is shown in the lower right corner (2.3–2.6 pc). The typical photometric and colour error are shown in each diagram as a function of object magnitude.

In the text
thumbnail Fig. 6

The derived extinction for sources in the ONC within the magnitude limits 10.5 < H < 18 mag, which corresponds to the range from where the ratio of stars to brown dwarfs is determined.

In the text
thumbnail Fig. 7

The ratio of low-mass stars (0.08–1.0 M) to brown dwarfs (0.03–0.08 M) for the ONC as a function of radius. The plus signs with error bars indicate the field star corrected measurements whereas the squares indicate the measurements before field star subtraction. The inner data point is from Slesnick et al. (2004) and the horizontal line is the average ratio with the 1σ error interval provided by the dashed lines.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.