## The low mass star and sub-stellar populations of the 25 Orionis group

Juan Jos´ e Downes 1 2 , /star , C´ esar Brice˜o, n 1 3 , † , Cecilia Mateu, 1 2 , Jes´s Hern´ndez, u a 1 Anna Katherina Vivas 3 Nuria Calvet 4 Lee Hartmann 4 Monika G. Petr-Gotzens 5 and Lori Allen 6

- 1 Centro de Investigaciones de Astronom´ ıa, AP 264, M´rida 5101-A, Venezuela e
- 2 Instituto de Astronom´ ıa, UNAM, Ensenada, C.P. 22860, Baja California, M´xico e
- 3 Cerro Tololo Interamerican Observatory Casilla 603, La Serena, Chile
- 4 Department of Astronomy, University of Michigan, 825 Dennison Building, 500 Church Street, Ann Arbor, MI 48109, USA
- 5 European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748, Garching bei M¨nchen, Germany u
- 6 National Optical Astronomy Observatoies, 950 N. Cherry Ave. Tucson, AZ 85719, USA

Accepted 2014 -- -. Received 2014 -- -; in original form -- -- -

## ABSTRACT

We present the results of a survey of the low mass star and brown dwarf population of the 25 Orionis group. Using optical photometry from the CIDA Deep Survey of Orion, near IR photometry from the Visible and Infrared Survey Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec at the MMT, we selected 1246 photometric candidates to low mass stars and brown dwarfs with estimated masses within 0 02 . &lt; M/M ∼ /circledot &lt; . ∼ 0 8 and spectroscopically confirmed a sample of 77 low mass stars as new members of the cluster with a mean age of ∼ 7 Myr. We have obtained a system initial mass function of the group that can be well described by either a Kroupa power-law function with indices α 3 = -1 73 . ± 0 31 and . α 2 = 0 68 . ± 0 41 . in the mass ranges 0 03 . /lessorequalslant M/M /circledot /lessorequalslant 0 08 . and 0 08 . /lessorequalslant M/M /circledot /lessorequalslant 0 5 . respectively, or a Scalo log-normal function with coefficients m c = 0 21 . +0 02 . -0 02 . and σ = 0 36 . ± 0 03 in the mass range 0 03 . . /lessorequalslant M/M /circledot /lessorequalslant 0 8. . From the analysis of the spatial distribution of this numerous candidate sample, we have confirmed the EastWest elongation of the 25 Orionis group observed in previous works, and rule out a possible southern extension of the group. We find that the spatial distributions of low mass stars and brown dwarfs in 25 Orionis are statistically indistinguishable. Finally, we found that the fraction of brown dwarfs showing IR excesses is higher than for low mass stars, supporting the scenario in which the evolution of circumstellar discs around the least massive objects could be more prolonged.

Key words: stars: low-mass, brown dwarf, open clusters and associations: individual (25 Orionis)

## 1 INTRODUCTION

Evidence collected over the past several years increasingly supports the general picture that low mass stars (LMS) and brown dwarfs (BDs), likely share a common formation mechanism. However, there are many open questions that remain to be addressed. For example, given typical temperatures and densities within molecular cloud cores (e.g. Padoan &amp; Nordlund 2004), the corresponding classi-

- /star E-mail: jdownes@cida.ve,jdownes@astrosen.unam.mx
- † E-mail: briceno@cida.ve,cbriceno@ctio.noao.edu

cal Jeans mass can be much higher than the sub-stellar mass limit ( ∼ 0 072 . M /circledot ; Baraffe et al. 1998), suggesting that forming BDs would be difficult under such conditions. However, we know that not only do LMS and BDs form, but they do so abundantly, being the most numerous objects in the Galaxy (e.g. Padoan &amp; Nordlund 2004; Bastian, Covey &amp; Meyer 2010). Therefore, their formation cannot be explained in terms of a simple-minded model of gravitational collapse alone, which has been an important driver of theoretical and observational efforts looking for the physical processes and conditions that control LMS and BD formation.

Theoretical efforts to explain the formation of LMS and BDs are usually classified in the following models, reviewed by Whitworth et al. (2007): (i) turbulent fragmentation, (ii) gravitational fragmentation (iii) disc fragmentation (iv) premature ejection of protostellar embryos and (v) photo-erosion of cores in HII regions. These processes are not mutually exclusive (Whitworth et al. 2007) in the sense that in a typical star forming region (SFR) all of them could take place simultaneously, resulting in a complex scenario. The relative importance of each model is currently explored in terms of its ability to reproduce observables such as the initial mass function (IMF), the properties of binary and multiple systems, the fraction of LMS and BDs showing circumstellar discs, magnetosferic accretion and/or outflows and studying the kinematic or spatial segregation as a function of stellar mass (e.g. Whitworth et al. 2007). Each model has proven its plausibility of forming LMS and BDs, though only a few of the available simulations offer results directly comparable with observational data, and several assumptions about the initial conditions are still a matter of debate.

Observationally, several SFRs and young clusters have been surveyed looking for their LMS and BD population (e.g. Luhman 2012; Luhman et al. 2010; Boudreault &amp; Bailer-Jones 2009; Bouvier et al. 2008; Lodieu et al. 2007; Caswell 2007; Luhman et al. 2007, and references therein). Some important results are: the continuity of the mass function across the substellar mass limit (e.g. Bastian, Covey &amp; Meyer 2010), the existence of circumstellar discs (e.g. Luhman et al. 2010), accretion phenomena (e.g. Muzerolle et al. 2000; Jayawardhana, Mohanty &amp; Basri 2003) and photometric variability (e.g. Scholz &amp; Eisl¨ffel o 2005) in BDs, which indicate the extension of the T-Tauri behaviour below the sub-stellar mass limit. There is also some initial insight on mass segregation (e.g. Kumar &amp; Schmeja 2007; Boudreault &amp; Bailer-Jones 2009), mass-dependent kinematics (e.g. Joergens 2006) and multiplicity (e.g. McCarthy &amp; Zuckerman 2004; Luhman 2004; Kraus, White &amp; Hillenbrand 2006; Luhman et al. 2009).

Here we present a new, deep optical/near infrared photometric survey of the 25 Orionis group, combined with limited follow-up spectroscopy, to search for and characterise the LMS (0 8 . &lt; M/M ∼ /circledot &lt; . ∼ 0 072) and BD (0 072 . &lt; M/M ∼ /circledot &lt; . ∼ 0 02) members of this population. The relevance of the 25 Orionis group for studies of star/disc formation and early evolution is given by its ∼ 7 Myr old age which roughly corresponds to the time when dissipation of primordial circumstellar discs is expected to have finished ( ∼ 10 Myr; Calvet et al. 2005), as evidenced by the decrease in disc fractions (Hern´ndez et al. 2006, 2007). Ada ditionally, the study of the spatial distribution of members and/or the IMF down to the BD regime at this age, is important to constraint the dynamical evolution of young clusters (e.g. Boudreault &amp; Bailer-Jones 2009).

Other associations of similar age are either too spatially extended because of their vicinity (e.g. TW Hydra, β Pictoris, /epsilon1 Cha; Torres et al. 2008) or too distant (e.g. NGC7160, NGC2169, NGC2362 Kun, Kiss &amp; Balog 2008; Jeffries et al. 2007; Dahm 2008). Lacerta OB1b at a distance of ∼ 400 pc and an estimated age within 12 to 16

Myr, is probably one of the best places to carry out new surveys for young LMS and BD, however, it also spans ∼ 70 deg 2 (Chen &amp; Lee 2008). Another interesting region, also spanning a large area of the sky, is the Upper Scorpius association at a distance between 125 and 165 pc (Luhman &amp; Mamajek 2012), whose last age determination is ∼ 11 Myr old (Pecaut, Mamajek &amp; Bubar 2012). Summarizing, the 25 Orionis group is the most numerous and spatially dense ∼ 7 Myr old population yet known within ∼ 500 pc from the Sun. Its modest angular extent ( &lt; deg ∼ 7 2 ) makes comprehensive studies of its entire population feasible, and being almost free of significant interstellar extinction, it is the ideal place to study slightly more evolved young LMS and BD.

The paper is organized as follows: In Section 2 we summarize previous results on the 25 Orionis group. In Section 3 we describe the observations, data reduction, photometric candidate selection, as well as the spectroscopic confirmation of memberships for a sub-sample of these candidates. In Section 4 we estimate and analyse the visual interstellar extinction affecting each photometric candidate and spectroscopically confirmed member, and estimate the contamination from field stars present in the sample of photometric candidates. In Section 5 we present the derivation of masses and ages, the analysis of the IMF, the analysis of the spatial distribution and infrared excesses as disc indicators and accretion signatures. In Section 6 we summarize our conclusions about the nature of the 25 Orionis group and explore possible implications on predictions from theoretical models for LMS and BD formation.

## 2 THE 25 ORIONIS GROUP

The 25 Orionis group was originally identified by Brice˜o et al. n (2005) as a ∼ 1 ◦ radius overdensity of ∼ 20 pre-main sequence (PMS) stars roughly centred on the Be star 25 Ori located at α J 2000 =05:24:44.8, δ J 2000 =+01:50:47.20. It was recognized on the basis of multi-epoch V and I-band photometry and follow-up optical spectroscopy to confirm memberships. Their initial member sample spanned magnitudes 14 5 . &lt; V ∼ &lt; ∼ 16 5, . corresponding to masses in the range 0 45 . &lt; M/M ∼ /circledot &lt; . ∼ 0 95 1 . Brice˜o et al. n (2005) reported a mean age for the Orion OB1a sub-association, which includes the 25 Orionis overdensity, of ∼ 7 -10 Myr and an inner-disc frequency of ∼ 11 per cent, smaller than the ∼ 23 per cent they had found in the younger ( ∼ 3 6 . Myr) Orion OB1b sub-association. These results were consistent with the decrease in disc fraction with age seen in other regions (Hern´ndez et al. 2008), a and supported the idea that inner discs around low-mass stars dissipate in less than ∼ 10 Myr (Hartmann 2001).

The newly identified 25 Orionis group was also recognized a few months later by Kharchenko et al. (2005) in their extensive survey of new galactic clusters, and listed as cluster ASCC-16. They derive central coordinates α J 2000 =5:24:36, δ J 2000 =+01:48:00, a cluster radius=0 62 . ◦ and core radius=0 25. They estimate a distance of 460 pc . ◦

1 Assuming models from Baraffe et al. (1998), an age of ∼ 7 Myr and distance modulus of 7.78 mag

(V -M V = 8 59), reddening E(B . -V) = 0 09, age . ∼ 8 5 . Myr, and a mean radial velocity of 0 75 . ± 8 75 km/s. How-. ever, the radial velocity comes from only one star, 25 Ori itself, and the age, distance and cluster radius are based on just a small number of stars earlier than K3.

The low-mass population of the 25 Orionis group was later considered by McGehee (2006). His variability-based survey of Orion OB1 using multi-epoch photometry from the Sloan Digital Sky Survey (Abazajian et al. 2005), within a 2 5 wide equatorial strip located roughly 0 6 south of the . ◦ . ◦ star 25 Ori. Based on the spatial distribution of his and Brice˜o et al. n (2005) samples, McGehee (2006) suggested that the 25 Orionis group could extend further south, with a total radius of ∼ 1 4, . ◦ greater than previously reported. McGehee (2006) suggested that the 25 Orionis group is an unbound association and derived a Classical T-Tauri star (CTTS) fraction of ∼ 10 per cent in what he called the 25 Orionis group.

Discs around the low-mass (0 12 . &lt; M/M ∼ /circledot &lt; . ∼ 1 2) members of 25 Orionis were searched for and characterised by Hern´ andez et al. (2007) in their Spitzer IRAC and MIPS imaging study of this group. They found that in the 25 Orionis group, discs around low-mass stars are less frequent ( ∼ 6 per cent) than in the younger ( ∼ 4 Myr) Orion OB1b ( ∼ 13 per cent), and that this frequency seemed to be a function of the spectral type with a maximum around M0. They also determined that disc dissipation takes place faster in the inner region of the discs; as suggested by several disc models (e.g. Weidenschilling 1997; Dullemond &amp; Dominik 2004).

After the initial discovery and first studies, Brice˜o et al. (2007) conducted the most extensive assessn ment so far of the low-mass population (0 3 . &lt; M/M ∼ /circledot &lt; ∼ 1) in the 25 Orionis group. They concluded that not only is 25 Orionis a kinematically distinct entity, but it is also different in velocity space from the widely spread stellar population of Orion OB1a, within which it is located. These results pointed to the idea of the 25 Orionis group as a physical stellar group, akin to structures like the σ Ori cluster. From the spatial distribution of the low-mass members they found a maximum surface density of ∼ 128 stars/ deg 2 with masses M &gt; 0 5 . M /circledot , and derived a cluster radius of ∼ 7 pc. In this work they also redetermined the fraction of Classical TTauri stars (CTTS) using a more numerous sample, and find disc fractions of ∼ 6 per cent in 25 Orionis and ∼ 13 per cent in Ori OB1b, confirming the factor of 2 × decline in the accretor fraction seen by Brice˜o et al. (2005). n

The 25 Orionis group had been also the target for several other works: Using the Spitzer IRAC and MIPS data, Hern´ andez et al. (2006) studied the debris disc fraction in stars with spectral types earlier than F5 in the 25 Orionis group. Biazzo et al. (2011) measured abundances for 8 T Tauri stars with spectral types from K7 to M0 and derived an average metallicity [Fe/H]=-0.05 ± 0.05. Ingleby et al. (2011), have studied the X-ray and far ultraviolet fluxes of 25 Orionis and determined accretion rates for a set of members with spectral types K5 to M5. Most recently, 25 Orionis has been a hunting ground for groups looking for transiting planets among young solar-like stars (van Eyken et al. 2011; Neuh¨user et al. 2011). In Table 1 we summarize the a derived parameters we assume here for the 25 Orionis group, based on the aforementioned studies and the present work. In Figure 1 we show the spatial distribution of candidates and confirmed members from previous surveys.

## 3 OBSERVATIONS, CANDIDATE SELECTION AND MEMBERSHIP DETERMINATION

## 3.1 Optical photometry

Multi-epoch optical V, R, I-band, and H α observations across the entire Orion OB1 association (spanning ∼ 180deg ) were obtained as part of the CVSO (Brice˜o et al. 2 n 2005), being conducted since 1998 with the J¨rgen Stock u 1.0/1.5 Schmidt-type telescope and the 8000 × 8000 pixel QUEST-I CCD Mosaic camera (Baltay et al. 2002), at the National Astronomical Observatory of Venezuela. The instrument is composed of sixteen 2048 × 2048 Loral CCD devices set in a 4 × 4 array covering most of the focal plane of the telescope, and yielding a field of view ∼ 2 3 wide . ◦ in declination. The chips are front illuminated with a pixel size of 15 µ m, which corresponds to a scale of 1.02' per pixel. Each row of four CCDs in the N-S direction can be fitted with a different filter. The system was optimized to operate in drift-scan mode, in which the telescope is fixed at a constant declination and hour angle, and an object moves across each filter as a consequence of sidereal motion, resulting in a survey rate of ∼ 34 deg 2 per hour per filter. The scans for the CVSO were performed along strips centred at declinations -5 , ◦ -3 , ◦ -1 , +1 , +3 ◦ ◦ ◦ and +5 ◦ (for more details of the survey strategy see Brice˜o et al. 2005, 2007; n Downes et al. 2008).

Because the exposure time in the drift-scan mode observations is fixed 2 , this defines the usable magnitude range for each individual observation. At the bright end, our data saturate around magnitude 13 in the V, R and I bands, which corresponds to pre-main sequence (PMS) stars with masses ∼ 1 M /circledot , at the distance and age of the 25 Orionis group (Baraffe et al. 1998). At the faint end, the 3 σ limiting magnitudes for individual scans are V lim ∼ 19 7, . R lim ∼ 19 7, I . lim ∼ 19 5, with completeness magnitude . 3 of V com ∼ 18 7, . R com ∼ 18 7, . I com ∼ 18 2. . Again assuming Baraffe et al. (1998) models, and the age and distance of the 25 Orionis group ( ∼ 7 Myr, ∼ 360 pc from Section 5.1), the substellar limit is located at V sub ∼ 19 8, . R sub ∼ 18 7 . and I sub ∼ 17, which clearly shows that individual scans are not sensitive enough for carrying out a complete study of the least massive BD in the Orion OB1 association.

In order to increase the signal to noise ratio (SNR) and therefore the limiting and completeness magnitudes, so that we could reach well into the substellar regime at the distance of Orion, we coadded the individual observations following the procedure explained in more detail in Downes et al. (2008). Summarizing, the process is based in

2 The exposure time is set by how long it takes the star to cross the entire CCD, which in our instrument is ∼ 140s near the equator.

3 The magnitude for which the logarithm of the number of sources as a function of magnitude departs from linear behaviour

Table 1. General properties of members and candidates for the 25 Orionis group from the present and previous works.

| Reference                 | Candidates [N]   | Members [N]   | Mass interval                         | LMS CTTS [per cent]   | LMS Discs [per cent]   |   D [pc] | Radius [ ◦ ]   | Age [Myr]     | ¯ AV [mag]   |
|---------------------------|------------------|---------------|---------------------------------------|-----------------------|------------------------|----------|----------------|---------------|--------------|
| Brice˜o n et al. (2005)   | · · ·            | 31            | 0 . 45 < ∼ M/M /circledot < ∼ 0 . 95  | ∼ 3                   | ∼ 3                    |      330 | · · ·          | 7-10          | 0.28         |
| Kharchenko et al. (2005)  | 90               | 29 a          | M > 1 . 4 M /circledot b              | · · ·                 | · · ·                  |      460 | 0.62           | ∼ 8 . 5       | 0.27         |
| McGehee (2006)            | 62               | · · ·         | 0 . 06 < ∼ M/M /circledot < ∼ 0 . 8 c | ∼ 10                  | · · ·                  |      330 | 1.4            | · · ·         | · · ·        |
| Hern´ndez a et al. (2007) | 41               | 115           | 0 . 12 < ∼ M/M /circledot < ∼ 1 . 2   | · · ·                 | ∼ 6                    |      330 | · · ·          | · · ·         | · · ·        |
| Brice˜o n et al. (2007)   | · · ·            | 124 d         | 0 . 3 < ∼ M/M /circledot < ∼ 1        | ∼ 6                   | · · ·                  |      330 | 1              | ∼ 7           | 0.29         |
| This work                 | 1246             | 77            | 0 . 02 < M/M /circledot < 0 . 8       | 3 . 8 ± 0 . 5         | 4 . 1 ± 0 . 4          |      360 | 0.5            | 6 . 1 ± 0 . 8 | 0.30         |

a Considered by Kharchenko et al. (2005) as high probability members.

b The most massive object predicted by Baraffe et al. (1998) models at such age is 1 4 . M /circledot , corresponding to a M V = 4 85 while . the fainter object in the Kharchenko et al. (2005) sample is M V = 4 05. .

c We interpolate into the Baraffe et al. (1998) models the Bessel magnitudes obtained from the SDSS photometry using Davenport et al. (2006) conversions.

d A total of 47 members have kinematic information.

the packages Offline and DQ developed by the QUESTI collaboration (Baltay et al. 2002). In a first step, the QUEST-I data pipeline ( Offline ) automatically processes every single drift-scan considering each CCD of the mosaic as an independent device, corrects the raw images by bias, dark current and flat field, performs the detection of point sources, aperture photometry, and determines detector coordinates for each object. Then the software solves the world coordinate system by computing the astrometric matrices for each CCD of the mosaic, based on the USNOA-2.0 astrometric catalog (Monet 1998).

After the basic processing for every individual scan has been done with the Offline package, the second step is to coadd these reduced single observations, using the package DQ . Using the transformation matrices output by Offline , the DQ package computes the offsets and rotations needed to place images of the same area of the sky, coming from different single drift-scans, in the same reference frame and add them pixel by pixel to produce the final coadded driftscan. The coadded drift-scans used here are listed in Table 2, together with the mean coordinates, filter set, resulting limiting and completeness magnitudes for each photometric band after coadding, average FWHM of point sources, and number of individual scans included in the coadd.

In a third step, photometry and astrometry on the resulting coadded images are performed using an IRAF 4 and IMWCS-based pipeline . The maximum number of detec5 tions at the 3 σ threshold in the I-band images, corresponds to a mean spatial density of ∼ 1 point source per ∼ 340 arcsec square, resulting in a mean distance between nearest sources of about ∼ 22'. Because of this, source crowding is not a problem, and we can safely perform aperture photometry; we used an aperture radius equal to the average FWHM of the coadded images ( ∼ 3'; Table 2).

The zero point calibration of the instrumental photometry in the Cousins system for each coadded drift-scan

4 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.

5 IMWCS is the Image World Coordinate Setting Program described in Mink (2006)

was done by normalizing the instrumental magnitudes to a master catalogue of ∼ 10 5 non-variable stars 6 with calibrated magnitudes. This large set of secondary standards was produced from a set of single drift-scans obtained with the same instrument and filters under photometric conditions, calibrated with Landolt (1992) standard star fields (Mateu et al. 2012). The calibration has an RMS of 0.021, 0.037 and 0.042 magnitudes in the V , R and I bands respectively. The mean positional uncertainty as inferred from the offsets in coordinates between the optical bands is 0.4 ′′ , with a standard deviation of 0.6 ′′ , similar to what we find for the offsets with respect to the infrared datasets we have used in this work ( § 3.2). Therefore, the photometric and astrometric calibration offer adequate accuracy for a reliable selection of LMS and BD candidates.

The catalogues with calibrated photometry for each photometric band from different coadded drift-scans, were finally compiled into a single deep catalog for each band, extracting repetitions and keeping, in such cases, always the measurement showing the highest SNR . We matched our 7 optical source lists using a search radius of 3 ′′ , always keeping the object with the minimum positional difference. We will refer to this collection of four deep catalogs (in the V, R, I passbands and in a narrow band - 10 nm wide - H α filter), as the CIDA Deep Survey of Orion (CDSO, Downes et al. 2014 in preparation). The sources in the CDSO span the region 65 ◦ /lessorequalslant α J2000 /lessorequalslant 100 ◦ and -6 ◦ /lessorequalslant δ J 2000 /lessorequalslant +6 , with ◦ positions and photometry for 8.6 million objects in the Iband, 8.1 million objects in the R-band, 6.9 million objects in the V-band, and 2.4 million objects in the H α filter.

In this contribution we focus on the LMS and BD populations of the 25 Orionis group, therefore we limit our analysis and discussion to drift-scans centred at declinations +1 ◦

6 Since we used individual scans spanning a time baseline of 7 years, we could determine which stars did not show hints of variability above the 0.02 mag level, and use these as secondary standards for our photometric calibration. Because of this, and the large number of stars used, we obtained a very robust calibration.

7 All positional matching between object catalogs throughout this paper was done using routines based on the the Starlink Tables Infrastructure Library Tool Set (STILTS, Taylor 2006) available at http://www.starlink.ac.uk/stilts/

Table 2. Observation log for the coadded drift-scans.

| Drift-Scan     | δc J 2000[ ◦ ]   |   α ini J 2000[ ◦ ] |   α end J 2000[ ◦ ] | Filters     |   FWHM a [ ′′ ] |   N b |   Vcom c | Rcom   |   Icom |   V lim | R lim   |   I lim |
|----------------|------------------|---------------------|---------------------|-------------|-----------------|-------|----------|--------|--------|---------|---------|---------|
| Single         | +1, +3           |               65    |               99    | All filters |             2.5 |     1 |     18.7 | 19.0   |   18.5 |    20.5 | 20.5    |    20   |
| Coadds 011-012 | +1               |               76.14 |               86.3  | VIIV        |             2.6 |    28 |     20.5 | . . .  |   20.3 |    22.3 | . . .   |    21.8 |
| Coadd 013      | +1               |               77.11 |               82.56 | VIRI        |             2.8 |    29 |     20.5 | 20.8   |   20.3 |    22.3 | 22.3    |    21.8 |
| Coadd 940      | +1               |               67    |               98.22 | RH α IV     |             2.9 |    13 |     20.1 | 20.4   |   19.9 |    21.9 | 21.9    |    21.4 |
| Coadd 950      | +3               |               71.15 |               89.41 | VRIV        |             2.6 |    10 |     19.9 | 20.2   |   19.7 |    21.7 | 21.7    |    21.2 |

- a The mean FWHM is from point sources in the I-band images.
- b The number of added single scans in the I-band.
- c The limiting and completeness magnitudes were computed for a 3 σ detection threshold.

and +3 , in particular within the region 0 35 ◦ . ◦ &lt; δ &lt; ∼ ∼ 3 35 and . ◦ 79 7 . ◦ &lt; α &lt; ∼ ∼ 82 7 shown in Figures 2, which encompasses the . ◦ area expected to be spanned by the 25 Orionis group and its surroundings. In this area the I-band catalog has a limiting I lim ∼ 22 and completeness I com ∼ 19 6, . which is enough to detect objects down to masses ∼ 30 M Jup , assuming Baraffe et al. (1998) models for the assumed age and distance of 25 Orionis. The full analysis of the entire CDSO will be presented in a forthcoming article.

0.6

## 3.2 Infrared photometry

## 3.2.1 Near Infrared Data

As will be explained in Section 3.3, our technique of searching for young BDs and LMS in the extended off-cloud, lowextinction regions of the Orion OB1 association hinges on the combination of deep optical and near-infrared photometry. Ever since its release, the Two Micron All-Sky Survey (2MASS Skrutskie et al. 2006) has been the tool of choice for all large scale studies that required access to near infrared photometry for large numbers of sources over extended areas of the sky, with uniform spatial coverage. However, with limiting magnitudes of J lim ∼ 17 4, H . lim ∼ 16 7, . Ks lim ∼ 16 2 . and completeness magnitudes J com ∼ 16 6, . H com ∼ 15 8, Ks . com ∼ 15 4, the 2MASS catalog is shallow . compared to the deep CDSO data, imposing restrictions on our search for BD in Orion, in particular for the faintest, lowest mass objects. Nevertheless, being at the time of our initial candidate selection the only spatially complete nearIR dataset available, we cross matched our I-band CDSO catalog with 2MASS using a search radius of 3 ′′ , resulting in a master catalog with IJHKs photometry for each source. With this dataset we created the optical-near-IR colourmagnitude diagrams used to select LMS and BD candidates (see § 3.3) for follow-up spectroscopy.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-b926dec3-f8ef-480e-b10a-581d8cd267f7.jpg)

Recently, a new, much deeper near-IR survey has become available. During 2009 a new dedicated 4m survey telescope, the Visible and Infrared Survey Telescope for Astronomy (VISTA), located at ESO's Paranal Observatory, was commissioned by the VISTA consortium (Emerson et al. 2004; Emerson &amp; Sutherland 2010). It is equipped with a 67 × 10 6 pixels camera with a scale of 0.34 ′′ /pix and a field of view of 1 65 in diameter. For the Galactic Science Verifi-. ◦ cation of VISTA, a ∼ 30deg 2 area of the Orion OB1 association, which included the Orion Belt region, part of the Orion A cloud, the 25 Orionis and σ Ori clusters, was imaged in the Z, Y, J, H and Ks filters, during October 16 to Novem-

Figure 1. Spatial distribution of photometric candidates and confirmed members of the 25 Orionis group from previous studies. Solid squares: LMS confirmed as members from Brice˜o et al. n (2005, 2007). The dotted horizontal line indicates the northern limit of the survey from Brice˜o et al. (2005). Small circles: n Photometric and astrometric candidates from Kharchenko et al. (2005). Crosses: members and candidates from Hern´ndez et al. a (2006). Triangles: Photometric candidates from McGehee (2006). The solid horizontal line indicates the northern limit of the survey of McGehee (2006). Diamonds: Eclipsing systems and/or CTTS candidates from van Eyken et al. (2011). Large circles indicate from the centre outward: the core and corona cluster radii (0 . ◦ 25 and 0 62 respectively) . ◦ computed by Kharchenko et al. (2005), the radius suggested by Brice˜o et al. n (2005) ( ∼ 1 ◦ ) and the radius proposed by McGehee (2006) ( ∼ 1 4); the dashed square . ◦ indicates the spatial limits of the IRAC survey. Background: Map of the extinction A V from Schlegel, Finkbeiner &amp; Davis (1998) with scale on the right. The big central asterisk represents the 25 Ori star and the black dots other clusters from Kharchenko et al. (2005).

ber 2, 2009, down to 5 σ limiting magnitudes Z ∼ 22 5, Y . ∼ 21 2, . J ∼ 20 4, . H ∼ 19 4, . Ks ∼ 18 6 . (Petr-Gotzens et al. 2011). The data reduction was performed using a dedicated pipeline, developed within the VISTA Data Flow System, and run by the Cambridge Astronomy Survey Unit in the UK. The pipeline delivers science-ready stacked images and mosaics, as well as photometrically calibrated source cat-

0.6

Figure 2. Spatial distribution of LMS and BD of the 25 Orionis group from the present work. Crosses: LMS and BD photometric candidates ( § 3.3). Circles: LMS confirmed as members ( 3.5). § The dashed circle indicates the field observed with Hectospec ( § 3.4). Other elements are as in Figure 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-cc1ffd41-50b6-46ed-8a8c-ba5daa7cf00f.jpg)

alogs 8 . The final, band-merged Orion OB1 VISTA catalog contains ∼ 3 million sources. The VISTA photometric calibration is based on (but different from) 2MASS. Instrumental magnitudes were converted into apparent calibrated magnitudes on the VISTA system and calibrated via colour equations. This also included the Z and Y bands, following the procedure described in Hodgkin et al. (2009). The point sources detected in the I-band CDSO survey catalog were cross-matched with the VISTA Orion OB1 master catalog, using a 3 ′′ match radius. The resulting IZYJHKs catalog was then matched with our V, R and H α catalogs to add additional measurements in these bands for those objects which had them; however, the VRH α data were not used in the process of selecting LMS and BD candidates. Because VISTA goes much deeper than 2MASS, and since for objects in common the VISTA and 2MASS magnitudes agree within the uncertainties, we adopted the CDSO-VISTA dataset as our final optical-near IR photometry catalog. For objects showing J &lt; 12, H &lt; 12 or Ks &lt; 11 in the VISTA catalogue we adopted J, H and Ks magnitudes from 2MASS in order to avoid possible saturation effects in VISTA data.

## 3.2.2 Mid-infrared Data

In order to identify disc-bearing sources among selected candidate LMS and BDs, we matched our multiband VRH IZYJHKs catalog with our existing Spitzer observaα tions within the 25 Orionis region obtained as part of GO3437 (Hern´ndez et al. 2007) and GO-50360 (Brice˜o et al. a n 2014, in preparation) and also with the recently released

8 More details on the VISTA data processing can be found at http://casu.ast.cam.ac.uk/surveys-projects/vista/technical.

Wide Infrared Survey Explorer (WISE Wright et al. 2010) database. First, a sub-sample of the photometric candidates and of the spectroscopically confirmed members (Sections 3.3 and 3.5), were identified in the object catalogs extracted from our IRAC images. The image reduction, point source detection and photometry for the datasets from Hern´ andez et al. (2007) and Brice˜o et al. (2014, in preparan tion) are analogous and were described in Hern´ndez et al. a (2007). The IRAC observations were performed in the passbands centred at 3.6 µ m, 4.5 µ m, 5.8 µ m and 8.0 µ m with the spatial coverage shown in Figure 1. As a second step, the photometric candidates and spectroscopically confirmed members were identified in the WISE database which offers measurements in four pass-bands, centred at 3.4 µ m, 4.6 µ m, 12 µ m and 22 µ m. Compared to the IRAC observations, WISE is a shallow survey with worse angular resolution, but its spatial coverage is uniform, providing also data for the surroundings of 25 Orionis, where we did not have Spitzer data. The limiting and completeness magnitudes for each VISTA, IRAC and WISE photometric bands are summarized in Table 3 together with the number of photometric candidates and confirmed members from this work with available photometry from these IR surveys.

## 3.3 Selection of photometric candidates

The selection of photometric candidates was performed according to their position in colour-magnitude diagrams which combine optical and near infrared magnitudes. This method has been extensively used and has proved highly successful in identifying LMS and BDs in young clusters and star-forming regions (e.g Brice˜o et al. 2002; Luhman et al. n 2003b; Downes et al. 2008).

We selected as candidates those objects located in what we call the membership locus , simultaneously in the I vs IJ, I vs I-H and I vs I-Ks colour-magnitude diagrams. These loci were defined as regions centred on the isochrone corresponding to the mean age from the Baraffe et al. (1998) models ( ∼ 7 Myr) reported by Brice˜o et al. n (2005, 2007) and a distance of ∼ 360 pc (Brice˜o et al. 2014, in preparan tion). The width of the loci were computed considering the 1 σ uncertainties in age and distance, the mean interstellar extinction throughout the region, unresolved binarity, the 1 σ photometric uncertainties and the expected mean intrinsic photometric variability. The resulting locus for the I vs. I-J diagram is shown in Figure 3. An object was considered as a photometric candidate if it was located within the membership locus in all three diagrams. Using the CDSOVISTA catalog, a total of 1246 photometric candidates were selected inside the region showed in Figure 2, spanning a ∼ 2 5 . ◦ × ∼ 2 5 area covering the previously known population . ◦ of the 25 Orionis group and its surroundings.

The estimated masses of the full candidate sample range from ∼ 0 8 . M /circledot down to ∼ 0 02 . M /circledot according to the models of Baraffe et al. (1998), with an age of ∼ 7 Myr and distance modulus of 7.78 mag. Due to the sensitivity and saturation of the I-band observations, the sample is complete from ∼ 0 8 . M /circledot ( I ∼ 13) down to ∼ 0 03 . M /circledot ( I ∼ 19 6). . In Section 5.1 we present the estimation of ages and masses for each candidate.

We emphasize that colour-magnitude diagrams constructed combining optical and near infrared magnitudes

Table 3. Limiting and completeness magnitudes and the number of candidates and members from VISTA, IRAC and WISE surveys.

| Survey   | Band     |   Completeness |   Limiting |   Candidates a |   Members |
|----------|----------|----------------|------------|----------------|-----------|
| VISTA    | Z        |           22.5 |       23.5 |           1246 |        77 |
|          | Y        |           20.8 |       22.5 |           1246 |        77 |
|          | J        |           20.2 |       21.7 |           1246 |        77 |
|          | H        |           19.2 |       20.5 |           1246 |        77 |
|          | Ks       |           18.4 |       19.6 |           1246 |        77 |
| IRAC     | 3 . 6 µm |           17.2 |       19.5 |            369 |        72 |
|          | 4 . 5 µm |           16.7 |       19   |            388 |        77 |
|          | 5 . 8 µm |           15.7 |       18.5 |            365 |        70 |
|          | 8 . 0 µm |           14.7 |       17.5 |            358 |        77 |
| WISE     | 3 . 4 µm |           16.2 |       17.2 |           1038 |        77 |
|          | 4 . 6 µm |           16   |       17   |           1038 |        77 |
|          | 12 µm    |           12.8 |       13.2 |           1038 |        77 |
|          | 22 µm    |            9.1 |        9.5 |           1038 |        77 |

- a VISTA and WISE surveys cover the entire region while IRAC survey covers the area shown in Figure 1.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-2d455b95-2dc8-4175-8221-272ee03b323a.jpg)

ness bias because BDs at this age and distance are red and very faint, so they drop out quickly in the bluer bands. Second, although our selection is not exempt from contamination by field stars, candidate selection performed in diagrams that include only near infrared data strongly increase the contamination by field dwarfs and by extragalactic sources, because of the limited leverage in colour; in these diagrams the isochrones are almost vertical and piled up, therefore, there is hardly any separation between the PMS locus and the field population. The contamination in our candidate sample by field stars and extragalactic sources is discussed in Section 4.2. The photometric variability in young LMS and BDs affects the optical and near-IR bands in a different way. Because our I-band data is the result from a coadd of images obtained in a period of 15 months, we considered that as the mean I-band magnitudes. The reported amplitudes of the variations in J, H and Ks-bands for LMS and BDs are within 0.2 to 0.5 mag (e.g. Scholz et al. 2009) We add to the width of the membership locus the corresponding change in colors that result from the maximum variation of 0.5 magnitudes in the near-IR bands.

Figure 3. Color-magnitude diagram used as part of the candidate selection, including the objects selected as candidates (empty circles) and those confirmed as new members (big dots). The thick horizontal line indicates the completeness limit. The solid curves indicates evolutionary tracks for 0.8 M /circledot , 0.7 M /circledot , 0.5 M , 0.4 M /circledot /circledot , 0.1 M /circledot , 0.06 M /circledot , 0.05 M /circledot , 0.03 M /circledot , 0.02 M /circledot , 0.01 M /circledot (thin) and 0.072 M /circledot (thick). The thick dashed line indicates the isochrone corresponding to 7 Myr and the thin doted lines enclose the membership locus. Theoretical isochrones and evolutionary tracks are from Baraffe et al. (1998) models, are ploted without extinction correction and assuming m M -= 7 78 magni-. tudes. The dereddening vector for the mean extinction throughout the region ( A V = 0 3) is also indicated. .

yield the best selection of photometric PMS stellar candidates and young sub-stellar objects mainly because of two reasons: first, diagrams based on optical data alone (V, R and I-band magnitudes) suffer from an important complete-

## 3.4 Optical spectroscopy

Even with a refined selection technique, photometric candidate samples suffer from some degree of contamination by non-members like foreground field dwarfs and also by background giant stars and extragalactic objects. Therefore, spectroscopic follow up is necessary to unambiguously confirm membership. Moreover, spectra also provide crucial information necessary to derive fundamental properties, like the object's effective temperature, A V , and hence its luminosity, whether it is accreting from a circumstellar disc, if it has a jet, and so on. Here we have extended in number and to much lower masses the member list discussed by Brice˜o et al. (2005, 2007), by obtaining low resolution n (R ∼ 1000) optical spectra of a sample of 90 photometric candidates of the 25 Orionis group.

The spectra of the sample were obtained during the

night of October 22 2007, during queue mode observations with the Hectospec multi-fiber spectrograph on the 6.5 MMT telescope (Fabricant et al. 2005). We observed a single field centred at α J 2000 =05:25:13, δ J 2000 =+01:39:25 (Figure 2), using one single fiber configuration for which we obtained five consecutive exposures, each 1800 s long. We used the 270 groove mm -1 grating, providing a spectral resolution of R ∼ 1000 with a spectral coverage from 3700 ˚ A to 9150 ˚ . Because of the poor seeing conditions ( A ∼ 2 2 . ′′ ) during the integrations and the diameter of the Hectospec fibers ( ∼ 1 5 . ′′ ), only the brightest candidates ( I &lt; ∼ 17) could be observed with reasonable SNR ( ∼ 15), allowing for an appropriate spectral classification.

The 25 Orionis group is placed far enough from the Orion molecular clouds, so that problems due to the high sky background caused by nebulosity are not present. All the Hectospec spectra were processed, extracted, corrected for sky lines, and wavelength calibrated by S. Tokarz at the CfA Telescope Data Center, using customized IRAF routines and scripts developed by the Hectospec team (see Fabricant et al. 2005). The five individual 1800 s exposures were coadded into a single, deep observation, with a total integration time of 9000 s. Light sources inside the Hectospec fiber positioner, contaminate ∼ 1 3 of the fibers at / wavelengths longer than ∼ 8500 A (N. Caldwell, ˚ private communication), therefore, during the analysis we ignored data at these wavelengths.

Of the 1246 photometric candidates, 277 are placed within our Hectospec field. Because of limitations with the positioning of fibers and the poor seeing conditions we observed a total of 90 candidates out of which 77 were confirmed as new members of the 25 Orionis group for a ∼ 86 per cent efficiency. The remaining 13 candidates were confirmed as field stars. An example of the spectra is shown in Figure 4.

## 3.5 Spectral classification and membership diagnosis

We performed the spectral classification by comparing the equivalent widths of sixteen spectral features typical of M type stars, such as TiO and VO bands, (see Table 3 of Downes et al. 2008), with their corresponding values in a library of standard spectra of field dwarfs from Kirkpatrick et al. (1999), and following the semi-automated scheme of Hern´ndez et al. (2004) a 9 . The differences in surface gravity between the young targets and the field dwarfs used as standards did not affect significantly the spectral classification for objects earlier than M6. We estimate the typical uncertainty resulting from our classification to be ± 0.5 spectral subclass.

We have used several criteria to assign 25 Orionis group membership to each spectroscopically observed candidate. We start by considering the detection of H α in emission. Strong H α in emission is a common feature of chromospherically active young objects, and of objects accreting from a circumstellar disc for the largest equivalent widths (e.g.

9 This procedure was done using the code SPTCLASS, available at http://www.astro.lsa.umich.edu/ ∼ hernandj/SPTclass/sptclass.html

Figure 4. Spectra of a LMS confirmed in this work as a new member of the 25 Orionis group. The top panel shows the complete spectrum obtained by Hectospec and the lower panels the H α emission line, the LiI λ 6708 absorption line and the NaI λ 8183 8195 absorption doublet. Labels indicate the member ID, , the spectral type and the equivalent width of H α . Vertical dashed lines in the lower panels indicate the central wavelength of the H , LiI α λ 6708 and NaI λ 8183 8195. , In the panel for the NaI λ 8183 8195 , absorption doublet we also show the corresponding spectra of a field dwarf (dashed red line; Kirkpatrick et al. (1999)), and a young dwarf from the Chamaeleon I region ( ∼ 2 Myr ; dotted line; Luhman et al. (2005)) of the same spectral type and smoothed to a same spectral resolution of 16 ˚ . The A spectra for the entire sample of new members (see Table 4) are available in the electronic edition.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-3401ea55-d08e-4b7a-a548-2229f339acb1.jpg)

Muzerolle et al. 2005). Figure 4 shows the H α emission line profile and its corresponding equivalent width for a confirmed member. At the age of the 25 Orionis group ( ∼ 7 Myr), most stars have stopped accreting (at least at significant levels), such that strong H α emission is not necessarily expected (Brice˜o et al. n 2007), and many members could exhibit equivalent widths similar to those measured in active dMe stars. Therefore, as an additional membership indicator we looked at the NaI λ 8183 8195 absorption doublet. , The NaI λ 8183 8195 feature is sensitive to the surface grav-, ity of the star, being very strong (deep) in old field dwarf stars, very weak in late type M giants, and intermediate in late type PMS stars and young BDs which are still contracting down to the main sequence or asymptotically approaching the main sequence in the case of BDs (Luhman et al. 2003a). The NaI λ 8183 8195 absorption doublet of one of , the LMS confirmed as member, superimposed on spectra of field dwarfs from Kirkpatrick et al. (1999), and of young M dwarfs of the same spectra type (from Luhman et al. 2003b) are shown in Figure 4. Given the typical I-J, I-H and I-Ks colours of the LMS and BD candidates considered here and the low visual extinction along the line of sight towards the 25 Orionis group, no reddened late giant stars are expected to contaminate our candidate sample, as we show in Section 4.2.

In addition to the above criteria, the presence of the Li I λ 6708 line in absorption was also considered. This spectral feature is a well-known indicator of youth in late K and Mtype objects (Strom et al. 1989; Brice˜o et al. 1998, 2001), n up to a spectral type ∼ M6.5, which at an age of ∼ 10 Myr corresponds to the mass limit for Li fusion ( M ∼ 0 06 . M /circledot Basri 1998). However, because of the low spectral resolution and SNR in most of our spectra, we could detect the LiI λ 6708 line in only ∼ 60 per cent of the confirmed members. Objects showing Li I λ 6708 line in absorption were selected as members; those without a detection were not discarded as members unless they also failed the Na I criteria explained above.

Summarizing, photometric candidates were classified as members if they showed H α in emission and NaI λ 8183 8195 , in absorption weaker than a field dwarf of the same spectral type. The detection of the Li I λ 6708 line in absorption was only considered as an additional youth indicator, in light of the bias imposed by the low spectral resolution and the SNR in most spectra. Additionally, as a sanity check, we looked at the A V of each star, computed from the observed colours and from the intrinsic colours we obtained from the spectral types and the relations from Luhman et al. (2003b), and noted whether it was consistent with the average A V for the 25 Orionis group (Brice˜o et al. 2007). n

Clearly, the spectral features we have considered above are youth indicators but do not, by themselves, guarantee membership to the 25 Orionis group. The reason is that these same criteria will select young members that belong to other regions of the Orion OB1 association, placed behind the 25 Orionis group. Indeed, we know that 25 Orionis is located in the extended OB1a sub-association, and the much younger OB1b region is a few degrees away. However, as we argue in Section 4.2 the expected number of such possible young contaminants is low and does not affect the general analysis presented here. Additionally, the Hectospec field is centred in the 25 Orionis spatial overdensity we explain in Section 5.2, increasing the probability of these objects being real members of 25 Orionis group.

We confirmed 77 new members of the group whose photometric information and resulting spectral types are listed in Tables 4 and 5. In Figure 5 we show the H-R diagram with the new members.

## 4 INTERSTELLAR EXTINCTION AND CONTAMINATION BY FIELD STARS

One of the problems encountered when carrying out photometric surveys over large areas of the sky is obtaining followup spectroscopy for every single source of interest, thereby making sure one has positively identified every member of the stellar population of interest. Even with the multiplexing capabilities of modern day multi-fiber spectrographs, this can be a daunting task, specially if one needs to use very long exposure times in order to get enough SNR on very faint objects. However, in our case, we can extend our study of the 25 Orionis group beyond the spatial and mass completeness of the spectroscopically confirmed member sample, by a judicious use of the photometric candidate member sample. At the expense of a larger degree of contamination from non-members, the photometric sample has the virtue

Figure 5. H-R diagram for the new members. The isochrones for 1, 3.2, 4.6, 8.5, 11.4 and 100 (dotted) Myr and the evolutionary tracks for 0.8 M /circledot , 0.6 M /circledot , 0.5 M /circledot , 0.3 M /circledot , 0.2 M /circledot , 0.1 M /circledot and 0.072 M /circledot (dashed) from the Baraffe et al. (1998) models are indicated. The top axis indicates the corresponding spectral types from Luhman et al. (2003b).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-bcf031b0-d400-435f-8b5c-ff65a7dff32d.jpg)

of providing a spatially complete rendering of the 25 Orionis aggregate, and it also allows us to sample large numbers of objects down to much lower masses. The disadvantage of this approach is that issues like contamination from nonmembers and variable extinction could potentially render any reliable analysis as useless. However, we will show that in the particular case of the 25 Orionis group, conditions are such that a photometric study of the group properties is not only feasible, but quite reliable.

## 4.1 Interstellar extinction

We determined the visual extinction A V for each of the 77 spectroscopically confirmed members, calculated by using the intrinsic I-J colour for the corresponding spectral type proposed by Luhman et al. (2003b), the observed I-J colour, and the Cardelli, Clayton &amp; Mathis (1989) extinction law assuming R V = 3.09. The resulting A V values are in excellent agreement with previous spectroscopic determinations by Brice˜o et al. (2005) for stellar members, as we n show in Figure 6.

The A V for the 1246 photometric candidates were derived using the Schlegel extinction maps (SEM; Schlegel, Finkbeiner &amp; Davis 1998). This approach is valid in this particular region because it has a spatially smooth extinction. Still, estimates from SEM have two main limitations: First, the A V reported in the SEM is the integrated extinction along the line of sight, limited by the sensitivity of the survey and the resulting A V for a nearby object will be overestimated, due to extinction sources behind the object. Second, the low spatial resolution of the SEM could result in an inadequate interpolation for objects placed in regions with abrupt variations of the extinction.

Table 4. Photometric catalog of the new spectroscopically confirmed LMS of the 25 Orionis group. The table is published in its entirety in the electronic edition.

|   ID |   α J2000 [ ◦ ] |   δ J2000 [ ◦ ] | V C   |   R C |   I C |     Z |     Y |     J |     H |    Ks | 3.6 µ m   | 4.5 µ m   | 5.8 µ m   | 8.0 µ m   |
|------|-----------------|-----------------|-------|-------|-------|-------|-------|-------|-------|-------|-----------|-----------|-----------|-----------|
| 1232 |         81.2583 |         1.62249 | 16.25 | 15.18 | 13.69 | 13.35 | 12.95 | 12.36 | 11.71 | 11.5  | 11.45     | 11.24     | 11.38     | 11.18     |
| 1242 |         81.0373 |         1.73113 | 16.11 | 15.09 | 13.66 | 13.18 | 12.68 | 12.07 | 11.57 | 11.39 | 11.17     | 11.11     | 11.07     | 11.10     |
| 1245 |         81.0633 |         1.90166 | 20.61 | 19.24 | 16.9  | 16.32 | 15.56 | 14.98 | 14.51 | 14.16 | 13.77     | 13.72     | 13.67     | 13.65     |
| 1246 |         81.2779 |         1.79876 | 19.16 | 17.83 | 15.65 | 15.12 | 14.46 | 13.89 | 13.4  | 13.11 | 12.69     | 12.64     | 12.57     | 12.54     |
| 1247 |         81.0647 |         1.931   | 17.22 | 16.03 | 14.28 | 13.85 | 13.28 | 12.75 | 12.25 | 11.96 |           |           |           |           |
| 1248 |         81.2433 |         1.97599 | 17.45 | 16.38 | 14.73 | 14.33 | 13.86 | 13.36 | 12.77 | 12.56 |           | 12.22     |           | 12.17     |
| 1252 |         81.2708 |         1.61483 | 16.92 | 15.81 | 14.1  | 13.64 | 13.17 | 12.61 | 12.04 | 11.81 | 11.49     | 11.49     | 11.48     | 11.47     |
| 1253 |         81.1664 |         1.3613  | 15.29 | 14.36 | 13.38 | 12.96 | 12.63 | 12.3  | 11.6  | 11.42 | 11.29     |           | 11.25     |           |
| 1258 |         81.2161 |         1.54132 |       | 18.58 | 16.29 | 15.72 | 14.98 | 14.38 | 13.91 | 13.56 | 13.17     | 13.10     | 13.11     | 13.14     |
| 1261 |         81.0869 |         1.4014  | 15.69 | 14.74 | 13.62 | 13.32 | 13    | 12.47 | 11.77 | 11.65 | 11.47     |           | 11.40     |           |

Table 5. Spectroscopic catalog of the new confirmed members of the 25 Orionis group. The table is published in its entirety in the electronic edition.

| ID                                                | ST           | W[H α ] [ ˚] A                       | LiI a         | Av [mag]               |
|---------------------------------------------------|--------------|--------------------------------------|---------------|------------------------|
| 1278 1273 1262 1277 1284 1269 1263 1266 1279 1242 | M4.5 M5 M4.5 | -5.8 -13 -8.4 -4.3 -9.9 -5 -8.9 -8.9 | 0 1 1 1 1 1 1 | 0 0.12 0 0 0 0.46 0.54 |
|                                                   | M3.5         |                                      | -1            |                        |
|                                                   | M4.5         |                                      |               |                        |
|                                                   | M4.5         | -16.7                                |               |                        |
|                                                   | M2.5         |                                      |               |                        |
|                                                   | M4           |                                      |               | 0.06                   |
|                                                   | M3.5         |                                      | 0             | 0.21                   |
|                                                   | M3.5         | -6.6                                 | 1             | 0.79                   |

a The flags -1, 0 and 1 indicate respectivelly non detection, uncertain detection and certain detection.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-b1a8d34b-9255-49e1-aa24-a9a42732eed8.jpg)

We consider this average ¯ A V as the mean excess extinction added by the SEM. Regarding the issue of limited spatial resolution of the SEM, from Figure 1 it is clear that the background towards the 25 Orionis group is quite smooth, showing no significant spatial variations ( σ A V ∼ 0 2). This . uniformity is also confirmed in our Spitzer-MIPS 24 µ m images of this area (Hern´ndez et al. 2007) and by preliminary a extinction maps performed on the basis of near IR photometry from VISTA (Alves, private communication).

Figure 6. Distribution of the visual extinction A V throughout the 25 Orionis group. Left panel: A V values measured from spectra as explained in the text, from Brice˜o et al. (2005, n 2007) (black), from this work (gray), and the complete sample (white). Right panel: Distribution of the residuals between visual extinctions computed from spectroscopy, and those from the Schlegel, Finkbeiner &amp; Davis (1998) maps, for the same complete sample of the left panel. Solid line indicates the normal distribution which fits the data, with ¯ A V = 0 16 and . σ = 0 39 magni-. tudes after removing the outliers.

In order to assess how much we would overestimate the extinction by using the SEM, we computed a mean excess in A V by comparing the distribution of the residuals between A V computed using spectroscopy for the members confirmed in this work and those from Brice˜o et al. (2005), n and A V from the SEM as we show in Figure 6. Except for a few outliers, the resulting distribution is essentially normal, with a mean ¯ A V = 0 16 and standard deviation . σ A V = 0 39. .

Thus, we estimated the A V for each photometric candidate by interpolating in the SEM and subtracting the ¯ A V estimated from the distribution of residuals. We consider the mean uncertainty of our estimation over the entire sample of candidates as σ A V / √ N = 0 03 magnitudes, where . N = 134 is the number of confirmed members used for the residual distribution.

## 4.2 Contamination of the photometric candidate sample by field stars and extragalactic sources

There are three main possible sources of contamination of a LMS and BD photometric candidate sample from a star forming region: field dwarf and giant stars, unresolved extragalactic sources and, in our case, young stars belonging to other sub-regions such as Orion OB1a. Such contamination could affect estimates of the IMF, the mean age and the spatial distribution of the 25 Orionis group when performed using only photometric candidates. We use a Galactic model to estimate the contamination by Orion non-members to-

wards our field, and since we obtained a large number of spectra, we can verify the predictions from such models.

The number of field stars, of all kinds, in the photometric candidate sample, was estimated from a simulation of the expected galactic stellar population in an area of ∼ 6 4 . deg 2 centred at the 25 Orionis group, performed with the Besan¸ con Galactic model (Robin et al. 2003). This model has been extensively tested and correctly reproduces the luminosity function within the sensitivity limits of the 2MASS survey (Robin et al. 2003). Second, we reproduce in the Besan¸ con data the typical photometric incompleteness from our survey according to the following procedure: First, we ran the model simulation in the range 13 &lt; &lt; I 22 and simulating the photometric errors as function of magnitude, using an exponential function with parameters obtained from fitting our survey data. Second, to estimate the fractional contamination in our photometric survey we need to remove objects from the Besan¸ con model in order to reproduce the fall-off as a function of magnitude observed in our survey. To do this we produced I vs. I-J Hess diagrams for our survey and the simulated sample using 0.1 × 0.1 magnitude bins and randomly discarded simulated objects in each bin until their number matched our observed number. After this correction was applied, a star from the Besan¸on simulation c was considered a contaminant of the photometric candidate sample if it satisfied the selection criteria explained in Section 3.3.

The resulting sample of contaminants is constituted exclusively by 212 low-mass main sequence stars that have an homogeneous spatial distribution with a mean spatial density of ∼ 33 objects per deg 2 , estimated spectral types from K7 to L0 and distances between ∼ 50 pc and ∼ 300 pc. In our spectroscopic sample, every single star rejected as a member turned out to be a field dwarf, several with some degree of emission in H α as seen in dMe stars, but all showing the deep Na I 8183,8195 absorption expected in older, low-mass disc stars. No late type giants were found among the rejected, non-member, spectroscopic sample as suggested by our results from Besan¸ con. In order to estimate the contamination as a function of member's mass, we assume such contaminant field dwarf stars mistaken for members of the 25 Orionis group, and we computed the masses and ages they would have, following the procedure we will explain in Section 5.1, resulting in the distributions we show in Figures 7 and 8.

We compared the observational J vs J-Ks diagram with those obtained from an appropiate simulation with the Besan¸ con Galactic model (Robin et al. 2003) of the expected galactic population in the area covered by the survey (see Section 4.2). We found that the point sources showing JKs¿1.5 are spatially unresolved extragalactic objects that are not predicted by the Besan¸ con galactic model. Then, we applyed to the selected extragalactic sources our criteria for the selection of photometric candidates to members of the 25 Ori group and found that only a negligible number of unresolved extragalactic sources ( ∼ 1 source per deg 2 ) do satisfy our criteria for photometric candidate selection; this means a total of ∼ 7 extragalactic sources is expected in the entire sample of photometric candidates.

Summarizing, we found 1246 member candidates in the 25 Orionis group in the mass range considered here. If we neglect the contamination due to extragalactic sources and

©

c

members from other regions of Orion OB1, from the Besan¸ con model we expect 316 field dwarfs within the sample for a mean density of ∼ 33 contaminants per deg 2 . The total number of candidates expected to be real members of the 25 Orionis group in the mass range considered here is 930 corresponding to a mean spatial density of ∼ 146 members per deg 2 . Thus, our photometric selection has a global efficiency of ∼ 82 per cent which is in good agreement with the 77 new confirmed members which corresponds to the ∼ 86 per cent of our spectroscopic survey because spectra were obtained in the overdensity of candidates explaines in Section 5.2. We emphasize that the contamination shows a dependence with photometric colour which results in the fiducial mass and age distributions shown in Figures 7 and 8, which would be the mass and age distributions resulting from wrongly considering such contaminants as members.

## 5 ANALYSIS AND RESULTS

## 5.1 The mean age and the initial mass function

The masses and ages we derived here for the 25 Orionis lowmass population were obtained through interpolation of the luminosity and effective temperatures in the Baraffe et al. (1998) models, adopting a distance of 360 pc from Brice˜o n et al. (2014, in preparation). For the photometric candidates we interpolated the luminosity and the effective temperature we obtained by interpolation of the I-band magnitude and the I-J photometric colour, both corrected by interstellar extinction as explained in Section 4.1, using the relationships from Luhman et al. (2003b) and Dahn et al. (2002), from which we obtained the effective temperatures and bolometric corrections. For the spectroscopically confirmed members we followed two procedures: First, as we did for the photometric candidates. Second, by interpolation of the measured spectral types in the same relationships. Within the uncertainties, no significant differences were found for masses and ages derived with either procedure for the spectroscopically confirmed member sample. This similarity is not surprising, given the low and smooth interstellar reddening towards this region, and the ∼ 0.5 spectral subclass mean uncertainty in the estimation of spectral types.

We derived the distribution of ages and the systemIMF for the photometric candidates and for the confirmed members in the entire area of our survey and in the spatial overdensity that will be defined in Section 5.2 as the 25 Orionis group overdensity. Figure 7 shows the resulting distribution of ages for photometric candidates and for spectroscopically confirmed members including also the members from Brice˜o et al. n (2005) within the same mass interval. The member sample is still largely incomplete and biased toward the brighter (more massive) objects, but as can be seen in Figure 7, this bias does not strongly affect the estimation of age. The distribution corresponding to photometric candidates was corrected by subtracting the age distribution for the expected contamination by field stars, i.e. the age distribution that results from mistaking the contaminants as members of the 25 Orionis group.

We found a good agreement between the distributions of ages computed from photometric candidates and from

Figure 7. Distribution of ages for the full sample of stellar and sub-stellar candidates and confirmed members in the complete area of the survey (left) and inside a circle with radius 0.5 ◦ centred at the 25 Orionis group overdensity defined in Section 5.2 (right). The distribution corresponding to photometric candidates (crosses) was corrected by subtracting the expected contamination by field stars (gray circles), resulting in the distribution indicated with black dots and the corresponding error bars. The triangles indicate the distribution for the spectroscopically confirmed members from this work and from Brice˜o et al. (2005). The distribution for members n is the same in both plots because the spectroscopic sample covers only the overdensity. Top labels and vertical lines indicate the mode and standard deviation of the distributions from confirmed members (left) and from photometric candidates (right).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-54c3f9a7-69de-49dd-b712-2261df5a3031.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-b5dd909e-cb50-45bc-99ed-810c8befdef4.jpg)

Figure 8. The system-IMF for the two samples considered in the Figure 7. The system-IMF from photometric candidates (crosses) was corrected by subtracting the expected contamination by field stars (gray circles), resulting in the system-IMF indicated with black dots and error bars. The empty circles indicate the system-IMF before the completeness correction was applied for objects with masses M &lt; 0 06 . M /circledot . The triangles indicate the system-IMF for the spectroscopically confirmed members, which is still largely incomplete. We compute all the distributions using bins of 0.2 dex. The dotted and solid curves indicate the best-fit to a log-normal distribution and a Kroupa power-law function.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-6ff75e6f-0d3f-4005-b8a3-0c785d6f2fc8.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-c97f1789-2868-4233-9142-8d80902f9d7f.jpg)

spectroscopically confirmed members. In the complete area of the survey, the derived age of the 25 Orionis group is 7 7 . ± 0 6 Myr from candidates and 6 1 . . ± 0 8 Myr from con-. firmed members. In the spatial overdensity, where the new members were detected, we found a mean age of 8 6 . ± 0 6 . Myr from candidates. All the ages reported correspond to the mode of each distribution. We adopt the age derived from the confirmed members (6 1 . ± 0 8 . Myr) as the age of the 25 Orionis group, although this value is statistically indistinguishable from that obtained from photometric candidates. Therefore we confirm the age derived in previous works by our group ( ∼ 7 Myr ; Brice˜o et al. 2007). In this n new sample we are including several new spectroscopically confirmed members and photometric candidates covering a much wider magnitude range which allows for a more robust estimation of the mean age. Additionally, our result shows that the population within the 25 Orionis group is not distinguishable from the surrounding population in terms of their ages alone.

Figure 8 shows the system-IMF in the LMS and BD regime for the photometric candidates of the 25 Orionis group. Although the photometric completeness of our survey (I=19.6) is equivalent to a completeness mass limit of 0.03 M /circledot assuming the Baraffe et al. (1998) models with an age of 7 Myr old and a distance modulus of m-M=7.78, the width of the membership loci defined in Section 3.3 and the dispersion of the candidates within it, suggest that some objects more massive than 0.03 M /circledot could be located below the photometric completeness. In order to account for such incompleteness in our sample we estimated, per magnitude bin, the number fraction of missing objects as the ratio between the expected and observed number of objects fainter than the completeness magnitude 10 . Figure 8 shows the system-IMF obtained after this completeness correction and the correction for the contamination by field stars were applied.

Because we did not apply any correction to account for unresolved binaries, what we report is the so called system-IMF of the 25 Orionis group (see e.g. Parravano, McKee &amp; Hollenbach 2010, and references therein). Although the system-IMF does not reveal by itself the real number of objects per mass bin, it is useful for comparison with system-IMFs from other SFR and young clusters, assuming that the samples have similar spatial resolution, and that the populations have similar multiplicity properties such as the number fraction of objects belonging to multiple systems, separation and mass ratio distributions.

We describe the derived system-IMF according to a Kroupa power-law function Ψ( m ) = dN/dm ∝ m -α i (Kroupa 2001, 2002), and a Scalo log-normal function (Scalo 1986) in the form Ψ(log m ) = dN/d (log m ) ∝ exp -((log m -log m c ) 2 / σ 2 2 ) (Chabrier et al. 2005). We fitted the mass distributions within the mass range 0 03 . &lt; M/M /circledot &lt; 0 8 . to the log-normal function and within the

10 The observed number counts (for the full catalogue) behave linearly on a log scale up to the completeness magnitude (by definition), after which these start to decline. We extrapolate this linear behaviour for magnitudes fainter than the completeness, and use this fit to compute the expected number counts.

Figure 9. Comparison of the system-IMF obtained for the 25 Orionis group with those from Chabrier (2003b) and Kroupa (2001, 2002). The black solid curve and straight lines correspond to the fit to a log-normal and power-law functions for the entire area of the survey respectively. The gridded area indicates, in the BD regime, the uncertainty in the fit to the log-normal function. The dotted curve and straight lines: fit to a log-normal and power-law functions for the spatial overdensity. The dashed area indicates, in the BD regime, the uncertainty in the fit to the lognormal function. Gray curve: IMF from Chabrier (2003b). Gray solid lines: IMF from Kroupa (2001, 2002).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-68ac9cbc-6b80-430d-a5e5-9c43016f37dc.jpg)

mass ranges 0 03 . &lt; M/M /circledot &lt; . 0 08 and 0 08 . &lt; M/M /circledot &lt; . 0 5 for the power-law description. As we did for the age estimate, we obtain the system-IMF for the candidates inside the spatial overdensity (see Section 5.2) and for the entire area of the survey. For the latter, the system-IMF for masses in the intervals 0 03 . &lt; M/M /circledot &lt; 0 08 . and 0 08 . &lt; M/M /circledot &lt; 0 5 . follows Kroupa-type functions with α 3 = -1 73 . ± 0 31 and . α 2 = 0 68 . ± 0 41 respectively, and . that in the complete mass interval 0 03 . &lt; M/M /circledot &lt; 0 8 . is well described by a log-normal function with coefficients m c = 0 21 . +0 02 . -0 02 . and σ = 0 36 . ± 0 02. . For the spatial overdensity we found α 3 = -1 97 . ± 0 02, . α 2 = 0 37 . ± 0 04, . m c = 0 22 . +0 02 . -0 02 . and σ = 0 42 . ± 0 05. In Figure 8 and Table . 6 we show the resulting fits and in Figure 9 we compare our derived IMF with other IMFs from the literature in which no correction for binarity was made. Additionally, we computed the distribution of spectral types, which is generally used as a proxy of the IMF which avoids most of the uncertainties related to the computation of masses (e.g. Preibisch 2012). Figure 10 shows the spectral type distributions for confirmed members and photometric candidates. For the later, the spectral types were estimated from the observed I -J colour, dereddened following the procedure explained in Section 4.1, and applying the colour-spectral-type relationship from Luhman et al. (2003b). The distribution of spectral types for candidates peaks around ∼ M . 4 5 and declines for earlier and latter objects, and is consistent with the spectral type distributions reported for confirmed mem-

Table 6. Fit parameters for the system-IMF.

| Population                                          | M interval                                                                       | α 2                                   | α 3                                       | m c                                                    | σ                                      | Reference                            |
|-----------------------------------------------------|----------------------------------------------------------------------------------|---------------------------------------|-------------------------------------------|--------------------------------------------------------|----------------------------------------|--------------------------------------|
| 25 Orionis overdensity 25 Orionis entire area Field | 0 . 03 < M/M /circledot < 0 . 8 0 . 03 < M/M /circledot < 0 . 8 M < M /circledot | 0 . 37 ± 0 . 04 0 . 68 ± 0 . 41 . . . | - 1 . 97 ± 0 . 02 - 1 . 73 ± 0 . 31 . . . | 0 . 22 +0 . 02 - 0 . 02 0 . 21 +0 . 02 - 0 . 02 0 . 22 | 0 . 42 ± 0 . 05 0 . 36 ± 0 . 02 0 . 57 | This work This work Chabrier (2003b) |

Figure 10. The distribution of spectral types in the 25 Orionis group. Triangles: confirmed members from this work and from Brice˜o et al. (2007). Black and empty dots indicate respectively n the candidates in the entire area of the survey and inside the overdensity previously corrected by the expected contamination by field stars.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-b40c7a74-9f72-4abf-812e-afc954b97faa.jpg)

bers in other regions such as IC348 and Chamaeleon I (e.g Luhman 2012).

The system-IMF, and the spectral type distribution are commonly used as observables to constrain which processes and conditions are preponderant during the formation of LMS and BDs. In the sub-stellar domain some variations have been observed in these quantities, Taurus probably being the most representative case (Brice˜o et al. n 2002). These variations are not well understood and more observations are needed in order to better constraint these phenomenon. What we have presented here is a new photometry-based system-IMF in the mass range 0 03 . &lt; M/M /circledot &lt; 0 8, . obtained for a ∼ 7 Myr old dispersed population in a low-density environment, using a numerous sample with 1246 objects. The peak of the spectral-types distribution and the characteristic mass m c of the log-normal representation of the system-IMF we found for the 25 Orionis group, is consistent with those reported for other SFR and young clusters (e.g. Oliveira, Jeffries &amp; van Loon 2009), which, as mentioned by Oliveira, Jeffries &amp; van Loon (2009), supports the idea of an IMF relatively independent of environmental properties such as density, temperature, metallicity and radiation field, as predicted by the models from Bonnell, Clarke &amp; Bate

(2006) and Elmegreen, Klessen &amp; Wilson (2008). On the other hand, the standard deviation σ from the log-normal fit results to be slightly smaller than is expected from Chabrier (2003b) which is also consistent with the small values for the α 3 coefficients obtained from the fit to a power-law IMF. This suggest that the 25 Orionis population could have a lower number of BDs. Because the fits are strongly dependent on small variations in the number of objects in the least massive bins, more observations are needed in order to better constraint the system-IMF in the sub-stellar domain of the 25 Orionis group and confirm such a trend. A sensitive photometric survey for the detection of a complete sample of BDs down to ∼ 0 01 . M /circledot , and the spectroscopic confirmation of members in the entire mass range (including also M &gt; 0 8 . M /circledot ) of this group are needed to confirm this low number of BD. Also, kinematic information is needed in order to evaluate whether the dynamical evolution of the group could be biasing the result. High spatial resolution imaging intended to detect binary or multiple systems are important in order to establish the effect that such a population of binaries could have on the observed system-IMF. Spectroscopic observations of our least massive candidates and much deeper photometric observations are currently under way.

## 5.2 The spatial distributions

Here we present the first comparative study of the spatial distribution of ∼ 1100 candidate LMS and ∼ 120 candidate BDs within the area covered by the 25 Orionis group and its surroundings. With this deep and spatially complete photometric sample we confirm the 25 Orionis group as the strong overdensity of stellar and sub-stellar candidate and members shown in Figure 11. We compute the number of objects on 3 ′ × 3 ′ bins and smoothed the resulting pattern with a gaussian function with standard deviation of 6 ′ in both spatial directions. The overdensity is centred at ( α, δ ) ∼ (81 2 1 7) . ◦ , . ◦ with a mean radius of ∼ 0 5 and . ◦ a slight elongation approximately in the East-West direction. Its mean density is of ∼ 0.5 objects per 3 ′ × 3 ′ bin ( ∼ 200 obj/deg 2 ), surrounded by a typical mean density of ∼ 0.25 objects per 3 ′ × 3 ′ square ( ∼ 100 obj/deg 2 ). The centre was computed as the position of the maximum density within the overdensity and the mean radius as the distance from the centre at which the density drops to the value of the mean density of the surroundings. Assuming a distance of ∼ 360 pc the angular radius corresponds to a mean linear radius of ∼ 3.14 pc. We note that the 25 Orionis overdensity appears to be superimposed on a low-density diagonal filament indicated in Figure 11, whose mean density is ∼ 0.3 objects per 3 ′ × 3 ′ square.

The overdensities could be intrinsic to the population, a consequence of observational biases or a combination of

both. We consider the detected overdensities as real because: (i) They are not a consequence of interstellar extinction, since the latter does not significantly affect the number of selected candidates and the SEM does not shows a spatial distribution similar with those from candidates, as can be seen in Figure 2. (ii) They cannot be explained as a consequence of contamination by field stars because these follow a smooth spatial distribution which depends on galactic latitude as predicted by the Besan¸ con model. (iii) They are not a consequence of a photometric bias because the completeness and limiting magnitudes in all bands used for the candidate selection are spatially homogeneous in the area where the overdensities were detected and the slight difference in the completeness magnitude within observations centred at δ = 1 ◦ and δ = 3 ◦ (Table 2) cannot produce the detected overdensities. (iv) They are not a consequence of any artifact in the images or inappropriate identification of sources from different catalogues as we concluded from visual inspection of the images.

The East-West elongation of the 25 Orionis overdensity is suggested from the spatial distribution of the members confirmed by Brice˜o et al. n (2005) and also by the spatial distribution of the high mass candidates of Kharchenko et al. (2005) with the highest probability of membership. In this work we confirm this trend with a numerous sample of well characterised candidates covering not only the 25 Orionis group but a reasonable area of its surroundings. Furthermore, our results rule out the idea of a possible extension of the 25 Orionis group to the South as suggested by McGehee (2006). With this new data set, deeper and covering the entire area surrounding the group, we find that the possible southern extension of the group reported by McGehee (2006) is probably part of the filament described above that extends also to the North-East direction, or the cluster ASCC18 from Kharchenko et al. (2005). This conclusion could not have been derived from the data considered by McGehee (2006), which covers only the southern part of the area shown in Figure 1. We note that the southern overdensity detected by McGehee (2006) at ( α, δ ) ∼ (81 2 0 8) is . ◦ , . ◦ also observed in Figures 11 and 12 as part of the southern extension of the filament (note the different scales between our plots and Figure 8 from McGehee (2006)), but its mean low-density and position do not suggest it is part of the overdensity we consider here as the 25 Orionis group.

Results from van den Bergh (2006) show that during their first ∼ 15 Myr, ∼ 20 per cent of galactic clusters and associations have diameters larger than 7 pc, being expanding associations rather than gravitationally bound clusters. Therefore the diameter obtained in this work for the 25 Orionis group (D ∼ 6.3 pc) suggests it could be a bound cluster given its age of about ∼ 7 Myr. Nevertheless, this computation is based on the radius found for the spatial overdensity, but we cannot rule out the cluster is spread over a larger area, since, as noted in Section 5.1, the age distribution and system-IMF are indistinguishable from its surroundings. A formal computation of the virial mass of this cluster is out of the scope of this work, since it depends almost entirely on its massive stars, which are not discussed here.

The panels of Figure 12 show separately the spatial density distribution of ∼ 1100 LMS and ∼ 120 BDs. In order to test whether the observed spatial distributions for

Figure 11. Spatial density distribution of the LMS and BD photometric candidates of the 25 Orionis group and its surroundings. Colored scales and isocontours indicate the mean density computed in bins of 3 ′ × 3 ′ and smoothed with a gaussian with standard deviation of 6 ′ in both directions. The central strong overdensity is the 25 Orionis group and the doted line indicates the filamentary overdensity explained in the text. The circles indicate the radii from the previous works shown also in Figure 1. The contamination from field stars present in the candidate sample is spatially smooth. The central big dot indicates the centre of the group computed in this work, and the big asterisk, the position of the star 25 Ori.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-dccf56ca-bf42-4c8a-8a5b-d7f29a072ac9.jpg)

BDs and LMS are drawn from the same parent distribution (null hypothesis), we have computed a two-dimensional twosample Kolmogorov-Smirnov (K-S) test (Press et al. 2007). For the full BD and LMS samples the K-S statistic results in a probability P = 0 30 of the null hypothesis being true, . which means the difference between both spatial distributions is not statistically significant (Press et al. 2007). In order to assess the robustness of this result, we generated 10 3 bootstrap resamples of both the BD and LMS distributions. A single bootstrap sample is made, e.g. for BDs, by randomly sampling with replacement from the BD catalogue, until a new set with the same number of objects (as the original catalogue) is generated. This is a standard technique that allows producing a set of random realizations drawn from the same distribution as the parent populations, which in this case are the BD and VLM samples (see e.g. Wall &amp; Jenkins 2012; Press et al. 2007, for more details). We computed the same K-S test on the 10 3 bootstrap BD and LMS samples and found that only in 12 per cent of the experiments the distributions turned out to be incompatible at the 99 per cent confidence level ( P &lt; 0 01). . This shows it is a robust result, independent on stochastic fluctuations due to sample size, that the spatial distributions observed for BDs and LMS in this 7 Myr old population are statistically indistinguishable.

Although this result is not consistent with the mass segregation predicted by the early premature ejection models from Reipurth &amp; Clarke (2001) subsequent work by

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-aa643891-67dc-4145-8f93-6cf0497d2b56.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-ea5635a7-0735-4c44-8579-966beec5223c.jpg)

)

°

(

δ

Figure 12. Comparison of the density distributions of LMS candidates ( M &gt; 0 08 . M /circledot , white dots, left) and of BD candidates ( M &lt; 0 08 . M /circledot , white dots, right). The K-S test probability between both distributions is 0.30 which indicates that both distributions are statistically indistinguishable. The big central asterisk represents the 25 Ori star and the black dots clusters from Kharchenko et al. (2005). Other elements are as in Figure 11.

Bate, Bonnell &amp; Bromm (2003) does not predict significantly different velocity dispersions for ejected BDs and stars. Thus, the similar spatial distributions do not necessarily reject the premature ejection scenarios. In any case, what we are presenting here constitutes a new observational evidence supporting the similarity of the spatial distribution of LMS and BDs for a ∼ 7 Myr old population.

## 5.3 Disc indicators and accretion signatures

We computed the number fraction of LMS and BDs showing near IR excesses as an indication of the presence of primordial circumstellar inner-discs. In the left panel of Figure 13 we show a J-H vs. H-Ks colour-colour diagram with the CTTS locus from Meyer, Calvet &amp; Hillenbrand (1997) and its extrapolation to M6 spectral type. We obtained that 3 8 . ± 0 4 . per cent of LMS candidates and 4 5 . ± 0 5 . per cent of spectroscopically confirmed LMS are placed inside the CTTS locus. Both results are in good agreement within the uncertainties and are consistent with those from Brice˜o et al. (2007) ( n ∼ 6 per cent). The fractions reported here for the photometric candidates were computed considering the expected contamination by field stars and the errors were computed according to the binomial distribution.

M9 are more separated, as we show in the right panel of Figure 13. Essentially we define a locus whose lower limit in the I-Ks vs J-H diagram is parallel to the photospheric locus from Luhman et al. (2003b), and which is separated from it by a distance equivalent to a reddening vector corresponding to A V = 0 5, which is slightly larger . than the mean value of the interstellar extinction across the 25 Orionis group, as we have shown in Section 4.1. We found that ∼ 5 6 . ± 0 5 per cent of the candidate LMS and . ∼ 4 1 . ± 0 4 per . cent of the confirmed LMS show excesses. For the substellar regime the number fraction of candidates showing such IR excess is 12 7 . ± 0 5 per cent. Again, the number fraction of . stars showing infrared excesses from inner discs agrees well with previous determinations by Hern´ andez et al. (2006) and Brice˜o et al. n (2007) using Spitzer data. In the substellar regime, our result is the first estimate obtained for the 25 Orionis group.

For objects later than M6, the J-H vs. H-Ks diagram is not a good choice for detecting infrared excesses because the photospheric locus between spectral types M6 and M9 is essentially colinear with the corresponding reddenning vectors, which prevents a reliable measurement of the fraction of objects showing IR excesses as a function of their spectral-types. Because of this limitation, we apply an empirical criterion to detect objects with infrared excesses, using the I-Ks vs J-H diagram, where the photospheric loci and reddening vectors for spectral types between M6 and

Emissions from warm dust in evolved discs or in the external regions of primordial discs, occurs at longer wavelengths. In Figures 14 and 15 we show colour-colour diagrams in the wavelength range from ∼ 2 2 . µ m to ∼ 12 µ m, for confirmed members and photometric candidates with magnitudes from VISTA, IRAC and WISE. The left panel of Figure 14 shows the 3.6 µ m-4.5 µ m vs. 4.5 µ m-5.8 µ m diagram from IRAC data with the excess region defined by Luhman et al. (2005) were 7 8 . ± 0 3 . per cent of the candidates to LMS, 7 5 . ± 0 4 . per cent of the LMS confirmed as members and 18 1 . ± 0 3 . per cent of the candidates to BD show IR excesses. The right panel of Figure 14 shows the IRAC 3.6 µ m-4.5 µ m vs. 5.8 µ m-8.0 µ m diagram showing the CTTS locus defined by Hartmann et al. (2005) and Luhman et al. (2005), were 12 0 . ± 0 4 per cent of the LMS . candidates, 9 4 . ± 0 8 per cent of the LMS confirmed as mem-. bers and 39 2 . ± 0 4 per cent of the BD candidates show IR . excesses.

Figure 13. Optical and near-infrared colour-colour diagrams of spectroscopically confirmed members and photometric candidates. Left panel: J-H vs. H-Ks. Right panel: I-Ks vs. J-H. The dashed lines indicates the photospheric loci from Luhman et al. (2003b) and the shaded areas indicate the excess loci explained in Section 5.3. Empty and filled circles indicate respectively candidates and confirmed LMS. Triangles indicates candidate BDs. Reddening vectors for several spectral types are also indicated for the typical mean visual extinction ( A V ∼ 0 5 mag) through the 25 Orionis group. Number fractions of objects showing IR-excesses are listed in Table 7. .

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-eb799403-b122-4ac9-ab77-fad6f8c07915.jpg)

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-a77239e9-c39d-4238-9f06-300ea94ab1a6.jpg)

Table 7. Number fraction of LMS and BDs showing infrared excesses.

| Diagram                           | LMS members [per cent]   | LMS candidates [per cent]   | BD candidates [per cent]   |
|-----------------------------------|--------------------------|-----------------------------|----------------------------|
| J-H vs. H-Ks                      | 4 . 5 ± 0 . 5            | 3 . 8 ± 0 . 4               | · · ·                      |
| I-Ks vs. J-H                      | 4 . 1 ± 0 . 4            | 5 . 6 ± 0 . 5               | 12 . 7 ± 0 . 5             |
| 3.6 µm -4.5 µm vs. 4.5 µm -5.8 µm | 7 . 5 ± 0 . 4            | 7 . 8 ± 0 . 3               | 18 . 1 ± 0 . 3             |
| 3.6 µm -4.5 µm vs. 5.8 µm -8.0 µm | 12 . 0 ± 0 . 4           | 9 . 4 ± 0 . 8               | 39 . 2 ± 0 . 4             |
| Ks-8.0 µm vs. Ks-5.8 µm           | 5 . 7 ± 0 . 4            | 7 . 9 ± 0 . 5               | 47 . 0 ± 0 . 5             |

Using the Ks-8 µ m vs. Ks-24 µ m and the Ks-5.8 µ m vs. Ks-24 µ m diagrams based on IRAC and MIPS photometry, Luhman et al. (2010) show the distinction between circumstellar discs at different evolutionary stages. We do not have MIPS observations for the candidates and members considered in this work but we used the Ks-8 µ m and Ks-5.8 µ m colour intervals defined by Luhman et al. (2010) in order to separate objects showing photospheric emissions from those showing primordial discs and those showing transitional or evolved discs. Following this procedure we obtain an estimation of the number fraction of such populations. The left panel of Figure 15 shows the Ks-5.8 µ m vs. Ks-8.0 µ m diagram, where we consider objects with excesses to be those with Ks-5.8 µ m &gt; 1 and Ks-8.0 µ m &gt; 1, resulting in a number fraction of 7 9 . ± 0 5 . per cent of the LMS candidates, 5 7 . ± 0 4 per cent for confirmed LMS and . ∼ 47 0 . ± 0 5 per . cent for the photometric candidates to BD. In the transitional and evolved disc locus, defined as the region between 0 75 . &lt; Ks -8 0 . µm &lt; 1 4 and 0 3 . . &lt; Ks -5 8 . µm &lt; 1, we find 4 LMS members, 5 LMS candidates and 2 BD candidates for number fractions of 7 4 . ± 0 4 per cent, 2 3 . . ± 0 5 per cent . and 15 6 . ± 0 5 per cent respectively. The objects showing IR . excesses indicative of possible transitional or evolved discs are listed in Table 8. In the right panel of Figure 15 we also show the excesses at longer wavelengths, using a Ks-4.6 µ m vs. Ks-12 µ m diagram from WISE data.

The samples of new confirmed members and candidates are not affected by the photometric completeness of the IRAC observations and the reported fractions are not biased towards the objects showing excesses. The number fraction of BD showing IR excesses turns out to be systematically higher in all colour-colour diagrams suggesting that at an age of ∼ 7 Myr the disc fraction among BDs is larger than for LMS which could be an indicator of a slow evolution of discs around BDs. This supports previous results from Luhman &amp; Mamajek (2012) for the Upper Scorpius association which with a estimated age of ∼ 11Myr (Pecaut, Mamajek &amp; Bubar 2012) shows ∼ 25 per cent of objects later than M5 harbouring discs. Our results at a slight younger age support the scenario in which discs around LMS and BD evolve in different time scales in the sense that a higher fraction of BDs could retain their discs for longer periods of time.

Finally, low-resolution spectra provide an appropriate means to look for CTTS emission signatures such as H , the α CaII triplet ( λ 8498, λ 8542 and λ 8662) and

Figure 14. IRAC colour-colour diagrams of spectroscopically confirmed members and photometric candidates. Symbols are as in Figure 13. Left panel: the 3.6 µ m-4.5 µ m vs. 5.8 µ m-8.0 µ m diagram showing the CTTS locus (dashed lines) defined by Hartmann (2005). Right panel: the 3.6 µ m-4.5 µ m vs. 4.5 µ m-5.8 µ m diagrams with the excess region (dashed lines) defined by Luhman et al. (2005). Number fractions of objects showing IR-excesses are listed in Table 7

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-9ededeae-5252-41cd-b7db-bf4d14a83019.jpg)

Figure 15. Color-colour diagram combining photometry from VISTA, IRAC and WISE. Symbols are as in Figure 13. Dashed polygons indicate disc loci and the solid polygon in the left panel encloses the possible transitional or evolved discs from the colour intervals defined by Luhman et al. (2010).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-44fd65cf-3858-4f80-ad47-5ecde8614c78.jpg)

HeI λ 6676. The dividing line between CTTS and weak (non-accreting) T-Tauri stars (WTTS) was first set in terms of the H α emission at W(H ) α = -10 ˚ A by Herbig &amp; Bell (1988). Afterwards, White &amp; Basri (2003) revised this classification and, based on high resolution optical spectra ( R ∼ 33000), redefined the classification considering both the H α equivalent width and the spectral type. Barrado y Navascu´ es &amp; Mart´ ın (2003) using low-resolution optical spectra developed an empirical classification method based on the H α equivalent width and the spectral type. Here we followed the criteria from Barrado y Navascu´ es &amp; Mart´ ın (2003) to classify the confirmed LMS as CTTS or WTTS. Figure 16 shows the resulting classification in the W(H α ) vs. spectral-type plot. We obtained a number fraction of CTTS to WTTS of 3 8 . ± 0 5 . per cent for objects earlier than M6. Our results are consistent with those from the colour-colour diagrams and with previous determinations by Brice˜o et al. n (2005). On the

Table 8. Members and candidates showing possible transitional or evolved discs.

|   α ( J 2000) [ ◦ ] |   δ ( J 2000) [ ◦ ] |     J |     I |   CH1 |   CH2 |   CH3 |   CH4 |    W1 |    W2 |    W3 |   W4 | ST    | WHa [ ˚] A   | Av Av   | Status Status   |
|---------------------|---------------------|-------|-------|-------|-------|-------|-------|-------|-------|-------|------|-------|--------------|---------|-----------------|
|             80.7118 |             1.51037 | 15.02 | 16.89 | 13.68 | 13.54 | 13.29 | 12.65 | 13.81 | 13.45 | 11.81 | 8.93 | · · · | · · ·        | · · ·   | Candidate       |
|             80.8809 |             2.22021 | 13.03 | 14.92 | 11.53 | 11.24 | 11.02 | 10.83 | 11.65 | 11.21 | 10.42 | 8.38 | · · · | · · ·        | · · ·   | Candidate       |
|             80.883  |             2.17874 | 13.88 | 15.46 | 12.6  | 12.49 | 12.31 | 11.87 | 12.73 | 12.43 | 11.13 | 9.12 | · · · | · · ·        | · · ·   | Candidate       |
|             80.9344 |             1.69221 | 16.26 | 18.54 | 15    | 14.77 | 14.91 | 14.22 | 15.1  | 14.73 | 11.96 | 8.83 | M1.0  | -3.5         | 6.1     | Member          |
|             81.0207 |             2.15769 | 14.48 | 16.47 | 13.1  | 12.87 | 12.59 | 12.06 | 13.33 | 12.85 | 11.54 | 8.31 | · · · | · · ·        | · · ·   | Candidate       |
|             81.1103 |             1.84858 | 14.81 | 16.55 | 13.59 | 13.5  | 13.39 | 12.82 | 13.76 | 13.57 | 11.72 | 8.76 | M4.0  | -2.1         | 1.1     | Member          |
|             81.2452 |             1.42179 | 14.08 | 15.81 | 12.92 | 12.74 | 12.51 | 11.94 | 13.09 | 12.72 | 11.03 | 8.95 | M5.0  | -7.4         | 0.0     | Member          |
|             81.3514 |             1.6854  | 13.78 | 15.61 | 12.56 | 12.44 | 12.32 | 11.88 | 12.77 | 12.47 | 10.99 | 8.72 | M5.5  | -15.1        | 0.0     | Member          |
|             81.4448 |             1.72511 | 13.15 | 14.54 | 11.97 | 11.9  | 11.72 | 11.08 | 12.1  | 11.84 |  9.91 | 6.74 | M3.5  | -10.1        | 0.0     | Member          |
|             81.4683 |             2.36819 | 13.16 | 14.52 | 11.27 | 10.97 | 10.82 | 10.58 | 11.38 | 10.96 |  9.88 | 8.1  | · · · | · · ·        | · · ·   | Candidate       |
|             81.5062 |             1.54283 | 14.71 | 16.42 | 13.31 | 13.12 | 12.91 | 12.45 | 13.64 | 13.31 | 11.47 | 9.16 | · · · | · · ·        | · · ·   | Candidate       |
|             81.6042 |             1.87741 | 15.28 | 17.08 | 13.99 | 13.91 | 13.71 | 13.5  | 14.18 | 13.91 | 12.18 | 8.47 | M2.5  | · · ·        | 2.8     | Member          |
|             81.7246 |             2.30153 | 15    | 16.67 | 13.76 | 13.55 | 13.36 | 12.6  | 13.86 | 13.52 | 12.31 | 8.62 | · · · | · · ·        | · · ·   | Candidate       |
|             81.8506 |             1.63384 | 17.09 | 20.32 | 15.18 | 15.04 | 15.34 | 14.32 | 15.66 | 15.44 | 12.41 | 9.14 | · · · | · · ·        | · · ·   | Candidate       |
|             82.3184 |             2.15333 | 13.76 | 15.46 | 12.4  | 12.16 | 11.84 | 11.32 | 12.52 | 12.12 | 10.92 | 8.91 | · · · | · · ·        | · · ·   | Candidate       |
|             82.4414 |             1.95883 | 14.02 | 15.73 | 12.72 | 12.53 | 12.24 | 11.69 | 12.86 | 12.46 | 11.04 | 8.71 | · · · | · · ·        | · · ·   | Candidate       |
|             82.5921 |             1.64195 | 15.15 | 17.32 | 13.52 | 13.31 | 13.05 | 12.92 | 13.71 | 13.3  | 12.3  | 8.59 | · · · | · · ·        | · · ·   | Candidate       |

Figure 16. Equivalent width of Hα versus spectral type and effective temperature (upper axis) for the new confirmed LMS members. Dotted lines indicate the limit between CTTS and WTTS according to Barrado y Navascu´ es &amp; Mart´ ın (2003). Circles indicates objects showing equivalent width of Hα consistent with CTTS. Arrows indicate objects with strong Hα emission whose WHα is labelled. We found 3.8 ± 0.5 per cent of the members of 25 Orionis group being CTTS (open circles).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0009v1-figure-00753133-c48e-47b5-81e5-37a11de1f5c6.jpg)

other hand the number fraction of CTTS presented here is larger than the fraction reported by McGehee (2006) ( ∼ 0 1 . per cent) that, as we explained in Section 5.2, most probably does not correspond to the 25 Orionis group.

## 6 SUMMARY AND CONCLUSIONS

We have presented the first extensive study of the lower mass population (0 02 . /lessorequalslant M/M /circledot /lessorequalslant 0 8) of the 25 Ori group . and its surrounding regions. We provide for this mass range the spatial distribution, mean age, IMF, and investigate the presence of circumstellar discs and accretion signatures. In the following we summarize our results:

- · We have obtained a sample of 1246 photometric candidates to LMS and BDs in the 25 Ori group and its surround-

ings with estimated masses within 0 02 . &lt; M/M ∼ /circledot &lt; . ∼ 0 8, as well as 77 new LMS spectroscopically confirmed as members with spectral types in the range ∼ M1 to ∼ M5.

- · We have characterised the possible contamination of our photometric candidate sample using the Besan¸ con Galactic model (Robin et al. 2003). The expected mean contamination was found to be ∼ 33 field dwarf stars from the foreground per deg 2 , which yields an expected average efficiency of ∼ 82 per cent for our photometric candidate selection.
- · From the age distributions of the photometric candidates and confirmed members we find mean ages of 7 7 . ± 0 6 . Myr and 6 1 . ± 0 8 Myr respectively. These results are in good . agreement with those reported by Brice˜o et al. (2005). n
- · We have derived the system-IMF for 25 Ori from a numerous sample of photometric candidate LMS and BDs previously corrected by the contamination from field stars. We have obtained a system-IMF which can be well described by either a Kroupa power-law function with indices α 3 = -1 73 . ± 0 31 and . α 2 = 0 68 . ± 0 41 in the mass . ranges 0 03 . &lt; M/M /circledot &lt; 0 08 and 0 08 . . &lt; M/M /circledot &lt; 0 5 re-. spectively, or a Scalo log-normal function with coefficients m c = 0 21 . +0 02 . -0 02 . and σ = 0 36 . ± 0 02 . in the mass range 0 03 . &lt; M/M /circledot &lt; 0 8. . For the log-normal representation of the IMF we obtained that the characteristic mass is consistent with an IMF independent of environmental properties, although the standard deviation is smaller than the system-IMF from Chabrier (2003a). A sensitive photometric survey for the detection of a complete sample of BD down to ∼ 0 01 . M /circledot , and their spectroscopic counterparts is needed to confirm the slightly lower value of the standard deviation and whether this trend of the system-IMF is a consequence of a real deficit of BDs.
- · From the analysis of the spatial distribution of the photometric candidate sample we find the peak of the 25 Orionis overdensity at ( α, δ ) ∼ (81 2 1 7) with a diameter . ◦ , . ◦ of ∼ 6 3 pc. We have confirmed the East-West elongation . of the 25 Orionis group, also observed in the works from Kharchenko et al. (2005) and Brice˜o et al. (2005), and rule n out the southern extension proposed by McGehee (2006). Additionally we found that the 25 Orionis group is projected over a filament distributed along the south-west to north-east direction. We have also found that the spatial distributions of LMS and BDs in 25 Orionis are statistically indistinguishable.
- · For the LMS we calculated the number fraction of stars

showing IR excess in the J-H vs H-Ks diagram according to the CTTS locus from Meyer, Calvet &amp; Hillenbrand (1997). We found that ∼ 4 per cent of the LMS show the levels of IR excess expected from CTTS with optically thick primordial discs. We found similar fractions for the empirical locus defined in the I -Ks vs J -H colour-colour diagram where we find also that ∼ 12 per cent of the BD show IR excesses indicative of primordial discs.

- · Using VISTA, IRAC and WISE photometry we have also found IR excess at longer wavelengths indicative of discs. Depending on the selected diagram we found that the number fraction of LMS showing IR excesses is between ∼ 7 to ∼ 10 per cent while the number fraction for BDs increases up to ∼ 20 to ∼ 50 per cent.
- · Additionally we found that 11 members and candidates show signs of possible transitional or evolved discs and that 3 8 . ± 0 5 per cent of the LMS members show . Hα emission consistent with a CTTS nature, confirming the CTTS number fraction found by Brice˜o et al. (2007). n
- · Our results show that the fraction of BD showing IR excesses is higher than for LMS, supporting the scenario in which the evolution of circumstellar discs around the least massive stars could occur at larger time scales.

## ACKNOWLEDGMENTS

This work has been supported in part by grant S1200101144 and 200400829 from FONACIT, Venezuela. J. J. Downes, acknowledges support from CIDA, Venezuela and project 152160 from CONACyT, M´ exico. C. Mateu acknowledges the support of the postdoctoral Fellowship of DGAPA-UNAM, M´ exico. We thank Kevin Luhman for useful comments regarding the spectral classification and Susan Tokarz, who is in charge of the reduction and processing of Hectospec spectra.

This work was based on observations collected at the J¨rgen Stock 1-m Schmidt telescope of the National Obseru vatory of Llano del Hato Venezuela (NOV), which is operated by Centro de Investigaciones de Astronom´ ıa (CIDA) for the Ministerio del Poder Popular para Ciencia y Tecnolog´ ıa, Venezuela, and at the MMT Observatory a joint facility of the Smithsonian Institution and the University of Arizona. This work is based [in part] on observations made with the Spitzer Space Telescope, which is operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with NASA. Support for this work was provided by NASA through an award issued by JPL/Caltech. This publication makes use of data products from the Wide-field Infrared Survey Explorer, which is a joint project of the University of California, Los Angeles, and the Jet Propulsion Laboratory/California Institute of Technology, funded by the National Aeronautics and Space Administration. We very much appreciate the great work done by the UK-based VISTA consortium who built and commissioned the VISTA telescope and camera. This work is based on observations made during VISTA science verification, under the program ID 60.A-9285(B).

We thank the assistance of the personnel, observers, telescope operators and technical staff at CIDA and FLWO, who made possible the observations at the J¨rgen u Stock telescope of the Venezuela National Astronomical Observa- tory (OAN) and at the MMT telescope at Fred Lawrence Whipple Observatory (FLWO) of the Smithsonian Institution, especially Perry Berlind, Daniel Cardozo, Orlando Contreras, Franco Della Prugna, Emilio Falco, Freddy Moreno, Hern´ an Ram´ ırez, Carmen Rodr´ ıguez, Richard Rojas, Gregore Rojas, Gerardo S´nchez, Gustavo S´nchez and a a Ubaldo S´ anchez.

This work makes an extensive use of R from the R Development Core Team (2011) available at http://www.R-project.org/ and described in R: A language and environment for statistical computing from R Foundation for Statistical Computing, Vienna, Austria, ISBN 3900051-07-0.

We thank the referee Phil Lucas for useful comments and suggestions to the manuscript.

## REFERENCES

Abazajian K. et al., 2005, AJ, 129, 1755 Baltay C. et al., 2002, PASP, 114, 780 Baraffe I., Chabrier G., Allard F., Hauschildt P. H., 1998, A&amp;A, 337, 403 Barrado y Navascu´ es D., Mart´ ın E. L., 2003, AJ, 126, 2997 Basri G., 1998, in Astronomical Society of the Pacific Conference Series, Vol. 134, Brown Dwarfs and Extrasolar Planets, R. Rebolo, E. L. Martin, &amp; M. R. Zapatero Osorio, pp. 394-+ Bastian N., Covey K. R., Meyer M. R., 2010, ARA&amp;A, 48, 339 Bate M. R., Bonnell I. A., Bromm V., 2003, MNRAS, 339, 577 Biazzo K., Randich S., Palla F., Brice˜o C., 2011, A&amp;A, n 530, A19+ Bonnell I. A., Clarke C. J., Bate M. R., 2006, MNRAS, 368, 1296 Boudreault S., Bailer-Jones C. A. L., 2009, ApJ, 706, 1484 Bouvier J. et al., 2008, A&amp;A, 481, 661 Brice˜o C., Calvet N., Hern´ndez J., Vivas A. K., Hartn a mann L., Downes J. J., Berlind P., 2005, AJ, 129, 907 Brice˜o C., Hartmann L., Hern´ndez J., Calvet N., Vivas n a A. K., Furesz G., Szentgyorgyi A., 2007, ApJ, 661, 1119 Brice˜o C., Hartmann L., Stauffer J., Mart´ ın E., 1998, AJ, n 115, 2074 Brice˜o C., Luhman K. L., Hartmann L., Stauffer J. R., n Kirkpatrick J. D., 2002, ApJ, 580, 317 Brice˜o C. et al., 2001, Science, 291, 93 n Calvet N., Brice˜o C., Hern´ndez J., Hoyer S., Hartmann n a L., Sicilia-Aguilar A., Megeath S. T., D'Alessio P., 2005, AJ, 129, 935 Cardelli J. A., Clayton G. C., Mathis J. S., 1989, ApJ, 345, 245 Caswell J. L., 2007, in IAU Symposium, Vol. 242, IAU Symposium, J. M. Chapman &amp; W. A. Baan, ed., pp. 194199 Chabrier G., 2003a, PASP, 115, 763 Chabrier G., 2003b, ApJL, 586, L133 Chabrier G., Baraffe I., Allard F., Hauschildt P. H., 2005, ArXiv Astrophysics e-prints

Chen W. P., Lee H. T., 2008, arXiv.org, astro-ph, to appear in the Handbook of Star Forming Regions, Vol. I, edited byBo Reipurth

Dahm S. E., 2008, The Young Cluster and Star Forming Region NGC 2264, Reipurth, B., ed., pp. 966-+ Dahn C. C. et al., 2002, AJ, 124, 1170

- Davenport J. R. A., West A. A., Matthiesen C. K., Schmieding M., Kobelski A., 2006, PASP, 118, 1679
- Downes J. J., Brice˜o C., Hern´ndez J., Calvet N., Hartn a mann L., Ponsot Balaguer E., 2008, AJ, 136, 51
- Dullemond C. P., Dominik C., 2004, A&amp;A, 421, 1075
- Elmegreen B. G., Klessen R. S., Wilson C. D., 2008, ApJ, 681, 365

Emerson J. P., Sutherland W. J., 2010, in Society of PhotoOptical Instrumentation Engineers (SPIE) Conference Series, Vol. 7733, Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series

- Emerson J. P., Sutherland W. J., McPherson A. M., Craig S. C., Dalton G. B., Ward A. K., 2004, The Messenger, 117, 27

Fabricant D. et al., 2005, PASP, 117, 1411

Hartmann L., 2001, AJ, 121, 1030

- Hartmann L., 2005, in Astronomical Society of the Pacific Conference Series, Vol. 341, Chondrites and the Protoplanetary Disk, A. N. Krot, E. R. D. Scott, &amp; B. Reipurth, ed., pp. 131-+
- Hartmann L., Megeath S. T., Allen L., Luhman K., Calvet N., D'Alessio P., Franco-Hernandez R., Fazio G., 2005, ApJ, 629, 881
- Herbig G. H., Bell K. R., 1988, Third Catalog of EmissionLine Stars of the Orion Population : 3 : 1988, Herbig, G. H. &amp; Bell, K. R., ed.
- Hern´ andez J., Brice˜o C., Calvet N., Hartmann L., Muzen rolle J., Quintero A., 2006, ApJ, 652, 472
- Hern´ andez J., Calvet N., Brice˜o C., Hartmann L., Berlind n P., 2004, AJ, 127, 1682
- Hern´ andez J. et al., 2007, ApJ, 671, 1784
- Hern´ andez J., Hartmann L., Calvet N., Jeffries R. D., Gutermuth R., Muzerolle J., Stauffer J., 2008, ApJ, 686, 1195
- Hodgkin S. T., Irwin M. J., Hewett P. C., Warren S. J., 2009, MNRAS, 394, 675
- Ingleby L., Calvet N., Hern´ndez J., Brice˜o C., Espaillat a n C., Miller J., Bergin E., Hartmann L., 2011, AJ, 141, 127 Jayawardhana R., Mohanty S., Basri G., 2003, ApJ, 592, 282
- Jeffries R. D., Oliveira J. M., Naylor T., Mayne N. J., Littlefair S. P., 2007, MNRAS, 376, 580

Joergens V., 2006, A&amp;A, 446, 1165

- Kharchenko N. V., Piskunov A. E., R¨ser S., Schilbach E.,
- o Scholz R., 2005, A&amp;A, 440, 403
- Kirkpatrick J. D. et al., 1999, ApJ, 519, 802
- Kraus A. L., White R. J., Hillenbrand L. A., 2006, ApJ, 649, 306
- Kroupa P., 2001, MNRAS, 322, 231
- Kroupa P., 2002, Science, 295, 82

Kumar M. S. N., Schmeja S., 2007, A&amp;A, 471, L33

- Kun M., Kiss Z. T., Balog Z., 2008, Star Forming Regions in Cepheus, Reipurth, B., ed., pp. 136-+
- Landolt A. U., 1992, AJ, 104, 340
- Lodieu N., Hambly N. C., Jameson R. F., Hodgkin S. T., Carraro G., Kendall T. R., 2007, MNRAS, 374, 372
- Luhman K. L., 2004, ApJ, 614, 398
- Luhman K. L., 2012, ARA&amp;A, 50, 65
- Luhman K. L., Allen P. R., Espaillat C., Hartmann L.,

Calvet N., 2010, ApJS, 186, 111

- Luhman K. L., Brice˜o C., Stauffer J. R., Hartmann L., n Barrado y Navascu´ es D., Caldwell N., 2003a, ApJ, 590, 348
- Luhman K. L., Joergens V., Lada C., Muzerolle J., Pascucci I., White R., 2007, Protostars and Planets V, 443 Luhman K. L. et al., 2005, ApJL, 631, L69
- Luhman K. L., Mamajek E. E., 2012, ApJ, 758, 31
- Luhman K. L., Mamajek E. E., Allen P. R., Muench A. A., Finkbeiner D. P., 2009, ApJ, 691, 1265

Luhman K. L., Stauffer J. R., Muench A. A., Rieke G. H., Lada E. A., Bouvier J., Lada C. J., 2003b, ApJ, 593, 1093 Mateu C., Vivas A. K., Downes J. J., Brice˜o C., Zinn R., n Cruz-Diaz G., 2012, MNRAS, 427, 3374

- McCarthy C., Zuckerman B., 2004, AJ, 127, 2871
- McGehee P. M., 2006, AJ, 131, 2959
- Meyer M. R., Calvet N., Hillenbrand L. A., 1997, AJ, 114, 288

Mink D., 2006, in Astronomical Society of the Pacific Conference Series, Vol. 351, Astronomical Data Analysis Software and Systems XV, C. Gabriel, C. Arviset, D. Ponz, &amp; S. Enrique, ed., p. 204

- Monet D. G., 1998, in Bulletin of the American Astronomical Society, Vol. 30, Bulletin of the American Astronomical Society, pp. 1427-+
- Muzerolle J., Brice˜o C., Calvet N., Hartmann L., Hillenn brand L., Gullbring E., 2000, ApJL, 545, L141
- Muzerolle J., Luhman K. L., Brice˜o C., Hartmann L., n Calvet N., 2005, ApJ, 625, 906
- Neuh¨ auser R. et al., 2011, Astronomische Nachrichten, 332, 547
- Oliveira J. M., Jeffries R. D., van Loon J. T., 2009, MNRAS, 392, 1034
- Padoan P., Nordlund ˚ ., 2004, ApJ, 617, 559 A

Parravano A., McKee C. F., Hollenbach D., 2010, 45

- Pecaut M. J., Mamajek E. E., Bubar E. J., 2012, ApJ, 746, 154

Petr-Gotzens M. et al., 2011, The Messenger, 145, 29 Preibisch T., 2012, Research in Astronomy and Astrophysics, 12, 1, provided by the SAO/NASA Astrophysics Data System

- Press W. H., Teukolsky S. A., Vetterling W. T., Flannery B. P., 2007, Numerical recipes in C, Cambridge University Press

Reipurth B., Clarke C., 2001, AJ, 122, 432

- Robin A. C., Reyl´ C., Derri`re S., Picaud S., 2003, A&amp;A, e e 409, 523

Scalo J. M., 1986, FCP, 11, 1

- Schlegel D. J., Finkbeiner D. P., Davis M., 1998, ApJ, 500, 525
- Scholz A., Eisl¨ffel J., 2005, A&amp;A, 429, 1007 o
- Scholz A., Xu X., Jayawardhana R., Wood K., Eisl¨ffel J., o Quinn C., 2009, MNRAS, 398, 873

Skrutskie M. F. et al., 2006, AJ, 131, 1163

- Strom K. M., Wilkin F. P., Strom S. E., Seaman R. L., 1989, AJ, 98, 1444

Taylor M. B., 2006, in Astronomical Society of the Pacific Conference Series, Vol. 351, Astronomical Data Analysis Software and Systems XV, C. Gabriel, C. Arviset, D. Ponz, &amp; S. Enrique, ed., p. 666

- Torres C. A. O., Quast G. R., Melo C. H. F., Sterzik M. F., 2008, Young Nearby Loose Associations, Reipurth, B.,

pp. 757-+ van den Bergh S., 2006, AJ, 131, 1559 van Eyken J. C. et al., 2011, AJ, 142, 60 Wall J. V., Jenkins C. R., 2012, Practical Statistics for Astronomers Weidenschilling S. J., 1997, ICARUS, 127, 290 White R. J., Basri G., 2003, ApJ, 582, 1109 Whitworth A., Bate M. R., Nordlund ˚ ., Reipurth B., ZinA necker H., 2007, Protostars and Planets V, 999, 459 Wright E. L. et al., 2010, AJ, 140, 1868

This paper has been typeset from a T E X/ LT E X file preA pared by the author.
Juan José Downes César Briceño Cecilia Mateu Jesús Hernández Anna Katherina Vivas Nuria Calvet Lee Hartmann Monika G. Petr-Gotzens Lori Allen
The low mass star and sub-stellar populations of the 25 Orionis group
We present the results of a survey of the low mass star and brown dwarf
population of the 25 Orionis group. Using optical photometry from the CIDA Deep
Survey of Orion, near IR photometry from the Visible and Infrared Survey
Telescope for Astronomy and low resolution spectroscopy obtained with Hectospec
at the MMT, we selected 1246 photometric candidates to low mass stars and brown
dwarfs with estimated masses within $0.02 \lesssim M/M_\odot \lesssim 0.8$ and
spectroscopically confirmed a sample of 77 low mass stars as new members of the
cluster with a mean age of $\sim$7 Myr. We have obtained a system initial mass
function of the group that can be well described by either a Kroupa power-law
function with indices $\alpha_3=-1.73\pm0.31$ and $\alpha_2=0.68\pm0.41$ in the
mass ranges $0.03\leq M/M_\odot\leq 0.08$ and $0.08\leq M/M_\odot\leq0.5$
respectively, or a Scalo log-normal function with coefficients
$m_c=0.21^{+0.02}_{-0.02}$ and $\sigma=0.36\pm0.03$ in the mass range $0.03\leq
M/M_\odot\leq0.8$. From the analysis of the spatial distribution of this
numerous candidate sample, we have confirmed the East-West elongation of the 25
Orionis group observed in previous works, and rule out a possible southern
extension of the group. We find that the spatial distributions of low mass
stars and brown dwarfs in 25 Orionis are statistically indistinguishable.
Finally, we found that the fraction of brown dwarfs showing IR excesses is
higher than for low mass stars, supporting the scenario in which the evolution
of circumstellar discs around the least massive objects could be more
prolonged.