## Mapping accretion and its variability in the young open cluster NGC 2264: a study based on u-band photometry ? ?? ,

L. Venuti 1 2 3 ; ; , J. Bouvier 1 2 ; , E. Flaccomio 4 , S. H. P. Alencar 5 , J. Irwin 6 , J. R. Stau GLYPH&lt;11&gt; er 7 , A. M. Cody 7 , P. S. Teixeira 8 , A. P. Sousa 5 , G. Micela 4 , J.-C. Cuillandre 9 , and G. Peres 3

- 1 Univ. Grenoble Alpes, IPAG, F-38000 Grenoble, France
- 2 CNRS, IPAG, F-38000 Grenoble, France
- e-mail: Laura.Venuti@obs.ujf-grenoble.fr
- 3 Dipartimento di Fisica e Chimica, Università degli Studi di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy
- 4 Istituto Nazionale di Astrofisica, Osservatorio Astronomico di Palermo G.S. Vaiana, Piazza del Parlamento 1, 90134 Palermo, Italy
- 5 Departamento de Física - ICEx - UFMG, Av. Antônio Carlos, 6627, 30270-901 Belo Horizonte, MG, Brazil
- 6 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA
- 7 Spitzer Science Center, California Institute of Technology, 1200 E California Blvd., Pasadena, CA 91125, USA
- 8 Universität Wien, Institut für Astrophysik, Türkenschanzstrasse 17, 1180, Vienna, Austria
- 9 Canada - France - Hawaii Telescope Corporation, 65-1238 Mamalahoa Highway, Kamuela, HI, 96743, USA

Received 7 March 2014 / Accepted 29 July 2014

## ABSTRACT

Context. The accretion process has a central role in the formation of stars and planets.

Aims. We aim at characterizing the accretion properties of several hundred members of the star-forming cluster NGC 2264 (3 Myr). Methods. We performed a deep ugri mapping as well as a simultaneous u -band + r -band monitoring of the star-forming region with CFHT MegaCam in order to directly probe the accretion process onto the star from UV excess measurements. Photometric properties / and stellar parameters are determined homogeneously for about 750 monitored young objects, spanning the mass range GLYPH&lt;24&gt; 0.1-2 M GLYPH&lt;12&gt; . About 40% of the sample are classical (accreting) T Tauri stars, based on various diagnostics (H GLYPH&lt;11&gt; , UV and IR excesses). The remaining non-accreting members define the (photospheric + chromospheric) reference UV emission level over which flux excess is detected and measured.

Results. We revise the membership status of cluster members based on UV accretion signatures, and report a new population of 50 CTTS candidates. A large range of UV excess is measured for the CTTS population, varying from a few times 0.1 to GLYPH&lt;24&gt; 3 mag. We convert these values to accretion luminosities and accretion rates, via a phenomenological description of the accretion shock emission. We thus obtain mass accretion rates ranging from a few 10 GLYPH&lt;0&gt; 10 to GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 7 M yr. Taking into account a mass-dependent GLYPH&lt;12&gt; / detection threshold for weakly accreting objects, we find a &gt; GLYPH&lt;27&gt; 6 correlation between mass accretion rate and stellar mass. A powerlaw fit, properly accounting for censored data (upper limits), yields ˙ Macc / M 1 4 : GLYPH&lt;6&gt; 0 3 : GLYPH&lt;3&gt; . At any given stellar mass, we find a large spread of accretion rates, extending over about 2 orders of magnitude. The monitoring of the UV excess on a timescale of a couple of weeks indicates that its variability typically amounts to 0.5 dex, i.e., much smaller than the observed spread in accretion rates. We suggest that a non-negligible age spread across the star-forming region may e GLYPH&lt;11&gt; ectively contribute to the observed spread in accretion rates at a given mass. In addition, di GLYPH&lt;11&gt; erent accretion mechanisms (like, e.g., short-lived accretion bursts vs. more stable funnel-flow accretion) may be associated to di GLYPH&lt;11&gt; erent ˙ Macc regimes.

Conclusions. A huge variety of accretion properties is observed for young stellar objects in the NGC 2264 cluster. While a definite correlation seems to hold between mass accretion rate and stellar mass over the mass range probed here, the origin of the large intrinsic spread observed in mass accretion rates at any given mass remains to be explored.

Key words. Accretion, accretion disks - Stars: formation - Stars: low-mass - Stars: pre-main sequence - open clusters and associations: individual: NGC 2264 - Ultraviolet: stars

## 1. Introduction

The disk accretion phase assumes a pivotal role in the scenarios of early stellar evolution and planetary formation. Circumstellar

- ? Based on observations obtained with MegaPrime MegaCam, a joint / project of CFHT and CEA DAPNIA, at the Canada-France-Hawaii / Telescope (CFHT) which is operated by the National Research Council (NRC) of Canada, the Institut National des Sciences de l'Univers of the Centre National de la Recherche Scientifique (CNRS) of France, and the University of Hawaii.
- ?? Tables 2, 3 and 4 are only available in electronic form at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http: // cdsweb.u-strasbg.fr cgi-bin qcat?J A / / / + A /

disks are the ubiquitous result of the earliest stages of star formation, surrounding the vast majority of solar-type protostars at an age of GLYPH&lt;24&gt; 1 Myr. Mass accretion from the inner disk to the central star regulates the star-disk interaction over the few subsequent million years; this process has a direct, long-lasting impact on the basic properties of the system, by providing at the same time mass and angular momentum. Angular momentum transport and mass infall have a most important part in dictating the dynamical evolution of protoplanetary disks, hence setting the environmental conditions which eventually lead to the formation of planetary systems (Armitage 2011).

The widely accepted paradigm for mass accretion in young stellar objects (YSOs) is that of magnetospheric accretion (Camenzind 1990; Königl 1991). A cavity of a few stellar radii extends from the stellar surface to the inner disk rim, as an e GLYPH&lt;11&gt; ect of the inner disk truncation by the stellar magnetosphere (B GLYPH&lt;24&gt; 2 kG). Accretion columns, channeled along the magnetic field lines, regulate the transport of material from the disk inner edge to the central star; as the accreting material impacts onto the stellar surface at near free-fall velocity, accretion shocks are produced near the magnetic poles. This picture is strongly supported by its predictive capability of the main observational features of accreting young solar-type stars (classical T Tauri stars, or CTTS): UV, optical and IR excess compared to the photospheric flux; spectral 'veiling'; broad emission lines; inverse P Cygni profiles; pronounced spectroscopic and photometric variability (e.g. Bouvier et al. 2007b).

Several studies (e.g. Muzerolle et al. 2003; Rigliaco et al. 2011) have shown that the rates ( ˙ Macc ) at which mass accretion occurs tend to scale with the mass of the central object (M ). GLYPH&lt;3&gt; Notwithstanding this general trend, a large spread in the ˙ Macc values at a given M GLYPH&lt;3&gt; appears to be a constant feature (e.g. Hartmann et al. 2006); this is indicative of a complex and dynamic picture, shaped by a variety of parameters and concurrent processes. Similarly, observations suggest a progressive decrease of ˙ Macc over the age of the systems (Hartmann et al. 1998; SiciliaAguilar et al. 2010), following the time evolution and dispersal of circumstellar disks (Haisch et al. 2001; Hillenbrand 2005). The large spread associated with the e GLYPH&lt;11&gt; ective relationship attests that a multiplicity of elements concur to the disk dissipation; moreover, the initial conditions in which stars are formed (e.g. Dullemond et al. 2006) and local environmental properties (e.g. Guarcello et al. 2007) are likely to significantly a GLYPH&lt;11&gt; ect the evolutionary pattern of individual star-disk systems.

Statistical surveys of the accretion properties of young stars provide essential insight into the dynamics of accretion. Large, homogeneously characterized samples of objects, spanning several ranges of parameters, are crucial to achieve a proper understanding of the physics governing the disk accretion history of forming stars and to explore di GLYPH&lt;11&gt; erent facets of the accretion mechanisms. To this respect, extensive mapping surveys of wellpopulated star-forming regions at UV wavelengths (e.g. Rigliaco et al. 2011; Manara et al. 2012) provide a most interesting standpoint, as they allow to directly probe accretion for large samples of stars from the excess emission produced in the shocked impact layer at the base of the accretion column.

Such single epoch surveys e GLYPH&lt;14&gt; ciently trace the instantaneous spread of accretion properties in a given star-forming region; however, they are unable to probe the nature of such a spread, i.e. to determine whether this is due to variability of individual objects or to an intrinsic spread amongst objects. And yet, a key aspect of the accretion dynamics is the intrinsic variability profile of the process (i.e., the amplitude and characteristic timescales of the variability in ˙ Macc ). The enhanced variability typical of young stars, cornerstone of the very definition of the T Tauri class (Joy 1945), is manifest over a broad range of wavelengths (X-rays, UV, optical, IR) and of time baselines (hours, days, years); a variety of processes concur to shape the variability pattern of individual objects, with a significant contribution coming from the geometry of the system. Understanding to which extent and on which timescales the accretion process is intrinsically variable, and characterizing its variability on a statistical basis compared to individual cases, is of utmost importance in order to accurately interpret the information provided by large accretion surveys and thus infer a detailed picture of the accretion process.

Multi-epoch surveys hence provide a more accurate description of the accretion process, as they allow to probe the time evolution of single snapshots of the accretion properties and, thus, to investigate how accretion occurs. This has been addressed in a few recent studies (e.g. Nguyen et al. 2009; Fang et al. 2013), which explored the variability of spectral accretion diagnostics for tens of objects with a few observing epochs over several months. Such e GLYPH&lt;11&gt; orts provided very valuable constraints on the average accretion variability shown by typical CTTS. However, a much tighter time sampling is needed in order to recover an estimate of the e GLYPH&lt;11&gt; ective intrinsic accretion variability, isolated from geometric contributions which are simply due to a varying visible portion of the accretion spots during stellar rotation.

We have recently conducted an extensive observational program specifically aimed at addressing the issue of YSOs variability over the full spectrum. The Coordinated Synoptic Investigation of NGC 2264 (CSI 2264; Cody et al. 2013) project was devised as a coordinated multi-wavelength exploration of the hours-to-weeks variability of the pre-Main Sequence (PMS) population of the star-forming region NGC 2264. The space telescopes CoRoT and Spitzer provided the backbone of the optical and IR investigations, with simultaneous continuous monitoring over 40 and 30 days respectively; the Flames multi-object spectrograph at the Very Large Telescope (VLT) provided a full set of spectra covering 20 di GLYPH&lt;11&gt; erent epochs for GLYPH&lt;24&gt; 100 young stars, while multi-band optical + UV photometric monitoring obtained at the Canada-France-Hawaii Telescope (CFHT) provided contemporaneous, synoptic measures of the accretion rates and their variability from the direct diagnostics of the UV excess, with an overall time coverage of 14 consecutive days and several measurements per night. The targeted region, NGC 2264, has long been a benchmark for star formation studies, thanks to its youth ( GLYPH&lt;24&gt; 3 Myr), its richness ( GLYPH&lt;24&gt; 1500 known members, both CTTS and non-accreting, weak-lined T Tauri stars, or WTTS), its relative proximity (d GLYPH&lt;24&gt; 800 pc), the low extinction on the line of sight and the association with a molecular cloud complex that significantly reduces the contamination from background stars (see Dahm 2008 for a recent review on the region).

As part of the CSI 2264 campaign, this paper specifically focuses on the issue of accretion properties and accretion variability in the PMS population of NGC 2264. Sect. 2 describes the photometric dataset obtained at CFHT, on which the study reported here is based. Sect. 3 reports on the characterization of photometric accretion signatures, which allowed us to perform a new, accretion-driven census of NGC 2264 members; an extensive, homogeneous investigation of individual stellar properties is also reported. In Sect. 4, we describe the procedure adopted to convert the UV excess measurement to a ˙ Macc estimate; we analyze the ˙ Macc over M GLYPH&lt;3&gt; distribution thus inferred and compare this picture with the amount of variability registered over a baseline of two weeks; we further probe the origin of this variability and infer a specific estimate of the intrinsic accretion variability as opposed to rotational modulation detected during the monitoring. Results are discussed in Sect. 5; our conclusions are synthetized in Sect. 6. This provides a complete picture of a whole star-forming region and its several hundred members at short wavelengths, from a specifically accretion-oriented perspective. A more detailed characterization of the variability properties monitored at CFHT and their relation to the underlying physics will be addressed in a forthcoming paper (Venuti et al., in preparation).

Fig. 1. CFHT MegaCam / u g r color image of NGC 2264.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-9f69f9e7-db36-45be-8a05-524bbe7ca0f0.jpg)

## 2. Observations and data reduction

2.1. CFHT MegaCam surveys of NGC 2264: data reduction and photometric calibration

Wesurveyed NGC 2264 in two multiwavelength observing campaigns at the CFHT, using the wide-field optical camera MegaCam; the camera has a field of view (FOV) of GLYPH&lt;24&gt; 1 deg. 2 (Boulade et al. 2003), thus fitting the whole region in a single telescope pointing. Five broad-band filters ( u GLYPH&lt;3&gt; ; g ; 0 r 0 ; i 0 ; z 0 ) are adopted at the camera; their design, close to the Sloan Digital Sky Survey (SDSS) photometric system (Fukugita et al. 1996), ensures high e GLYPH&lt;14&gt; ciency for faint object detection and deep sky mapping.

## 2.1.1. Mapping survey

The first MegaCam survey of NGC 2264, held on Dec. 12, 2010, consisted of a deep u GLYPH&lt;3&gt; g 0 r 0 mapping of the region. In the r 0 -band, five short exposures (10 s) were obtained, using a dithering pattern in order to compensate for the presence of bad pixels and small gaps on the CCD mosaic. In the u GLYPH&lt;3&gt; -band and g 0 -band, the same 5-step dithering procedure was repeated in two di GLYPH&lt;11&gt; erent exposure modes, short (10 s) and long (60 s), in order to obtain a good signal noise (S N) for the whole sample of objects with-/ / out saturating the brightest sources. All images were obtained over a continuous temporal sequence within the night. On Feb. 28, 2012, concurrently with the monitoring survey reported in Sect. 2.1.2, we additionally obtained a deep i 0 -band mapping, in order to reconstruct a full 4-band picture of NGC 2264 together with the u GLYPH&lt;3&gt; g 0 r 0 survey of Dec. 2010. Similarly to previous observations, a dithering pattern was used, along with two di GLYPH&lt;11&gt; erent exposure times: 180 s, to detect all sources in the FOV, and 5 s, to recover the bright sources saturated in the long exposure.

Raw data have been preliminarily processed at the CFHT, in order to correct for instrumental signatures across the whole CCD mosaic (i.e., bad pixel clusters or columns, dark current, bias, non-uniformity in the response), and a first astrometric and photometric calibration has been performed. The subsequent photometric processing of the images obtained during the survey has been handled as described below.

For each spectral band and each exposure time, single dithered images have been combined, in order to discard spurious signal detections (e.g., cosmic rays) and correct for small CCD gaps. The photometric images thus assembled (one in the r 0 -band and two for each other filter, corresponding to the two di GLYPH&lt;11&gt; erent exposure times) have then been processed using the SExtractor (Bertin &amp; Arnouts 1996) and PSFEx (Bertin 2011) tools, in order to identify true sources in the field of view and extract the relevant photometry. At first, we run a preliminary aperture photometry procedure with a detection threshold of 10 GLYPH&lt;27&gt; above the local background, in order to extract only brighter sources; this selection of objects has been used as input to run a pre-PSF fitting routine, aimed at testing the optimal parameters for the PSF function to be used for analyzing the images. A final, complete single-band catalog has eventually been produced from each image processed with SExtractor by performing a new PSF-fitting run with a lower detection threshold of 3 GLYPH&lt;27&gt; , thus recovering the fainter component of the population.

A first filtering of each catalog population, based on the properties of the extracted sources and on the relevant extraction flags, allowed us to discard the stellar component potentially a GLYPH&lt;11&gt; ected by saturation or other detection problems. Long and short exposure catalogs corresponding to the same spectral band, when present, have then been merged, in order to obtain a unique photometric catalog for each spectral band and avoid duplicates. The sample of objects common to both the longexposure and short-exposure catalogs has been used to evaluate the photometric accuracy, to correct the short-exposure photometry for small magnitude o GLYPH&lt;11&gt; sets with photometry obtained from the long exposures and to define a confidence magnitude threshold below which photometry could begin to be a GLYPH&lt;11&gt; ected by saturation. Objects have been preferentially taken from the long-exposure catalog; sources brighter than the long-exposure threshold have been retrieved from the short-exposure catalog, when available; sources brighter than the short-exposure threshold have been discarded.

A complete calibration of the CFHT instrumental magnitudes to the SDSS photometric system, including a correction for both o GLYPH&lt;11&gt; set and color e GLYPH&lt;11&gt; ect, has then been performed. The calibration procedure has not been anchored to a sample of standard stars in the field of view; instead, we used an ensemble of around 3000 stars, present in the CFHT FOV, for which SDSS photometry was available from the seventh SDSS Data Release (Abazajian et al. 2009), to statistically calibrate the CFHT magnitudes to the SDSS system. The matching of sources between the SDSS and the CFHT catalogs has been performed using the TOPCAT (Taylor 2005) tool, fixing a matching radius of 1 arcsec and retaining only closest matches. The statistical ensemble of stars used for the calibration has been selected among the objects with complete CFHT u GLYPH&lt;3&gt; g 0 r i 0 0 photometry and that matched the magnitude limits for 95% detection repeatability for point sources defined for the SDSS photometric survey (Abazajian et al. 2004). The calibration procedure adopted is based on the results of the photometric calibration study of the CFHT MegaCam photom-/ etry to the SDSS system performed by Regnault et al. (2009). We compared the distribution of SDSS colors over the ensemble of stars with the SDSS dwarf color loci tabulated in literature (Covey et al. 2007; Lenz et al. 1998) and retained, for calibration purposes, only objects whose colors matched the tabulated ranges. After checking their consistency with the observed trends, Regnault's linear calibration relations have then been used to fit the ( m CFHT GLYPH&lt;21&gt; GLYPH&lt;0&gt; m S DSS GLYPH&lt;21&gt; ) vs. ( m S DSS GLYPH&lt;21&gt; GLYPH&lt;0&gt; m S DSS GLYPH&lt;21&gt; + 1 ) distri-

butions in u GLYPH&lt;3&gt; , g 0 and r 0 , by keeping the angular coe GLYPH&lt;14&gt; cient fixed and adapting the intercept value. In the i 0 -band, for which the filter in use at MegaPrime has been changed, due to accidental breakage, after the epoch of acquisition of data used in Regnault et al. (2009)'s study, no evident color dependence of the o GLYPH&lt;11&gt; set between the two photometric systems has been observed; the i 0 -band photometry has therefore been uniquely corrected for a constant o GLYPH&lt;11&gt; set. The system of 4 calibration equations thus obtained (one for each spectral band; Eq. 1-4) has then been solved to obtain the final conversion from CFHT instrumental magnitudes to SDSS photometry.

$$u _ { C F H T } - u _ { S D S S } = - 0. 2 1 1 \, ( u _ { S D S S } - g _ { S D S S } ) - 0. 6 3 \, \quad \ \ ( 1 )$$

$$g _ { C F H T } - g _ { S D S S } = - 0. 1 5 5 \, ( g _ { S D S S } - r _ { S D S S } ) - 0. 5 1 \, \quad \ \ ( 2 ) \, \quad \ w e \, \text{exam}$$

$$r _ { C F H T } - r _ { S D S S } & = - 0. 0 3 \left ( r _ { S D S S } - i _ { S D S S } \right ) & ( 3 ) \underset { \text{sequenc} } { \text{$\text{$\text{$max\infty$} } } }$$

$$i _ { C F H T } - i _ { S D S S } & = - 0. 0 2 & ( 4 ) \quad \text{to the } \mathfrak { m } _ { \mathfrak { a. c. d. } }$$

Acomplete u GLYPH&lt;3&gt; g 0 r i 0 0 catalog of NGC 2264, containing GLYPH&lt;24&gt; 9000 sources in the field of view, has been built by cross-correlating the single-band catalogs on TOPCAT. A matching radius of 1 arcsec has been fixed for identifying common sources detected in the u GLYPH&lt;3&gt; ; g ; 0 r 0 concurrent observations, while a larger matching radius of 2 arcsecs has been introduced to match u GLYPH&lt;3&gt; g 0 r 0 sources with their i 0 -band counterpart, to account for a lower astrometric accuracy (lower angular resolution) characterizing the i 0 -band field, due to poorer seeing conditions at the time of acquisition of the i 0 -band images. Sources have been identified by retaining only best matches; we address potential problems in the derived colors associated with misidentifications and multiple matches in Sect. 2.2. Final SDSS u 0 g 0 r i 0 0 (hereafter u g ri ) photometry has then been associated to each source in the catalog using the derived calibration relations.

## 2.1.2. Monitoring survey

The second MegaCam survey of NGC 2264, performed in Feb. 2012, consisted in monitoring the u GLYPH&lt;3&gt; -band and r 0 -band variability of the entire stellar population of the region over a mid-term timescale (Feb.14-28, i.e., 2 weeks), as a part of the CSI 2264 project. On each observing night during the run, the region was imaged repeatedly, with a temporal cadence varying from 20 min to 1.5 h. Each u r GLYPH&lt;3&gt; 0 observing block was performed using a 5-step dithering pattern, with single exposures of 3 s in the r 0 -band and 60 s in the u GLYPH&lt;3&gt; -band.

The procedure adopted for the photometric processing of the images obtained during the u r GLYPH&lt;3&gt; 0 monitoring survey is similar to that applied to the mapping survey images, but each exposure has been individually processed in order to retrieve the luminosity variations over di GLYPH&lt;11&gt; erent timescales. The i -band image (Sect. 2.1.1) has been used as the astrometric reference exposure, as the nebulosity component was the least conspicuous in this image compared to the other filters; for each source, an aperture has been placed down on each u GLYPH&lt;3&gt; and r 0 image at the position predicted from the i -band source detection and the amount of flux within this aperture has been extracted to build the light curve.

A preliminary quality check on the CFHT light curves consisted in closely examining a sample of NGC 2264 members with known periodicities, in order to attempt to recover their periods from the phase-folded CFHT light curves. This allowed us to ascertain the global accuracy of the u r GLYPH&lt;3&gt; 0 time series photometry, albeit occasionally a GLYPH&lt;11&gt; ected by isolated discrepant points and non-negligible scatter due to poorer observing conditions. In order to locate the potentially problematic observing sequences,

Table 1. Detection limits and completeness in the CFHT u g ri survey of NGC 2264.

| (mag)            | u       | g       | r         | i       |
|------------------|---------|---------|-----------|---------|
| mag range *      | 15-23.5 | 14-21.5 | 13.5-20.5 | 13-19.5 |
| Saturation start | < 12.5  | < 13    | < 13.5    | 12.5    |
| Detection limit  | 23.5    | 21.5    | 20.5      | 19.5    |
| Completeness     | 21.5    | 19.5    | 18.5      | 17.5    |

Notes. Values reported are conservative estimates.

( * ) Values referred to the main body of the magnitude distribution of the catalog population.

we examined the time series photometry for the whole sample of field stars in the CFHT FOV and derived, for each observing sequence, the distribution of median zero-point o GLYPH&lt;11&gt; sets referred to the master frame. The zero-point o GLYPH&lt;11&gt; set, for a given object, is defined as the mag di GLYPH&lt;11&gt; erence between the frame of interest and the master frame; in case of large deviation from zero, the zeropoint o GLYPH&lt;11&gt; set distribution associated to an observing sequence thus signals that the observation has been carried out through clouds. This translates to significantly lower than average S N for a given / source, which implies that, even for the brightest stars, no accurate measurements can be inferred from that specific observation; hence, epochs matching this description have not been used for any subsequent analysis.

As done earlier for the u g ri snapshot survey, the instrumental light curve photometry has been calibrated to the SDSS system, referring the procedure to the calibrated photometry from the first CFHT survey. No o GLYPH&lt;11&gt; sets and color e GLYPH&lt;11&gt; ects have been noticed in the r 0 -band (hereafter r -band) photometry, while a correction for both contributions has been applied to the u GLYPH&lt;3&gt; -band photometry from the ( u GLYPH&lt;3&gt; CFHT lc \_ GLYPH&lt;0&gt; u ) vs. ( u GLYPH&lt;0&gt; r ) plane:

$$u _ { C F H T _ { - l c } } ^ { * } - u = - 0. 1 3 4 3 \ ( u - r ) - 0. 1 0 0 9$$

## 2.2. The CFHT catalog of NGC 2264

Fig. 1 shows a MegaCam color picture of the field imaged at CFHT during the NGC 2264 campaign. The star-forming region extends over the central part of the field, with the most active sites of star formation located at the center of the image, northward of the Cone Nebula (cf. Dahm 2008).

We obtained a complete u g ri dataset and photometric monitoring for GLYPH&lt;24&gt; 9000 sources in the field, exploring the area projected onto the star-forming region as well as the periphery of the cluster and background foreground sky. Table 1 provides some de-/ tails on the photometric properties of the population of the CFHT catalog. The survey is complete down to u GLYPH&lt;24&gt; 21.5 and r GLYPH&lt;24&gt; 18.5 and the monitored objects span a range of GLYPH&lt;24&gt; 7 mags. The photometry is globally accurate up to a few GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 2 mag and the relative accuracy rises up to order of 10 GLYPH&lt;0&gt; 3 mag for brighter objects.

In order to evaluate the probability of misidentifications as a result of the procedure adopted for matching sources in different catalogs, and thus evaluate the potential impact on the photometric properties inferred for individual sources, we reconsidered all multiple identifications within the allowed matching radii and closely examined their spatial distribution. The sources with potentially multiple identifications amount to GLYPH&lt;24&gt; 2.9% of the population of the catalog. The di GLYPH&lt;11&gt; erence in o GLYPH&lt;11&gt; set distance between the best match and the next best match is less than 0.5 arcsec for only 23% of these cases. Less than 20% of the multiply matched sources are located within the region occupied by

Fig. 2. ( r GLYPH&lt;0&gt; i , g GLYPH&lt;0&gt; r ) color-color diagram for objects monitored at CFHT. Empirical dwarf color sequences from Covey et al. (2007) and Kraus &amp; Hillenbrand (2007) are shown as black crosses and blue triangles, respectively. The reddening vector is traced based on the reddening parameters in the SDSS system provided by the Asiago Database on Photometric Systems (ADPS; Fiorucci &amp; Munari 2003).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-ad740bb6-cd0a-421b-b1c4-b427ece9f709.jpg)

most NGC 2264 members. For the statistical purposes of this study, this component is assumed to be negligible.

## 3. The PMS population of NGC 2264

A global picture of the photometric properties of the stellar population can be achieved from the analysis of color-color diagrams, as shown in Figs. 2-3. Empirical dwarf SDSS color sequences (Covey et al. 2007; Kraus &amp; Hillenbrand 2007) have been overplotted to the diagrams in order to follow the observed color distribution and investigate the relationship between color loci and stellar properties.

In Fig. 2, earlier-type field stars are located along the diagonal ellipse, with foreground or slightly reddened objects on the lowermost part of the distribution and background stars at the uppermost end; M-type stars are predominantly located on the horizontal branch, whose g GLYPH&lt;0&gt; r position is in good agreement with the value of GLYPH&lt;24&gt; 1.4 mentioned in previous studies (Covey et al. 2007; Finlator et al. 2000; Fukugita et al. 1996; Annis et al. 2011). A number of objects can be located to the left (i.e., blueward) of the standard locus in r GLYPH&lt;0&gt; i ; only for a part of them this displacement may be ascribed to photometric inaccuracies. The most likely explanation for these objects is that they are highly variable, and the inferred r GLYPH&lt;0&gt; i color is spurious due to the large epoch di GLYPH&lt;11&gt; erence between the r -band and i -band photometry (the over one year time lapse between the u g r and the i -band mappings of the region, as discussed in Sect. 2).

Similarly, a progression of increasing spectral types and reddening e GLYPH&lt;11&gt; ects can be observed from left to right along the ellipse in Fig. 3; an interesting locus is defined by the point dispersion at 0.5 GLYPH&lt;20&gt; g GLYPH&lt;0&gt; r GLYPH&lt;20&gt; 2.0, below the main body of the color distribution, as it shows u GLYPH&lt;0&gt; g color excesses presumably linked with the YSO accretion activity, manifest at UV wavelengths with a flux ex-

Fig. 3. ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) color-color diagram for objects monitored at CFHT. Empirical dwarf color sequences from Covey et al. (2007) and Kraus &amp; Hillenbrand (2007) are shown as black crosses and blue triangles, respectively. The reddening vector is traced based on the reddening parameters in the SDSS system provided by the ADPS (Fiorucci &amp; Munari 2003).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-965f81c1-40bb-4747-be41-8d8d07c54dae.jpg)

cess compared to the photospheric emission. An o GLYPH&lt;11&gt; set of a few tenths of a mag is observed between the mean locus of field stars in our survey and the sequences of Covey et al. (2007) and Kraus &amp; Hillenbrand (2007). This o GLYPH&lt;11&gt; set cannot be explained with imprecisions in our calibration. A possible explanation may lie in a lower than solar standards metallicity of the field population probed here; at any rate, this disagreement has no impact on our analysis, as we uniquely used our internal reference sequence to probe the photometric accretion signatures of members, as illustrated in the following sections.

We further pursued the color signatures of di GLYPH&lt;11&gt; erent YSO types in the SDSS filters by analyzing the color properties of known members of the region with respect to the loci traced on the diagrams by the large population of field stars. Known members were identified amidst the population of the region based on a wide collection of photometric spectroscopic data available / in the literature from various surveys of the cluster; membership and further information on the WTTS (no disk evidence) vs. CTTS (disk-accreting) nature of selected members have been inferred based on one or more of the following criteria: i) H GLYPH&lt;11&gt; emitter from narrow-band photometry and variability from the data of Lamm et al. (2004) and following their criteria; ii) X-ray detection (Ramírez et al. 2004; Flaccomio et al. 2006) and location on the cluster sequence in the (I, R-I) diagram when R,I photometry is available; iii) H GLYPH&lt;11&gt; emitter from spectroscopy (H GLYPH&lt;11&gt; EW &gt; 10 Å, H GLYPH&lt;11&gt; width at 10% intensity &gt; 270 km s GLYPH&lt;0&gt; 1 ); iv) radial velocity member according to F˝ urész et al. (2006); v) member according to Sung et al. (2008) (cf. criteria enumerated earlier) and Sung et al. (2009) (Spitzer Class I Class II). It is worth remarking that / this preliminary analysis of membership and TTS classification is not based on direct accretion criteria.

Around 700 known members have been matched in the CFHT catalog within a matching radius of 1 arcsec; among

Fig. 4. ( r GLYPH&lt;0&gt; i , g GLYPH&lt;0&gt; r ) color-color diagram for NGC 2264 members monitored at CFHT. Field stars, non-accreting cluster members (WTTS) and accreting cluster members (CTTS) are depicted as grey dots, blue triangles and red circles, respectively. The reddening vector is traced based on the reddening parameters reported in the ADPS (Fiorucci &amp; Munari 2003).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-ca719500-3e9d-44a8-8219-e22e43d95f6e.jpg)

these, GLYPH&lt;24&gt; 60% are classified as WTTS. Fig. 4 and 5 show the color properties of these groups of objects respectively in the ( r GLYPH&lt;0&gt; i , g GLYPH&lt;0&gt; r ) and ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) diagrams.

In redder filters, as shown in Fig. 4, both CTTS and WTTS define a single color sequence following the one traced by field stars. A few CTTS appear well blueward of this sequence, reflecting multi-year variability that is evidently more common in the CTTS than the WTTS. Conversely, accreting and nonaccreting stars define two clear distinct distributions on the ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) diagram in Fig. 5, with the latter displaying colors consistent with the photometric properties of field stars and the former located at smaller (i.e., bluer) u GLYPH&lt;0&gt; g values compared to the main color locus defined by dwarfs.

These distinctive features are well observed in Fig. 6, which shows how members distribute on the ( u GLYPH&lt;0&gt; r , r ) color-magnitude diagram. There, the cluster sequence is clearly traced by WTTS, while CTTS appear broadly scattered to the left as a result of their u -band excess, indicative of active accretion activity, compared to WTTS. A number of presumed non-accreting members show a non-negligible displacement blueward redward of / the main locus; these two groups of objects are discussed further in Sect. 3.2.2.

Two main inferences can be drawn from these diagrams:

- 1. CTTS and WTTS members do di GLYPH&lt;11&gt; erentiate in the SDSS colors on the basis of photometric accretion signatures at short wavelengths (UV flux excess);
- 2. the UV excess linked with accretion, detected in the u -band, does not a GLYPH&lt;11&gt; ect significantly the observations in filters at longer wavelengths, as emphasized by the substantial agreement of the CTTS and WTTS color sequences in Fig. 4 and the color saturation at g GLYPH&lt;0&gt; r GLYPH&lt;24&gt; 1.4, a behavior common to field stars.

Fig. 5. ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) color-color diagram for NGC 2264 members monitored at CFHT. Field stars, WTTS and CTTS are depicted as grey dots, blue triangles and red circles, respectively. The reddening vector is traced based on the reddening parameters reported in the ADPS (Fiorucci &amp; Munari 2003).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-bd6727b8-eb90-4263-9830-91da9f71b29a.jpg)

These photometric properties, observed for di GLYPH&lt;11&gt; erent ensembles of stars, enable us to perform a new membership and population study on an individual basis, from an accretion-driven per-

Fig. 6. ( u GLYPH&lt;0&gt; r , r ) color-magnitude diagram for NGC 2264 members monitored at CFHT. Field stars, WTTS and CTTS are depicted as grey dots, blue triangles and red circles, respectively. The reddening vector is traced based on the reddening parameters reported in the ADPS (Fiorucci &amp; Munari 2003).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-ac921baa-6c51-480a-80e1-9566ba321eba.jpg)

spective. This analysis has the twofold purposes of investigating the presence, in the CFHT FOV, of additional objects that are new candidate members, and of re-examining the classification of known members, looking in both cases for classification outliers in the diagram loci dominated by accretion.

## 3.1. UV excess vs. different accretion diagnostics

In order to ascertain the coherence of the classification based on the UV excess detection on these diagrams and define some confidence limits for a reliable identification of accreting members, we compared the u -band excess information drawn from CFHT photometry with di GLYPH&lt;11&gt; erent diagnostics commonly used to investigate accretion, such as the equivalent width and the width at 10% intensity in the H GLYPH&lt;11&gt; emission line, probing the accretion funnel, or the mid-infrared emission, that signals the presence of material in the inner disk when a flux excess compared to the photospheric level is detected.

Fig. 7 shows the same color-color diagram as in Fig. 5 with additional information on the H GLYPH&lt;11&gt; EW values measured for the subsample of objects having VLT Flames spectra from the / CSI 2264 campaign (A. Sousa, UFMG). A threshold value commonly adopted to identify accreting YSOs from H GLYPH&lt;11&gt; emission is H GLYPH&lt;11&gt; EW GLYPH&lt;21&gt; 10 Å (Herbig &amp; Bell 1988), while a more robust scheme relating the threshold value to the stellar spectral type has been introduced by White &amp; Basri (2003). As can be observed in Fig. 7, the UV excess and the H GLYPH&lt;11&gt; EW diagnostics are fully consistent, with the objects having the least H GLYPH&lt;11&gt; EWs located on the color locus traced by WTTS and field stars, while a region of the diagram dominated by accretion can clearly be delimited below the field stars distribution and identified from the detection of a UV excess.

A useful tool for characterizing presence and properties of a circumstellar disk is the examination of the infrared spectral

Fig. 7. u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r colors for NGC 2264 members are compared to the H GLYPH&lt;11&gt; EW. Asterisks mark members for which H GLYPH&lt;11&gt; EW is available from the CSI 2264 campaign (see text). Asterisk colors and sizes are scaled according to the value of H GLYPH&lt;11&gt; EW.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-e6779e80-d92a-4a60-a202-367abae9dbbb.jpg)

energy distribution (SED) of the system. A particularly interesting indicator of the SED morphology is the measure of its slope in the spectral range of interest (i.e., the spectral index GLYPH&lt;11&gt; first introduced by Lada 1987): smaller amounts of material in the inner disk result in smaller contributions, i.e. smaller flux excesses, at mid-IR wavelengths and hence more negative values of the GLYPH&lt;11&gt; -index, which approaches the slope of a reddened blackbody, while positive values of the mid-IR SED slope are specific to deeply embedded, Class I sources. In Lada et al. (2006), the GLYPH&lt;11&gt; -index diagnostics for disks is applied to the mid-IR emission sampling between 3.6 and 8 GLYPH&lt;22&gt; m performed with the Spitzer instrument IRAC for the young cluster IC 348; the same tool has been adopted by Teixeira et al. (2012) to characterize the evolutionary state of NGC 2264 members based on the inner disk properties. We compare their results to the colors of NGC 2264 YSOs in the ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) diagram in Fig. 8. Following the classification proposed in Table 2 of Teixeira et al. (2012), values of GLYPH&lt;11&gt; IRAC &lt; GLYPH&lt;24&gt; GLYPH&lt;0&gt; 2 are indicative of objects with anaemic disks or naked photospheres, while objects with thick disks (hence likely to be actively accreting) are characterized by GLYPH&lt;11&gt; IRAC &gt; GLYPH&lt;0&gt; 2. Indeed, as can be observed in Fig. 8, objects in the first GLYPH&lt;11&gt; IRAC group are mainly located on the WTTS dwarf locus in / u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r , while objects showing a UV excess (and thus actively accreting stars) have the largest values of GLYPH&lt;11&gt; IRAC , as globally expected.

## 3.2. UV census of the NGC 2264 population

## 3.2.1. New CTTS candidates in NGC 2264

A most interesting locus to identify accreting members is the data points dispersion below the main body of the distribution in Fig. 5. In fact, this allows easy detection and measurement of the u -band excess with respect to WTTS. The inferred color excess estimate is independent of distance and has no significant dependence on visual extinction A V , as WTTS define a color

Fig. 8. u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r colors for NGC 2264 members are compared to the GLYPH&lt;11&gt; IRAC diagnostics, which probes the inner disk properties (see text). Colors and sizes of symbols are scaled according to the value of GLYPH&lt;11&gt; IRAC from the data of Teixeira et al. (2012).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-a2b59523-4c83-4aad-9790-70d07f61d441.jpg)

Fig. 9. Selection of new CTTS candidates from the ( g GLYPH&lt;0&gt; r , u GLYPH&lt;0&gt; g ) diagram. The grey dashed line traces a conservative boundary separating the accretion-dominated region of the diagram from the color locus of non-accreting stars (see text). Field stars, WTTS and known CTTS are depicted as grey dots, blue triangles and red circles, respectively; the new CTTS candidates selected from this diagram are marked as green crosses. The purple line fitting the upper envelope of the WTTS distribution traces the reference sequence of photospheric colors that has been adopted for measuring the UV excess of accreting stars (see Sect. 4).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-809c00f8-f0b5-4477-a92a-ef58081c5eb0.jpg)

sequence nearly parallel to the reddening vector down to GLYPH&lt;24&gt; M0 spectral type. For later-type stars, the color sequence traced by WTTS turns down, saturating at g GLYPH&lt;0&gt; r GLYPH&lt;24&gt; 1.4 (cf. Fig. 4); the implications of this for the UV excess measurement are discussed in Sect. 4. In order to minimize the contamination, we traced a straight line, parallel to the reddening vector, bordering the accretion-dominated area below the lower envelope of the dwarf color locus (Fig. 9); we then selected as new CTTS candidates all the objects lying in this region and marked as field stars. For all these objects, a preliminary check of the photometry on the CFHT field images has been performed, in order to discard the sources whose flux measurements could be a GLYPH&lt;11&gt; ected by detection issues (e.g., too faint source, presence of two close sources, partial projection of the source on CCD gaps or bad pixel rows).

A second group of new CTTS candidates has been extracted from the color distribution in Fig. 4, namely among the color outliers located to the left and, to a smaller extent, below the main color distribution. Similarly to the procedure adopted for the selection of objects on the diagram in Fig. 5, a check on the photometry quality for all the objects of interest has preceded the identification of new candidates.

A third color-color diagram of interest for the analysis of membership is the one depicting u GLYPH&lt;0&gt; g vs. r GLYPH&lt;0&gt; i (Fig. 10); as a 'transitional' diagram between u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r (Fig. 5; direct UV excess detection) and g GLYPH&lt;0&gt; r vs. r GLYPH&lt;0&gt; i (Fig. 4; nearly horizontal photospheric color sequence for late-type stars, with no color saturation on r GLYPH&lt;0&gt; i ), this diagram provides a greater sensitivity to the UV excess for spectral types GLYPH&lt;21&gt; M0. Additional, later-type CTTS candidates have then been identified on this diagram, and

Fig. 10. ( r GLYPH&lt;0&gt; i , u GLYPH&lt;0&gt; g ) diagram for cluster members and field stars in the NGC 2264 field. Field stars, WTTS and known CTTS are depicted as grey dots, blue triangles and red circles, respectively. Grey dots encircled in red correspond to new CTTS candidates selected on di GLYPH&lt;11&gt; erent diagrams (Fig. 9 or Fig. 4.) Green crosses mark the additional CTTS candidates selected on this diagram.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-7a741280-10fd-40af-8d41-cb8c578fbcd6.jpg)

their photometric properties inspected similarly to what done for the candidate members extracted earlier.

A second aspect of our selection of new candidate, actively accreting YSOs concerns the group of already known members which had previously been considered WTTS (that is, YSOs without evidence of significant, ongoing accretion at the time of the previous surveys). We examined the locations of these putative WTTS in our several diagrams (Figures 4 through 9) and identified a significant number whose colors fell generally within the UV-bright regions - for example, the small group of blue points located below the normal dwarf locus in Fig. 5. This procedure allowed us to re-classify 19 WTTS members as CTTS, which add to the 50 newly accretion-identified candidates; this varied status does not necessarily reflect an erroneous previous classification, but may as well be the result of a long-term ( GLYPH&lt;24&gt; years) variability in the accretion activity of the objects.

The results of our UV-based membership investigation are provided in Table A.1 (new CTTS candidates) and Table A.2 (re-classified CTTS). This revised census of accreting members we infer for NGC 2264 likely encompasses all objects showing a strong UV excess, hence a significant amount of accretion. However, our sample might not be exhaustive regarding the weakly-accreting CTTS component. Indeed, a non-negligible overlap between the CTTS and WTTS populations can be seen in Figs. 5-6, which cannot be e GLYPH&lt;14&gt; ciently probed with our diagnostics. Hence, we do identify some new CTTS, but we might be missing a number of lower accretors, showing a smaller UV excess that cannot be easily separated from the color properties of WTTS field stars. /

The wide FOV of CFHT MegaCam allowed us to probe a / quite extended region around the cluster, encompassing most of the fields of view adopted in previous surveys. It is thus of interest to examine the spatial distribution of the new CTTS candi-

Fig. 11. Spatial distribution of the new CTTS candidates from the CFHT u -band survey. Field stars, known cluster members and new CTTS candidates are depicted as small grey dots, blue dots and black crosses, respectively.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-bca3b9f6-de19-4f62-bc2b-5a2131cb9f10.jpg)

dates, selected as described in Sect. 3.2.1, compared to the distribution of known NGC 2264 members on the CFHT field of view. This comparison is shown in Fig. 11. A fraction of the object positions are projected onto the main cluster region, while a significant part are located on the periphery or beyond. This provides evidence for a wider spatial spread of the PMS population associated with NGC 2264, and attests to the presence of actively accreting members far from the original sites of star formation, currently known to be located towards the north and the center of the field imaged at CFHT (Dahm 2008).

## 3.2.2. Field stars contamination in the NGC 2264 sample

Based on the discussion in Sect. 3.2.1, another potentially interesting location for our accreting members census is the data point dispersion to the left of the cluster sequence (WTTS distribution) in Fig. 6. A direct search for candidate members in this color-magnitude locus would hardly be straightforward, due to the distance-dependence and to the vast preponderance of field stars; on the other hand, this diagram provides an additional test bench for investigating the nature of individual groups of objects.

Similarly to Fig. 5, we can identify on Fig. 6 a group of putative non-accreting members with inferred u GLYPH&lt;0&gt; r colors bluer than expected. This feature, suggestive of ongoing accretion activity for a young cluster member, is consistently observed on di GLYPH&lt;11&gt; erent color-magnitude diagrams ( g vs. u GLYPH&lt;0&gt; g , i vs. u GLYPH&lt;0&gt; i ) for several of these objects; interestingly, these objects populate the thin lower branch of the WTTS distribution, below the main WTTS locus, on the lower part of the dwarf locus in Fig. 5 ( u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r diagram). Little spectroscopic information is available for this group of stars; contrary to what is expected for objects clearly showing accretion activity, the level of variability exhibited by many of these objects during our monitoring is not significant over the noise level defined, as a function of magnitude, by field stars in our survey (see Venuti et al., in prep.). Additionally, several of these putative members would be located close to or onto the MS-turno GLYPH&lt;11&gt; (or below) in the HR diagram of the region

(Fig. 12). We therefore suspect at least a part of these objects to possibly be misclassified dwarfs instead. A list of objects with questioned membership, of the relevant CFHT photometry and of the original selection criteria is provided in Table A.3.

Another interesting group of objects in Fig. 6 is the sample of WTTS displaced redward of the main cluster sequence. Most of these objects consistently appear significantly reddened through all CFHT photometric diagrams; in several cases, they display little variability. While contamination by field stars might be partly reflected in the observed properties, most of the objects in this stellar group would actually look like very young, or equivalently overluminous, objects on the HR diagram of the region (i.e., they are located on the upper right envelope of the data point distribution in Fig. 12). This would rather suggest binarity and or earlier (i.e., more embedded) evolutionary status as pos-/ sible explications of the distinctive properties of these objects.

## 3.3. NGC 2264 members in the CFHT survey: individual extinction and stellar parameters

The sample of PMS stars retrieved or newly identified in the CFHT survey amounts to about 750 objects. Based on the multiband photometry obtained during our campaign, as well as on the collection of previous literature data in di GLYPH&lt;11&gt; erent spectral bands, we performed an extensive investigation of individual stellar parameters (like extinction, mass and radius); this is aimed at characterizing the properties of members in a uniform and selfconsistent manner.

Extinction towards NGC 2264 is known to be low (cf. Dahm 2008, for a review of literature results), with a typical value of &lt; GLYPH&lt;24&gt; 0.4 mag for the visual extinction A V . Indeed, the good agreement of the cluster sequence with the color locus traced by foreground stars on Fig. 4 and the color saturation at GLYPH&lt;24&gt; 1.4 in g GLYPH&lt;0&gt; r (cf. discussion at the beginning of Sect. 3) is consistent with a low value of the average A V . We use the photometric properties displayed in this diagram to derive an estimate of the individual A V for late-type stars, i.e. the group of members located on the horizontal branch in the color locus. For this group of objects, A V is determined by displacing each point along the reddening direction until the lowermost envelope of the branch distribution defined by field stars is reached and computing the corresponding amount of reddening. This procedure has the advantage of relying solely on observed quantities, without involving the comparison with an external reference sequence that would intrinsically introduce an additional source of uncertainty on the result.

For earlier-type stars ( &lt; GLYPH&lt;24&gt; M0), the cluster sequence is nearly parallel to the reddening direction in Fig. 4, thus preventing a direct A V determination from the object location on the diagram. For this group of objects, we derive a first estimate of the individual A V from the analysis of the infrared photometry in the JHK S filters, as measured in the 2MASS survey (Skrutskie et al. 2006). This consists in comparing their infrared colors with the JHK S color sequence for dwarfs (Covey et al. 2007) and with the CTTS infrared locus defined in Meyer et al. (1997) and converted to the 2MASS photometric system as in Covey et al. (2010) (K. Covey, private communication). Additional A V estimates for a part of these sources are available from previous, similar studies (e.g., Rebull et al. 2002, where A V is investigated from the color excess on R-I, or Cauley et al. 2012, where A V is computed by fitting stellar spectra with spectral templates) and or derived / from our optical photometry collection (from the measurement of the color excess on R-I, V-I or B-V in order of preference). In the absence of a direct A V estimate from CFHT photometry, all available A V derivations have been examined and their con-

sistency checked on the color properties displayed in the CFHT diagrams, in order to discard discrepant estimates; an average value has then been adopted as a final estimate.

In order to evaluate the uncertainty on our individual A V estimates, we collected all available A V determinations for any object in our sample, and measured the amount of scatter observed among values inferred from di GLYPH&lt;11&gt; erent authors methods. / This comparison allowed us to conclude that di GLYPH&lt;11&gt; erent A V estimates are typically consistent within a radius of a few 0.1 mag, which then provides an order-of-magnitude uncertainty on the A V values adopted in this study.

Of the sample of CFHT members, spectral type is known for around 50% and has been retrieved, in order of preference, from the studies of Dahm &amp; Simon (2005), Rebull et al. (2002) or Walker (1956). For the remaining objects, the spectral type has been derived from the comparison of the dereddened colors with the empirical optical color sequence of Covey et al. (2007). We used the dereddened photometry of members with known spectral type, earlier than M1, to recalibrate the empirical spectral type-color sequence onto CFHT photometry. Spectral types were converted to e GLYPH&lt;11&gt; ective temperatures T e f f following the scale of Cohen &amp; Kuhi (1979), checked for consistency against the scale of Luhman et al. (2003).

Bolometric luminosities have been derived from the dereddened J-band photometry retrieved from the 2MASS catalog, adopting T e f f -dependent J-band bolometric corrections (BC ) J obtained as a fit to the T e f f -BC J scales of Pecaut &amp; Mamajek (2013) and Bessell et al. (1998). A cluster distance value of 760 pc (Sung et al. 1997) has been adopted for the conversion to absolute magnitudes; this commonly adopted estimate of distance has been recently strengthened by the results of Gillen et al. (2014), who performed a detailed characterization of a newly identified, well sampled low-mass PMS eclipsing binary in NGC 2264, with strong evidence of membership and an inferred distance of 756 GLYPH&lt;6&gt; 96 pc. For a few tens of objects, for which J-band photometry from the 2MASS survey was not available, L bol has been derived using an empirical calibration relationship between the dereddened r -band photometry and L bol , which has been inferred from the distribution of dereddened r -band magnitudes vs. L bol values derived from the J-band photometry over the whole sample.

Stellar masses have been determined by placing each object on the Hertzsprung-Russell diagram and comparing their position with the mass tracks of the PMS model grid of Siess et al. (2000) (Fig. 12). A track-fitting tool provided by L. Siess on his webpage has been used to interpolate the mass value between two mass tracks. CFHT members span a wide range of masses, varying from 0.1 to 2 M GLYPH&lt;12&gt; . Stellar radii have been determined from the estimates of L bol and T e f f as

$$R _ { * } = \sqrt { \frac { L _ { b o l } } { 4 \pi \sigma T _ { e f f } ^ { 4 } } }$$

where GLYPH&lt;27&gt; is the Stefan-Boltzmann constant.

Photometric data and stellar parameters for the whole population of NGC 2264 members monitored at CFHT are reported in Table 2 and Table 3, respectively.

## 4. UV excess and accretion in NGC 2264

The investigation of YSOs at short wavelengths o GLYPH&lt;11&gt; ers a unique window on the accretion mechanisms: it provides one of the most direct probes to the hot emission from the impact layer

Fig. 12. HR diagram for NGC 2264 members. Mass tracks and isochrones from Siess et al. (2000)'s models are plotted over the point distribution.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-da33d759-4d62-45f0-9a85-a402e07f6738.jpg)

where the accreting material hits the stellar surface at near freefall velocities, producing a shocked area up to several GLYPH&lt;2&gt; 10 3 K hotter than the stellar photosphere. This enhanced luminosity at UV wavelengths is reflected in the characteristic color properties of accretion-dominated objects compared to non-accreting young stars (cf. Figs. 5 and 6) and hence provides a direct proxy to the accretion luminosity, which in turn enables the investigation of the rates of mass accretion in individual systems.

The u -band flux excess is trivially given by

$$\underset { \text{and} \, \text{mag} } { \text{mag} } \underset { \text{and} } { \text{mag} } \underset { \text{and} } { \text{mag} } F _ { u } ^ { e x c } = F _ { u } ^ { o b s } - F _ { u } ^ { p h o t },$$

where F obs u is the measured flux in the u -band and F phot u is the amount of flux that would be expected from a purely photospheric emission. It is important to remark that, in order to isolate the component of additional flux produced in the accretion shock and thus probe accretion onto the star, the flux excess (compared to dwarfs) due to the enhanced chromospheric activity of T Tauri stars has to be included in the photospheric emission. This issue is addressed by taking the photometric properties observed for the WTTS population as the reference for the definition of the UV excess. The color excess will thus correspond to the di GLYPH&lt;11&gt; erence between the observed stellar colors and the expected photospheric colors based on the properties of the non-accreting counterparts of the objects of interest:

$$E ( u - m ) = ( u - m ) _ { o b s } - ( u - m ) _ { p h o t + c h r o m }$$

As discussed in Sect. 3, no significant impact of the accretion luminosity is observed on the stellar flux detected at filters redder than the u -band; hence, the color excess of Eq. 7 basically corresponds to the u -band excess revealing the presence of accretion, E u ( ) ' E u ( GLYPH&lt;0&gt; m ).

We measured the UV color excess of accreting stars in two di GLYPH&lt;11&gt; erent ways, from the properties displayed on the u GLYPH&lt;0&gt; g vs. g GLYPH&lt;0&gt; r and r vs. u GLYPH&lt;0&gt; r diagrams (Fig. 5 and Fig. 6, respectively), as detailed below.

Table 2. Photometry for NGC 2264 members monitored at CFHT.

| Mon-ID 1   |      RA |      Dec |      u |      g |      r |      i |
|------------|---------|----------|--------|--------|--------|--------|
| Mon-000007 | 100.471 |  9.96755 | 15.85  | 14.807 | 14.674 | 14.192 |
| Mon-000008 | 100.452 |  9.90322 | 19.656 | 17.083 | 15.839 | 19.887 |
| Mon-000009 | 100.538 |  9.80134 | 18.208 | 16.043 | 15.011 | 14.532 |
| Mon-000011 | 100.322 |  9.909   | 18.548 | 16.808 | 15.531 | 14.535 |
| Mon-000014 | 100.528 |  9.69214 | 18.895 | 15.879 | 14.28  | 13.445 |
| Mon-000015 | 100.538 |  9.9841  | 19.427 | 16.172 | 14.589 | 13.786 |
| Mon-000017 | 100.383 | 10.0068  | 18.402 | 15.972 | 14.93  | 14.496 |
| Mon-000018 | 100.305 |  9.91909 | 19.311 | 16.526 | 15.366 | 14.613 |
| Mon-000020 | 100.538 |  9.73427 | 19.964 | 17.213 | 15.891 | 15.267 |
| Mon-000021 | 100.248 |  9.99594 | 18.174 | 15.835 | 14.81  | 14.452 |

Notes. A full version of the Table is available in electronic form at the CDS. A portion is shown here for guidance regarding its form and content. Uncertainties on the photometric measurements are estimated to range from GLYPH&lt;24&gt; 0.15 to GLYPH&lt;24&gt; 0.004 from the faintest to the brightest objects in u and g , from GLYPH&lt;24&gt; 0.015 to GLYPH&lt;24&gt; 0.0015 in r and 5-10 times smaller in i .

( 1 ) Object identifiers adopted within the CSI 2264 project (see Cody et al. 2014).

Table 3. Spectral type, extinction and stellar parameters for NGC 2264 members monitored at CFHT.

| Mon-ID     | Status 1   | SpT     | SpT_ref. 2   | A V   |   L bol (L GLYPH<12> ) |   M GLYPH<3> (M GLYPH<12> ) |   R GLYPH<3> (R GLYPH<12> ) | log t(yr)   |
|------------|------------|---------|--------------|-------|------------------------|-----------------------------|-----------------------------|-------------|
| Mon-000007 | c          | K7      | s            |       |                  1.28  |                        0.69 |                        2.36 | 6.09        |
| Mon-000008 | w          | K5      | s            | 0.3   |                  0.241 |                        0.8  |                        0.85 |             |
| Mon-000009 | w          | F5      | p            | 2.3   |                  2.31  |                        1.3  |                        1.24 |             |
| Mon-000011 | c          | K7      | s            | 0.3   |                  0.79  |                        0.7  |                        1.86 | 6.34        |
| Mon-000014 | w          | K7:M0   | p            | 1.2   |                  3.9   |                        0.66 |                        4.22 | 5.68        |
| Mon-000015 | w          | K7:M0   | p            | 1.0   |                  2.7   |                        0.65 |                        3.5  |             |
| Mon-000017 | c          | K5      | s            | 0.2   |                  0.71  |                        1.13 |                        1.45 |             |
| Mon-000018 | w          | K3:K4.5 | p            | 0.6   |                  1.7   |                        1.47 |                        2.03 | 6.60        |
| Mon-000020 | w          | K7      | s            | 0.4   |                  0.44  |                        0.71 |                        1.39 | 6.70        |
| Mon-000021 | c          | K5      | s            | 0.4   |                  0.96  |                        1.2  |                        1.69 | 6.69        |

A full version of the Table is available in electronic form at the CDS. A portion is shown here for guidance regarding its form and content.

Notes. ( 1 ) 'c' = CTTS; 'w' = WTTS; 'cc' = CTTS candidate.

( 2 ) 's' = spectroscopic SpT estimate (retrieved from literature); 'p' = photometric SpT estimate (from CFHT optical colors).

- -In the first case (Fig. 5), we referred the color excess measurement to a straight line nearly parallel to the reddening vector and following the upper part of the WTTS distribution, as shown in Fig. 9. E u ( ) is thus measured as

$$E ( u ) & = ( u - g ) _ { o b s } - ( u - g ) _ { r e f } \,, & ( 8 ) & \quad \text{the } \varsigma$$

where ( u GLYPH&lt;0&gt; g ) re f is the reference (i.e., expected) color at the observed g GLYPH&lt;0&gt; r level, ( u GLYPH&lt;0&gt; g ) re f = 1 312 ( : g GLYPH&lt;0&gt; r ) obs + 1 03. The : WTTS color sequence is in good agreement with the reference line until spectral type GLYPH&lt;24&gt; M0, and hence the adoption of this reference sequence allows us to derive a UV excess estimate essentially una GLYPH&lt;11&gt; ected by reddening. For later-type stars, however, the real WTTS color sequence turns down on the diagram; this implies that the extrapolation of the WTTS trend at spectral types earlier than M0 to the whole sample yields an excess overestimation (i.e., more negative values of E u ( )) for M-type stars. As the bulk of M-type WTTS lies about 0.2 mag below the reference line, we corrected for this bias by adding a constant o GLYPH&lt;11&gt; set of 0.2 mag to the E u ( ) estimates derived for this group of objects.

- -In the second case, as can be observed on Fig. 6, WTTS nicely trace the cluster sequence over the whole range of magnitude. We thus dereddened the u r ; photometry and defined the sequence of reference colors as the least-squares fit polynomial to the resulting WTTS distribution. The UV

excess is thus provided by

$$\underset { \text{ibu-} } { \text{mag} } \quad E ( u ) = ( u - r ) _ { o b s } - ( u - r ) _ { r e f } \,,$$

where ( u GLYPH&lt;0&gt; r ) re f is the reference color at brightness robs . Due to the small number of point and the large scatter a GLYPH&lt;11&gt; ecting the definition of the sequence at the brighter end, we decided to restrict our UV excess analysis to objects fainter than 14.5 in r .

This second method, compared to the first, has the disadvantage of being intrinsically more uncertain, as it is subject to errors in A V determination (while the first method is essentially A V -independent); on the other hand, the application of this method to u - and r -band observations taken during the monitoring survey allows us to derive both a global picture of accretion in NGC 2264 from the average photometry, and a variability range for the UV excess on a timescale of a few weeks (cf. Sect. 4.2).

An average rms error of GLYPH&lt;24&gt; 0.16 mag has been associated to the UV excess determinations, in order to account for the scatter of the WTTS distribution around the reference sequence.

From Eq. 6 and

$$\underset { \text{range of} } { \text{$E(u)=u_{obs}-u_{phot}$} }$$

## we obtain

$$\tilde { \text{The UV} } \ \ F _ { u } ^ { e x c } = F _ { u } ^ { o b s } \left ( 1 - 1 0 ^ { + 0. 4 E ( u ) } \right ) = F _ { u } ^ { 0 } \, 1 0 ^ { - 0. 4 u } \left ( 1 - 1 0 ^ { + 0. 4 E ( u ) } \right ), \ \ ( 1 1 )$$

Article number, page 11 of 24

Fig. 13. UV excess luminosity measurements obtained in this study are compared to Rebull et al. (2002)'s for objects in common to both samples (L. Rebull, private communication). On the x-axis, dots correspond to the median L exc u detected during the 2 week-long CFHT monitoring, while variability bars depict the actual range in L exc u values measured at di GLYPH&lt;11&gt; erent observing epochs during the survey. The equality line is traced as a solid line to guide the eye.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-d0fbb11e-5220-4b6b-aa36-9087c5532327.jpg)

where u is the (dereddened) apparent u -band magnitude of the object, E u ( ) is the UV excess computed as in Eqs. 8-9 and Fu 0 is the SDSS u -band zero-point flux ( Fu 0 = 3767.2 Jy).

We used Eq. 11 (cf. also Rigliaco et al. 2011) to measure the u -band flux excess over the whole sample and then converted F exc u to L exc u using a distance of 760 pc.

The results of a similar, single epoch photometric survey to detect and measure mass accretion in NGC 2264 from the UV excess diagnostics have been reported in Rebull et al. (2002). The authors measured the UV excess displayed by accreting stars from U-band and V-band photometry, by comparing the observed colors with colors expected based on the spectral type of individual objects. In Fig. 13, we compare their measured UVexcess luminosities (L. Rebull, private communication) with both the median L exc u values and the relevant variability ranges we detect here, on a timescale of a couple of weeks, for CTTS targeted in both surveys (about 100). As can be observed, data points result to be well distributed around the equality line; the rms scatter measured about the line, amounting to 0.5 dex, is quite consistent with the average amount of variability on L exc u (about 0.6 dex) detected across the sample from CFHT monitoring.

## 4.1. Accretion luminosity and mass accretion rates

The u -band excess luminosity represents a single component of the total excess emission over the whole spectral range. In order to characterize the properties of the accretion process, it is therefore necessary to also take into account the flux excess at unobserved wavelengths, i.e., to apply a bolometric correction.

Gullbring et al. (1998) showed that a direct correlation exists between the total accretion luminosity L acc and the excess luminosity observed in the U-band. This is indeed expected, since the excess emission likely dominates the observed flux at wavelengths &lt; 4500 Å (i.e., below the g -band filter). In order to achieve this result, they measured the excess emission from spectra covering the range 3200-5400 Å, for a sample of 26 CTTS. They modeled the excess spectrum from the accretion spot assuming a region of constant temperature and density (slab model) and used this model fit to estimate the amount of excess emission at unobserved wavelengths. The authors eventually resolved to adopt a uniform correction factor, across the whole sample, for the ratio of total-to-observed excess flux, corresponding to a T GLYPH&lt;24&gt; 10 000 K. For each object, they subsequently synthetized a U-band magnitude from the observed spectrum and computed the corresponding flux excess. The comparison of the two quantities (L acc and L exc U ) enabled the authors to yield a well-determined calibration equation, log( Lacc = L GLYPH&lt;12&gt; ) = 1 09 : + 0 04 : GLYPH&lt;0&gt; 0 18 : log( L exc U = L GLYPH&lt;12&gt; ) + 0 98 : + 0 02 : GLYPH&lt;0&gt; 0 07 : , that was ultimately consistent with the model predictions of Calvet &amp; Gullbring (1998), who introduced a detailed description of the accretion shock region.

Here we adopt a simple phenomenological approach to derive a direct relationship between the u -band excess luminosity and the total accretion luminosity L acc . In a first approximation, we describe the stellar emission as a blackbody emission at the photospheric temperature T phot . Similarly, the emission from the accretion spot is approximated by a blackbody at T = T spot . The excess emission linked with accretion is thus given by the difference between the flux emitted by the distribution of accretion spots and the flux emitted by a corresponding area, of the same extension, on the stellar photosphere. In terms of blackbodies, this becomes the di GLYPH&lt;11&gt; erence between a blackbody flux density at T = T spot and a blackbody flux density at T = T phot , both integrated over the same area. In this picture, the ratio r of total-to-observed flux excess in the u -band, i.e. the corrective factor to be introduced to derive the total L acc from L exc u , can be estimated as

$$& \text{sure the} \\ & \text{onverted} \quad r = \frac { \sigma T _ { s p o t } ^ { 4 } - \sigma T _ { p h o t } ^ { 4 } } { \pi \int _ { u } ( B _ { T _ { s p o t } } ( \lambda ) - B _ { T _ { p h o t } } ( \lambda ) ) \, d \lambda } \,, & ( 1 2 ) \\ & \text{urvey to} \quad \int _ { u } ( B _ { T _ { s p o t } } ( \lambda ) - B _ { T _ { p h o t } } ( \lambda ) ) \, d \lambda \\ & \text{under the UV}$$

where the integration is performed over the u -band window. No assumptions on the filling factor are needed in this description.

In order to investigate the relationship between L exc u and L acc , we selected a subsample of 44 CTTS, covering the spectral type range M3.5:K2 (or equivalently the mass range GLYPH&lt;24&gt; 0.2-1.8 M GLYPH&lt;12&gt; ) and whose variability, monitored in the u -band ( GLYPH&lt;21&gt; e f f = 3557 Å) and r -band ( GLYPH&lt;21&gt; e f f = 6261 Å) at CFHT, is likely accretion-dominated. An extensive spot modeling (cf. Bouvier et al. 1993) of the simultaneous u + r variability amplitudes has been performed in order to assess the dominating features of the variability displayed by individual objects; accretion-dominated sources have been identified as well-represented variables in terms of a surface spot distribution GLYPH&lt;24&gt; 10 3 to several GLYPH&lt;2&gt; 10 3 Khotter than the stellar photosphere. A full description of the variability analysis for the CFHT NGC 2264 members will be provided in a forthcoming paper (Venuti et al., in prep.).

For each of the selected objects, we adopted the individual T spot estimate inferred from spot models, that corresponds, for a given object, to the color temperature of the spot distribution that best reproduces the simultaneous flux variations observed at di GLYPH&lt;11&gt; erent wavelengths; T phot is derived from the spectral type as described in Sect. 3.3. We then derived, for each object, the ratio r of Eq. 12, through a numerical integration of B( GLYPH&lt;21&gt; ) in the spectral range 3257-3857 Å (cf. Fukugita et al. 1996); an uncertainty

on the value of r has been inferred correspondingly to the uncertainty on the spot temperature T spot . Based on these estimates, the percentage of the total excess flux detectable in the u -band is on average GLYPH&lt;24&gt; 10%. The accretion luminosity L acc is then computed as Lacc = r L exc u .

A least-squares fit to the log GLYPH&lt;0&gt; L exc u = L GLYPH&lt;12&gt; GLYPH&lt;1&gt; vs. log ( Lacc = L GLYPH&lt;12&gt; ) distribution inferred as previously described provided us with the following calibration relationship:

$$\log \left ( \frac { L _ { a c c } } { L _ { \odot } } \right ) = ( 0. 9 7 \pm 0. 0 3 ) \, \log \left ( \frac { L _ { u } ^ { e x c } } { L _ { \odot } } \right ) + ( 1. 0 9 \pm 0. 0 7 ) \quad \ ( 1 3 ) \, \begin{array} { c } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{\quad \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{  \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{ } \text{$$

Fig. 14 compares our result with the calibration inferred from the study of Gullbring et al. (1998); the parameters of the two calibrations are fairly consistent within the error bars. The locus traced by our point distribution is globally consistent with the prediction of Gullbring et al.'s study, albeit showing a systematic shift above their calibration line. Our estimates result to be larger by a factor of . 2 to 4, corresponding to a di GLYPH&lt;11&gt; erence from GLYPH&lt;24&gt; 0.2 dex at the highest ˙ Macc regimes ( GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 7 M yr) to GLYPH&lt;12&gt; / GLYPH&lt;24&gt; 0.6 dex at the lowest ˙ Macc regimes ( GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 10 M yr). GLYPH&lt;12&gt; /

This di GLYPH&lt;11&gt; erence between the two conversions is significantly reduced if we adopt, in our approach, a constant spot temperature of 10 000 K, an assumption closer to the uniform correction, at T GLYPH&lt;24&gt; 10 000 K, that was introduced in Gullbring et al.'s method. Indeed, individual T spot estimates adopted in our study cover a wide range of values, with a mean around GLYPH&lt;24&gt; 6500 K, a spread of GLYPH&lt;24&gt; 2000-3000 K at a given stellar temperature and nominal T spot that can reach up to GLYPH&lt;24&gt; 9000-11 000 K and down to GLYPH&lt;24&gt; 5000 K. As described earlier, spot temperatures, in our approach, are inferred case by case by fitting simultaneous optical ( r -band) and UV ( -band) variability amplitudes with spot models. Hence, u these correspond to the color temperatures that best reproduce

Fig. 14. L acc -L exc u calibration relationship. The black solid line traces the relationship derived in this study, while the orange dashed line represents the relationship by Gullbring et al. (1998). The data points represent the u -band excess luminosity and the accretion luminosity estimates inferred on an individual basis for a sample of 44 CTTS in our survey (see text). The error bars reflect the uncertainty on the individual T spot estimates. The black dashed-dotted line shows the calibration line that would be inferred, with the procedure followed in this study, in the assumption T spot = 10 000 K.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-5378c729-a145-48a3-bfa0-a7cb69250896.jpg)

the observed, GLYPH&lt;21&gt; -dependent contrast of the spot distribution with the stellar reference emission level (which is the e GLYPH&lt;11&gt; ect we aim at accounting for when introducing a bolometric correction). As the peak of the emission of a blackbody in the u -band occurs around 10 000 K, lower T spot estimates naturally imply that a larger correction factor is needed in order to recover the total flux emitted by the spot. While, for a typical photosphere (T GLYPH&lt;24&gt; 4000 K), 12% of the excess emission of a 10 000 K spot is observed in the u -band, this percentage reduces to GLYPH&lt;24&gt; 10% for a 7000 K spot and to less than 8% for a 6000 K spot. In order to evaluate how much of the discrepancy could actually be due to this e GLYPH&lt;11&gt; ect, we repeated our procedure assuming a constant T spot = 10 000 K over the whole subsample of objects; as shown in Fig. 14, this indeed would allow us to mostly cancel the o GLYPH&lt;11&gt; set between the two calibrations.

We are aware of the simplified nature of the blackbody approximation for stellar and spot fluxes; at the same time, we stress that the results inferred from this description are not inconsistent with those obtained from more sophisticated pictures. The discrepancy, albeit systematic, between L acc estimates from Eq. 13 and those from Gullbring et al. (1998)'s calibration is smaller than, or comparable to, the usual estimates of the uncertainty a GLYPH&lt;11&gt; ecting the accretion rate measurements ( GLYPH&lt;24&gt; 0.5 dex), and a non-negligible contribution to this o GLYPH&lt;11&gt; set may actually derive from the di GLYPH&lt;11&gt; erent assumption on the correction factors adopted across the sample of objects considered in each study. Fig. 10 of the theoretical study performed by Calvet &amp; Gullbring (1998), though it provides full endorsement of Gullbring et al.'s relationship, yet depicts the complexity of achieving a detailed physical representation of the accretion shock. The range of predictions from model configurations they explored would lie, on Fig. 14 of our work, across the region delimited by our (Eq. 13) and Gullbring et al. (1998)'s relationship.

Considering this, and the fact that each model has its own limitations, we believe it beneficial to explore whether, and to which extent, information deduced from data might actually be sensitive to di GLYPH&lt;11&gt; erent model assumptions. We therefore adopt the calibration in Eq. 13 to estimate the total L acc from the measured u -band luminosity excess, and will refer at the same time to Gullbring et al.'s, where relevant, for comparison purposes.

If we assume that the accretion energy is reprocessed entirely into the accretion continuum (cf. Herczeg &amp; Hillenbrand 2008), the mass accretion rate ˙ Macc can be derived from L acc as

$$\begin{array} { c } \Big | & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \\ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &  & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & && & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &
 & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & - & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ; & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & } & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & _ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & \ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & + & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & = & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &, & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & ^ & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & b & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & | & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &.. & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & & &$$

where M GLYPH&lt;3&gt; and R GLYPH&lt;3&gt; are the parameters of the central star, R in is the inner disk radius and we are assuming that the accretion funnel originates at the truncation radius of the disk, at R in GLYPH&lt;24&gt; 5 R GLYPH&lt;3&gt; (Gullbring et al. 1998).

Fig. 15 compares the ˙ Macc distributions inferred from each of the two UV derivations described at the beginning of Sect. 4. The diagram reports the ˙ Macc determinations for every accreting object, in the sample, that displays a UV excess according to the identification procedure previously described. The two distributions are on average consistent, within a small mass-dependent o GLYPH&lt;11&gt; set likely reflecting a residual overestimation of the ( u GLYPH&lt;0&gt; g ) color excess from Fig. 9 (see discussion around Eq. 8). Variability is likely to a GLYPH&lt;11&gt; ect, at least partly, the scatter observed within GLYPH&lt;6&gt; 0.25 dex around the fitting line (see Sect. 4.2).

Fig. 16 shows the mass accretion rate distribution inferred for a population of 236 high-confidence NGC 2264 members, as a function of stellar mass, whose range extends down to GLYPH&lt;24&gt; 0.1 M GLYPH&lt;12&gt;

Fig. 15. ˙ Macc values obtained with the two methods described in the text. ˙ Macc values on the y-axis are inferred from u g r snapshot photometry, while ˙ Macc values on the x-axis correspond to the median UV excess detected during the ur monitoring. Colors and sizes of symbols are scaled according to stellar mass M . The solid line marks the equality GLYPH&lt;3&gt; line; the linear regression to the point distribution is traced as a dashed line.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-006acdce-d6a6-43a1-87c8-e7fddfb5a08d.jpg)

and up to GLYPH&lt;24&gt; 1.5 M GLYPH&lt;12&gt; . These objects have been classified as CTTS, i.e., accreting objects, as opposed to WTTS among the PMS population analyzed so far, based on spectroscopic criteria (H GLYPH&lt;11&gt; EW, H GLYPH&lt;11&gt; width at 10% intensity) and or photometric criteria (signifi-/ cant UV excess from the CFHT sample, GLYPH&lt;11&gt; IRAC value consistent with Class II sources from the data of Teixeira et al. 2012, large variability from the CFHT photometry; see Venuti et al., in prep., for this last point). The accretion parameters for these objects are reported in Table 4. The nominal values of the accretion rates are computed from the photometry corresponding to the median UV excess detected during the CFHT MegaCam 2-week long mon-/ itoring and are estimated to be accurate within a factor of GLYPH&lt;24&gt; 3 typically.

A confidence detection threshold for ˙ Macc has been computed by evaluating the average scatter of the WTTS around the reference color sequence and testing the corresponding level of 'noise' on the UV excess - and consequently accretion rate - determination. Namely, the magnitude-dependent rms scatter of WTTS around the reference sequence has been measured to compute a (magnitude-dependent) minimum UV excess level to be detected in order to be confident that accretion is indeed the dominant source of excess. This allowed us to trace a conservative boundary separating the objects with an unambiguous ˙ Macc detection from the objects for which it is not possible, a priori, to evaluate the relevance of spurious contributions (i.e., chromospheric activity) to the apparent UV excess. For the accreting objects for which the actually measured UV excess is less than the corresponding boundary value, only an upper limit at the detection threshold has been reported for ˙ Macc .

Little extensive investigation of mass accretion rates in NGC 2264 has been reported in the literature prior to the survey we are currently performing. Here we refer again to the study of Rebull et al. (2002). Those authors reported ˙ Macc estimates for about 75 CTTS in the region; their accretion rates are derived from the measured UV excess luminosities by adopting Gullbring et al. (1998)'s relationship to obtain the total L acc values. Results of the comparison of our ˙ Macc estimates with theirs, for the few tens of sources in common (that cover nonetheless about 2 orders of magnitude in ˙ Macc ), are quite consistent with what deduced from the comparison of the measured L exc UV (Fig. 13). When using the same L exc u -L acc calibration, ˙ Macc estimates from the two studies distribute around the equality line, with an rms scatter of about 0.5 dex, pretty consistent with the average amount of ˙ Macc variability we detect on week timescales during our monitoring (see Sect. 4.2).

Some important features can be observed on Fig. 16. The ˙ Macc distribution globally traces a mass-dependent trend, showing an average increase of the mass accretion rates with stellar mass. The existence of such a trend has been investigated from observations over the last decade, albeit with a large uncertainty on the quantitative estimate of the actual relationship, which is likely to be dependent, among other factors, upon the probed stellar mass range and the mean age and age spread across the cluster (cf., e.g., Hartmann et al. 2006 and references therein, Rigliaco et al. 2011, Manara et al. 2012). Scattering of objects at a given stellar mass (see discussion later in the Section) also contributes to increasing the uncertainty on the actual relationship. Another important point to this respect is to understand the actual role played by detection biases in the observed mass-˙ Macc trend. As previously discussed, the statistical sample of accretion rates is limited at the lower ˙ Macc end (Fig. 16), due to the impossibility of an accurate detection below a certain threshold; for this group of objects, only an upper limit to the actual ˙ Macc value can be inferred. Contrast e GLYPH&lt;11&gt; ects favor the detection of smaller accretion rates on cooler, i.e. lower-mass stars over higher-mass stars (as can be see in Fig. 16), thus unavoidably introducing a mass dependence. It is thus of interest to assess the degree of disentanglement between the observed mass-˙ Macc dependence and the e GLYPH&lt;11&gt; ects of censoring in the data.

Astatistical approach that allows us to infer a robust measure of correlation in a given sample, in presence of multiple censoring (both ties, here points at the same x-value but with di GLYPH&lt;11&gt; erent y-values, and non-detections, i.e. upper limits), is a generalization of Kendall's GLYPH&lt;28&gt; test for correlation (see Feigelson &amp; Babu 2012; Helsel 2012). The test consists in pairing the data points (measurements and upper limits) in all possible configurations and, for each pair, measuring the slope (positive or negative) of the line connecting the two points. A pair of data points having the same x-value will be counted as an indeterminate relationship, likewise a pair of upper limits; upper limits will uniquely contribute a definite relationship when paired with actual detections at higher y-value. The number of times positive slopes occur is then compared to the number of negative slopes, and their di GLYPH&lt;11&gt; erence weighted against the total number of available pairs. This defines the statistical correlation coe GLYPH&lt;14&gt; cient, GLYPH&lt;28&gt; . ln this picture, upper limits have the e GLYPH&lt;11&gt; ect of lowering the likelihood of a correlation, if this is present in the data, as they introduce a number of indeterminate relationships among the tested pairs. Conversely, no correlation will be found in a pure sample of upper limits, even though these trace a well-defined trend, since the number of pairs with definite relationship will be zero. In the null hypothesis of no correlation, GLYPH&lt;28&gt; is expected to follow a normal distribution centered around zero; by measuring the deviation from zero of the actually found value, in terms of GLYPH&lt;27&gt; , it is thus possible to assess a confidence level for the presence of correlation in a given sample.

Fig. 16. ˙ Macc distribution as a function of stellar mass for the population of NGC 2264 accreting members observed at CFHT. Upper limits are shown for those objects that fall below the estimated confidence threshold for accretion detection (see text). For some objects, no significant ˙ Macc has been detected at the median state of the system, but they do show a significant accretion activity (i.e., above the detection level) at the brightest state; in these cases, an upper limit is placed at the detection threshold, while a cross marks the maximum ˙ Macc actually detected. Orange dots and green triangles mark two subgroups of objects dominated respectively by stochastic accretion bursts and variable extinction from a rotating inner disk warp (Stau GLYPH&lt;11&gt; er et al. 2014).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-52a3f7da-0dad-4f12-afc0-f9a4cb379dcf.jpg)

Table 4. UV excess, u -band excess luminosity, ˙ Macc and ˙ Macc variability estimates for NGC 2264 accreting members.

| Mon-ID     |   ( u GLYPH<0> g ) exc. |   ( u GLYPH<0> r ) exc. 1 | log L exc u / L GLYPH<12> 1   | log ˙ M ( u GLYPH<0> g ) acc   | log ˙ M ( u GLYPH<0> r ) acc 1   | log ˙ M max = ˙ M med 2   | log ˙ M med = ˙ M min 3   |
|------------|-------------------------|---------------------------|-------------------------------|--------------------------------|----------------------------------|---------------------------|---------------------------|
| Mon-000007 |                   -0.16 |                     -0.65 | -1.50                         | -7.17                          | -7.22                            | 0.32                      |                           |
| Mon-000011 |                   -0.97 |                     -1.95 | -0.88                         | -7.53                          | -6.73                            | 0.24                      | 0.19                      |
| Mon-000017 |                    0.03 |                     -0.27 | -2.14                         |                                | -8.24 *                          | 0.02                      |                           |
| Mon-000028 |                   -1.79 |                     -1.4  | -3.38                         | -8.90                          | -9.13                            | 0.40                      | 0.59                      |
| Mon-000040 |                   -0.46 |                     -0.34 | -3.89                         | -9.24                          | -9.37                            | 0.20                      | 0.55                      |
| Mon-000042 |                   -0.17 |                     -0.01 | -5.08                         | -9.34                          | -9.48 *                          |                           |                           |
| Mon-000053 |                   -0.35 |                     -0.11 | -4.31                         | -9.06                          | -9.46                            | 0.33                      |                           |
| Mon-000056 |                    0.04 |                      0.01 |                               |                                | -8.24 *                          |                           |                           |
| Mon-000059 |                   -1.91 |                     -1.11 | -3.32                         | -8.36                          | -8.91                            | 0.26                      | 0.57                      |
| Mon-000063 |                   -0.53 |                     -1.03 | -2.56                         | -8.62                          | -8.21                            | 0.15                      |                           |

Notes. A full version of the Table is available in electronic form at the CDS. A portion is shown here for guidance regarding its form and content. ( 1 ) Value measured at the median luminosity state of the system.

( 2,3 ) Variability range detected around the median ˙ Macc level during the 2 week-long ur monitoring. If no accretion activity is detected above the detection limit for the object at the minimum median luminosity level, but some significant accretion is detected at the brightest state, only the / upper variability bar is reported (i.e., log ˙ Mmax = ˙ Mmed upper ( ) ); if the object is a non-detection in accretion, no variability measurement is reported. ( * ) Upper limit.

The application of this statistical tool to our data allowed us to establish the presence of a correlation ( GLYPH&lt;28&gt; GLYPH&lt;24&gt; 0 28) that appears : to hold across the whole mass range studied here with a confidence of more than 6 GLYPH&lt;27&gt; . The robustness of this result was tested against a similar search for correlation in randomly generated distributions of ˙ Macc in the log range [-10.5 , -6.5] for the same stellar population with the same individual detection thresholds. The GLYPH&lt;28&gt; distribution inferred from 100 iterations is peaked around GLYPH&lt;28&gt; = 0.00 (no correlation) with a GLYPH&lt;27&gt; of 0.04; this suggests that indeed an intrinsic correlation is observed in our data to a significance of &gt; 6 GLYPH&lt;27&gt; . We then inferred a statistical estimate for the slope following the approach of Akritas-Theil-Sen nonparametric regression (see Feigelson &amp; Babu 2012):

- 1. we derived a first guess for the slope from a least-squares fit to the data distribution and explored a range of GLYPH&lt;6&gt; 0.5 around this first estimate with a step of 0.005;
- 2. for each value x , we subtracted the x M trend from the initial GLYPH&lt;3&gt; ˙ Macc distribution and repeated the Kendall's GLYPH&lt;28&gt; computation procedure on the distribution of residues thus inferred.

The best fitting trend is the one that, subtracted to the point distribution, produces GLYPH&lt;28&gt; = 0, while an error bar can be inferred corresponding to the interval of values that produce GLYPH&lt;28&gt; &lt; j GLYPH&lt;27&gt; j . We applied this procedure to both the mass-˙ Macc distribution inferred from the u GLYPH&lt;0&gt; r excess distribution (see discussion around Eq. 9) and the one inferred from the u GLYPH&lt;0&gt; g excess distribution (see discussion around Eq. 8), obtaining respectively:

$$\dot { M } _ { a c c } & \in M _ { * } ^ { 1. 5 \pm 0. 2 } & ( 1 5 ) & \text{negligibl} \\ & \text{cretion } t$$

$$\dot { M } _ { a c c } \in M _ { * } ^ { 1. 3 \pm 0. 2 } & & ( 1 6 ) \quad \text{hand, } m _ { \substack { m } { \shortmid } \in \lambda \end{m}$$

The two trends are consistent within 1 GLYPH&lt;27&gt; ; the smaller exponent in Eq. 16 compared to Eq. 15 may reflect the small, mass-dependent o GLYPH&lt;11&gt; set between the ˙ Macc distributions inferred from the two methods (cf. Fig. 15). A slope of 1.7 GLYPH&lt;6&gt; 0.2, instead of the estimate in Eq. 15, would be inferred if we converted u GLYPH&lt;0&gt; r excesses to L acc values following Gullbring et al. (1998). The two estimates are still consistent within 1 GLYPH&lt;27&gt; ; this comparison illustrates that a (small) range in values around the true slope may be induced by external factors such as the L exc u -L acc conversion.

The analysis of correlation described above refers to the homogeneous sample of Fig. 16 as a whole. While the relationship we infer seems to be robust across the whole sample, this likely provides a general, average information on the global trend of ˙ Macc with mass, not a punctual description of the trend. We notice, for instance, that little mass dependence is observed in the upper envelope of the ˙ Macc distribution above M GLYPH&lt;3&gt; GLYPH&lt;24&gt; 0.25 M GLYPH&lt;12&gt; . We attempted a more detailed investigation of the observed ˙ Macc vs. M relationship by probing separately two di GLYPH&lt;3&gt; GLYPH&lt;11&gt; erent, similarly populated mass bins, M GLYPH&lt;3&gt; GLYPH&lt;20&gt; 0.3 M GLYPH&lt;12&gt; and M GLYPH&lt;3&gt; &gt; 0.3 M GLYPH&lt;12&gt; . In both subgroups, we found some evidence of correlation, albeit with different values and significance levels: a correlation with a slope of 3.3 GLYPH&lt;6&gt; 0.7 was inferred, to the 4 GLYPH&lt;27&gt; level, in the lower-mass group, while a correlation with a slope of 0.9 GLYPH&lt;6&gt; 0.4 was obtained for the higher-mass group with a significance of 2.4 GLYPH&lt;27&gt; . These results might suggest that the dependence of ˙ Macc on M GLYPH&lt;3&gt; is not uniform across the whole M - ˙ GLYPH&lt;3&gt; Macc range, but shallower at higher masses and higher accretion rates; on the other hand, a closer investigation of the M GLYPH&lt;3&gt; ˙ Macc relationship from separate, small subsamples is necessarily more uncertain than the assessment of a general correlation trend observed across the whole sample.

Another remarkable feature on Fig. 16 is the large dispersion in the ˙ Macc values observed at each stellar mass, covering a range of GLYPH&lt;24&gt; 2 to 3 dex. Similar results have been inferred from previous statistical surveys of accretion in a given cluster (see, e.g., the recent studies of Rigliaco et al. 2011 on GLYPH&lt;27&gt; Ori and of Manara et al. 2012 on the ONC). This implies that only the average behavior of the actual trend of ˙ Macc with M GLYPH&lt;3&gt; can be deduced, and suggests that di GLYPH&lt;11&gt; erent accretion regimes may co-exist within the same young stellar population. To this respect, we highlight in Fig. 16 the accretion properties of two, morphologically well-distinct, classes of YSOs in NGC 2264, both actively accreting. The first corresponds to the newly identified (Stau GLYPH&lt;11&gt; er et al. 2014) class of objects whose light curves are dominated by accretion bursts. The second class includes objects whose light curves show periodic or quasi-periodic flux dips likely resulting from a rotating inner disk warp partly occulting the stellar emission (Bouvier et al. 2007a; Alencar et al. 2010). As can be seen in Fig. 16, both classes of objects show a level of accretion on average larger than the mean level of accretion activity detected across the population of NGC 2264; however, the two classes are well distinguished, with the first showing a mean accretion rate about three times larger than the second.

## 4.2. Mid-term ˙ Macc variability

The scatter in Fig. 16 could in principle partly arise from the huge photometric variability exhibited by T Tauri stars, both on short-term (hours) and mid-term (days) timescales (e.g., Ménard &amp; Bertout 1999). Short-term variability can a GLYPH&lt;11&gt; ect the ˙ Macc estimates inferred for single objects, but its contribution is likely negligible for our purposes (i.e., for a statistical mapping of accretion throughout the population of the cluster). On the other hand, mid-term variability, comprising the geometric e GLYPH&lt;11&gt; ects linked with stellar rotation and the intrinsic accretion variability (timescales relevant for hot spots), might e GLYPH&lt;11&gt; ectively blur the ˙ Macc distribution inferred from the mapping of the region, since each object is observed at an arbitrary phase.

In order to assess the contribution of mid-term variability to the ˙ Macc spread, we probed the variability of the detected ˙ Macc , over timescales relevant to stellar rotation, by measuring the UV excess and correspondingly computing the mass accretion rate from each observing epoch during the CFHT u -band and r -band monitoring. On average, 16-17 points distributed over the 2-week long survey have been considered to probe variability. This procedure allowed us to robustly associate a range of variability to the individual average ˙ Macc , inferred for each object as described in Sect. 4.1.

Fig. 17 shows the results inferred from this analysis; each variability bar encompasses all values of ˙ Macc detected for the corresponding object during the monitoring (i.e., it traces the di GLYPH&lt;11&gt; erence between the highest and the lowest detected ˙ Macc ). Objects for which only an upper limit to the average ˙ Macc had been derived have not been considered for this step of the analysis; for objects whose nominal ˙ Macc value would fall below the detection threshold (Sect. 4.1) at certain epochs, an upper limit to the lowest apparent ˙ Macc has been drawn.

As can be observed on Fig. 17, variability properties are not uniform across the sample. No clear trends emerge with stellar mass or mass accretion rates; the geometry of individual systems might have a more important contribution in determining the apparent ˙ Macc variability on week timescales, as suggested by the very few points on the diagram, maybe observed in a pole-on configuration, that show a negligible variation in the measured mass accretion rate. The (a)symmetry of individual bars, albeit potentially a GLYPH&lt;11&gt; ected by the specific temporal sampling, may reflect the way accretion proceeds onto the central object, as well

Fig. 17. Variability of the mass accretion rates observed for NGC 2264 members over the 2-week long CFHT monitoring. Black dots correspond to the median ˙ Macc measured for each object (see Fig. 16). For each object, the relevant variability bar is displayed, corresponding to the di GLYPH&lt;11&gt; erence between the maximum and the minimum ˙ Macc detected during the monitoring. A typical variability bar is shown in red.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-6072b8ce-8578-46c5-88a5-be8a6ae163f4.jpg)

as the dominant features of the system: a strongly asymmetric bar with the median ˙ Macc point close to the lower end can identify an object showing sudden, short-lived accretion bursts, while the opposite case well describes a system whose typical luminosity state is interspersed by short duration flux dips. These properties are indeed observed for some of the objects marked in Fig. 16.

The detected amplitudes of variability in ˙ Macc range from small fractions to GLYPH&lt;24&gt; 1 order of magnitude. In order to better visualize the average e GLYPH&lt;11&gt; ects, we subdivided the sample of objects in Fig. 17 in three broad mass intervals (M GLYPH&lt;3&gt; &lt; 0 4 M : GLYPH&lt;12&gt; ; 0 4 M : GLYPH&lt;12&gt; GLYPH&lt;20&gt; M GLYPH&lt;3&gt; GLYPH&lt;20&gt; 1 M ; M GLYPH&lt;12&gt; GLYPH&lt;3&gt; &gt; 1 M ) and for each of them comGLYPH&lt;12&gt; puted the mean ˙ Macc and amplitude of variability. These results are reported in Table 5, while a typical variability bar is represented in red in the lower right part of Fig. 17. On average, a variability of GLYPH&lt;24&gt; 0.5 dex in ˙ Macc is observed over a couple of weeks (which correspond to GLYPH&lt;24&gt; 2-3 rotational cycles for most of the objects; see A GLYPH&lt;11&gt; er et al. 2013), against a typical ˙ Macc spread of GLYPH&lt;24&gt; 2 dex. This comparison clearly shows that the observed ˙ Macc spread at each mass is significantly larger than what can be accounted for by variability on week timescales, be it intrinsic (non-steady accretion) or geometric (rotational modulation).

## 4.3. Accretion variability: intrinsic vs. modulated

As discussed in Sect. 4.2 (see also Sect. 5.2), an average variability of GLYPH&lt;24&gt; 0.5 dex is derived when monitoring the observed ˙ Macc for individual objects over days-to-weeks timescales. This amount of variability is expected to incorporate both the intrinsic variability associated with the accretion process and the geometric e GLYPH&lt;11&gt; ect of varying accretion shock projection along the line of sight as the star rotates.

An illustration of this is shown in Fig. 18. The underlying

Fig. 18. Phased diagram of the variability of the UV excess monitored over a baseline varying from a few hours to a few rotational cycles for an NGC 2264 accreting member (M GLYPH&lt;3&gt; = 1.13 M ; GLYPH&lt;12&gt; ˙ Macc = 5.98 GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 9 M yr; GLYPH&lt;12&gt; / P rot = 4.767 days, from A GLYPH&lt;11&gt; er et al. 2013). The phase is computed starting from the epoch of the first measurement; the UV excess at each epoch is computed as described in Eq. 9. Di GLYPH&lt;11&gt; erent symbols correspond to di GLYPH&lt;11&gt; erent rotational cycles.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-f335c99e-14d6-4c60-9bd4-0761b789eeff.jpg)

sinusoidal profile describes the geometric modulation of the accretion spot linked with stellar rotation, which is reflected in the measured UV excess and, consequently, in the value of ˙ Macc inferred at a given phase. Two additional e GLYPH&lt;11&gt; ects contribute to the variability picture: the intrinsic variability of accretion on days timescales, reflected in the variation of the UV excess measured at the same phase in two di GLYPH&lt;11&gt; erent cycles, and the intrinsic variability of accretion on hours timescales, reflected in the rms variation of the UV excess measured within the same rotational cycle around a given phase (i.e., measurements obtained a few hours apart during the same night).

In order to probe the relative significance of these di GLYPH&lt;11&gt; erent contributions to the monitored ˙ Macc variability and characterize the intrinsic accretion variability over hours-to-days timescales, we selected a subsample of 35 accreting members of NGC 2264, subject to the following conditions:

- 1. being an actual detection in ˙ Macc (i.e., not an upper limit);
- 2. displaying a detectable periodicity, known from the study of A er et al. (2013); GLYPH&lt;11&gt;
- 3. having rotational period small enough to be able to detect at least a partial superposition of di GLYPH&lt;11&gt; erent cycles in the phased diagram.

A er et al. (2013)'s study of rotation in NGC 2264 is based on GLYPH&lt;11&gt; a 23 day-long optical monitoring of the star-forming region per-

Table 5. Average variability in ˙ Macc displayed, on week timescales, by di GLYPH&lt;11&gt; erent, mass-sorted, groups of accreting members of NGC 2264.

| M GLYPH<3> group (M GLYPH<12> ) :         | < 0.4             | 0.4-1            | > 1              |
|-------------------------------------------|-------------------|------------------|------------------|
| < M GLYPH<3> > (M GLYPH<12> )             | 0.26              | 0.60             | 1.20             |
| < log( ˙ M med acc ) > (M GLYPH<12> / yr) | -8.36             | -7.92            | -7.58            |
| < GLYPH<1> log ˙ M acc 1 > (dex)          | 0.50              | 0.46             | 0.52             |
| < log bar + bar GLYPH<0> 2 >              | -0.1 GLYPH<6> 0.3 | 0.0 GLYPH<6> 0.3 | 0.1 GLYPH<6> 0.2 |

Notes. Estimates only take into account objects for which a detection, and not an upper limit, is available for the lowest ˙ Macc phase. ( 1 ) log ˙ M max GLYPH&lt;0&gt; log ˙ M min acc acc ( 2 ) Index probing the symmetry of the variability bars around ˙ M med acc : bar + = log ˙ M max acc GLYPH&lt;0&gt; log ˙ M med acc ; bar GLYPH&lt;0&gt; = log ˙ M med acc GLYPH&lt;0&gt; log ˙ M min acc .

formed with the CoRoT satellite in March 2008. Accreting stars with detected rotational modulation amount to about 16% of the CTTS in our sample. These are distributed over the whole ˙ Macc range at all masses in Fig. 16, except at the lowest ˙ Macc regimes ( . 5 GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 10 M yr), typically populated by objects fainter than GLYPH&lt;12&gt; / 17 in r and hence not attainable with CoRoT.

For each of these objects selected as described above, we phased the light curve of the UV excess and measured independently the average rms variation of detections obtained within the same observing night and the average rms variation between the average UV excess measured in a given phase bin of 0.05 at di GLYPH&lt;11&gt; erent rotational cycles.

The average peak-to-peak variation in the color excess throughout the sample amounts to GLYPH&lt;24&gt; 0.47 mag, which corresponds to a variation of GLYPH&lt;24&gt; 0.52 dex in ˙ Macc , consistent with what we found in Sect. 4.2. 50% of the objects show an average intrinsic rms variation of &lt; 0.1 mag at a given phase from cycle to cycle (80% at GLYPH&lt;20&gt; 0.15 mag and 90% at GLYPH&lt;20&gt; 0.25 mag). Hence, on average a ratio of GLYPH&lt;24&gt; 0.25 between the intrinsic rms variation on day timescales and the peak-to-peak variation can be inferred. The average rms variation on hour timescales is smaller, amounting to GLYPH&lt;20&gt; 0.05 mag. On the assumption that the two contributions (accretion variability on hours vs. days timescales) are mutually independent, we conclude that these two sources of variability together contribute on average just GLYPH&lt;24&gt; 28% of the observed peakto-peak variation in the UV excess. For a typical GLYPH&lt;1&gt; ˙ Macc of 0.52 dex, intrinsic accretion variability is thus expected to contribute for about 0.15 dex, while a major contribution arises from the geometric modulation of the accretion shock.

Interestingly, in the picture of modulation-dominated variability and assuming a homogeneous phase sampling, we would expect to statistically observe symmetric variability bars around the average value of ˙ Macc ; this is indeed consistent with the results reported in Sect. 4.2 and synthetized in Table 5.

## 5. Discussion

Accretion is a key ingredient in early stellar evolution: it regulates the star-disk interaction over the first few million years, with a major impact on the dynamical evolution of both the central object and the protoplanetary disk. Inferring a complete picture of accretion within large, well-characterized samples of young stars allows us to explore the accretion properties over a variety of parameters (e.g., stellar mass, age, disk properties), and hence represents a crucial issue for understanding how accretion evolves on an individual vs. statistical basis.

## 5.1. The M GLYPH&lt;3&gt; GLYPH&lt;0&gt; ˙ Macc relationship

Understanding what relationship links disk accretion rates and the mass of central objects has been a major focus for studies of accretion in YSOs over the last years. A prime motivation for this lies in that the M GLYPH&lt;3&gt; - ˙ Macc relationship likely reflects the actual mechanisms governing the accretion process.

Hartmann et al. (2006) reviewed a number of di GLYPH&lt;11&gt; erent mechanisms which may be relevant to disk accretion and its variability and to the M GLYPH&lt;3&gt; ˙ Macc relationship. As noted by those authors, several unrelated mechanisms may contribute to establishing this relationship, and those mechanisms may act di GLYPH&lt;11&gt; erently in di GLYPH&lt;11&gt; erent mass regimes. A number of previous studies have attempted to better understand the mass dependence of the M GLYPH&lt;3&gt; ˙ Macc relationship (Vorobyov &amp; Basu 2008; Rigliaco et al. 2011; Fang et al. 2013); these analyses suggest a diminution of the steepness of the ˙ Macc -M dependence as M GLYPH&lt;3&gt; GLYPH&lt;3&gt; increases from a few tenths or hundredths of solar mass to a few solar masses. This is consistent with what we observe here for NGC 2264; on the other hand, the overall trend observed in the subsolar-to-solar mass regimes appears to be significantly shallower than that observed in the intermediate (1-10 M GLYPH&lt;12&gt; ) mass regime (Ercolano et al. 2014).

In this study, we examined the distribution of mass accretion rates over mass for a homogeneous population of 236 confirmed, accreting members in NGC 2264, covering the mass range GLYPH&lt;24&gt; 0.1-1.5 M GLYPH&lt;12&gt; . We inferred the presence of a correlation between ˙ Macc and M GLYPH&lt;3&gt; and further showed that this result is not significantly a GLYPH&lt;11&gt; ected by a mass-dependence in the ˙ Macc detection limits. The same fraction of accreting stars over the total young population ( GLYPH&lt;24&gt; 40%) is roughly observed at lower ( &lt; 0.4 M GLYPH&lt;12&gt; ) and higher ( &gt; 1 M ) masses; unless the true fraction GLYPH&lt;12&gt; of accreting objects is considerably larger at higher masses than at lower masses, this attests that no significant selection biases towards lower-mass objects a GLYPH&lt;11&gt; ect our accretion survey. Hence, we obtain ˙ Macc / M , with GLYPH&lt;11&gt; GLYPH&lt;3&gt; GLYPH&lt;11&gt; GLYPH&lt;24&gt; 1.4 GLYPH&lt;6&gt; 0.3 (Eqs. 15-16). This exponent is smaller than the GLYPH&lt;11&gt; = 2 dependence inferred by Muzerolle et al. (2003) over the mass range 0.04-1 M GLYPH&lt;12&gt; ; the trend they obtained, however, is strongly influenced by the substellar component of the PMS population, which is not probed in our study. Indeed, results similar to Muzerolle et al.'s findings have been subsequently inferred from studies probing a comparable mass range in the substellar down to planetary mass regimes (e.g., Mohanty et al. 2005 for a composite sample of members from di GLYPH&lt;11&gt; erent star-forming regions, Natta et al. 2006 in GLYPH&lt;26&gt; Ohiuchi, Herczeg &amp; Hillenbrand 2008 in Taurus, Alcalá et al. 2014 in Lupus). On the other end, studies addressing more specifically the T Tauri mass regime have recently proposed somewhat lower estimates for the quantitative dependence of ˙ Macc on M . Rigliaco et al. (2011) GLYPH&lt;3&gt; probed a mass range, in GLYPH&lt;27&gt; Ori, similar to the one investigated in Muzerolle et al.'s study (0.06-1 M GLYPH&lt;12&gt; ) and inferred a value of GLYPH&lt;11&gt; = 1 6 : GLYPH&lt;6&gt; 0 4. Barentsen et al. (2011) studied a population of 158 : accreting stars, covering the mass range 0.2-2 M GLYPH&lt;12&gt; , in IC 1396, and inferred a M - ˙ GLYPH&lt;3&gt; Macc relationship with GLYPH&lt;11&gt; = 1 1 : GLYPH&lt;6&gt; 0 2. : All these results are quite consistent with the estimate of slope here inferred. The range of GLYPH&lt;11&gt; from these studies may in part reflect the inhomogeneous nature of the methodologies and stellar samples; however, taken at face value, the results suggest a possible dependence of GLYPH&lt;11&gt; on mass, with the higher mass regime showing a shallower slope.

## 5.2. Accretion variability: what contribution to the ˙ Macc spread at a given mass?

A ubiquitous result of the studies of ˙ Macc distributions for large samples of young stars is the significant spread observed in ˙ Macc at a given mass, enveloping the average trend over 2 to 3 orders of magnitude. CTTS are known to display a strong variability over di GLYPH&lt;11&gt; erent timescales; this takes contribution both from geometric e GLYPH&lt;11&gt; ects linked with the rotational modulation of surface starspots and disk features (Herbst et al. 1994) and from the intrinsic variability of the accretion process. Variability could thus be partly responsible for the observed scatter.

In order to probe this e GLYPH&lt;11&gt; ect, we monitored the UV excess and ˙ Macc variability for about 150 NGC 2264 accreting members over two full weeks, with typically 16-17 epochs for each object; this allowed us to show that variability on week timescales determines on average a variation of the detected ˙ Macc over GLYPH&lt;24&gt; 0.5 dex, significantly smaller than the extent of the average ˙ Macc scatter.

This study complements the results inferred from recent studies that explored the variability in mass accretion rates

on longer timescales. Nguyen et al. (2009) analyzed the ˙ Macc variability of a sample of 40 members in Taurus-Auriga and Chamaeleon I, from multi-epoch optical spectra distributed over GLYPH&lt;24&gt; 10 months, with typically 4 epochs per object on a baseline varying from hours to months. They probed the accretion properties from H GLYPH&lt;11&gt; 10% width and Ca II flux and inferred a median variability estimate amounting respectively to 0.65 dex and 0.35 dex. Fang et al. (2013) probed the ˙ Macc variability of several tens of members in L1641, monitored on timescales of 1, 10 and 22 months; measuring the amount of accretion from the full width of the H GLYPH&lt;11&gt; line at 10% intensity, the authors inferred a typical ˙ Macc variability of GLYPH&lt;24&gt; 0.6 dex.

A recent spectroscopic survey of the ˙ Macc variability in T Tauri and Herbig Ae stars, addressing shorter ( . hour to days) timescales, has been performed by Costigan et al. (2014). The authors probed accretion properties and their variability from H GLYPH&lt;11&gt; EW measurements for a sample of 14 objects; targets have been monitored over . 1 hour blocks, in some cases iterated within a single night, and repeated over a few nights. Combining the results of this short-term variability monitoring with multiyear information for a few objects, as well as with results from di GLYPH&lt;11&gt; erent surveys, the authors found that the dominant contribution to the observed variability in ˙ Macc arises from the days timescale, with a spread in accretion rates ranging from 0.04 to 0.4 dex.

All these results robustly show that variability on day-tomonth timescales does not impact significantly on the large spread observed in ˙ Macc around the average ˙ Macc -M GLYPH&lt;3&gt; trend; hence, the spread is real and possibly linked to the di GLYPH&lt;11&gt; erent evolutionary stages of individual objects (see Sect. 5.3). Additionally, a strong overlap in the age distributions for the populations of accreting (CTTS) and non-accreting (WTTS) members seems to hold; modulo the issue of an accurate age determination for young stars, this suggests that an intrinsically large range in ˙ Macc exists at a given mass and age, which may reflect various other parameters (like, e.g., the properties of the stellar magnetic field or the initial mass of the disk).

Remarkably, the order-of-magnitude estimates for the ˙ Macc variability inferred on day-to-week timescales (this study; Costigan et al. 2014) and on month-to-year timescales (Nguyen et al. 2009; Fang et al. 2013) are quite consistent; this suggests that the same mechanisms dominate variability over these baselines, with a major contribution arising from the mid-term variability component. This picture is supported by the results of the comparison, shown in Fig. 13, between the L exc UV measurements obtained from ours and Rebull et al. (2002)'s survey of NGC 2264, respectively; what we observe in this case is that, indeed, UV excess luminosities measured at the two epochs are globally consistent, with an rms scatter, about the equality line, of the same order of magnitude as the average amount of variability detected during our 2 week-long CFHT monitoring.

## 5.3. What is the origin of the intra-cluster spread in ˙ Macc ?

A star-to-star variety in the evolutionary stage and evolutionary scenario may provide an important contribution to the interpretation of the large range of accretion properties observed within the region. The cluster shows a hierarchical structure (see Dahm 2008 and references therein), with regions of active, current star formation embedded in a more widespread population of somewhat older stars (Sung et al. 2009). Several episodes of star formation may have occurred, as proposed by Adams et al. (1983). Despite the uncertainty on a quantitative estimate, there is a general consensus about an age spread of up to several Myr within

Fig. 19. Evolution of the accretion properties as a function of stellar age for NGC 2264 accreting members. Three di GLYPH&lt;11&gt; erent mass groups are compared: M GLYPH&lt;3&gt; &lt; 0 4 : M GLYPH&lt;12&gt; (filled circles); 0 4 : M GLYPH&lt;12&gt; GLYPH&lt;20&gt; M GLYPH&lt;3&gt; GLYPH&lt;20&gt; 1 M GLYPH&lt;12&gt; (open squares); M GLYPH&lt;3&gt; &gt; 1 M GLYPH&lt;12&gt; (crosses).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-71b77d40-d590-4e3c-aa69-69aee2dc3a09.jpg)

the population of NGC 2264. It is thus of interest to attempt a characterization of the accretion properties relative to age indicators and disk evolutionary stages.

Looking at how points distribute on the HR diagram of the cluster (Fig. 12) and comparing this to model isochrones provides a debated indication on nominal ages. In order to probe the evolution of accretion properties with stellar age, we inferred an age estimate for each source from isochrone fitting, adopting the models of Siess et al. (2000). The results for accreting sources are displayed vs. the ˙ Macc estimates in Fig. 19. Only objects which are actual detections in ˙ Macc and whose age estimate from models is between GLYPH&lt;24&gt; 10 6 yr and 10 7 yr have been retained for the comparison.

An overall decreasing trend of the average ˙ Macc with increasing stellar ages is traced in Fig. 19, consistently with previous results regarding the average evolution of accretion with stellar age (e.g. Hartmann et al. 1998; Sicilia-Aguilar et al. 2010; Manara et al. 2012). The anticorrelation trend is well described in the M GLYPH&lt;3&gt; &lt; 1 M GLYPH&lt;12&gt; regime, reaching a significance of over 5 GLYPH&lt;27&gt; according to a Spearman's rank order test (Press et al. 1992). A linear regression to this group of objects yields

$$\dot { M } _ { a c c } \infty \, t ^ { - \eta }, \, \eta \simeq 1. 5 1,$$

quite consistent with the GLYPH&lt;17&gt; = 1 5 estimate inferred by Hartmann : et al. (1998). Little information can be achieved on the accretion properties of the higher-mass group (M GLYPH&lt;3&gt; &gt; 1 M ) from our UV GLYPH&lt;12&gt; excess analysis.

The actual ˙ Macc -age relationship may vary depending on the mass range probed, as suggested by the looser correlation that can be inferred from the least massive group of objects (M GLYPH&lt;3&gt; &lt; 0.4 M GLYPH&lt;12&gt; ) compared to the 0.4 - 1 M GLYPH&lt;12&gt; group. This might reflect truly di GLYPH&lt;11&gt; erent evolutionary timescales, but also be affected by uncertainties on the age estimates inferred from models. Indeed, a clear-cut decreasing trend is observed on the lowerleft envelope of the data point distribution in Fig. 19, with the

youngest group being among the most luminous sources in the sample. While limits on the detection of small accretion rates in luminous stars might somewhat a GLYPH&lt;11&gt; ect the observed trend, a more significant bias here is likely introduced by the model age estimates. We notice that lower-mass stars appear to be systematically older than higher-mass stars, based on model ages. If this is not a physical e GLYPH&lt;11&gt; ect, but rather an artifact driven by systematics in the model, the actual anticorrelation trend that can be derived from Fig. 19 is, at best, doubtful.

The reliability of such time evolution analyses of accretion is subject to the assumption that the interpretation of the luminosity spread for a given PMS population on the HR diagram in terms of a true age spread is overall correct. However, this interpretation is subject to a severe controversy. Indeed, recent studies have questioned the accuracy of individual age estimates from the stellar luminosity. Bara GLYPH&lt;11&gt; e et al. (2009) claimed that the luminosity spread observed on the HR diagram of a given starforming region does not correspond to an actual age spread, but merely reflects the individual episodic accretion history, whose dynamics may mimic an age spread of up to a few Myr for coeval stars. While no consensus has obviously been achieved on this issue, an intrinsic age spread of the order of Myr, if present, is likely to be reflected in di GLYPH&lt;11&gt; erent properties of the region, such as the spatial distribution of members. This point is discussed in Sect. 5.4.

As mentioned in Sect. 3.1, information on the circumstellar disk morphology, hence on the evolutionary status, is conveyed by the mid-IR slope of the SED of the system. In Fig. 20, the accretion properties for three di GLYPH&lt;11&gt; erent mass groups of members are compared to the GLYPH&lt;11&gt; IRAC classification from Teixeira et al. (2012). Smaller values of GLYPH&lt;11&gt; IRAC correspond to more evolved disks; the average accretion rates decrease as disks become more evolved, albeit within a significant spread for each mass group.

These comparisons show that several parameters, beyond stellar mass, regulate the accretion evolution of individual systems and eventually contribute to shape the complex picture of a star-forming region such as we may infer from an extensive accretion survey.

## 5.4. A spatial mapping of accretion properties in NGC 2264

Characterizing the spatial distribution of accreting members according to their accretion properties is a useful tool to investigate the evolution dynamics of the region.

Fig. 21 shows a spatial mapping of accretion in NGC 2264. Most accretors are projected onto the cluster; a part of the population of accreting or mildly accreting members are instead located on the peripheral regions. Some interesting features can be observed: the strongest accretors tend to be concentrated around the two known active sites of star formation within the cloud, a few parsecs apart, while the population of more moderate accretors tends to be more evenly distributed throughout the whole cluster. The halo of members is spread over a few parsecs around the star-forming sites.

This di GLYPH&lt;11&gt; erential distribution of mass accretion rates may be the product of a sequential star formation process: the youngest stars, likely the most actively accreting, are indeed expected to be found close to their birth sites, while dynamical evolution over a few Myr (the average age of the cluster) can explain displacements up to a few parsecs relative to the original location where stars were formed. If this is the case, the distribution of ˙ Macc also reflects the presence of an intrinsic age spread within the region, hence corroborating the interpretation in terms of a time evolution over di GLYPH&lt;11&gt; erent ages.

Fig. 20. ˙ Macc -GLYPH&lt;11&gt; IRAC relationship. The GLYPH&lt;11&gt; IRAC -based disk morphology classification follows Table 2 of Teixeira et al. (2012). Three di GLYPH&lt;11&gt; erent mass groups are compared: M GLYPH&lt;3&gt; &lt; 0 4 : M GLYPH&lt;12&gt; (filled circles); 0 4 : M GLYPH&lt;12&gt; GLYPH&lt;20&gt; M GLYPH&lt;3&gt; GLYPH&lt;20&gt; 1 M GLYPH&lt;12&gt; (open squares); M GLYPH&lt;3&gt; &gt; 1 M GLYPH&lt;12&gt; (crosses).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-91a5f76b-7f57-4b03-964f-b6ffdacacfe0.jpg)

We attempted to pursue further this scenario by investigating the spatial distributions of distinct, equal-sized groups of objects at the two age extremities in Fig. 19. If our interpretation is correct, we would expect the presumed younger member group to be spatially more compact than the older one. A statistical method to probe spatial compactness is the minimum spanning tree (MST) method (Prim 1957). This (see also Allison et al. 2009; Parker et al. 2011) consists in computing the minimum path-length required to connect all points of a given sample in a network of direct point-to-point links with no closed loops. Shorter MST lengths will thus identify, to a certain degree of confidence, more compact populations. The significance of the result can be tested by randomly selecting two groups of objects of the same size and belonging to the same population as the two groups of objects under study, and computing the MST of the first random group, the MST of the second random group and the ratio between the two. The iteration of this procedure allows us to reconstruct the normal distribution of MST1 / MST2 values when no significant di GLYPH&lt;11&gt; erence exists between the spatial distributions of two given groups of objects selected among the population of interest.

For our population of accreting objects, sampled in random groups of 20 to 30 elements, we derived a normal distribution of MST ratios centered on 1.0 (homogeneous spatial distribution) with a GLYPH&lt;27&gt; of 0.2. Little evidence for an actual di GLYPH&lt;11&gt; erence in the spatial distributions of 'younger' vs. 'older' members in Fig. 19 can be inferred when comparing the result of their MST test to this distribution. Somewhat smaller values than the mean of 1.0 are obtained from the MST y oun g er / MST older ratio, oscillating around a significance of 1.5 GLYPH&lt;27&gt; that this di GLYPH&lt;11&gt; erence is truly meaningful. This might suggest that the younger population is to some extent more compact than the older population; however, no conclusive evidence can be drawn from this analysis, and the small statistics at lower and higher ages in Fig. 19 yields some dependence of the specific result on the specific selection

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0432v1-figure-575cad55-2772-47e0-9508-4a6aa4735fb2.jpg)

a complete picture of accretion in the cluster and monitored its variability over a few rotational cycles. A broad interval of ˙ Macc values has been inferred within the region, ranging from GLYPH&lt;24&gt; 10 GLYPH&lt;0&gt; 7 to a few GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 10 M yr; at any given stellar mass M , a large GLYPH&lt;12&gt; / GLYPH&lt;3&gt; spread in ˙ Macc values, typically extending over 2 orders of magnitude, is observed.

The average ˙ Macc shows a statistically significant correlation with M ; more uncertain is the functional form of this depenGLYPH&lt;3&gt; dence, which in turn appears to be mass-dependent and likely time-evolving, as highlighted by the results inferred from previous studies for various young stellar populations. A statistical slope of 1.4 GLYPH&lt;6&gt; 0.3 for the log ˙ Macc over log M GLYPH&lt;3&gt; relationship is here derived for NGC 2264 in the mass sample probed; while the absence of the lower-mass ( &lt; 0.1 M GLYPH&lt;12&gt; ) stellar component may contribute to explain this less steep dependence than the value of GLYPH&lt;24&gt; 2 proposed in several studies, the large spread of ˙ Macc values observed at a given mass attests the mere indicative value of such a relationship, which weighs the contributions of several trends as well as of undetected regimes.

-10.0

Fig. 21. Spatial mapping of accretion properties in NGC 2264. Small grey dots indicate the distribution of field stars, large dots indicate accreting cluster members. Only objects with an actual ˙ Macc detection are shown. Symbol colors and sizes are scaled according to the value of ˙ Macc . The black boxes mark the regions of maximum stellar density in the cluster (from Lamm et al. 2004). Double arrows mark a physical distance of 1.5 parsecs.

of objects for the two groups. On the other hand, the complexity of the cluster, with its clumpy structure, the presence of at least two separated poles of star formation, the association with several massive OB stars, intrinsically translates to a non-linearity in the cluster evolution and, consequently, to a certain degree of ambiguity in the information provided by the spatial distribution of more evolved members.

## 6. Conclusions

Wehave performed an extensive photometric investigation of accretion in the young open cluster NGC 2264. Our sample, homogeneous in the dataset and in the characterization of individual properties, includes GLYPH&lt;24&gt; 750 members, spanning in mass the range GLYPH&lt;24&gt; 0.1-2 M GLYPH&lt;12&gt; ; of these, about 40% show ongoing accretion. Accretion activity has been traced from the detection and monitoring of the u -band excess emission relative to non-accreting young stars; we showed that this excess is a reliable indicator of accretion, matching the properties observed at di GLYPH&lt;11&gt; erent wavelengths (H GLYPH&lt;11&gt; emission, IR excess), with the advantage of providing alone a direct proxy to the accretion luminosity.

We re-examined the membership of the cluster from an accretion-related perspective, exploring a wider area around the spatial distribution of currently known members within the starforming region. This allowed us to expand the census of accreting members, with the identification of 50 new candidates based on their accretion signatures at short wavelengths. We detected UV color excess ranging from a few tenths of a mag to GLYPH&lt;24&gt; 3 mag, with a sensitivity of GLYPH&lt;24&gt; 0.2 mag deriving from the definition of the reference color locus from the properties of non-accreting young stars. Following a new empirical prescription for the conversion of the UV excess into accretion luminosity, we inferred

Little contribution to the spread in the ˙ Macc vs. M GLYPH&lt;3&gt; distribution arises from ˙ Macc variability, which is found to amount to GLYPH&lt;24&gt; 0.5 dex on the average, at least on days to months timescales; hence, this spread is physically relevant and provides valuable insight on the complexity of the accretion and evolutionary picture for a given stellar population. A di GLYPH&lt;11&gt; erential distribution of evolutionary stages among members, such as we may deduce from the HR diagram of the cluster and or from disk evolution / tracers ( GLYPH&lt;11&gt; IRAC ), would naturally impact the global accretion picture we may infer for the cluster as a whole. Indeed, ˙ Macc is shown to anticorrelate with such intrinsic evolutionary age in-/ dicators, although these trends are in turn not exempt from a large amount of scatter. The detected spread in ˙ Macc may thus reflect the dynamics of cluster evolution and the presence of a non-negligible age spread among its members, be it the result of sequential star formation episodes and or of a multiplicity / of evolutionary patterns. Moreover, a diversity in the accretion regimes at play, from episodic, short-lived accretion bursts to steadier, funnel-flow accretion, may concur to shape the features of di GLYPH&lt;11&gt; erent stellar groups and result in a range of observed accretion rates, from the upper envelope of the ˙ Macc vs. M GLYPH&lt;3&gt; distribution to lower accretion rate regimes.

An important contribution toward the achievement of a complete physical picture of the accretion process is provided by a detailed analysis of the variability properties. Geometric e GLYPH&lt;11&gt; ects constitute an important filter to the e GLYPH&lt;11&gt; ective variability profile for a given system, contributing up to 75%, on the average, of the detected ˙ Macc variability; however, much information on the underlying physics is conveyed by the specific photometric signatures of this variability, and a comparative analysis over intrinsically di GLYPH&lt;11&gt; erent stellar groups may be of remarkable interest to discern di GLYPH&lt;11&gt; erent scenarios. This aspect will be fully addressed for the young stellar population of NGC 2264 in a forthcoming paper (Venuti et al., in prep.); this is intended to provide a complete characterization of the accretion process and its impact onto the dynamics of the cluster.

Acknowledgements. We would like to thank Christian Veillet, former director of CFHT, for granting discretionary time to perform the mapping survey in December 2010, and the Terapix center at Institut d'Astrophysique de Paris, and in particular Yannick Mellier, for the prompt processing of the MegaCam images obtained during this run. We also thank Nadine Manset and Jim Thomas at CFHT for e GLYPH&lt;14&gt; cient run scheduling and data retrieval procedures. We warmly thank Luisa Rebull for retrieving for us the L acc values she derived for her 2002 paper on disk-bearing objects in NGC 2264. We thank Kevin Covey for discussions on SDSS dwarf sequences and Lynne Hillenbrand for discussions on bolometric corrections scales; we also acknowledge useful discussion on the ˙ Macc -

M relationship with Beate Stelzer. This publication makes use of data products GLYPH&lt;3&gt; from the Sloan Digital Sky Survey and the Two Micron All Sky Survey. This project was in part supported by the grant ANR 2011 Blanc SIMI5-6 020 01. SHPA and APS acknowledge support from CNPq, CAPES and Fapemig.

## References

Abazajian, K., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2004, AJ, 128, 502

- Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543
- Adams, M. T., Strom, K. M., &amp; Strom, S. E. 1983, ApJS, 53, 893
- A er, L., Micela, G., Favata, F., Flaccomio, E., &amp; Bouvier, J. 2013, MNRAS, GLYPH&lt;11&gt; 430, 1433
- Alcalá, J. M., Natta, A., Manara, C. F., et al. 2014, A&amp;A, 561, A2

Alencar, S. H. P., Teixeira, P. S., Guimarães, M. M., et al. 2010, A&amp;A, 519, A88

Allison, R. J., Goodwin, S. P., Parker, R. J., et al. 2009, MNRAS, 395, 1449

- Annis, J., Soares-Santos, M., Strauss, M. A., et al. 2011, ArXiv e-print 1111.6619 Armitage, P. J. 2011, ARA&amp;A, 49, 195
- Bara GLYPH&lt;11&gt; e, I., Chabrier, G., &amp; Gallardo, J. 2009, ApJ, 702, L27

Barentsen, G., Vink, J. S., Drew, J. E., et al. 2011, MNRAS, 415, 103

Bertin, E. 2011, in Astronomical Society of the Pacific Conference Series, Vol.

442, Astronomical Data Analysis Software and Systems XX, ed. I. N. Evans,

A. Accomazzi, D. J. Mink, &amp; A. H. Rots, 435

Bertin, E. &amp; Arnouts, S. 1996, A&amp;AS, 117, 393

Bessell, M. S., Castelli, F., &amp; Plez, B. 1998, A&amp;A, 333, 231

- Boulade, O., Charlot, X., Abbon, P., et al. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. M. Iye &amp; A. F. M. Moorwood, Vol. 4841, 72-81

Bouvier, J., Alencar, S. H. P., Boutelier, T., et al. 2007a, A&amp;A, 463, 1017

- Bouvier, J., Alencar, S. H. P., Harries, T. J., Johns-Krull, C. M., &amp; Romanova, M. M. 2007b, Protostars and Planets V, 479
- Bouvier, J., Cabrit, S., Fernandez, M., Martin, E. L., &amp; Matthews, J. M. 1993, A&amp;A, 272, 176
- Calvet, N. &amp; Gullbring, E. 1998, ApJ, 509, 802
- Camenzind, M. 1990, in Reviews in Modern Astronomy, ed. G. Klare, Vol. 3, 234-265
- Cauley, P. W., Johns-Krull, C. M., Hamilton, C. M., &amp; Lockhart, K. 2012, ApJ, 756, 68
- Cody, A. M., Stau GLYPH&lt;11&gt; er, J., Baglin, A., et al. 2014, AJ, 147, 82
- Cody, A. M., Stau GLYPH&lt;11&gt; er, J. R., Micela, G., Baglin, A., &amp; CSI 2264 Team. 2013, Astronomische Nachrichten, 334, 63

Cohen, M. &amp; Kuhi, L. V. 1979, ApJS, 41, 743

- Costigan, G., Vink, J. S., Scholz, A., Ray, T., &amp; Testi, L. 2014, MNRAS, 440, 3444

Covey, K. R., Ivezi´ c, Ž., Schlegel, D., et al. 2007, AJ, 134, 2398

Covey, K. R., Lada, C. J., Román-Zúñiga, C., et al. 2010, ApJ, 722, 971

- Dahm, S. E. 2008, in Handbook of Star Forming Regions, Volume I, ed. B. Reipurth, 966

Dahm, S. E. &amp; Simon, T. 2005, AJ, 129, 829

Dullemond, C. P., Natta, A., &amp; Testi, L. 2006, ApJ, 645, L69

- Ercolano, B., Mayr, D., Owen, J. E., Rosotti, G., &amp; Manara, C. F. 2014, MNRAS Fang, M., Kim, J. S., van Boekel, R., et al. 2013, ApJS, 207, 5

Feigelson, E. D. &amp; Babu, G. J. 2012, Modern Statistical Methods for Astronomy, Cambridge University Press

F˝ urész, G., Hartmann, L. W., Szentgyorgyi, A. H., et al. 2006, ApJ, 648, 1090

Finlator, K., Ivezi´, Ž., Fan, X., et al. 2000, AJ, 120, 2615 c

Fiorucci, M. &amp; Munari, U. 2003, A&amp;A, 401, 781

Flaccomio, E., Micela, G., &amp; Sciortino, S. 2006, A&amp;A, 455, 903

Fukugita, M., Ichikawa, T., Gunn, J. E., et al. 1996, AJ, 111, 1748

Gillen, E., Aigrain, S., McQuillan, A., et al. 2014, A&amp;A, 562, A50

Guarcello, M. G., Prisinzano, L., Micela, G., et al. 2007, A&amp;A, 462, 245

Gullbring, E., Hartmann, L., Briceno, C., &amp; Calvet, N. 1998, ApJ, 492, 323

Haisch, Jr., K. E., Lada, E. A., &amp; Lada, C. J. 2001, ApJ, 553, L153

Hartmann, L., Calvet, N., Gullbring, E., &amp; D'Alessio, P. 1998, ApJ, 495, 385

Hartmann, L., D'Alessio, P., Calvet, N., &amp; Muzerolle, J. 2006, ApJ, 648, 484

Helsel, D. R. 2012, Statistics for Censored Environmental Data, Wiley

- Herbig, G. H. &amp; Bell, K. R. 1988, Third Catalog of Emission-Line Stars of the Orion Population : 3 : 1988

Herbst, W., Herbst, D. K., Grossman, E. J., &amp; Weinstein, D. 1994, AJ, 108, 1906 Herczeg, G. J. &amp; Hillenbrand, L. A. 2008, ApJ, 681, 594

Hillenbrand, L. A. 2005, ArXiv astro-ph 0511083

/

Joy, A. H. 1945, ApJ, 102, 168

Königl, A. 1991, ApJ, 370, L39

Kraus, A. L. &amp; Hillenbrand, L. A. 2007, AJ, 134, 2340

Lada, C. J. 1987, in IAU Symposium, Vol. 115, Star Forming Regions, ed. M. Peimbert &amp; J. Jugaku, 1-17

Lada, C. J., Muench, A. A., Luhman, K. L., et al. 2006, AJ, 131, 1574

Article number, page 22 of 24

- Lamm, M. H., Bailer-Jones, C. A. L., Mundt, R., Herbst, W., &amp; Scholz, A. 2004,
- A&amp;A, 417, 557 Lenz, D. D., Newberg, J., Rosner, R., Richards, G. T., &amp; Stoughton, C. 1998, ApJS, 119, 121 Luhman, K. L., Stau GLYPH&lt;11&gt; er, J. R., Muench, A. A., et al. 2003, ApJ, 593, 1093 Manara, C. F., Robberto, M., Da Rio, N., et al. 2012, ApJ, 755, 154 Ménard, F. &amp; Bertout, C. 1999, in NATO ASIC Proc. 540: The Origin of Stars and Planetary Systems, ed. C. J. Lada &amp; N. D. Kylafis, 341 Meyer, M. R., Calvet, N., &amp; Hillenbrand, L. A. 1997, AJ, 114, 288 Mohanty, S., Jayawardhana, R., &amp; Basri, G. 2005, ApJ, 626, 498 Muzerolle, J., Hillenbrand, L., Calvet, N., Briceño, C., &amp; Hartmann, L. 2003, ApJ, 592, 266 Natta, A., Testi, L., &amp; Randich, S. 2006, A&amp;A, 452, 245 Nguyen, D. C., Scholz, A., van Kerkwijk, M. H., Jayawardhana, R., &amp; Brandeker, A. 2009, ApJ, 694, L153 Parker, R. J., Bouvier, J., Goodwin, S. P., et al. 2011, MNRAS, 412, 2489 Pecaut, M. J. &amp; Mamajek, E. E. 2013, ApJS, 208, 9 Press, W. H., Teukolsky, S. A., Vetterling, W. T., &amp; Flannery, B. P. 1992, Numerical recipes in C. The art of scientific computing Prim, R. C. 1957, Bell Syst. Tech. J., 36, 1389 Ramírez, S. V., Rebull, L., Stau GLYPH&lt;11&gt; er, J., et al. 2004, AJ, 127, 2659 Rebull, L. M., Makidon, R. B., Strom, S. E., et al. 2002, AJ, 123, 1528 Regnault, N., Conley, A., Guy, J., et al. 2009, A&amp;A, 506, 999 Rigliaco, E., Natta, A., Randich, S., Testi, L., &amp; Biazzo, K. 2011, A&amp;A, 525, A47 Sicilia-Aguilar, A., Henning, T., &amp; Hartmann, L. W. 2010, ApJ, 710, 597 Siess, L., Dufour, E., &amp; Forestini, M. 2000, A&amp;A, 358, 593 Skrutskie, M. F., Cutri, R. M., Stiening, R., et al. 2006, AJ, 131, 1163 Stau GLYPH&lt;11&gt; er, J., Cody, A. M., Baglin, A., et al. 2014, AJ, 147, 83 Sung, H., Bessell, M. S., Chun, M.-Y., Karimov, R., &amp; Ibrahimov, M. 2008, AJ, 135, 441 Sung, H., Bessell, M. S., &amp; Lee, S.-W. 1997, AJ, 114, 2644 Sung, H., Stau GLYPH&lt;11&gt; er, J. R., &amp; Bessell, M. S. 2009, AJ, 138, 1116 Taylor, M. B. 2005, in Astronomical Society of the Pacific Conference Series, Vol. 347, Astronomical Data Analysis Software and Systems XIV, ed. P. Shopbell, M. Britton, &amp; R. Ebert, 29 Teixeira, P. S., Lada, C. J., Marengo, M., &amp; Lada, E. A. 2012, A&amp;A, 540, A83 Vorobyov, E. I. &amp; Basu, S. 2008, ApJ, 676, L139 Walker, M. F. 1956, ApJS, 2, 365

White, R. J. &amp; Basri, G. 2003, ApJ, 582, 1109

Appendix A: Accretion-oriented revision of the census of NGC 2264 members: new CTTS candidate members and possible field contaminants

Table A.1. List of new candidate members from the accretion-based census of the population of NGC 2264.

| Mon-ID       |       RA |      Dec | Selection criteria    | Additional information 1                          |
|--------------|----------|----------|-----------------------|---------------------------------------------------|
| Mon-001428   | 100.193  |  9.71302 | UV, var I             | likely CTTS based on IRAC light curve             |
| Mon-005009   | 100.541  |  9.79834 | var II, lc            | irregular IRAC / CoRoT light curves               |
| Mon-005278   | 100.648  |  9.43207 | UV, var II            | IR excess                                         |
| Mon-005326   | 100.621  |  9.4455  | UV                    |                                                   |
| Mon-005385   | 100.354  |  9.75632 | UV, var I             | flat IRAC light curve                             |
| Mon-005455   | 100.215  |  9.37313 | UV                    | flat IRAC light curve                             |
| Mon-005596   | 100.516  |  9.45916 | var I                 | flat IRAC light curve                             |
| Mon-005664   | 100.227  |  9.15885 | UV                    | spotted-like CoRoT light curve                    |
| Mon-005745   | 100.518  |  9.16134 | UV                    |                                                   |
| Mon-005807   | 100.282  |  9.23241 | UV, var I             |                                                   |
| Mon-005836   | 100.371  |  9.3043  | UV                    | flat IRAC light curve                             |
| Mon-006037 * |  99.8723 |  9.72772 | var II, lc            | spotted-like CoRoT light curve                    |
| Mon-006079   |  99.8713 |  9.7107  | UV, var I, var II     | eclipsing binary / very short period binary       |
| Mon-006144   |  99.9284 |  9.55441 | UV                    | flat IRAC light curve                             |
| Mon-006183   |  99.972  |  9.82939 | UV, var II            | flat IRAC light curve                             |
| Mon-006286   | 100.126  |  9.19233 | UV                    |                                                   |
| Mon-006324   |  99.9531 |  9.29308 | UV                    |                                                   |
| Mon-006325   | 100.06   |  9.22705 | UV, var II            |                                                   |
| Mon-006369   |  99.8757 |  9.29094 | varI, var II          |                                                   |
| Mon-006398   |  99.8103 |  9.61444 | UV, var II            | flat IRAC light curve                             |
| Mon-006409   |  99.8134 |  9.64136 | UV                    | flat IRAC light curve                             |
| Mon-006429   |  99.7954 |  9.53839 | UV, var I, var II     | strong IR excess; IRAC light curve of YSO         |
| Mon-006465   |  99.8548 |  9.54391 | UV, var I, var II     | eclipsing binary                                  |
| Mon-006491   |  99.8563 |  9.5276  | UV, var II, lc        | IR excess; CoRoT light curve of CTTS              |
| Mon-006498   |  99.8275 |  9.39888 | UV                    | flat IRAC light curve                             |
| Mon-006515   |  99.7764 |  9.66813 | UV                    |                                                   |
| Mon-006873   |  99.7924 |  9.15089 | var I, var II         |                                                   |
| Mon-006902   |  99.7511 |  9.23571 | var I, var II         |                                                   |
| Mon-006930   |  99.767  |  9.27055 | var II, lc            | likely CTTS based on CoRoT light curve            |
| Mon-006950   |  99.7792 |  9.85195 | UV, var II            |                                                   |
| Mon-006985   |  99.7794 | 10.1191  | UV, var II            |                                                   |
| Mon-006986   |  99.7915 | 10.0869  | UV, var I, var II, lc | likely CTTS based on CoRoT light curve            |
| Mon-006991   |  99.8417 | 10.1065  | var II, lc            | likely CTTS based on CoRoT light curve; IR excess |
| Mon-006999   |  99.8665 | 10.0635  | var I                 |                                                   |
| Mon-007004   |  99.8581 | 10.0854  | UV                    |                                                   |
| Mon-007018   |  99.9725 | 10.1093  | var I, var II         |                                                   |
| Mon-007402   | 100.625  |  9.94183 | var I                 |                                                   |
| Mon-007418   | 100.576  | 10.0519  | UV                    |                                                   |
| Mon-007461   | 100.241  | 10.1179  | UV, var II            | probable MIPS 24 mexcess                          |
| Mon-007496   | 100.515  | 10.0884  | var I, var II         | GLYPH<22>                                         |
| Mon-007845   |  99.8919 |  9.40008 | UV                    |                                                   |
| Mon-008139   |  99.7565 | 10.0797  | UV                    |                                                   |
| Mon-008140   |  99.8043 |  9.61337 | UV, var II            |                                                   |
| Mon-008141   | 100.665  |  9.67704 | UV                    |                                                   |
| Mon-008142   |  99.7682 |  9.69662 | UV, var I, var II     |                                                   |
| Mon-014132   |  99.7648 |  9.27109 | var II, lc            | AA Tau-type CoRoT light curve; IR excess          |
| Mon-020662   | 100.211  | 10.0384  | UV                    |                                                   |
| Mon-021771   | 100.518  |  9.57561 | var I, var II         | flat IRAC light curve                             |
| Mon-025015   |  99.967  |  9.40013 | UV                    |                                                   |
| Mon-025101   |  99.8522 |  9.92528 | UV, var II            |                                                   |
| Mon-124054   | 100.391  |  9.92262 | UV, var I             |                                                   |

Notes. Selection criteria: UV = UV excess; var I = multi-year variability; var II = mid-term variability (i.e., on week timescales); lc = light curve morphology.

- ( 1 ) From Spitzer IRAC and or CoRoT monitoring obtained during the CSI 2264 campaign (and or additional IR photometry), when available. / /

( * ) WTTS candidate member

Table A.2. List of NGC 2264 members, with no apparent accretion activity from previous surveys, reclassified as CTTS from this accretion-based census of the population of NGC 2264.

| Mon-ID     |       RA |      Dec | Selection criteria   | Additional information 1                                                                                                                                      |
|------------|----------|----------|----------------------|---------------------------------------------------------------------------------------------------------------------------------------------------------------|
| Mon-000293 | 100.252  |  9.68758 | UV                   |                                                                                                                                                               |
| Mon-000316 | 100.232  |  9.64025 | UV                   | IR excess; likely CTTS based on IRAC light curve flat IRAC light curve flat IRAC light curve; possible 8 mexcess                                              |
| Mon-000340 | 100.483  |  9.67968 | UV                   | IR excess; likely CTTS based on IRAC light curve flat IRAC light curve flat IRAC light curve; possible 8 mexcess                                              |
| Mon-000507 | 100.261  |  9.85702 | UV                   | GLYPH<22>                                                                                                                                                     |
| Mon-000685 | 100.189  |  9.8195  | UV                   | flat IRAC light curve                                                                                                                                         |
| Mon-000718 | 100.267  |  9.34567 | UV                   | flat IRAC light curve                                                                                                                                         |
| Mon-000750 | 100.327  |  9.52487 | UV                   | IRAC light curve; possible 8 GLYPH<22>                                                                                                                        |
| Mon-000823 | 100.209  |  9.95111 | UV                   | IR excess; likely CTTS based on IRAC light curve flat IRAC light curve flat IRAC light curve; possible 8 mexcess                                              |
| Mon-000914 | 100.216  |  9.6322  | UV                   | very strong IR excess                                                                                                                                         |
| Mon-000957 | 100.19   |  9.20354 | UV                   | IR excess                                                                                                                                                     |
| Mon-000959 | 100.188  |  9.2685  | var                  | IR excess                                                                                                                                                     |
| Mon-001060 | 100.128  | 10.0008  | UV                   | IR excess; likely CTTS based on IRAC light curve IR excess                                                                                                    |
| Mon-001211 |  99.9598 |  9.52143 | UV                   | IR excess; likely CTTS based on IRAC light curve IR excess                                                                                                    |
| Mon-001224 | 100.07   |  9.67295 | UV                   | IR excess; likely CTTS based on IRAC light curve IR excess                                                                                                    |
| Mon-001241 | 100.141  |  9.81216 | UV                   | IR excess; likely CTTS based on IRAC light curve strong IR excess; IRAC light curve of YSO strong IR excess strong MIPS 24 GLYPH<22> mexcess strong IR excess |
| Mon-001267 | 100.143  |  9.42142 | var                  | IR excess; likely CTTS based on IRAC light curve strong IR excess; IRAC light curve of YSO strong IR excess strong MIPS 24 GLYPH<22> mexcess strong IR excess |
| Mon-001280 |  99.9177 |  9.30604 | UV                   | IR excess; likely CTTS based on IRAC light curve strong IR excess; IRAC light curve of YSO strong IR excess strong MIPS 24 GLYPH<22> mexcess strong IR excess |
| Mon-001576 | 100.133  | 10.119   | UV                   | IR excess; likely CTTS based on IRAC light curve strong IR excess; IRAC light curve of YSO strong IR excess strong MIPS 24 GLYPH<22> mexcess strong IR excess |
| Mon-001582 | 100.499  | 10.0638  | UV                   | IR excess; likely CTTS based on IRAC light curve strong IR excess; IRAC light curve of YSO strong IR excess strong MIPS 24 GLYPH<22> mexcess strong IR excess |

Notes. Selection criteria: UV = UV excess; var = multi-year variability.

( 1 ) From Spitzer IRAC monitoring obtained during the CSI 2264 campaign (and or additional IR photometry), when available. / /

Table A.3. Objects with questioned membership from the CFHT UV-census.

| Mon-ID       |       RA |      Dec |      u |      g |      r |      i | Membership criteria   |
|--------------|----------|----------|--------|--------|--------|--------|-----------------------|
| Mon-000003   | 100.439  |  9.78583 | 20.46  | 18.204 | 16.924 | 16.206 | P, X                  |
| Mon-000004   | 100.561  |  9.84104 | 17.808 | 15.998 | 15.075 | 14.594 | P, RV                 |
| Mon-000025   | 100.511  |  9.97441 | 18.025 | 16.134 | 15.193 | 14.69  | P, RV                 |
| Mon-000031   | 100.555  |  9.76763 | 17.801 | 16.118 | 15.307 | 14.888 | P, RV                 |
| Mon-000387   | 100.455  |  9.68496 | 20.53  | 18.408 | 17.273 | 16.681 | X                     |
| Mon-000414   | 100.559  |  9.59579 | 18.446 | 16.44  | 15.431 | 14.944 | P, RV                 |
| Mon-000593 * | 100.257  |  9.35157 | 17.793 | 15.822 | 15.077 | 14.836 | P, RV                 |
| Mon-000600 * | 100.205  |  9.98095 | 22.567 | 19.684 | 18.106 | 17.212 | X                     |
| Mon-000696 * | 100.337  |  9.5717  | 20.309 | 19.512 | 19.589 | 18.516 | h                     |
| Mon-000820   | 100.481  |  9.50528 | 19.538 | 17.513 | 16.49  | 15.902 | X                     |
| Mon-000944   | 100.458  |  9.68479 | 20.258 | 18.274 | 17.185 | 16.629 | X                     |
| Mon-000946 * | 100.202  |  9.99424 | 21.252 | 18.81  | 17.472 | 16.715 | II                    |
| Mon-000966   | 100.196  |  9.30923 | 17.729 | 15.882 | 14.942 | 14.415 | P, RV                 |
| Mon-000976   | 100.433  |  9.15588 | 21.325 | 19.108 | 17.798 | 17.058 | X                     |
| Mon-000983   | 100.242  |  9.30085 | 16.047 | 14.663 | 14.206 | 14.042 | RV                    |
| Mon-000988   | 100.088  |  9.98401 | 20.823 | 18.249 | 16.758 | 15.93  | P, RV                 |
| Mon-001043 * | 100.064  |  9.71678 | 21.746 | 19.149 | 17.695 | 16.854 | P, II / III           |
| Mon-001046   |  99.9849 |  9.72542 | 15.811 | 14.455 | 14.048 | 13.881 | RV                    |
| Mon-001086   |  99.8899 |  9.86127 | 18.279 | 16.524 | 15.639 | 15.144 | P, RV                 |
| Mon-001192   | 100.162  |  9.36066 | 17.954 | 15.866 | 15.079 | 14.799 | RV                    |
| Mon-001207   |  99.9804 |  9.98644 | 21.113 | 18.345 | 16.911 | 16.009 | h                     |
| Mon-001283 * |  99.861  |  9.51055 | 21.364 | 18.75  | 17.59  | 17.012 | RV                    |
| Mon-001289   |  99.8527 |  9.37588 | 16.778 | 15.082 | 14.379 | 14.084 | P, RV                 |
| Mon-001293   |  99.8972 |  9.54233 | 18.331 | 16.315 | 15.239 | 14.681 | P, RV                 |
| Mon-001311   |  99.8151 |  9.65357 | 18.014 | 16.134 | 15.181 | 14.684 | P, RV                 |
| Mon-001393   | 100.552  |  9.98531 | 17.072 | 15.373 | 14.609 | 14.195 | P, RV                 |
| Mon-006019 * | 100.142  |  9.42022 | 19.511 | 17.507 | 16.217 | 15.507 | P                     |
| Mon-007454   | 100.129  | 10.0801  | 17.799 | 16.099 | 15.494 | 15.14  | P                     |

Notes. Membership criteria, in the literature, for the listed objects: P = photometric member (i.e., located on a cluster sequence in the R vs. R-I diagram); X = X-ray emitter, from Sung et al. (2009); RV = radial velocity member, from F˝ urész et al. (2006); h = H GLYPH&lt;11&gt; emission candidate, from Sung et al. (2009); II = Class II object, from Sung et al. (2009); II III / = Class II III object, from Sung et al. (2009). /

( * ) More uncertain cases
Laura Venuti Jérôme Bouvier Ettore Flaccomio Silvia H. P. Alencar Jonathan Irwin John R. Stauffer Ann Marie Cody Paula S. Teixeira Alana P. Sousa Giuseppina Micela Jean-Charles Cuillandre Giovanni Peres
Mapping accretion and its variability in the young open cluster NGC 2264: a study based on u-band photometry
We aim at characterizing the accretion properties of several hundred members
of the star-forming cluster NGC 2264 (3 Myr). We performed a deep u,g,r,i
mapping and a simultaneous u+r monitoring of the region with CFHT/MegaCam in
order to directly probe the accretion process from UV excess measurements.
Photometric properties and stellar parameters are determined homogeneously for
about 750 monitored young objects, spanning the mass range 0.1-2 Mo. About 40%
are classical (accreting) T Tauri stars, based on various diagnostics (H_alpha,
UV and IR excesses). The remaining non-accreting members define the
(photospheric+chromospheric) reference UV emission level over which flux excess
is detected and measured. We revise the membership status of cluster members
based on UV accretion signatures and report a new population of 50 CTTS
candidates. A large range of UV excess is measured for the CTTS population,
varying from a few 0.1 to 3 mag. We convert these values to accretion
luminosities and obtain mass accretion rates ranging from 1e-10 to 1e-7 Mo/yr.
Taking into account a mass-dependent detection threshold for weakly accreting
objects, we find a >6sigma correlation between mass accretion rate and stellar
mass. A power-law fit, properly accounting for upper limits, yields M_acc
$\propto$ M^{1.4+/-0.3}. At any given stellar mass, we find a large spread of
accretion rates, extending over about 2 orders of magnitude. The monitoring of
the UV excess on a timescale of a couple of weeks indicates that its
variability typically amounts to 0.5 dex, much smaller than the observed
spread. We suggest that a non-negligible age spread across the cluster may
effectively contribute to the observed spread in accretion rates at a given
mass. In addition, different accretion mechanisms (like, e.g., short-lived
accretion bursts vs. more stable funnel-flow accretion) may be associated to
different M_acc regimes.