# CoRoT 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk

E. Gillen, S. Aigrain, A. Mcquillan, J. Bouvier, S. Hodgkin, S. H. P. Alencar, C. Terquem, J. Southworth, N. P. Gibson, A. Cody, et al.

► To cite this version:

E. Gillen, S. Aigrain, A. Mcquillan, J. Bouvier, S. Hodgkin, et al. CoRoT 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk. Astronomy and Astrophysics - A&A, 2014, 562, 10.1051/0004-6361/201322493. insu-03618607

HAL Id: insu-03618607

https://insu.hal.science/insu-03618607

Submitted on 24 Mar 2022

**HAL** is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire **HAL**, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# CoRoT* 223992193: A new, low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk

E. Gillen<fnref n="1" />, S. Aigrain<fnref n="1" />, A. McQuillan<fnref n="1" />,<fnref n="2" />, J. Bouvier<fnref n="3" />, S. Hodgkin<fnref n="4" />, S. H. P. Alencar<fnref n="5" />, C. Terquem<fnref n="1" />, J. Southworth<fnref n="6" />, N. P. Gibson<fnref n="7" />, A. Cody<fnref n="8" />, M. Lendl<fnref n="9" />, M. Morales-Calderón<fnref n="10" />, F. Favata<fnref n="11" />, J. Stauffer<fnref n="8" />, and G. Micela<fnref n="12" />

Received 15 August 2013 / Accepted 10 November 2013

## ABSTRACT

We present the discovery of CoRoT 223992193, a double-lined, detached eclipsing binary, comprising two pre-main sequence M dwarfs, discovered by the CoRoT space mission during a 23-day observation of the 3 Myr old NGC 2264 star-forming region. Using multi-epoch optical and near-IR follow-up spectroscopy with FLAMES on the Very Large Telescope and ISIS on the *William Herschel* Telescope we obtain a full orbital solution and derive the fundamental parameters of both stars by modelling the light curve and radial velocity data. The orbit is circular and has a period of 3.8745745 ± 0.0000014 days. The masses and radii of the two stars are 0.67 ± 0.01 and 0.495 ± 0.007 $M_\odot$ and 1.30 ± 0.04 and 1.11$^{+0.04}_{-0.05} R_\odot$, respectively. This system is a useful test of evolutionary models of young low-mass stars, as it lies in a region of parameter space where observational constraints are scarce; comparison with these models indicates an apparent age of ~3.5–6 Myr. The systemic velocity is within $1\sigma$ of the cluster value which, along with the presence of lithium absorption, strongly indicates cluster membership. The CoRoT light curve also contains large-amplitude, rapidly evolving out-of-eclipse variations, which are difficult to explain using starspots alone. The system’s spectral energy distribution reveals a mid-infrared excess, which we model as thermal emission from a small amount of dust located in the inner cavity of a circumbinary disk. In turn, this opens up the possibility that some of the out-of-eclipse variability could be due to occultations of the central stars by material located at the inner edge or in the central cavity of the circumbinary disk.

**Key words.** binaries: eclipsing – stars: pre-main sequence – binaries: spectroscopic – stars: individual: CoRoT 223992193 – open clusters and associations: individual: NGC 2264 – protoplanetary disks

## 1. Introduction

Detached, double-lined eclipsing binaries (EBs) are extremely valuable objects because their masses, radii, effective temperatures and luminosities can be determined in a model-independent manner from the light and radial velocity curves of the system. When these reach a precision of a few percent or less, they provide one of the most powerful tests of stellar evolution models available (Andersen 1991; Torres et al. 2010). As these models underpin much of astrophysics, it is vital that they are tested as rigorously as possible. The two components of a given EB can generally be assumed to share the same age and metallicity, which adds to the tightness of the constraints.

Figure 1 shows the existing mass and radius measurements for low-mass, detached EBs<fnref n="1" />. While there are now many well-characterised systems on the main sequence, there are very few on the pre-main sequence (PMS): to the best of our knowledge, there are only eight published low-mass EBs (with known stellar parameters), where both components are PMS objects with masses below 1.5 $M_\odot$ (see Table 1), and all but one of these are located in Orion. The PMS is an important region of parameter space because it corresponds to a period of very rapid evolution, when the models are still sensitive to their initial conditions, motivating efforts to discover and characterise more of these systems.

Perhaps the most important discrepancy between current observations and models is the fact that low-mass stars are

<footnote n="1">¹ Sub-department of Astrophysics, Department of Physics, University of Oxford, Keble Road, Oxford, OX1 3RH, UK
e-mail: ed.gillen@astro.ox.ac.uk</footnote>

<footnote n="2">² School of Physics and Astronomy, Raymond and Beverly Sackler, Faculty of Exact Sciences, Tel Aviv University, 69978 Tel Aviv, Israel</footnote>

<footnote n="3">³ UJF-Grenoble 1/CNRS-INSU, Institut de Planétologie et d'Astrophysique de Grenoble (IPAG) UMR 5274, 38041 Grenoble, France</footnote>

<footnote n="4">⁴ Institute of Astronomy, Madingley Road, Cambridge, CB3 0HA, UK</footnote>

<footnote n="5">⁵ Departamento de Física – ICEx – UFMG, Av. Antônio Carlos, 6627, 30270-901 Belo Horizonte, MG, Brazil</footnote>

<footnote n="6">⁶ Astrophysics Group, Keele University, Staffordshire, ST5 5BG, UK</footnote>

<footnote n="7">⁷ European Southern Observatory, Karl-Schwarzschild-Str. 2, 85748 Garching bei München, Germany</footnote>

<footnote n="8">⁸ Spitzer Science Center, California Institute of Technology, 1200 E California Blvd., Pasadena, CA 91125, USA</footnote>

<footnote n="9">⁹ Observatoire de Genève, Université de Genève, Chemin des Maillettes 51, 1290 Sauverny, Switzerland</footnote>

<footnote n="10">¹⁰ Centro de Astrobiología (INTA-CSIC); ESAC Campus, PO Box 78, 28691 Villanueva de la Canada, Spain</footnote>

<footnote n="11">¹¹ European Space Agency, 8–10 rue Mario Nikis, 75738 Paris Cedex 15, France</footnote>

<footnote n="12">¹² INAF – Osservatorio Astronomico di Palermo, Piazza del Parlamento 1, 90134 Palermo, Italy</footnote>

<footnote>* The CoRoT space mission was developed and is operated by the French space agency CNES, with participation of ESAs RSSD and Science Programmes, Austria, Belgium, Brazil, Germany, and Spain.</footnote>

<footnote n="1">¹ Data from John Southworth’s catalogue, http://www.astro.keele.ac.uk/~jkt/debdata/debs.html</footnote>
<figure><img src="image_2.png" /></figure>

**Fig. 1.** Mass–radius relation for low-mass EBs. The black points show measurements for stars with masses <1.5 $M_\odot$ in detached EBs¹, and the lines show, from top to bottom, the theoretical isochrones of Baraffe et al. (1998, BCAH98) for 1–5, 6.3, 8, 10 Myr and 1 Gyr (brown, blue, purple, lilac, green, grey, ochre, cyan and black, respectively. $Y = 0.282$, [M/H] = 0, mixing length $\alpha = 1.9$). The components of the new system presented in this paper are shown by the red diamonds. Note that it lies in a very sparsely populated region of the diagram, making it a valuable test of PMS stellar evolution models. For comparison, we have also labelled the three lowest mass systems known in the Orion Nebula Cluster (see Table 1 for details. For clarity, the higher mass systems are not labelled). Colour versions of all figures are available in the online article.

<table><caption>**Table 1.** Known low-mass, pre-main sequence eclipsing binary systems.</caption><thead><tr><th>Name</th><th><i>M</i><sub>pri</sub><br>(<i>M</i><sub>☉</sub>)</th><th><i>M</i><sub>sec</sub><br>(<i>M</i><sub>☉</sub>)</th><th><i>R</i><sub>pri</sub><br>(<i>R</i><sub>☉</sub>)</th><th><i>R</i><sub>sec</sub><br>(<i>R</i><sub>☉</sub>)</th><th>Cluster*</th><th>Age<br>(Myr)</th><th>Reference(s)</th></tr></thead><tbody><tr><td>2MJ0535-05</td><td>0.06</td><td>0.04</td><td>0.68</td><td>0.49</td><td>ONC</td><td>~1</td><td>Stassun et al. (2006, 2007)</td></tr><tr><td>JW 380</td><td>0.26</td><td>0.15</td><td>1.19</td><td>0.90</td><td>ONC</td><td>~1</td><td>Irwin et al. (2007)</td></tr><tr><td>Par 1802</td><td>0.41</td><td>0.41</td><td>1.82</td><td>1.69</td><td>ONC</td><td>~1</td><td>Cargile et al. (2008); Stassun et al. (2008)</td></tr><tr><td>ISOY J0535-0447</td><td>0.83<sup>a</sup></td><td>0.05<sup>a</sup></td><td></td><td></td><td>ONC</td><td>~1</td><td>Morales-Calderón et al. (2012)</td></tr><tr><td>V1174 Ori</td><td>1.00</td><td>0.73</td><td>1.34</td><td>1.07</td><td>Ori OB 1c</td><td>~5–10</td><td>Stassun et al. (2004)</td></tr><tr><td>RX J0529.4+0041A</td><td>1.27</td><td>0.93</td><td>1.44</td><td>1.35</td><td>Ori OB 1a</td><td>~7–13</td><td>Covino et al. (2000, 2001, 2004)</td></tr><tr><td>ASAS J0528+03</td><td>1.38</td><td>1.33</td><td>1.83</td><td>1.73</td><td>Ori OB 1a</td><td>~7–13</td><td>Stempels et al. (2008)</td></tr><tr><td>MML 53</td><td>≥0.97</td><td>≥0.84</td><td></td><td></td><td>UCL</td><td>~15</td><td>Hebb et al. (2010, 2011)</td></tr></tbody></table>

**Notes.** (a) Preliminary estimates. (b) ONC = Orion Nebula Cluster and UCL = Upper Centaurus Lupus.

observed to be cooler and larger than predicted by the models: theoretical models under-predict radii by up to 20% (Covino et al. 2004; Stassun et al. 2007; Coughlin et al. 2011) and masses by 10 to 30% (the greater discrepancy corresponding to lower masses; Hillenbrand & White 2004; Mathieu et al. 2007), for both PMS and main sequence objects. Testing stellar evolution models with EBs assumes that the stars can be considered to evolve independently. This has been questioned for close-separation binaries, which are spun up by tidal interactions, and are more magnetically active than single stars of comparable masses and ages: Chabrier et al. (2007) studied the impact of this enhanced magnetic field on convective processes within low-mass stars, finding that it could explain the observed cool temperatures and large radii. This is especially relevant to young binaries, which are particularly active. Coughlin et al. (2011)

find that the radius excess seen in members of EBs decreases as the orbital period increases, suggesting that the radius excess may be a by-product of interactions between the two components of a close binary system. Indeed, long period main sequence EBs, with component masses $M \ge 0.7 M_\odot$, appear not to exhibit a radius excess (e.g. Lacy et al. 2005; Bass et al. 2012; Southworth 2013). However, it is worth noting that at lower masses a radius excess is observed in some main sequence systems, even at long periods (e.g. Irwin et al. 2011).

Focussing on binary stars for model evaluation makes it difficult to distinguish between the effects of binarity and other physics which may be missing in the models. An alternative approach, which does not suffer from this problem, is to combine interferometric angular diameter measurements and distance estimates from parallaxes to obtain radii, luminosities and effective
temperatures for single stars. However, this is feasible only for
very nearby, bright stars, and their mass remains unknown (un-
less the target is a member of a visual binary with a well char-
acterised orbit). Asteroseismology is also an exquisite probe of
stellar interiors, but while it can provide model-independent den-
sity estimates, it only gives model-dependent masses and radii.
Therefore, detached, double-lined EBs remain an important ob-
servational test of stellar evolutionary models, and significant
resources are dedicated to discovering and characterising them
across as wide a range of mass, age and metallicity as possible.

The lack of PMS EBs is widely recognised, and numerous
observational programmes have been set up to detect and charac-
terise more systems; see for example the Monitor and YSOVAR
projects (Aigrain et al. 2007; Morales-Calderón et al. 2011).
With their wide field of view, excellent photometric precision
and continuous monitoring capability over weeks or even years,
the CoRoT (Baglin 2003) and Kepler (Borucki et al. 2010) space
missions are extremely efficient at detecting EBs (see e.g. Prša
et al. 2011), but they normally target main sequence field stars.
Fortunately, the young open cluster NGC 2264 falls within the
visibility zone of CoRoT, which enabled it to be observed con-
tinuously for over three weeks in March 2008. The resulting
dataset offers an unprecedented insight into the variability of
PMS stars, and was used to study accretion (Alencar et al. 2010),
rotation (Favata et al. 2010; Affer et al. 2013) and PMS pulsa-
tions (Zwintz et al. 2011), as well as to search for EBs.

NGC 2264 is a very well-studied young star forming re-
gion due to its relative proximity, well-defined membership list
and low foreground extinction (Dahm 2008). It is the domi-
nant component of the Mon OB1 association in the Monoceros
constellation, situated in the local spiral arm at a distance of
~700–900 pc (e.g. Sung et al. 1997; Baxter et al. 2009; Sung &
Bessell 2010). The cluster age has been estimated to be ~3 Myr
old from isochrone fitting to the low-mass stellar population
(e.g. Sung et al. 1997; Park et al. 2000; Rebull et al. 2002) but
there is an apparent dispersion of ~5 Myr from the broadened
sequence of suspected members, and star formation is ongo-
ing (e.g. Reipurth et al. 2004; Young et al. 2006; Dahm 2008).
Literature estimates for the total stellar population are in the
region of ~1000 members (Dahm 2008). Many new candidate
members have been identified from the CoRoT observation it-
self and the deep CFHT ugr imaging (Venuti et al., in prep.)
bringing the total number of confident members to at least 1500
(Cody & Flaccomio, priv. comm.). The cluster has a recessional
velocity of V = 22 ± 3.5 km s⁻¹ (Fúrész et al. 2006), and prelim-
inary metallicity estimates are [Fe/H] ≈ −0.15 and near-solar
for other elements (King et al. 2000). The mean reddening along
the line of sight to the cluster is E(B − V) = 0.06–0.15 mag
(Walker 1956; Park et al. 2000; Rebull et al. 2002; Sung et al.
2004; Mayne & Naylor 2008).

The 2008 CoRoT observation of NGC 2264 revealed sev-
eral tens of new EBs, over a dozen of which may be cluster
members. These will be described in more detail in future pub-
lications; here we present the discovery and characterisation of
one particular system, CoRoT 223992193. This system is also
referred to as 2MASS J06414422+0925024, W6712 (Sung et al.
2008) and has a CSI 2264 (Coordinated Synoptic Investigation
of NGC 2264) identifier of Mon-000256. Its light curve dis-
plays deep (~30%) eclipses with a period of 3.87 days, as
well as significant (~15% peak-to-peak) irregular out-of-eclipse
(OOE) variability. Its optical and near-IR colours are compatible
with cluster membership, and preliminary light curve modelling
and follow-up radial velocity measurements showed it to be
a double-lined, near-equal mass system with a total mass of

~1.2 M<sub>☉</sub> and individual radii >1 R<sub>☉</sub>, as expected for a PMS sys-
tem. This motivated additional follow-up spectroscopy and mod-
elling, enabling us to refine the fundamental parameters of both
stars and to confirm that the systemic velocity is compatible with
that of the cluster.

The system also shows indirect evidence for a proto-
planetary disk, most likely circumbinary. If confirmed, it will
be the first PMS EB system found to harbour such a disk. This
is particularly interesting in the light of the recent discovery of a
number of circumbinary exoplanets (e.g. Kepler-16 Doyle et al.
2011; Kepler-34 and 35 Welsh et al. 2012; Kepler-38 Orosz et al.
2012a; Kepler-47 Orosz et al. 2012b; and Kepler-64 Schwamb
et al. 2013). In all of these systems, the orbital angular momen-
tum of the planets is closely aligned with that of the host binary,
strongly suggesting in situ formation within a circumbinary disk.
Initial theoretical modelling suggests that circumbinary planets
around close separation binaries should be a common occurrence
(Alexander 2012).

CoRoT 223992193 is a very rich system. This paper focuses on the first step in its characterisation, namely the parameters of the two stars based on the 2008 CoRoT data and follow-up radial velocities. The detailed study of the interaction of the stars with their environment requires additional, multi-band photometry and spectroscopy, and will be the subject of a later paper. In Sect. 2, we give details of the CoRoT observations and follow-up spectroscopy. In Sect. 3 we model the light curve and radial velocity data, deriving fundamental parameters. We discuss the properties of the stars in Sect. 4, where we also present a preliminary analysis of the spectral energy distribution (SED) of the system, before concluding in Sect. 5.

2. Observations

2.1. Photometry

CoRoT observed NCG 2264 continuously for 23.4 days between
7 and 31 March 2008 (run SRa01). CoRoT observations are con-
ducted in a broad 300–1000 nm bandpass with a standard ca-
dence of 512 s, giving 3936 photometric data points for this run.
A total of 8150 stars with magnitudes 9.5 < R < 17 were mon-
itored in the 1.3 × 2.6 degrees field-of-view of CoRoT’s exo-
planet channel. About 1000 of these were previously known or
suspected members of NGC 2264. These observations<fnref n="2" /> represent
an unprecedented photometric dataset for a young cluster, both
in terms of sampling and precision. Basic aperture photometry is
carried out on board the satellite, and the data were further pro-
cessed by the CoRoT pipeline (Auvergne et al. 2009) to correct
them for known instrumental effects (including background and
pointing jitter).

To search for eclipses, we first applied a short-baseline run-
ning median filter to exclude outlying data points and a 1 day
baseline iterative non-linear filter (see Aigrain & Irwin 2004)
to remove long-term variations such as stellar activity. An au-
tomated, least-squares search for trapezoidal eclipses was then
performed in the light curves of all targets, followed by visual
examination of all candidates with an eclipse signal-to-noise ra-
tio (S/N) above 30. Systems with only one detected eclipse were
discarded, along with spurious detections induced by hot pixels
(identified from the shape of the “eclipses”, and the three-colour
light curves, which are available for stars with R < 15).

This process yielded 103 eclipsing systems, which were
cross-matched with 2MASS and other published data on

<footnote n="2">² The data are publicly available from the IAS CoRoT archive:
http://idoc-corot.ias.u-psud.fr/</footnote>
<table><caption>**Table 2.** Coordinates and photometric properties of CoRoT 223992193.</caption><thead><tr><th>RA</th><th>Dec</th></tr><tr><td>06 41 44.22</td><td>+09 25 02.398</td></tr><tr><th>Passband</th><th>(ref.)</th><th>Magnitude</th></tr></thead><tbody><tr><td>u</td><td>a</td><td>20.021 ± 0.048</td></tr><tr><td>g</td><td>a</td><td>17.543 ± 0.005</td></tr><tr><td>r</td><td>a</td><td>16.090 ± 0.004</td></tr><tr><td>i</td><td>a</td><td>15.069 ± 0.004</td></tr><tr><td>z</td><td>a</td><td>14.479 ± 0.004</td></tr><tr><td>J</td><td>b</td><td>13.329 ± 0.029</td></tr><tr><td>H</td><td>b</td><td>12.614 ± 0.022</td></tr><tr><td>K</td><td>b</td><td>12.331 ± 0.029</td></tr><tr><td>[3.6]</td><td>c</td><td>11.731 ± 0.031</td></tr><tr><td>[4.5]</td><td>c</td><td>11.533 ± 0.049</td></tr><tr><td>[5.8]</td><td>c</td><td>11.336 ± 0.065</td></tr><tr><td>[8.0]</td><td>c</td><td>10.951 ± 0.034</td></tr></tbody></table>

**Notes.** The photometric uncertainties listed here are the formal measurement errors; they do not account for the intrinsic variability of the system. (a) SDSS (Abazajian et al. 2009; Adelman-McCarthy et al. 2009), AB magnitudes; (b) 2MASS (Cutri et al. 2003), Vega magnitudes; (c) Spitzer/IRAC (Sung et al. 2009), Vega magnitudes.

NGC 2264, covering: X-rays (Flaccomio et al. 1997, 2006; Ramírez et al. 2004); disks (Rebull et al. 2002); optical and infrared variability (Makidon et al. 2004; Lamm et al. 2004) and spectroscopic and proper motion studies (Walker 1956). A 2MASS *J* vs. *J*−*K* colour–magnitude diagram was then used in conjunction with other membership information (where available), the spatial distribution of the targets, and a magnitude limit of *R* = 15.5, to select 12 possible cluster member systems for spectroscopic follow-up.

CoRoT 223992193 was selected, despite being fainter than our nominal cut-off (*R* = 15.74), because it fulfilled all the other criteria, and its light curve shows interesting OOE variations. The photometric properties of the system are shown in Table 2. Its magnitude and colours are consistent with a fairly low mass member of the cluster. Its light curve is shown in Fig. 2. The large-amplitude, rapidly evolving variability seen outside the eclipses could be due to spots, but would require unusually rapid spot evolution (on timescales of hours, this is discussed in more detail in Sect. 4.4). An alternative hypothesis is that the variability might be caused by partial occultation of the stars by circumstellar or circumbinary material.

A second set of CoRoT observations, lasting 40 days, took place in December 2011/January 2012 (run SRa05, see Fig. 3). We were not able to combine the two CoRoT light curves directly, because the background appears to be significantly more spatially variable during the second run, and the current CoRoT pipeline performs only a globally averaged background correction. Improvements to the pipeline are being tested, which will hopefully resolve this problem. In the mean time, we opted to work mainly with the 2008 light curve, using the 2011/2012 data only to refine the ephemeris of the system (see Sect. 3.1.2)

## 2.2. Spectroscopy

We performed low resolution spectroscopy to infer the combined spectral type and medium resolution spectroscopy to extract radial velocities.

We obtained a low-resolution (~7 Å) optical spectrum of CoRoT 223992193 on 13 April 2011 (JD 2455 665.3334) using

<figure><img src="image_3.png" /></figure>

**Fig. 2.** Light curve of CoRoT 223992193 obtained by CoRoT in 2008. The original light curve is shown in black in the top panel. The cyan line and light cyan shaded area show the mean and 95% confidence interval of the predictive distribution of the Gaussian Process used to model the out-of-eclipse variations. The bottom panel shows the light curve after subtracting the mean of the predictive distribution and re-normalising to unity. The detrended light curve, which was used to model the eclipses, has a typical photometric precision of ~0.26% per 512 s exposure.

<figure><img src="image_4.png" /></figure>

**Fig. 3.** Light curve of CoRoT 223992193 obtained by CoRoT in December 2011/January 2012.

the CAFOS focal reducer (Meisenheimer 1994) on the 2.2 m Calar Alto telescope equipped with the G-100 grism and the SITe#1d_15 CCD camera, over the 4600–7700 Å wavelength range. The spectrum was acquired at an airmass of 1.4 using a slit width of 1.5 arcsec, and an exposure time of 1800 s. A spectrum of the spectrophotometric standard, Feige 34, was obtained immediately after, at an airmass of 1.0, with the same set-up. The spectra of the two objects were reduced in the same way, namely: the 2D images were bias and flat-field corrected, and the stellar and neighbouring sky spectra were extracted using the IRAF/TWODSPEC package<fnref n="3" /> (Tody 1993). The wavelength scale was calibrated using the HgHeRb spectral lamp, and owing to the small number of lines available in the arc spectrum over the observed wavelength range, the wavelength calibration is accurate to only about 20 Å. Finally, the response of the instrument was corrected for in the target spectrum using that of the spectrophotometric standard. The resulting spectrum

<footnote n="3">³ IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the Association of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation.</footnote>
<figure><img src="image_5.png" /></figure>

**Fig. 4.** Low-resolution optical spectrum of CoRoT 223992193 taken with the 2.2 m CAHA telescope. Note the strong molecular bandheads (from TiO) characteristic of the early-M combined spectral type of the target, as well as strong Hα and Hβ emission.

is shown in Fig. 4. It exhibits clear emission at Hβ and Hα, with equivalent widths of 3.2 and 5.6 Å, respectively. The strong TiO absorption bands are indicative of a late spectral type. Comparing CoRoT 223992193’s photospheric spectrum with a grid of young spectral standards from Alves de Oliveira et al. (2012), we derive a spectral type of M2 and negligible visual extinction ($A_V \approx 0$).

We then obtained moderate resolution, near-IR spectra from the Intermediate dispersion Spectrograph and Imaging System (ISIS) on the 4 m William Herschel Telescope (WHT) at Roque de Los Muchachos observatory on La Palma, Spain. We observed the target at seven epochs with WHT/ISIS between 3 and 5 December 2011, covering the wavelength range, ~7850–8900 Å, with a spectral resolution of R ~ 12 000 and using the R1200R grating. Six observations were at quadrature and one close to eclipse. At each epoch we took three successive spectra, with exposure times of 300 or 600 s, directly followed by CuAr and CuNe arc lamp spectra. An example ISIS spectrum is shown in the top panel of Fig. 5. The ISIS spectra were processed with IRAF, using the ccdproc packages to perform the standard CCD reduction steps (trimming and overscan subtraction, bias subtraction and flat fielding) and specred.doslit to extract and wavelength calibrate the spectra using the arc lamp spectra. The three exposures taken at a given epoch were then combined to maximise the S/N and perform cosmic ray rejection. The resulting spectra typically have S/N ~ 20 per pixel (~25 per spectral resolution element; see Table 3).

A further 21 medium resolution optical spectra were obtained with the FLAMES instrument on the Very Large Telescope (VLT) at Paranal, Chile, as part of a multi-object monitoring programme primarily designed to study accretion (GO program 088.C-0239(A), PI Alencar). The spectra were obtained over a ~3 month period (4 December 2011–24 February 2012), with both sparsely and densely sampled time intervals. Telescope pointing issues resulted in significant loss of flux through the fibre during 6 of the 21 epochs, so that only 15 spectra could be used in the present analysis, but the phase coverage of the orbit is nonetheless fairly complete. Observations were taken using the standard setting HR15N,

<figure><img src="image_6.png" /></figure>

**Fig. 5.** Example continuum-corrected spectra from ISIS (top), and FLAMES (bottom). Grey shaded regions indicate the wavelength ranges used when cross-correlating the spectra with theoretical templates to measure radial velocities. Note that some of the emission lines visible in the FLAMES spectra are caused by the background nebula.

yielding a resolving power R ~ 17 000 over the wavelength range ~6440–6820 Å. The FLAMES spectra were reduced, extracted and wavelength-calibrated using the standard ESO pipeline<fnref n="4" />. An example FLAMES spectrum is shown in the bottom panel of Fig. 5. The most prominent absorption feature is the Li line at 6707.8 Å, which is consistent with the youth of the system. In this work, the FLAMES spectra were used only to derive radial velocities, but they contain a wealth of additional information, particularly in the broad, resolved and highly variable Hα emission line, which will be studied in more detail in a future paper.

Finally we also obtained three spectra with the Intermediate Dispersion Spectrograph (IDS) on the 2.5 m Isaac Newton Telescope (INT), also on La Palma. The IDS spectra cover the wavelength range ~7650–9300 Å at R ~ 9300, and were obtained between 7 and 12 March 2012. They were reduced, extracted and wavelength calibrated in the same way as the ISIS spectra. The ISIS and FLAMES spectra already gave good phase coverage of the orbit. With IDS, we obtained one spectrum at quadrature and one during each eclipse in the hope of disentangling the spectra of the two components. However, the S/N of the IDS spectra turned out to be insufficient for this purpose, so we did not include them in the rest of the analysis.

# 3. Analysis

## 3.1. Light curve modelling

### 3.1.1. Out of eclipse variability removal

The target displays significant out-of-eclipse (OOE) flux variations, which must be accounted for when modelling the CoRoT

<footnote n="4">⁴ http://www.eso.org/sci/software/pipelines/giraffe/giraf-pipe-RECIPES.html</footnote>
light curve to determine the fundamental parameters of the two stars. Given the complex nature of the OOE variations, it is not practical to model them at the same time as the eclipses. We therefore opted to remove the OOE variations before modelling the residuals with standard eclipse modelling software. We first tried fitting cubic splines to the OOE data (masking the in-eclipse regions) and interpolating across the eclipses. This gives fairly satisfactory results, but there is no principled way to propagate uncertainties arising from the OOE removal process through to the residuals. We therefore implemented a second approach, based on Gaussian process (GP) regression, which provides both a better fit to the OOE data, and a more robust means of interpolating across the eclipses.

A detailed description of GP regression is beyond the scope of the present paper; we refer the interested reader to Rasmussen & Williams (2006) for a textbook-level introduction, Gibson et al. (2012) for a relatively detailed description in the context of astrophysical time-series data, and Aigrain et al. (2012) for specific examples of GP regression applied to stellar light curves. For the present purpose, it is sufficient to think of GPs as a means of modelling the light curve by parameterising the covariance between pairs of flux measurements, rather than writing down an explicit expression for the fluxes themselves. The joint distribution of the observed fluxes is then taken to be a multi-variate Gaussian, with a covariance matrix whose elements depend on the observation times through a covariance function. After experimenting with a range of widely-used covariance functions, we opted for a member of the Matérn family with smoothness parameter $\nu = 3/2$. This is appropriate for data displaying a relatively rough behaviour (Rasmussen & Williams 2006), such as we observe in the OOE light curve (see Fig. 2). The covariance $k$ between the fluxes observed at times $t$ and $t'$ is modelled as:

$$k_{3/2}(r) = A^2 \left( 1 + \frac{\sqrt{3}r}{l} \right) \exp \left( -\frac{\sqrt{3}r}{l} \right) + \sigma^2 \delta(r) \quad (1)$$

where $r = |t-t'|$ is the time-interval between the observations and $\delta(x)$ is the Kronecker delta function. The first term represents the OOE variations, with amplitude $A$ and characteristic time scale $l$, and the second term represents white noise with standard deviation $\sigma$. For a given set of parameters, the likelihood of the model, marginalised over all the possible flux vectors which share the same covariance matrix, can be estimated analytically. It requires an inversion of the covariance matrix, which is computationally expensive, but nonetheless feasible for up to a few thousand data points. In this analysis, we opt to subsample the data to speed up this inversion. We bin sets of ~8 measurements, i.e. retaining information on much shorter timescales than the OOE variations, and take the median and mean values for the flux and time respectively. Starting with an initial guess for the OOE and white noise parameters (obtained from visual inspection of the light curve, and from photon counting statistics, respectively), we use a Nelder-Mead optimiser to find the values which maximise the likelihood. These are $A = 11.6\%$, $l = 1.25$ days and $\sigma = 1.5$ mmag. We note that the white noise estimate is comparable to the estimated photon noise of $\sigma = 1.1$ mmag for the residuals (subsampled to the same frequency as used in the OOE modelling). Once the parameters are fixed to these values, we can then compute a predictive distribution for the fluxes at any given set of times. We do this for the times of all observations (both in and out of eclipse), and obtain a light curve corrected for the OOE variations by subtracting the mean of that distribution from the observed values (Fig. 2, bottom panel). The standard

deviation of the predictive distribution provides an estimate of the uncertainty on the corrected flux.

Once we have a satisfactory model for the OOE variations, we must decide whether to subtract it from the original, or divide the latter by the OOE model, before fitting the eclipses. The appropriate course of action depends on how the OOE variability affects an eclipse. For example, spots on the background star will tend to have a multiplicative effect across an eclipse, while those on the foreground star will have an additive effect. With eclipses on both stars, and with both stars likely to be spotted, the effect across all eclipses will not be the same. In addition, if some of the OOE variability arises from obscuration of one or both stars by material located outside the binary orbit (see Sect. 4.5), then this complicates matters further. Neither simply subtracting nor dividing will fully account for all the effects caused by the OOE variations. It may be feasible to disentangle the different effects using simultaneous light curves in widely separated bandpasses, but this is not possible using the CoRoT light curve alone. We therefore simply tested both approaches to remove the OOE variability. Subtracting gives slightly better results, in the sense that the depths of both primary and secondary eclipses are more consistent from one orbit to the next. We use the OOE-subtracted light curve when subsequently fitting the eclipses.

### 3.1.2. JKTEBOP

Analysis of EB light curves yields a wealth of physical information, from orbital geometry to relative stellar parameters. Light curve modelling was performed with JKTEBOP (Southworth et al. 2004, 2007), an extension of the EBOP code (Eclipsing Binary Orbit Program; e.g. Popper & Etzel 1981; Etzel 1981; Nelson & Davis 1972). JKTEBOP models each star as a sphere and computes light curves through numerical integration of concentric circles over each star, employing a Levenberg-Marquardt optimisation algorithm in finding the best fit. The approximation of modelling each star as a sphere is valid in the case of well-detached systems with modest tidal distortion.

The parameters adjusted during modelling are: the central surface brightness ratio, $J = SB_{sec}/SB_{pri}$; the sum of the radii as a fraction of the orbital separation, $(R_{pri}+R_{sec})/a$; the radius ratio, $R_{sec}/R_{pri}$; the orbital inclination, $i$; the orbital period, $P$; the time of primary eclipse centre, $T_{prim}$; the orbital eccentricity and the longitude of periastron, $e$ and $\omega$, in the form of the combination terms, $e \cos \omega$ and $e \sin \omega$. The initial guesses for these parameters were derived from the output of the eclipse search, from visual estimates of the ratio of eclipse depths, their durations and relative phases, and from blackbody approximations relating radius, effective temperature, luminosity and surface brightness. We checked that the final results were not sensitive to the initial guesses (provided these were sufficiently close to the final values for JKTEBOP to converge). Our GP treatment of the OOE variability removes reflection and gravity darkening effects from the light curve, and we therefore do not incorporate these into our modelling.

JKTEBOP models single-band light curves and unlike other techniques to determine eclipsing binary parameters, e.g. PHOEBE (Physics of Eclipsing Binaries; Prša & Zwitter 2005), based on the Wilson-Devinney approach (Wilson & Devinney 1971), does not require the use of model atmospheres. Effective temperatures and surface gravities therefore are only introduced in the determination of the limb darkening coefficients. We used a quadratic limb-darkening law with coefficients from ATLAS model spectra for the CoRoT bandpass (Sing 2010). The coefficients were specified from estimates of $T_{eff}$, $\log g$
and [M/H]: effective temperatures were estimated from tem-
plate spectra giving the highest stellar peaks in cross-correlation
(see Sect. 3.2), log *g* was set at 4.0 (cgs) for both stars (reason-
able for low-mass PMS stars) and metallicity was taken to be
[M/H] = -0.1, the closest value available to the cluster metallic-
ity (King et al. 2000). We initially allowed the limb darkening
parameters to vary but found that this yielded unphysical values
for the parameters, while it did not significantly change the re-
sults for the other parameters. This is because the eclipses are
grazing and the light curve simply does not contain enough in-
formation to constrain the limb darkening parameters. We there-
fore fixed them to the theoretical values in the final fit.

The fact that the eclipses are grazing also gives rise to a de-
generacy between the surface brightness ratio and the radius ra-
tio. To break this degeneracy, we added an additional constraint
in the form of a light ratio, estimated from the relative cross-
correlation function (CCF) peak heights. We used a FLAMES
spectrum for this purpose, because the wavelength coverage of
FLAMES is more similar to that of CoRoT than either ISIS
or IDS, and so the relative CCF peak heights are expected to
be more representative of the light ratio of the two stars in the
CoRoT bandpass. We chose a single CCF (HJD = 2 455 940.64,
see Fig. 8, left panel) with strong, well-separated, stellar peaks,
and obtained a peak height ratio of 0.641 (see Sect. 3.2 for de-
tails of the CCF fitting procedure).

As explained in Sect. 2.1, we rely primarily on the 2008
light curve to model the eclipses, because the spatially depen-
dent background in the 2011/2012 observation affects the eclipse
depths (the differences between the 2008 and the 2011/2012
depths is significantly larger than the small variation in the ap-
parent eclipse depths during either run). However, we obtain the
best-fit period and time of primary eclipse centre ($T_{prim}$) from
a combined fit to both light curves (after removing the OOE
variations in the second run, in same way as for the first). The
background correction error should not affect these parameters
and the extended time coverage gives significantly improved
accuracy.

Figure 6 shows the phase folded light curve along with the best-fit JKTEBOP model. The scatter of the residuals is significantly larger in eclipse than out of eclipse. This is partly due to occultation by the foreground star of starspots on the background star, as seen in numerous EB and transiting planet systems. In the present case, it is also due to increased uncertainty in the GP prediction across the eclipses, and to the residual variations in the depths of individual eclipses, already noted in Sect. 3.1.1, and which could be caused by changes in the global spot coverage of either star between one eclipse and the next. To account for this increased in-eclipse scatter, we performed a preliminary JKTEBOP run on the detrended, 2008 light curve (Fig. 2, bottom panel), and used the resulting fit to compute the reduced $\chi^2$ for the primary, secondary and out-of eclipse data separately. The photometric uncertainties on the fluxes were then rescaled so as to give a reduced $\chi^2$ of unity in each subset; the rescaling factors were 1.55 (primary), 2.68 (secondary) and 1.07 (out-of eclipse). JKTEBOP was then re-run (with the same initial guesses), using the scaled photometric errors.

The resulting best-fit parameters are given in Table 4, along
with uncertainties derived from a Monte Carlo analysis. This
involves generating a model light curve from the best-fit pa-
rameters, adding Gaussian white noise (matching the obser-
vational errors), fitting the result, and repeating the proce-
dure 10 000 times to obtain distributions for each parameter,
which are shown in Fig. 7. This procedure highlights the de-
generacies between some parameters, most notably the radii,

<figure><img src="image_7.png" /><figcaption>**Fig. 6.** Top: phase folded, detrended CoRoT light curve (black points, 2008 observation) with the JKTEBOP best-fit model shown in cyan. The residuals of the best-fit model are also shown, with a vertical offset for clarity. Phase zero marks the centre of the primary eclipse. The bottom panels show zooms on the primary and secondary eclipses (left and right, respectively).</figcaption></figure>

luminosity ratio and inclination, as expected for a grazing sys-
tem, and ensures they are accounted for in the reported uncer-
tainties. Interestingly, the use of the light ratio constraint from
the FLAMES spectra enables us to break the usual degener-
acy between the radius and surface brightness ratios, which is
often the limiting factor in the final radius estimates for near
equal-mass EBs (see e.g. the case of JW 380, Irwin et al. 2007).
We note that the sum of the radii is roughly a fifth of the or-
bital separation and that the expected oblateness of each star is
0.0025 (primary) and 0.0016 (secondary), well within the allow-
able range for JKTEBOP (the EBOP model breaks down around
an oblateness of 0.04; Popper & Etzel 1981). This validates our
initial assumption that the stars are well-detached and can be
treated as spherical. The best-fit model has an eccentricity con-
sistent with zero at <1.5σ level.

3.2. Radial velocities

We measured radial velocities by cross-correlating the spectra
described in Sect. 2.2 with MARCS theoretical model spectra
(Gustafsson et al. 2008).

Before performing the cross-correlation, we masked out re-
gions contaminated by telluric emission. We also masked the
calcium infrared triplet in the ISIS spectra, as the lines are very
broad and tend to smear out the resulting CCF. The remain-
ing wavelength intervals used were: 8380–8392, 8455–8484,
8510–8525, 8575–8644 and 8685–8745 Å (see Fig. 5). There
were many emission lines present in the FLAMES spectra, not
only from the system but also the nebula. These emission lines
were not present in the model spectra and inhibited clean cross-
correlation. We therefore clipped and interpolated a few data
points before and after each emission line to try and fully remove
<table><caption>**Table 3.** Radial velocities derived from the WHT/ISIS and VLT/FLAMES spectra.</caption><thead><tr><th rowspan="2">HJD</th><th rowspan="2">Orbital phase</th><th rowspan="2">S/N<br>(/pixel)</th><th colspan="3">Primary</th><th colspan="3">Secondary</th></tr><tr><th>C</th><th>W (km s<sup>-1</sup>)</th><th>RV (km s<sup>-1</sup>)</th><th>C</th><th>W (km s<sup>-1</sup>)</th><th>RV (km s<sup>-1</sup>)</th></tr></thead><tbody><tr><td colspan="9" style="text-align:center; font-weight:bold;">ISIS</td></tr><tr><td>2455 899.51482</td><td>0.716</td><td>~14</td><td>0.36</td><td>33.3 ± 1.8</td><td>+73.6 ± 1.9</td><td>0.28</td><td>35.4 ± 2.5</td><td>-68.7 ± 2.7</td></tr><tr><td>2455 899.73483</td><td>0.773</td><td>~13</td><td>0.25</td><td>36.9<sup>+2.5</sup><sub>-3.9</sub></td><td>(+83.4 ± 3.9)</td><td>0.08</td><td>94<sup>+537</sup><sub>-78</sub></td><td>(-57<sup>+73</sup><sub>-236</sub>)</td></tr><tr><td>2455 900.55487*</td><td>0.985</td><td>~19</td><td></td><td></td><td></td><td>0.69</td><td>30.9 ± 1.0</td><td>(+16.5 ± 1.0)</td></tr><tr><td>2455 901.50491</td><td>0.230</td><td>~18</td><td>0.37</td><td>32.7 ± 1.2</td><td>(-40.6 ± 1.4)</td><td>0.28</td><td>31.2 ± 1.6</td><td>(+96.2 ± 1.8)</td></tr><tr><td>2455 901.57491</td><td>0.248</td><td>~25</td><td>0.60</td><td>32.7 ± 1.1</td><td>-50.0 ± 1.1</td><td>0.42</td><td>35.0 ± 2.0</td><td>+95.0 ± 1.8</td></tr><tr><td>2455 901.69492</td><td>0.279</td><td>~24</td><td>0.56</td><td>32.8 ± 1.3</td><td>-49.2 ± 1.2</td><td>0.42</td><td>33.3 ± 1.7</td><td>+93.2 ± 1.6</td></tr><tr><td>2455 901.77492</td><td>0.300</td><td>~21</td><td>0.51</td><td>34.8 ± 1.3</td><td>-49.0 ± 1.2</td><td>0.43</td><td>34.8 ± 1.6</td><td>+88.7 ± 1.4</td></tr><tr><td colspan="9" style="text-align:center; font-weight:bold;">FLAMES</td></tr><tr><td>2455 915.65584</td><td>0.882</td><td>~26</td><td>0.36</td><td>20.8 ± 1.8</td><td>+60.5 ± 1.6</td><td>0.21</td><td>16.1 ± 1.7</td><td>-38.1 ± 1.7</td></tr><tr><td>2455 917.70423</td><td>0.411</td><td>~22</td><td>0.37</td><td>19.9 ± 1.3</td><td>-13.2 ± 1.3</td><td>0.17</td><td>34.2<sup>+6.4</sup><sub>-5.4</sub></td><td>+66.5 ± 5.3</td></tr><tr><td>2455 918.68222</td><td>0.663</td><td>~22</td><td>0.38</td><td>17.6 ± 1.1</td><td>+74.9 ± 1.2</td><td>0.15</td><td>16.1<sup>+2.1</sup><sub>-2.8</sub></td><td>-46.4 ± 2.7</td></tr><tr><td>2455 922.69870</td><td>0.700</td><td>~22</td><td>0.40</td><td>18.2 ± 1.1</td><td>+75.3 ± 1.2</td><td>0.20</td><td>17.9<sup>+2.1</sup><sub>-2.5</sub></td><td>-53.1 ± 2.4</td></tr><tr><td>2455 939.65778</td><td>0.077</td><td>~22</td><td>0.37</td><td>19.0 ± 1.6</td><td>-8.1 ± 1.2</td><td>0.18</td><td>16.7 ± 1.6</td><td>(+52.6 ± 2.2)</td></tr><tr><td>2455 940.64342</td><td>0.331</td><td>~34</td><td>0.38</td><td>17.6 ± 1.1</td><td>-33.2 ± 1.1</td><td>0.24</td><td>16.1 ± 1.4</td><td>+89.9 ± 1.5</td></tr><tr><td>2455 940.67976</td><td>0.341</td><td>~34</td><td>0.41</td><td>18.4 ± 1.1</td><td>-31.1 ± 1.1</td><td>0.27</td><td>16.2 ± 1.2</td><td>+88.2 ± 1.4</td></tr><tr><td>2455 941.66163</td><td>0.594</td><td>~35</td><td>0.38</td><td>18.4 ± 1.2</td><td>+50.9 ± 1.0</td><td>0.22</td><td>24.2 ± 3.0</td><td>(-19.4 ± 2.6)</td></tr><tr><td>2455 943.68750</td><td>0.117</td><td>~23</td><td>0.38</td><td>18.5 ± 1.3</td><td>-21.6 ± 1.2</td><td>0.19</td><td>21.4 ± 4.0</td><td>+76.3 ± 3.0</td></tr><tr><td>2455 944.63515</td><td>0.362</td><td>~32</td><td>0.36</td><td>19.4 ± 1.8</td><td>-24.9 ± 1.5</td><td>0.22</td><td>17.1 ± 1.9</td><td>+84.4 ± 2.1</td></tr><tr><td>2455 945.66720</td><td>0.628</td><td>~36</td><td>0.34</td><td>16.7 ± 1.2</td><td>+62.0 ± 1.2</td><td>0.18</td><td>16.8 ± 1.9</td><td>-39.8 ± 2.2</td></tr><tr><td>2455 946.68154</td><td>0.890</td><td>~22</td><td>0.36</td><td>20.7 ± 1.6</td><td>+56.9 ± 1.4</td><td>0.23</td><td>19.3 ± 2.3</td><td>-29.1 ± 1.9</td></tr><tr><td>2455 977.63971</td><td>0.880</td><td>~30</td><td>0.33</td><td>19.4 ± 1.8</td><td>+62.1 ± 1.4</td><td>0.19</td><td>20.4<sup>+3.8</sup><sub>-6.0</sub></td><td>-34.7<sup>+3.1</sup><sub>-4.2</sub></td></tr><tr><td>2455 979.57946</td><td>0.380</td><td>~31</td><td>0.40</td><td>18.7 ± 1.0</td><td>-21.9 ± 1.0</td><td>0.17</td><td>26.0<sup>+6.0</sup><sub>-4.8</sub></td><td>(+81.4 ± 4.8)</td></tr><tr><td>2455 981.61509</td><td>0.906</td><td>~28</td><td>0.40</td><td>18.1 ± 1.0</td><td>+50.5 ± 1.0</td><td>0.24</td><td>24.2 ± 2.6</td><td>-21.1 ± 2.5</td></tr></tbody></table>

Notes. RV values in brackets were not used in deriving the orbital solution, for reasons detailed in the text. C and W are the height and standard deviation of each peak, respectively. (*) Only one stellar peak was resolved in the CCF because the spectra was taken close to primary eclipse.

their effect. We also masked the $H\alpha$ emission region (6545–6587 Å). The regions used for cross-correlation of the FLAMES spectra were: 6521–6545, 6587–6598, 6610–6642, 6648–6700 and 6702–6713 Å. After experimenting with a range of parameters for the MARCS model spectra, we opted to use $T = 4000$ K and $\log g = 4$, which gave the best CCF contrast. We did not broaden the template, as that results in broader CCF features which makes it more difficult to separate the peaks corresponding to the two stars. The CCFs were computed and fitted using purpose written Python code.

We initially performed a simple least-squares fit to the largest peak(s) in the CCFs, using a model consisting of either one or two Gaussians plus a constant offset. However, we found that the resulting RVs were significantly affected by the pronounced correlated noise in the CCFs (examples of which are shown in Fig. 8). This is a major issue when measuring RVs via cross-correlation, particularly for late-type stars. This correlated noise arises in part from noise in the object spectrum, but also – more critically – from mismatch between the object and template spectra. A number of approaches have been developed to account for it when evaluating the uncertainties. For example, Tonry & Davis (1979) decompose the CCF into components which are symmetric and anti-symmetric with respect to the peak, and use the root mean square (RMS) scatter

of the antisymmetric component to estimate the CCF noise, on the basis that any “real” CCF signal should be symmetric about the peak. This is the most widely used approach, and it is implemented in the IRAF cross-correlation package FCOR. Working with VLT/FLAMES spectra of moderately faint stars from the OGLE survey, Bouchy et al. (2005) used the relation $\sigma_{RV} = 3\sqrt{FWHM/SC}$, where $\sigma_{RV}$ is the RV uncertainty, C and FWHM are the height and full-width at half-maximum of the CCF peak, respectively, and S is the signal-to-noise of the object spectrum. The form of this relation was deduced from photon-noise considerations, and it was calibrated on multiple observations of a large sample of stars, most of which were presumed to be non-variable. In practice, both of these methods give good results for the objects they were initially designed for, but they tend to systematically under-predict the uncertainties for late-type, rapidly rotating stars.

We therefore decided to model the CCF noise at the same time as fitting for the CCF peaks. Each CCF was modelled as the sum of two Gaussians plus a stochastic noise term, which is described by a GP. The choice of GP kernel was based on a careful examination of the CCFs such as those shown in Fig. 8. Away from the main peak(s), the CCFs display variations on both moderate velocity scales (few tens of m/s) and long velocity scales (few 100 m/s). These cannot be adequately described by a single
<table><caption>**Table 4.** Fitted and derived parameters of the system.</caption><thead><tr><th>Parameter</th><th>Symbol</th><th>Unit</th><th>Value</th></tr></thead><tbody><tr><td colspan="4" style="text-align: center;">JKTEBOP light curve fit</td></tr><tr><td>Central surface brightness ratio</td><td>J</td><td></td><td>0.871<sup>+0.037</sup><sub>-0.035</sub></td></tr><tr><td>Sum of radii</td><td>(R<sub>pri</sub> + R<sub>sec</sub>)/a</td><td></td><td>0.2198<sup>+0.0017</sup><sub>-0.0018</sub></td></tr><tr><td>Radius ratio</td><td>R<sub>sec</sub>/R<sub>pri</sub></td><td></td><td>0.854<sup>+0.058</sup><sub>-0.061</sub></td></tr><tr><td>Fractional primary radius</td><td>R<sub>pri</sub>/a</td><td></td><td>0.1186<sup>+0.0036</sup><sub>-0.0033</sub></td></tr><tr><td>Fractional secondary radius</td><td>R<sub>sec</sub>/a</td><td></td><td>0.1013<sup>+0.0039</sup><sub>-0.0045</sub></td></tr><tr><td>Orbital inclination</td><td>i</td><td>(°)</td><td>85.09<sup>+0.16</sup><sub>-0.11</sub></td></tr><tr><td>Orbital period</td><td>P</td><td>(days)</td><td>3.8745746 ± 0.0000014</td></tr><tr><td>Time of primary eclipse centre</td><td>T<sub>prim</sub></td><td>(HJD)</td><td>2 454 536.76357 ± 0.00043</td></tr><tr><td></td><td>e cos ω</td><td></td><td>0.00050<sup>+0.00029</sup><sub>-0.00028</sub></td></tr><tr><td></td><td>e sin ω</td><td></td><td>-0.0049<sup>+0.0077</sup><sub>-0.0075</sub></td></tr><tr><td colspan="4" style="text-align: center;">Photometrically constrained RV fit</td></tr><tr><td>Primary semi-amplitude</td><td>K<sub>pri</sub></td><td>(km s<sup>-1</sup>)</td><td>60.49 ± 0.39</td></tr><tr><td>Secondary semi-amplitude</td><td>K<sub>sec</sub></td><td>(km s<sup>-1</sup>)</td><td>81.56 ± 0.62</td></tr><tr><td>Systemic velocity</td><td>V<sub>sys</sub></td><td>(km s<sup>-1</sup>)</td><td>19.42 ± 0.26</td></tr><tr><td>Orbital period</td><td>P</td><td>(days)</td><td>3.8745745 ± 0.0000014</td></tr><tr><td>Time of primary eclipse centre</td><td>T<sub>prim</sub></td><td>(HJD)</td><td>2 454 536.76355<sup>+0.00042</sup><sub>-0.00043</sub></td></tr><tr><td></td><td>e cos ω</td><td></td><td>0.00049<sup>+0.00028</sup><sub>-0.00027</sub></td></tr><tr><td></td><td>e sin ω</td><td></td><td>-0.0033 ± 0.0040</td></tr><tr><td colspan="4" style="text-align: center;">Derived parameters</td></tr><tr><td>Luminosity ratio</td><td>L<sub>sec</sub>/L<sub>pri</sub></td><td></td><td>0.642<sup>+0.092</sup><sub>-0.091</sub></td></tr><tr><td>Orbital eccentricity</td><td>e</td><td></td><td>0.0037<sup>+0.0036</sup><sub>-0.0025</sub></td></tr><tr><td>Semi-major axis</td><td>a</td><td>(R<sub>☉</sub>)</td><td>10.921 ± 0.056</td></tr><tr><td>Primary mass</td><td>M<sub>pri</sub></td><td>(M<sub>☉</sub>)</td><td>0.668<sup>+0.012</sup><sub>-0.011</sub></td></tr><tr><td>Secondary mass</td><td>M<sub>sec</sub></td><td>(M<sub>☉</sub>)</td><td>0.4953<sup>+0.0073</sup><sub>-0.0072</sub></td></tr><tr><td>Primary radius</td><td>R<sub>pri</sub></td><td>(R<sub>☉</sub>)</td><td>1.295<sup>+0.040</sup><sub>-0.037</sub></td></tr><tr><td>Secondary radius</td><td>R<sub>sec</sub></td><td>(R<sub>☉</sub>)</td><td>1.107<sup>+0.044</sup><sub>-0.050</sub></td></tr><tr><td>Primary surface gravity</td><td>(log g)<sub>pri</sub></td><td>(cm s<sup>-2</sup>)</td><td>4.038<sup>+0.025</sup><sub>-0.026</sub></td></tr><tr><td>Secondary surface gravity</td><td>(log g)<sub>sec</sub></td><td>(cm s<sup>-2</sup>)</td><td>4.045<sup>+0.040</sup><sub>-0.033</sub></td></tr><tr><td>Primary semi-major axis<sup>a</sup></td><td>a<sub>pri</sub></td><td>(R<sub>☉</sub>)</td><td>4.651 ± 0.030</td></tr><tr><td>Secondary semi-major axis<sup>a</sup></td><td>a<sub>sec</sub></td><td>(R<sub>☉</sub>)</td><td>6.270<sup>+0.047</sup><sub>-0.048</sub></td></tr></tbody></table>

**Notes.** (a) Relative to centre of mass.

squared exponential covariance function, which only allows for variations on a single characteristic length scale:

$$k_{\text{SE}}(r) = A^2 \exp\left(-\frac{r^2}{2l^2}\right). \qquad (2)$$

Instead, we use the rational quadratic:

$$k_{\text{RQ}}(r) = A^2 \left(1 + \frac{r^2}{2\alpha l^2}\right)^{-\alpha}. \qquad (3)$$

This can be seen as a squared exponential with a certain amount of additional covariance on large scales, controlled by the

parameter $\alpha$ (when $\alpha$ is $\gg 1$, it reduces to the squared exponential), or alternatively as a scale mixture of squared exponential covariance functions with different characteristic length-scales (Rasmussen & Williams 2006). We also incorporate a very small white noise term on the diagonal of the covariance matrix to aid convergence. For the sake of computational efficiency, we model the CCFs only in the range $-580 \le v \le 620$ km s$^{-1}$, and use only every other data point in the fit. We used a Metropolis-Hastings Markov chain Monte Carlo (MCMC) with Gaussian proposal distributions to marginalise over the parameters of the Gaussian terms and of the covariance function. We performed five relatively short-chains of 15 000 steps, which were sufficient
<figure><img src="image_8.png" /><figcaption>**Fig. 7.** 2D contours and 1D histograms of the Monte Carlo chains for selected JKTEBOP parameters (bottom left) and of the MCMC chains for the photometrically constrained radial velocity solution (top right; see Sect. 3.2). On the 2D plots, red crosses show the median values and black contours represent 1, 2 and 3σ confidence intervals. On the histograms, solid and dotted grey vertical lines show the median and ±1σ intervals, respectively.</figcaption></figure>

to achieve convergence in all cases; the first 5000 steps were discarded to minimise sensitivity to the initial guesses, and the parameter distributions were derived from the remaining 10 000 steps. Figure 8 shows example CCFs with the best-fit model obtained in this way.

Not all the CCFs contain well-defined peaks that clearly stand out above the noise. To avoid inferring erroneous RV measurements, we visually inspected every CCF and fit, and discarded any epochs where the contrast of stellar peaks was below ~0.1, or where the properties of the GP component of the fit appeared to change markedly in the region containing the stellar peaks. In a few cases, the best fit to the secondary peak is low compared to the noise, has an unusually large width, and has an RV slightly offset from the expected value (given a preliminary

orbital solution based on the rest of the observations). An example of such a case is shown in the right panel of Fig. 8. In such cases, we discarded the secondary measurement, but retained the primary RV. We note that it might be possible to obtain RV measurements from a larger fraction of our spectra by fixing the widths of the two stellar peaks, but we opted not to do this, as we already have enough useful RVs for a good orbital solution, and this would only add data points with large error bars.

### 3.2.1. Orbital solution

We fit the RV measurements using Keplerian orbits, with free parameters $K_{\text{pri}}$ and $K_{\text{sec}}$ (the semi-amplitudes of the primary and secondary), $V_{\text{sys}}$ (the systemic velocity), $\Delta V_{\text{ISIS}}$ (the offset for
<figure><img src="image_9.png" /><figcaption>**Fig. 8.** Example CCFs (black crosses) obtained from FLAMES spectra taken at HJD = 2 455 940.6 (**left**) and 2 455 979.6 (**right**), showing the best-fit model (cyan) together with the associated 68.5% confidence interval (light cyan shaded region). The Gaussian and GP terms are also shown separately by the red dot-dashed and blue dashed lines respectively, with a vertical offset added for clarity. The **top panel** in each case shows the full range used in the fit and the **bottom panel** a zoom on the region around the peaks.</figcaption></figure>

the ISIS RV points), as well as the period, *P*, the time of primary eclipse centre, $T_{\text{prim}}$, $e \cos \omega$ and $e \sin \omega$. The values and uncertainties obtained from the light curve modelling for the last four of these parameters were used as priors in the RV fit. We used a Metropolis-Hastings MCMC with Gaussian proposal distributions to estimate the uncertainties on the parameters of the RV fit and performed 5 chains of $10^5$ steps each to ensure convergence. Figure 9 shows the resulting orbital fit, whose parameters are listed in Table 4. The distributions from the MCMC are shown in Fig. 7 and display negligible degeneracy.

We find that our systemic velocity of $V_{\text{sys}} = 19.42 \pm 0.26 \text{ km s}^{-1}$ is consistent with the literature value of the cluster's recessional velocity, $V = 22 \pm 3.5 \text{ km s}^{-1}$ (Fűrész et al. 2006), suggesting that the system is kinematically associated with the cluster.

### 3.3. Lithium absorption

The presence of lithium (6707.8 Å) absorption in both stars (Fig. 10) indicates youth and therefore membership. We determine equivalent widths of both lines by first fitting a cubic spline through the continuum data (masking out all absorption and emission features) and interpolating across each lithium line. We then determine relative flux contributions from each star (using the CCF peak heights for that spectrum) and scale both the data and the spline, before calculating the equivalent widths. We select three spectra in which the lines are well separated and calculate approximate equivalent widths of 0.56 and 0.55 Å for the primary and secondary, respectively (each individual measurement is mutually consistent with these values to within 0.01 and 0.02 Å, respectively). The equivalent widths are very similar, tentatively suggesting that there has not been significant lithium depletion, as the rate of depletion is expected to be faster in the

<figure><img src="image_10.png" /><figcaption>**Fig. 9. Top:** phase-folded RV data (points with error bars) and best-fit orbital solution for the primary (red solid line) and secondary (blue dashed line). Crosses and circles represent RVs derived from ISIS and FLAMES spectra, respectively. The grey horizontal dotted line shows the systemic velocity. **Bottom:** residuals of the best-fit model.</figcaption></figure>

less massive secondary. They are also consistent with other cluster members (e.g. Sergison et al. 2013), including early/mid-K spectral types, which are expected to still have their initial abundances at this age.
<figure><img src="image_11.png" /><figcaption>**Fig. 10.** Examples of the lithium 6707.8 Å absorption feature for both primary and secondary stars from FLAMES spectra. From top to bottom, HJD = 2 455 915.7, 2 455 922.7, 2 455 944.6 and 2 455 981.6. In each case the predicted positions of the two stars are indicated.</figcaption></figure>

### 3.4. Fundamental parameters

From our photometric and radial velocity analysis, we can calculate masses and radii. We propagate our distributions for the period, semi-amplitudes of the velocities, radius sum, radius ratio and inclination to derive masses and radii. These are presented in the bottom part of Table 4, along with the surface gravities and orbital separations. The masses of $0.668_{-0.011}^{+0.012}$ and $0.4953_{-0.0072}^{+0.0073} M_{\odot}$ are consistent with the combined spectral type derived from the CAHA spectrum. The large radii ($1.295_{-0.037}^{+0.040}$ and $1.107_{-0.050}^{+0.044} R_{\odot}$) are what we would expect for a PMS system.

## 4. Discussion

### 4.1. Comparison to pre-main sequence stellar evolution models

We are now in a position to place the components of CoRoT 223992193 on the mass-radius diagram (Fig. 1, red diamonds). We note that it lies in a very sparsely populated region

of the diagram, highlighting its value in testing and constraining PMS stellar evolution models. In Fig. 11, we compare the masses and radii of both stars to five sets of widely used PMS isochrones: BCAH98 (Baraffe et al. 1998), DM98 (D’Antona & Mazzitelli 1998), SDF00 (Siess et al. 2000), Dartmouth (Dotter et al. 2008) and Pisa (Tognelli et al. 2011). The parameters of the models shown are: helium mass fraction $Y = 0.282$, metallicity [M/H] = 0 relative to solar, and mixing length parameter $\alpha_{ML} = 1.9$ for the BCAH98 models; $Y = 0.28$, metal mass fraction $Z = 0.02$ and deuterium abundance $X_D = 2 \times 10^{-5}$ for the DM98 models; $Y = 0.277$ and $Z = 0.02$ for the SDF00 models; $Y = 0.274$, $Z = 0.0189$, [Fe/H] = [α/Fe] = 0 and $\alpha_{ML} = 1.938$ for the Dartmouth models and $Y = 0.268$, $Z = 0.01$, $X_D = 2 \times 10^{-5}$ and $\alpha_{ML} = 1.2$ for the Pisa models. These were chosen to match the metallicity of NGC 2264. For BCAH98 we selected models with a solar-calibrated mixing length $\alpha_{ML} = 1.9$, which have been extended down to $0.1 M_{\odot}$ (Baraffe, priv. comm.<fnref n="5" />). For the Pisa models we chose $\alpha_{ML} = 1.2$, as models with lower convective efficiency are known to be more consistent with existing observations of low-mass PMS objects (e.g. Stassun et al. 2004; Mathieu et al. 2007).

BCAH98, SDF00, Dartmouth and Pisa all manage to fit both components of CoRoT 223992193 simultaneously. All models under-predict the radius of the primary slightly, compared to that of the secondary. This is a feature which is also seen in other close binaries, and Chabrier et al. (2007) proposed that it could be due to enhanced magnetic activity on the primary. However, in the present case the discrepancy is not statistically significant ($<1.5\sigma$). More precise masses and radii for this system are needed to distinguish between the different models.

The best match to the BCAH98, SDF00, Dartmouth and Pisa isochrones is obtained for ages ~3.5–6 Myr. Ages from all isochrones are consistent with the literature estimates of the age of NGC 2264, given the wide dispersion of the latter. For example, Dahm (2008) report a median age of 3 Myr, but infer an apparent age spread of 5 Myr from the broadened sequence of cluster members, and Naylor (2009) determine a slightly older age of 5.5 Myr. Various other studies have reported similar ages and large dispersions, all of which depend strongly on the choice of models used to fit the low mass stellar population (Park et al. 2000). Rebull et al. (2002) found systematic differences in model-derived ages of up to half an order in magnitude for a spectroscopically classified sample of variable stars. Sequential star formation has also been suggested (Adams et al. 1983), with the peak rate of low-mass ($M < 0.5 M_{\odot}$) star formation preceding that of the higher mass population. On the other hand, Baraffe et al. (2009, 2012) showed that the broadened cluster sequences observed in a number of star forming regions, which are usually interpreted as arising from an age spread, could also be caused by episodic accretion early on in the star formation process. Importantly, they predict that, if most of the accreted energy is radiated away (which is expected for the modest initial core masses that will eventually become low-mass stars), accreting stars should have smaller radii than non-accreting stars of the same age and mass. Differential accretion history could therefore contribute not only to the apparent age spread within star forming regions such as NGC 2264, but also to the discrepancies between the parameters of individual components of binary systems such as CoRoT 223992193.

In addition, it is worth noting that a non-steady accretion history is also expected to cause enhanced lithium depletion, due

<footnote n="5">⁵ The extended models are available at http://perso.ens-lyon.fr/isabelle.baraffe/</footnote>
<figure><img src="image_12.png" /><figcaption>**Fig. 11.** Comparison of measured masses and radii for the two components of CoRoT 223992193 to the predictions of five sets of theoretical PMS models (see text for details). In each subplot the ages of the isochrones increase from top to bottom. The different values obtained when varying the limb-darkening coefficients (see Sect. 4.2) are shown in grey and the final, adopted values by red diamonds.</figcaption></figure>

to higher temperatures at the base of the convective envelope (Baraffe & Chabrier 2010). This is not seen in the present case.

## 4.2. Stellar temperatures and the effect of limb darkening

When modelling the eclipses, we fixed the limb darkening (LD) coefficients to theoretical values from Sing (2010), given an estimate of the effective temperature of each star. We initially chose temperatures of 4000 and 3750 K, which maximised the peak of the CCF for the primary and secondary, respectively. However, our best-fit masses and radii correspond to somewhat lower effective temperatures, of 3670 and 3570 K respectively, according to the BCAH98 isochrones. To check whether this has a significant effect on our analysis, we repeated the eclipse modelling using every combination of temperatures for which Sing (2010) tabulated LD coefficients, such that the $4250$ K $\ge T_{\text{eff,pri}} \ge T_{\text{eff,sec}}$ (LD coefficients were available for 3500, 3750, 4000 and 4250 K). The results are shown in Fig. 11: the values obtained using different combinations of temperatures are shown by the grey symbols, while the adopted values are shown by red diamonds. These correspond to effective temperatures of 3750 K and 3500 K for the primary and secondary respectively, and provide the best internal consistency between the effective temperatures used to fix the LD coefficients and those derived from the measured masses and radii according to the BCAH98 isochrones. We note, however, that the masses and radii themselves are essentially insensitive to the choice of LD coefficients, within the range we examined.

## 4.3. Distance to CoRoT 223992193

Having measured radii and obtained approximate effective temperatures for the two stars, we can evaluate their luminosities and hence the distance to the system. This can then be compared

to published estimates of the distance to NGC 2264, as a further consistency check. We performed this check using both the V and K band magnitudes of the system. The V-band magnitude was obtained by averaging the *ugriz* magnitudes’ approximation (according to the prescription of Jester et al. 2005), and the quoted V-band magnitude from Dahm et al. (2007, V = 16.81). This gives $V = 16.74 \pm 0.1$, where the uncertainty has been inflated to account for the system’s variability.

To derive distances to the system, we used the theoretical bolometric corrections (BCs) of Bessell et al. (1998), as they provide the best sampling in our region of the $T_{\text{eff}}-\log g$ parameter space. Empirical BCs are also available in the literature, but these have been shown to give comparable results to the theoretical BCs (e.g. Southworth et al. 2005). We adopt the values tabulated for solar metallicity, the closest match available to the metallicity of the cluster, and for $\log g = 4.0$, which again provides the best match to our estimate of the surface gravities of both stars. The BCs are tabulated for $T_{\text{eff}} = 3500$ and 3750 K, and are very sensitive to temperature, particularly in the V-band, so we computed a combined bolometric magnitude for the system using either value for each star (always ensuring $T_{\text{eff,pri}} \ge T_{\text{eff,sec}}$). Finally, we carried out the calculation for different amounts of extinction along the line of sight, both to the cluster, corresponding to the different determinations available in the literature (Sung et al. 1997; Rebull et al. 2002; Dahm & Simon 2005), and for CoRoT 223992193, from the CAHA spectrum (Fig. 4). The results are reported in Table 5. Note that although we have matched $T_{\text{eff}}$ and $\log g$ to our estimates for the system, giant stars with the same values would possess different photospheres and so a systematic uncertainty will be present, the effect of which is not well understood.

The V-band extinction estimates in the literature vary widely, mainly because of the different samples of stars used to measure
<table><caption>**Table 5.** Estimates of the distance to CoRoT 223992193 using V and K-band magnitudes, bolometric corrections (BCs; Bessell et al. 1998) for different effective temperature estimates, and a range of extinction estimates (from both the literature and as indicated from the CAHA spectrum, Fig. 4).</caption><thead><tr><th>Extinction<br>(mag)</th><th colspan="2">T<sub>eff</sub> (K)</th><th>Distance<br>(pc)</th></tr><tr><th></th><th>(pri)</th><th>(sec)</th><th></th></tr></thead><tbody><tr><td rowspan="3">A<sub>V</sub> = 0.71<br>(Dahm & Simon 2005)</td><td>3500</td><td>3500</td><td>561 ± 61</td></tr><tr><td>3750</td><td>3500</td><td>627 ± 71</td></tr><tr><td>3750</td><td>3750</td><td>668 ± 73</td></tr><tr><td rowspan="3">A<sub>V</sub> = 0.45<br>(Rebull et al. 2002)</td><td>3500</td><td>3500</td><td>633 ± 69</td></tr><tr><td>3750</td><td>3500</td><td>707 ± 80</td></tr><tr><td>3750</td><td>3750</td><td>754 ± 83</td></tr><tr><td rowspan="3">A<sub>V</sub> = 0.25<br>(Sung et al. 1997)</td><td>3500</td><td>3500</td><td>694 ± 76</td></tr><tr><td>3750</td><td>3500</td><td>775 ± 88</td></tr><tr><td>3750</td><td>3750</td><td>826 ± 91</td></tr><tr><td rowspan="3">A<sub>V</sub> = 0.0<br>(CAHA spectrum)</td><td>3500</td><td>3500</td><td>776 ± 85</td></tr><tr><td>3750</td><td>3500</td><td>868 ± 98</td></tr><tr><td>3750</td><td>3750</td><td>924 ± 101</td></tr><tr><td rowspan="3">A<sub>K</sub> = 0.073<br>(Rebull et al. 2002)</td><td>3500</td><td>3500</td><td>766 ± 85</td></tr><tr><td>3750</td><td>3500</td><td>743 ± 82</td></tr><tr><td>3750</td><td>3750</td><td>728 ± 81</td></tr><tr><td rowspan="3">A<sub>K</sub> = 0.0<br>(CAHA spectrum)</td><td>3500</td><td>3500</td><td>791 ± 87</td></tr><tr><td>3750</td><td>3500</td><td>768 ± 84</td></tr><tr><td>3750</td><td>3750</td><td>751 ± 83</td></tr></tbody></table>

them: the lowest value of $A_V = 0.25$ was obtained by Sung et al. (1997) using O and B main sequence stars, and the highest ($A_V = 0.71$) by Dahm & Simon (2005) using lower mass Hα emitters which may well have dustier environments. Rebull et al. (2002) measured an intermediate value of $A_V = 0.45$, as well as $A_K = 0.073$, using a spectroscopically selected sample of K and M stars. As the K-band distance estimates are less sensitive to reddening, and the BCs in that band less temperature-dependent, we consider the distance estimates we obtain using the K-band magnitude of the system and effective temperatures of 3750 K for the primary and 3500 K for the secondary, to be the most reliable. We consider the full distance range given by the Rebull et al. (2002) determination of $A_K$ and that inferred from the CAHA spectrum.

This places the system at a distance of 756 ± 96 pc, in good agreement with the distances reported for the cluster by Sung et al. (1997) (760 ± 90 pc from isochrone fitting to main sequence B stars) and Sung & Bessell (2010) (815 ± 95 pc, from fitting the spectral energy distribution of individual members using the models of Robitaille et al. 2007). Baxter et al. (2009) obtained a somewhat larger estimate of 913 ± 40 ± 110 pc (sampling and systematic errors respectively) by comparing projected rotational velocities and rotational periods, assuming an isotropic inclination distribution. This estimate relies heavily on $v \sin i$ measurements which are notoriously difficult using low signal-to-noise spectra, but is still within $2\sigma$ of our value for CoRoT 223992193. Thus we can confirm that the system's magnitude and fundamental parameters are consistent with cluster membership.

## 4.4. Spot modelling

In this section, we investigate the possibility that the OOE variations in the 2008 CoRoT light curve could be caused by starspots. Young low-mass stars tend to be heavily spotted, and

<figure><img src="image_13.png" /><figcaption>**Fig. 12.** Top: CoRoT 2008 light curve (black crosses) with eclipses removed, showing three versions of our toy spot model (fit to the second half of the light curve only). Red solid line: $T_{spot} = 3000$ K, angular radius =55° (fraction of stellar surface covered =0.21) and latitude =70°. Green dotted line: $T_{spot} = 2500$ K, angular radius =50° (fraction of stellar surface covered =0.18) and latitude =70°. Blue dashed line: $T_{spot} = 3250$ K, angular radius =60° (fraction of stellar surface covered =0.25) and latitude =65°. Bottom: the red spot model with the period set to the binary orbital period (3.875 days).</figcaption></figure>

the OOE variations seen in the second half of the light curve, although unusually large in amplitude, have a sinusoidal shape that is not atypical of spotted stars. At the age of NGC 2264, low-mass stars in close binary systems are expected to be synchronised (due to tidal effects) up to orbital periods of ~7 days (e.g. Mazeh 2008), and indeed the OOE variations in the second half of the light curve appear to be in phase with the eclipses. A Lomb-Scargle periodogram (Horne & Baliunas 1986) analysis of the second half of the OOE light curve yields a best period of 3.644 days. The associated formal $1\sigma$ error is 0.4 days, measured as being the $1\sigma$ width of a Gaussian fit to the periodogram peak. The OOE period thus appears to be consistent, to within $1\sigma$, with the orbital period of 3.875 days. However, fitting spot models to the OOE light curve, as shown in Fig. 12, reveals that the former period yields much better results than the latter. While the phase of successive minima and maxima are well reproduced using a period of 3.644 days, there is a clear phase drift between the model and the observations when using $P = 3.875$ days.

As a preliminary test of whether starspots could explain the OOE variability, we set up a toy model consisting of a single, large spot located on the primary, facing the secondary, such that flux drops are observed around primary eclipse, as seen in the light curve. The secondary is assumed to be unspotted, and we simply subtract its contribution (39%) from the total flux. We adopt linear limb-darkening coefficients from Sing (2010), assuming the same temperature for the primary's photosphere ($T_{pri} = 3750$ K) as in the eclipse modelling. We fit three models, with different spot temperatures (2500 to 3250 K), sizes (60 to 50°) and latitudes (70 to 65°), to the second half of the light curve. As can be seen in Fig. 12, each model is able to reproduce the large scale structure of the second half of the light curve.
<figure><img src="image_14.png" /><figcaption>**Fig. 13.** Spectral energy distribution of CoRoT 223992193 (cyan points) based on the magnitudes listed in Table 2, plus upper limits in the far-IR (cyan arrows). Left: the grey line and magenta triangles show the best-fit two naked photospheres model. Right: the black line and magenta triangles show the best-fit two naked photospheres model with a small amount of hot dust in the inner cavity of the circumbinary disk (see text for details). The stellar and hot dust terms are shown separately as the grey solid and dashed green lines, respectively. Also shown for completeness, but not used in the fit, is the expected emission from a razor-thin circumbinary disk extending down to $22 R_\odot$ (brown dot-dashed line), which is illuminated by the central star and heated by the gravitational potential energy released from accretion with $\dot{M} = 10^{-11} M_\odot \text{yr}^{-1}$.</figcaption></figure>

The spots are very large, typically ~20% of the primary surface, but this is not unusual for young active stars; this active region should more realistically be interpreted as an assembly of smaller spots covering a large fraction of the stellar surface near the pole.

Although the bulk of the large-scale modulations seen in the second half of the light curve can be reproduced by large, cold spots, rapid spot evolution would be required between the first and second halves of the light curve: the active region would essentially have to appear in less than a day. This seems implausible given its large size (e.g. Grankin et al. 2008). In addition, the small scale, short-timescale variability seen throughout the light curve cannot be accounted for by starspots, unless again they evolve very rapidly. Comparing to weak-lined T Tauri stars in the 2008 CoRoT observation that show sinusoidal variability, i.e. which should only have spot variability, it is rare to see spot evolution on timescales less than a few weeks. In a few cases we do see evolution, but it is in the form of slow amplitude changes, not erratic day timescale changes. While we cannot rule out spots as a source of variability, it is clear that they cannot be solely responsible for all the variability observed.

## 4.5. Spectral energy distribution

As a likely member of NGC 2264, the system is presumably no more than a few-to-several Myr old. It is therefore natural to ask whether there is any evidence for circumstellar or circumbinary material in the spectral energy distribution (SED) of the system. This material could conceivably contribute to the large-amplitude, rapidly evolving, out-of-eclipse variability in the CoRoT light curve, particularly as we have seen in Sect. 4.4 that it is difficult to explain this variability using starspots alone. Indeed, the OOE variability of CoRoT 223992193 is similar to that of classical T Tauri stars in the same region (Alencar et al. 2010). It is also, to a certain extent, reminiscent of the

quasi-periodic flux variations seen in AA Tau (e.g. Bouvier et al. 1999, 2003, 2007), which are attributed to occultations of the central star by a warped circumstellar disk. We therefore constructed an SED for CoRoT 223992193 using the magnitudes listed in Table 2. We also sought evidence of emission at longer wavelengths from archival data, namely *Spitzer*/MIPS 24 µm<fnref n="6" /> and *Herschel*/PACS 70 and 160 µm<fnref n="7" />. The nebula is prominent at these wavelengths (increasing with wavelength) and displays significant structure. Unfortunately, this system lies on the edge of a nebula filament making detection and analysis more difficult. We performed aperture photometry with APT (Aperture Photometry Tool)<fnref n="8" /> and found no clear evidence of emission associated with this system. At the PACS wavelengths, the nebula is very bright, and while there is emission above the sky level at CoRoT 223992193’s location, its profile is not Gaussian, it is not centred on the system and its structure largely follows the nebula. At MIPS 24 µm, the nebula is less prominent but the emission observed follows the structure of the nebula, and while there is a hint of a slightly larger flux at the system’s location, it is not statistically significant given the scatter in the sky flux. A more detailed analysis, accounting for the structured nebula emission would be necessary to determine whether this small excess is real. We therefore attribute the emission at the system’s location to the nebula in all cases and compute an upper limit for the flux of CoRoT 223992193 at each wavelength by taking the “nebula” emission (above that of the median sky level) and quoting the 3σ upper limit on this value. These upper limits are: 4.7, 35.2 and 1916.3 mJy at 24, 70 and 160 µm respectively. The resulting SED is shown in Fig. 13.

<footnote n="6">⁶ *Spitzer* Heritage Archive, http://sha.ipac.caltech.edu/applications/Spitzer/SHA/</footnote>

<footnote n="7">⁷ *Herschel* Science Archive, http://herschel.esac.esa.int/Science_Archive.shtml</footnote>

<footnote n="8">⁸ http://www.aperturephotometry.org/</footnote>
<figure><img src="image_15.png" /><figcaption>Fig. 14. Schematic representation of the proposed system geometry, showing distance from the centre of mass against height above the system plane (z). The primary and secondary stars, along with their circumstellar disks (truncated to a third of the binary separation), are shown in red and blue respectively (the circumstellar disks are shown for completeness but there is no direct evidence for their presence). The sizes and separations of the two stars are to scale. The circumbinary disk (brown) has its inner radius truncated at twice the binary separation. The green dots indicate the general location of dust lying within the inner cavity of the circumbinary disk, such as one could expect to find in accretions streams.</figcaption></figure>

A preliminary inspection of the SED suggests that there is a moderate amount of excess flux in the mid-IR. We confirmed this by fitting a stellar photosphere model to the Sloan and 2MASS fluxes, and comparing the best fit model to the observed Spitzer fluxes. To do this we constructed a grid of two-photosphere models from pairs of MARCS spectra with $2500 \le T_{\text{eff}} < 8000$ K, $\log g = 4.0$ and solar metallicity, ensuring that $T_{\text{eff,pri}} \ge T_{\text{eff,sec}}$, and that the radii of the two stars were within $5\sigma$ of the measured values. We then convolved the model spectra with the bandpass of each filter and optimised the $\chi^2$ of the fit with respect to the temperatures and radii of the two stars, and the distance and amount of extinction to the system. We constrain the temperatures such that $T \ge 3300$ K to be consistent with the combined spectral type and temperature ratio from JKTEBOP. The best fit parameters are: $T_{\text{pri}} = 3700$ K, $T_{\text{sec}} = 3600$ K, $R_{\text{pri}} = 1.42 R_{\odot}$ (+3$\sigma$), $R_{\text{sec}} = 1.15 R_{\odot}$ (+1$\sigma$), distance = 830 pc and $A_V = 0.1$. The results are shown in the left panel of Fig. 13, showing that there is clear excess emission in the mid-IR compared to any reasonable combination of naked stellar photospheres.

We therefore tested whether this excess could be due to extended dust emission in the vicinity of the two stars. To do this, we must first set out the basic geometry of the system. As the system is young, each star could be surrounded by a circumstellar disk, and both by a circumbinary disk. Tidal truncation of the circumstellar disks by the other star would limit their outer radii to about one third of the binary separation (Paczynski 1977; Papaloizou & Pringle 1977, i.e. ~0.017 AU or 3.6 $R_{\odot}$). The circumbinary disk would be centred on the centre of mass of the binary, and its inner radius is expected to be roughly twice the binary separation, i.e. ~0.1 AU or 22 $R_{\odot}$, because the tidal torque exerted by the binary prevents the disk from extending further in (Lin & Papaloizou 1979). However, material may stream through the inner cavity to be accreted onto the circumstellar disks (Artymowicz & Lubow 1996). A schematic representation of this geometry is shown in Fig. 14.

There is no evidence for excess emission in the near-IR, which is where we would expect any emission from material in the circumstellar disks to peak. If the primary were a single star, its sublimation radius, i.e. the distance beyond which amorphous dust grains typically survive ($T_{\text{dust grain}} \le 1500$ K), would be 3.9 $R_{\odot}$. This assumes that the temperature of a dust grain at a distance $r$ from a star of luminosity $L$ is given by $T_{\text{dust}} = [L/(16\pi\sigma r^2)]^{1/4}$, i.e. that the grains are spherical blackbodies, and are located in an optically thin environment. This close to the

primary star, the effect of the secondary's flux on the dust temperature is not significant. Similarly, if the secondary star were single, its sublimation radius would be $3.2 R_{\odot}$, at which distance the flux of the primary star has a negligible effect. Therefore, the circumprimary disk would not be expected to contain dust and the circumsecondary disk would only be expected to contain a small amount of very hot dust near the sublimation temperature. Added to this, the CCF peaks are not obviously rotationally broadened, which given the resolution of the FLAMES spectrograph, implies that the two stars have rotational velocities of order 17 km s$^{-1}$ or less. Taking this limiting value gives corotation radii (the radii in the circumstellar disks where the Keplerian period is equal to the stellar rotation period) of $\ge 9.1 R_{\odot}$ and $\ge 7.4 R_{\odot}$ for the primary and secondary stars, respectively. The strong magnetic fields of young, low-mass stars are usually assumed to lead to truncation of the inner disk at or near the corotation radius, which makes it less likely that circumstellar disks exist in this system. Note that the allowed rotational velocities of both stars are consistent with synchronisation; the corresponding rotation periods are $\ge 3.87$ and $\ge 3.31$ days for the primary and secondary respectively.

We therefore attempted to model the mid-IR excess (essentially the residuals of the photosphere-only model shown in the left panel of Fig. 13 but with the radii of the two stars constrained to within $1\sigma$) as emission from dust in a circumbinary disk. As the inner edge of the circumbinary disk is $> 20 R_{\odot}$ from the centre of mass of the system, we can approximate the incident flux by that of a single star with luminosity $L_{\star} = L_1 + L_2$, temperature $T_{\star}$ such that $\sigma T_{\star}^4 = L_1/(4\pi R_1^2) + L_2/(4\pi R_2^2)$ and radius $R_{\star}$ such that $\sigma T_{\star}^4 = L_{\star}/(4\pi R_{\star}^2)$.

We consider a razor thin disk that is aligned with the plane of the binary's orbit and which radiates away energy from both incident stellar flux and the gravitational potential energy released by the accretion of gas. We estimate an order of magnitude accretion rate of $\dot{M} = 10^{-11} M_{\odot} \text{yr}^{-1}$ from the multicomponent H$\alpha$ emission, using the relation described in Fang et al. (2009). This is also consistent with the UV excess upper limit derived from fitting the SED models of Robitaille et al. (2007). We compute the disk temperature as

$$T_{\text{disk}} = (T_{\text{irr}}^4 + T_{\text{acc}}^4)^{1/4} \quad (4)$$

where, following Armitage (2010), $T_{\text{irr}}$ and $T_{\text{acc}}$ are the temperatures arising from stellar irradiation and accretion respectively, and are given by

$$T_{\text{irr}}^4 = \frac{T_{\star}^4}{\pi} \left( \sin^{-1}\left(\frac{R_{\star}}{r}\right) - \frac{R_{\star}}{r} \sqrt{1 - \left(\frac{R_{\star}}{r}\right)^2} \right) \quad (5)$$

$$T_{\text{acc}}^4 = \frac{3GM_{\star}\dot{M}}{8\pi\sigma r^3} \left(1 - \sqrt{\frac{R_{\star}}{r}}\right). \quad (6)$$

The disk SED is shown by the brown dot-dashed line in the right panel of Fig. 13, which is two orders of magnitude smaller than that observed between 3.6 and 8 $\mu$m. We need an accretion rate of $\sim 10^{-7} M_{\odot} \text{yr}^{-1}$ to reproduce the observed fluxes. Such a high accretion rate would lead to much stronger H$\alpha$ emission, and a large UV excess, three orders of magnitude larger than observed (from the u-band magnitude of the system). We conclude that the mid-IR excess cannot be explained by emission from a circumbinary disk.
We now consider the possibility that the mid-IR excess is due to optically thin dust located in the inner cavity of the circumbinary disk, following Jensen & Mathieu (1997), who found evidence for this in the SEDs of a number of young spectroscopic binaries. This could arise from ongoing low-level accretion from the circumbinary disk onto the two stars, as predicted by recent numerical simulations (Shi et al. 2012). The geometry and kinematics of any dust located within the cavity would, of course, be complex. However, for the purpose of estimating the approximate temperature of the dust and the amount of it needed to reproduce the observed SED, a very simplified toy model is sufficient. We assume that the dust surface mass density varies as $r^{-1/2}$, and that the gas-to-dust ratio is 100. Again, we approximate the incident radiation field by using a single star with the parameters given above. We are now considering dust located much closer to the two stars, so the validity of this approximation is more doubtful, but it is acceptable for a simple order of magnitude calculation. We model the line-of-sight optical depth through the disk as $\tau_{\lambda}(r) = \kappa_{\lambda}\Sigma(r)/\cos i$, where $\Sigma$ is the gas surface density, $i$ is the angle between the normal to the disk plane and the line of sight, and $\kappa_{\lambda} = 0.1(\lambda/250\mu\text{m})^{-1}\text{cm}^2\text{g}^{-1}$ is the opacity (Jensen & Mathieu 1997). The emitted flux is then given by:

$$ \lambda F_{\lambda} = \frac{\cos i}{D^2} \int_{r_{\text{in}}}^{r_{\text{out}}} \lambda B_{\lambda}[T_{\text{dust}}(r)] (1 - e^{-\tau_{\lambda}(r)}) 2\pi r \, dr. \quad (7) $$

where the inner and outer radii of the cavity, $r_{\text{in}}$ and $r_{\text{out}}$, should be larger than the sublimation radius of the single star we consider ($5 R_{\odot}$) and similar to the circumbinary disk's inner radius, respectively. The flux received from the cavity increases with the mass of dust as long as $\tau_{\lambda} < 1$ and becomes insensitive to this mass when the optical depth reaches unity<fnref n="9" />. We find that we can explain the observed mid-IR excess with as little as $\sim 1 \times 10^{-13} M_{\odot}$, which gives $\tau_{\lambda} \sim 1$ at the outer edge of the cavity for $\lambda$ between 3.6 and $8 \mu\text{m}$. In this model, the cavity extends from $\sim 5$ to $32 R_{\odot}$ with corresponding dust temperatures of $\sim 1450$ to $600 \text{ K}$. The resulting fit to the SED is shown in the right panel of Fig. 13. Although the outer radius of the dust model is slightly larger than the expected inner edge of the circumbinary disk, they are of the same order of magnitude, which is sufficient given the very simplified nature of our model. Note that the mass of dust required and the extent of the cavity are insensitive to interstellar reddening (up to $A_{\text{V}} = 1.0$). As a consistency check, we compared the mass of dust we calculate to that required to reproduce the estimated mass accretion rate, and find them in agreement to within an order of magnitude.

Our model is very crude, but it does demonstrate that the observed SED can be reproduced by invoking a very small amount of dust within the inner cavity of a circumbinary disk, such as one could expect to find in accretion streams. Unfortunately, the flux measurements beyond $10 \mu\text{m}$ are only rather weak upper limits, because of emission by interstellar dust filaments (nebulosity) superimposed on the target, so we cannot place any constraints on the circumbinary disk itself. However, the presence of such a disk would certainly be required to replenish the dust in the cavity. This is because when dust grains enter the inner cavity, their motion is essentially ballistic; the accretion timescale is therefore close to the free fall time.

## 5. Conclusions and future work

We report a new double-lined, detached eclipsing binary, which is comprised of two pre-main sequence M-dwarfs, discovered by the CoRoT space telescope in the NGC 2264 star forming region. We have measured the fundamental parameters of both stars using the continuous 23.4 day light curve obtained by CoRoT in March 2008, as well as 22-epochs of radial velocity data obtained almost four years later with VLT/FLAMES and WHT/ISIS. The orbit is consistent with circular ($e = 0.0037^{+0.0036}_{-0.0025}$) and has a period of $3.8745745 \pm 0.0000014$ days and a separation of $10.921 \pm 0.056 R_{\odot}$. The primary and secondary stars have masses of $0.668^{+0.012}_{-0.011}$ and $0.4953^{+0.0073}_{-0.0072} M_{\odot}$, and radii of $1.295^{+0.040}_{-0.037}$ and $1.107^{+0.044}_{-0.050} R_{\odot}$, respectively. The systemic velocity is within $1\sigma$ of the cluster median which, along with the presence of lithium absorption, strongly indicates cluster membership.

This is only the ninth PMS EB system with component masses below $1.5 M_{\odot}$, and it lies in a region of the mass-radius plane where existing observational constraints are very scarce. Within the current uncertainties, the parameters of the two stars are essentially consistent with the predictions of PMS stellar evolution models for ages of $\sim 3.5-6$ Myr. Although we have broken the degeneracy between the radius and surface brightness ratios, which can be a severe limitation in grazing EB systems, the final uncertainties on the component masses and radii are still a few percent. As highlighted by Torres et al. (2010), sub-percent accuracies are needed to place truly significant constraints on evolutionary models. We hope to improve the constraints on this system in the future, by incorporating additional data and refining our treatment of the out-of-eclipse variability.

Modelling of the system's broadband optical and infrared SED reveals a mid-IR excess that cannot be explained by emission from two stellar photospheres alone; additional, cool material is required. Dynamical arguments and simple disk models indicate that neither circumstellar nor circumbinary disks can explain this excess. We find that the excess can be reproduced by dust emission from within the inner cavity of a circumbinary disk and show that only a very small dust mass of $\sim 1 \times 10^{-13} M_{\odot}$ is required. Such small amounts of dust could be found in accretion streams from a circumbinary disk. This opens up the possibility that some of the OOE variations may be due to obsolescence of the central stars by dust located at the inner edge, or in the central cavity, of the circumbinary disk.

NGC 2264 was re-observed by CoRoT in December 2011 and January 2012, as part of a co-ordinated multi-observatory campaign involving *Spitzer*, *Chandra*, MOST<fnref n="10" /> (Microvariability and Oscillations of STars) and a host of ground-based observatories, including CFHT *ugr* photometry and VLT/FLAMES spectroscopy<fnref n="11" />. Although the sampling of the observations by other telescopes is significantly sparser than that of the CoRoT data, the widely separated bandpasses of the CFHT, CoRoT and warm *Spitzer* data should help to disentangle the contributions of starspots and any occultations by dusty material in the out-of-eclipse variability. We will also model the time-dependent $H\alpha$ emission seen in the FLAMES spectra, which may shed light on any ongoing accretion in the system.

A number of other EBs were discovered in the field of NGC 2264 during the same CoRoT observation. Once we have

<footnote n="9">⁹ Even if the dust is optically thin when considering light from the star incident on a particle within the cavity, it could be optically thick when viewed (near) edge-on, i.e. from the point of view of an observer on Earth.</footnote>

<footnote n="10">¹⁰ Note that CoRoT 223992193 fell outside the *Chandra* and MOST fields of view.</footnote>

<footnote n="11">¹¹ The FLAMES spectra used to derive RVs in the present paper were obtained as part of this campaign.</footnote>
characterised those which belong to NGC 2264 in more detail, they will form a unique sample of near-coeval EB systems formed from the same molecular cloud, and may shed further light on the range of ages and accretion histories in this region.

**Acknowledgements.** We thank Isabelle Baraffe for extending her models of stellar evolution (with mixing length $\alpha = 1.9$) down to lower masses allowing us to compare the source's masses and radii to these model predictions. We thank D. Barrado, the Director of Calar Alto Observatory, for granting Discretionary Time to this project. We also thank Catarina Alves de Oliveira for providing her automatic fitting method to derive the spectral type of the source from the Calar Alto low resolution spectrum. This research is based on data collected by the CoRoT satellite, which is publicly available via the IAS Data Center at http://idoc-corot.ias.u-psud.fr/. Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 088.C-0239(A). The WHT and INT are operated on the island of La Palma by the Isaac Newton Group in the Spanish Observatory del Roque de los Muchachos of the Instituto de Astrofísica de Canarias. Based on observations collected at the Centro Astronómico Hispano Alemán (CAHA) at Calar Alto, operated jointly by the Max-Planck Institut für Astronomie and the Instituto de Astrofísica de Andalucía (CSIC). This work is based in part on observations made with the *Spitzer* Space Telescope, obtained from the NASA/IPAC Infrared Science Archive, both of which are operated by the Jet Propulsion Laboratory, California Institute of Technology under a contract with the National Aeronautics and Space Administration. This work is based in part on observations made with Herschel, which is an ESA space observatory with science instruments provided by European-led Principal Investigator consortia and with important participation from NASA. This publication makes use of data products from the 2MASS, which is a joint project of the University of Massachusetts and the Infrared Processing and Analysis Center/California Institute of Technology, funded by the National Aeronautics and Space Administration and the NSF. This research has made use of the SIMBAD data base, operated at DCS, Strasbourg, France. Funding for the SDSS and SDSS-II has been provided by the Alfred P. Sloan Foundation, the Participating Institutions, the National Science Foundation, the US Department of Energy, the National Aeronautics and Space Administration, the Japanese Monbukagakusho, the Max Planck Society, and the Higher Education Funding Council for England. The SDSS Web Site is http://www.sdss.org/, where the list of the funding organisations and collaborating institutions can be found. This work was supported in part by the UK Science and Technology Facilities Council through grant ST/G002266 (SA) and studentship ST/J500641/1 (EG). A.M. acknowledges support from a Science and Technology Facilities Council studentship. J.B. thanks CNES for the partial funding of this project. S.H.P.A. acknowledges financial support from CNPq, CAPES and Fapemig. J.S. acknowledges funding from STFC in the form of an Advanced Fellowship. Finally, we would like to thank the anonymous reference for his/her careful reading of the manuscript and helpful suggestions for improvement.

## References

Abazajian, K. N., Adelman-McCarthy, J. K., Agüeros, M. A., et al. 2009, ApJS, 182, 543

Adams, M. T., Strom, K. M., & Strom, S. E. 1983, ApJS, 53, 893

Adelman-McCarthy, J. K., et al. 2009, VizieR Online Data Catalog: II/294

Affer, L., Micela, G., Favata, F., Flaccomio, E., & Bouvier, J. 2013, MNRAS, 430, 1433

Aigrain, S., & Irwin, M. 2004, MNRAS, 350, 331

Aigrain, S., Hodgkin, S., Irwin, J., et al. 2007, MNRAS, 375, 29

Aigrain, S., Pont, F., & Zucker, S. 2012, MNRAS, 419, 3147

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

Alexander, R. 2012, ApJ, 757, L29

Alves de Oliveira, C., Moraux, E., Bouvier, J., & Bouy, H. 2012, A&A, 539, A151

Andersen, J. 1991, A&ARv, 3, 91

Armitage, P. J. 2010, Astrophysics of Planet Formation (Cambridge: Cambridge University Press)

Artymowicz, P., & Lubow, S. H. 1996, ApJ, 467, L77

Auvergne, M., Bodin, P., Boisnard, L., et al. 2009, A&A, 506, 411

Baglin, A. 2003, Adv. Space Res., 31, 345

Baraffe, I., & Chabrier, G. 2010, A&A, 521, A44

Baraffe, I., Chabrier, G., Allard, F., & Hauschildt, P. H. 1998, A&A, 337, 403

Baraffe, I., Chabrier, G., & Gallardo, J. 2009, ApJ, 702, L27

Baraffe, I., Vorobyov, E., & Chabrier, G. 2012, ApJ, 756, 118

Bass, G., Orosz, J. A., Welsh, W. F., et al. 2012, ApJ, 761, 157

Baxter, E. J., Covey, K. R., Muench, A. A., et al. 2009, AJ, 138, 963

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

Borucki, W. J., Koch, D., Basri, G., et al. 2010, Science, 327, 977

Bouchy, F., Pont, F., Melo, C., et al. 2005, A&A, 431, 1105

Bouvier, J., Chelli, A., Allain, S., et al. 1999, A&A, 349, 619

Bouvier, J., Grankin, K. N., Alencar, S. H. P., et al. 2003, A&A, 409, 169

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

Cargile, P. A., Stassun, K. G., & Mathieu, R. D. 2008, ApJ, 674, 329

Chabrier, G., Gallardo, J., & Baraffe, I. 2007, A&A, 472, L17

Coughlin, J. L., López-Morales, M., Harrison, T. E., Ule, N., & Hoffman, D. I. 2011, AJ, 141, 78

Covino, E., Catalano, S., Frasca, A., et al. 2000, A&A, 361, L49

Covino, E., Melo, C., Alcalá, J. M., et al. 2001, A&A, 375, 130

Covino, E., Frasca, A., Alcalá, J. M., Paladino, R., & Sterzik, M. F. 2004, A&A, 427, 637

Cutri, R. M., Skrutskie, M. F., van Dyk, S., et al. 2003, 2MASS All Sky Catalog of point sources

Dahm, S. E. 2008, in The Young Cluster and Star Forming Region NGC 2264, ed. B. Reipurth, 966

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

Dahm, S. E., Simon, T., Proszkow, E. M., & Patten, B. M. 2007, AJ, 134, 999

D'Antona, F., & Mazzitelli, I. 1998, in Brown Dwarfs and Extrasolar Planets, eds. R. Rebolo, E. L. Martin, & M. R. Zapatero Osorio, ASP Conf. Ser., 134, 442

Dotter, A., Chaboyer, B., Jevremović, D., et al. 2008, ApJS, 178, 89

Doyle, L. R., Carter, J. A., Fabrycky, D. C., et al. 2011, Science, 333, 1602

Etzel, P. B. 1981, in Photometric and Spectroscopic Binary Systems, eds. E. B. Carling, & Z. Kopal, 111

Fang, M., van Boekel, R., Wang, W., et al. 2009, A&A, 504, 461

Favata, F., Micela, G., Alencar, S., Aigrain, S., & Zwintz, K. 2010, Highlights of Astronomy, 15, 752

Fúrez, G., Hartmann, L. W., Szentgyorgyi, A. H., et al. 2006, ApJ, 648, 1090

Flaccomio, E., Sciortino, S., Micela, G., et al. 1997, Mem. Soc. Astron. It., 68, 1073

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

Gibson, N. P., Aigrain, S., Roberts, S., et al. 2012, MNRAS, 419, 2683

Grankin, K. N., Bouvier, J., Herbst, W., & Melnikov, S. Y. 2008, A&A, 479, 827

Gustafsson, B., Edvardsson, B., Eriksson, K., et al. 2008, A&A, 486, 951

Hebb, L., Stempels, H. C., Aigrain, S., et al. 2010, A&A, 522, A37

Hebb, L., Cegla, H. M., Stassun, K. G., et al. 2011, A&A, 531, A61

Hillenbrand, L. A., & White, R. J. 2004, ApJ, 604, 741

Horne, J. H., & Baliunas, S. L. 1986, ApJ, 302, 757

Irwin, J., Aigrain, S., Hodgkin, S., et al. 2007, MNRAS, 380, 541

Irwin, J. M., Quinn, S. N., Berta, Z. K., et al. 2011, ApJ, 742, 123

Jensen, E. L. N., & Mathieu, R. D. 1997, AJ, 114, 301

Jester, S., Schneider, D. P., Richards, G. T., et al. 2005, AJ, 130, 873

King, J. R., Soderblom, D. R., Fischer, D., & Jones, B. F. 2000, ApJ, 533, 944

Lacy, C. H. S., Torres, G., Claret, A., & Vaz, L. P. R. 2005, AJ, 130, 2838

Lamm, M. H., Bailer-Jones, C. A. L., Mundt, R., Herbst, W., & Scholz, A. 2004, A&A, 417, 557

Lin, D. N. C., & Papaloizou, J. 1979, MNRAS, 188, 191

Makidon, R. B., Rebull, L. M., Strom, S. E., Adams, M. T., & Patten, B. M. 2004, AJ, 127, 2228

Mathieu, R. D., Baraffe, I., Simon, M., Stassun, K. G., & White, R. 2007, in Protostars and Planets V (Tuscon: University of Arizona Press), 411

Mayne, N. J., & Naylor, T. 2008, MNRAS, 386, 261

Mazeh, T. 2008, in EAS Pub. Ser. 29, eds. M.-J. Goupil, & J.-P. Zahn, 1

Meisenheimer, K. 1994, Sterne und Weltraum, 33, 516

Morales-Calderón, M., Stauffer, J. R., Hillenbrand, L. A., et al. 2011, ApJ, 733, 50

Morales-Calderón, M., Stauffer, J. R., Stassun, K. G., et al. 2012, ApJ, 753, 149

Naylor, T. 2009, MNRAS, 399, 432

Nelson, B., & Davis, W. D. 1972, ApJ, 174, 617

Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012a, ApJ, 758, 87

Orosz, J. A., Welsh, W. F., Carter, J. A., et al. 2012b, Science, 337, 1511

Paczynski, B. 1977, ApJ, 216, 822

Papaloizou, J., & Pringle, J. E. 1977, MNRAS, 181, 441

Park, B.-G., Sung, H., Bessell, M. S., & Kang, Y. H. 2000, AJ, 120, 894

Popper, D. M., & Etzel, P. B. 1981, AJ, 86, 102

Prša, A., & Zwitter, T. 2005, ApJ, 628, 426

Prša, A., Batalha, N., Slawson, R. W., et al. 2011, AJ, 141, 83

Ramírez, S. V., Rebull, L., Stauffer, J., et al. 2004, AJ, 127, 2659

Rasmussen, C. E., & Williams, C. K. I. 2006, Gaussian processes for machine learning (Cambridge: MIT Press)

Rebull, L. M., Makidon, R. B., Strom, S. E., et al. 2002, AJ, 123, 1528

Reipurth, B., Yu, K. C., Moriarty-Schieven, G., et al. 2004, AJ, 127, 1069

Robitaille, T. P., Whitney, B. A., Indebetouw, R., & Wood, K. 2007, ApJ, 169, 328

Schwamb, M. E., Orosz, J. A., Carter, J. A., et al. 2013, ApJ, 768, 127

Sergison, D. J., Mayne, N. J., Naylor, T., Jeffries, R. D., & Bell, C. P. M. 2013, MNRAS, 434, 966
E. Gillen et al.: A low-mass, pre-main sequence eclipsing binary with evidence of a circumbinary disk

Shi, J.-M., Krolik, J. H., Lubow, S. H., & Hawley, J. F. 2012, ApJ, 749, 118

Siess, L., Dufour, E., & Forestini, M. 2000, A&A, 358, 593

Sing, D. K. 2010, A&A, 510, A21

Southworth, J. 2013, A&A, 557, A119

Southworth, J., Maxted, P. F. L., & Smalley, B. 2004, MNRAS, 351, 1277

Southworth, J., Maxted, P. F. L., & Smalley, B. 2005, A&A, 429, 645

Southworth, J., Bruntt, H., & Buzasi, D. L. 2007, A&A, 467, 1215

Stassun, K. G., Mathieu, R. D., Vaz, L. P. R., Stroud, N., & Vrba, F. J. 2004, ApJS, 151, 357

Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2006, Nature, 440, 311

Stassun, K. G., Mathieu, R. D., & Valenti, J. A. 2007, ApJ, 664, 1154

Stassun, K. G., Mathieu, R. D., Cargile, P. A., et al. 2008, Nature, 453, 1079

Stempels, H. C., Hebb, L., Stassun, K. G., et al. 2008, A&A, 481, 747

Sung, H., & Bessell, M. S. 2010, AJ, 140, 2070

Sung, H., Bessell, M. S., & Lee, S.-W. 1997, AJ, 114, 2644

Sung, H., Bessell, M. S., & Chun, M.-Y. 2004, AJ, 128, 1684

Sung, H., Bessell, M. S., Chun, M.-Y., Karimov, R., & Ibrahimov, M. 2008, AJ, 135, 441

Sung, H., Stauffer, J. R., & Bessell, M. S. 2009, AJ, 138, 1116

Tody, D. 1993, in Astronomical Data Analysis Software and Systems II, eds. R. J. Hanisch, R. J. V. Brissenden, & J. Barnes, ASP Conf. Ser., 52, 173

Tognelli, E., Prada Moroni, P. G., & Degl'Innocenti, S. 2011, A&A, 533, A109

Tonry, J., & Davis, M. 1979, AJ, 84, 1511

Torres, G., Andersen, J., & Giménez, A. 2010, A&ARv, 18, 67

Walker, M. F. 1956, ApJS, 2, 365

Welsh, W. F., Orosz, J. A., Carter, J. A., et al. 2012, Nature, 481, 475

Wilson, R. E., & Devinney, E. J. 1971, ApJ, 166, 605

Young, E. T., Teixeira, P. S., Lada, C. J., et al. 2006, ApJ, 642, 972

Zwintz, K., Kallinger, T., Guenther, D. B., et al. 2011, ApJ, 729, 20