the question of the origin of stellar clusters is an old one .
as early as 1785 , william herschel first considered this problem in the pages of these _ transactions _ when he speculated on the origin of the clusters messier 80 and messier 4 in ophiuchus .
more than two centuries later , despite profound advances in astronomical science , we find that the question of the physical origin of stellar clusters remains largely a mystery .
this is , in part , due to the fact that cluster formation is a complex physical process that is intimately linked to the process of star formation , for which there is as yet no complete theory .
the development of a theoretical understanding of star and cluster formation therefore depends on first acquiring detailed empirical knowledge of these phenomena .
this can be a very difficult task .
consider , for example , the globular clusters , the most massive stellar clusters in our own galaxy , the milky way .
these stellar systems are more than 12 billion years old and are no longer formed in the milky way .
consequently , direct empirical study of their formation process is not possible .
the situation is considerably better for galactic open clusters .
these typically span ages from 1 myr to 1 gyr and thus must be continually forming in the milky way , making direct observational study of their formation process possible , at least in principle .
however , such studies have been seriously hindered by the fact that open clusters are born in molecular clouds and , during their formation and early evolution , are completely embedded in molecular gas and dust .
they are thus obscured from view at optical wavelengths , where the traditional astronomical techniques are most effective .
fortunately , molecular clouds are considerably less opaque at infrared wavelengths and the development and deployment of infrared imaging cameras on optical- and infrared - optimized telescopes during the past two decades has provided astronomers with the ability to detect , survey and systematically study the extremely young embedded stellar clusters within nearby molecular clouds .
such studies indicated that embedded clusters are quite numerous and that they account for a significant fraction , if not the majority , of all star formation presently taking place in the galaxy . similarly , as discussed in accompanying articles by de grijs ( 2010 ) and larsen ( 2010 ) , young clusters also appear to account for a significant fraction of star formation in other galaxies , such as in the very active and luminous starburst galaxies and in the vigorous star - formation episodes which accompany galaxy mergers and close interactions . in our galaxy ,
the embedded - cluster phase lasts only 24 myr and the vast majority of embedded clusters which form in molecular clouds dissolve within 10 myr or less of their birth .
high cluster infant mortality has also been inferred for clusters in other galaxies ( see , e.g. , de grijs 2010 ; larsen 2010 ) .
this early mortality of embedded clusters is likely a result of the low star - formation efficiency that characterizes the massive molecular - cloud cores within which the clusters form .
the physical reason for the low formation efficiencies is not well understood and the origin of the massive cores themselves is one of the many mysteries of modern astrophysical research .
the embedded clusters are the primary laboratory for research into the question of the physical origin of stellar clusters . at the present time , the answer to this question is shrouded by the dusty veils of the dark molecular clouds in which stars and clusters are born . the ultimate goal of a theoretical understanding of the physical process of cluster formation is far from being realized .
the first step towards achieving this goal is to construct a solid empirical foundation upon which a physical theory can eventually be built .
this then requires a detailed understanding of the physical nature of embedded clusters and the environments in which they form . in this article
, i will review the empirical knowledge that is setting the stage for the eventual development of a theory of cluster formation .
i will describe existing knowledge concerning many of the basic physical properties of the embedded - cluster population , including the cluster mass function , cluster birthrate , structural properties , etc .
, as well as the nature of the gaseous cloud cores in which they form .
i will ( at the end ) briefly speculate on the implications of these results for understanding the origin of stellar clusters . a more theoretical discussion of possible physical mechanisms for forming such clusters is contained in the accompanying article by clarke ( 2010 ) .
it is useful to begin any discussion of embedded clusters with the definition of what constitutes such objects . following lada & lada ( 2003 ) ,
we first define a cluster to be a physically related group of stars that satisfy the following two requirements .
first , the group must have a stellar - mass volume density sufficiently large to render it stable against tidal disruption by the galaxy ( i.e. , @xmath1 pc@xmath2 ; bok 1934 ) and by passing interstellar clouds ( i.e. , @xmath3 pc@xmath2 ; spitzer 1958 ) .
a group of only 810 , 0.5 stars within a radius of 1 pc would satisfy this first requirement .
second , the group should contain a sufficient number of members to ensure that , if in a state of virial equilibrium , its evaporation time ( i.e. , the time it takes for internal , gravitational stellar encounters to eject all of its members ) is greater than the typical lifetime of galactic open clusters of @xmath4 yr . for a cluster in virial equilibrium ,
this number is about 35 stars .
an embedded cluster is then defined as a cluster that is fully or partially embedded in interstellar gas and dust .
the above definition of a cluster includes both gravitationally bound and unbound stellar systems .
bound stellar clusters are those whose total energy ( kinetic and potential ) is negative . when determining the total energy of a cluster , we take into account contributions from all forms of mass within the boundary of the cluster , including interstellar material .
this is crucial in evaluating the dynamical states of embedded clusters , since generally the mass of interstellar material exceeds the stellar mass within the cluster boundaries .      because most of their members are often optically obscured by interstellar dust , infrared imaging is generally required to identify and investigate embedded clusters .
the existence of a cluster is established observationally by an excess density of stars over the background .
consequently , the detectability of a cluster depends on its richness and compactness , the brightness of the cluster members and the density of background objects .
figure 1 shows optical and infrared images of the famous orion nebula and its embedded cluster .
it clearly illustrates the advantage of infrared imaging for revealing the stellar content of embedded clusters .
a number of authors have compiled catalogues of embedded clusters based on infrared observations .
( 2003 ) compiled a list of 34 clusters , which is likely nearly complete for clusters within 1 kpc of the sun .
lada & lada ( 2003 ) list 76 embedded clusters with well - determined properties within 2.5 kpc of the sun .
they estimate that their list is representative but incomplete and that there are likely a total of about 200 such clusters within this volume of space , most of which are unobserved to date .
( 2003 ) compiled a highly incomplete list containing more than 300 infrared clusters , stellar groups and candidate clusters , most within about 5 kpc of the sun , some fraction of which are likely embedded .
even cursory consideration of these studies suggests that embedded clusters are quite numerous in the galactic disc .
although the existence of a cluster can be established by the increase of source density over the background , identification of the individual members is considerably more difficult , especially for intrinsically faint members whose numbers are only comparable or significantly lower than those of the background / foreground stellar field population . the size of the cluster s membership can be determined statistically through comparison with the background / foreground field - star surface density over the cluster area . however , determining whether or not a specific star is a cluster member requires independent information .
since embedded clusters are very young stellar systems , independent indicators of stellar youth , such as the presence of circumstellar discs , variable emission lines , x - ray emission , etc .
, can be employed to ascertain membership of individual stars .    for typical embedded - cluster ages of 13 myr , at least half of the members will be surrounded by circumstellar discs , remnants of the star - formation process ( e.g. , haisch _ et al .
_ 2001 ; hernandez _ et al .
many of these discs will be accretion discs and exhibit strong , variable hydrogen recombination lines at optical and infrared wavelengths ( e.g. , green & lada 1996 ; white _ et al .
infrared spectroscopic surveys , especially those done with multi - object spectrographs , can be used to identify accreting objects in a cluster field . however , such observations can be time consuming and even prohibitive for faint members , which can account for a significant portion of the cluster population .
circumstellar discs can be detected directly with infrared broad - band photometry as a readily measurable excess above the expected emission from a normal stellar atmosphere .
excesses due to circumstellar discs and envelopes are most robustly measured at longer infrared wavelengths ( i.e. , @xmath5 m ) . in recent years ,
the _ spitzer space telescope _
, operating at wavelengths between 3 and 70 @xmath6 m , has provided the capability to survey and measure infrared excesses and circumstellar - disc emission with unprecedented sensitivity across the complete mass range of stars in embedded - cluster regions ( e.g. , allen _ et al .
these observations are providing the first statistically complete accountings of disc populations in individual young clusters ( e.g. , lada _ et al .
_ 2006 ; sicilia  aguilar _ et al .
_ 2006 ) and important new insights into their structural properties ( gutermuth _ et al .
_ 2009 ) .
although infrared observations are sufficiently sensitive to detect the presence of infrared excess from stars across the entire spectrum of stellar mass , there are typically many cluster members which have no discs or detectable infrared excess , presumably because these stars have passed quickly through the disc phase of early stellar evolution .
as many as 50% of all members of a young cluster may have lost their discs and any other accretion signatures at an age of 23 myr .
observations made from space with x - ray observatories such as the _ chandra x - ray observatory _ have shown that young stars emit x - rays at levels @xmath7 times that of normal stars , particularly during the first 10 myr of their lives ( feigelson & prebisch 2005 ) .
such emission appears to be detectable across almost the entire stellar - mass spectrum , whether or not the stars have discs ( e.g. , telleschi _ et al .
_ 2007 ) . in principle
, x - ray surveys can be used to identify the discless membership of a young stellar population and thus , in combination with infrared surveys , lead to identification of all members of an embedded cluster .
for example , a deep _ chandra observatory _ observation of the trapezium cluster in the orion nebula resulted in the x - ray image shown in figure 2 , in which all cluster members ( with and without discs ) are detected as x - ray sources ( getman _ et al .
_ 2005 ) . however , these _ chandra _ observations were very deep , and in general the degree to which x - ray surveys can sample the membership of a cluster depends on the sensitivity and depth of the observations . in typical situations , infrared observations are more sensitive to the faintest members of a cluster population ( e.g. , winston _ et al . _ 2007 ) . nonetheless , to accurately obtain a census of the membership of an embedded cluster
requires a varied set of observations and such information is , at the present time , only available for a limited set of clusters .
the frequency distribution of embedded - cluster masses is a fundamental property of these systems which any physical theory of cluster formation must explain .
lada & lada ( 2003 ) used the data from their cluster catalogue to uniformly derive masses for the known embedded clusters within 2.5 kpc .
they then constructed the embedded - cluster mass spectrum ( ecms ) and found it to be given by a power law of the form @xmath8 , where @xmath92.0 for masses in excess of about 50 . for a distribution with @xmath10 , the total mass of clusters in any logarithmic interval ( bin ) is constant .
thus , even though there are many more 50 than 1000 clusters , the very rare 1000 clusters contribute the same fraction of the total mass in clusters as the more numerous 50 clusters .
this , in turn , implies that if a star in the galaxy is born in a cluster , it is just as likely to form in a 1000 as in a 50 cluster .
the ecms very likely represents the primordial mass function of galactic open clusters . in this context , it is interesting to note that a recent attempt to determine the mass function of classical open clusters also produced a power - law relation with an index of @xmath11 ( piskunov _ et al .
_ 2008 ; see also discussion by larsen 2010 ) .
observations also suggest that the ecms peaks near 50 and then falls off towards lower masses .
indeed , in a detailed study of all known stellar groups and clusters within 1 kpc of the sun , porras _
( 2003 ) found that although small stellar ( @xmath12 30 ) groups greatly outnumber clusters , clusters with 100 or more members account for 80% of the total number of stars contained in both groups and clusters .
this is consistent with the notion that most stars that form in the galaxy form or have formed in embedded clusters with 100 or more members , and these are the clusters with masses of 50 or greater .
is shown as the dotted histogram .
massive embedded clusters from the list of ascenso ( 2008 ) are represented by the dashed trace .
( adapted from ascenso 2008 . ) ]    the most massive cluster in the lada & lada ( 2003 ) sample was the orion nebula cluster , weighing in at a little above @xmath13 .
more massive systems are known to exist , so it seems natural to assume that the lack of massive clusters in the lada & lada mass function is not caused by a truncation at high masses but is instead a result of the distance restriction of 2.5 kpc coupled with small - number statistics at higher cluster mass .
recently , ascenso ( 2009 ) compiled a list of known clusters independent of distance in the galaxy with masses estimated as in excess of 10@xmath14 .
the largest of these appears to be westerlund 1 , whose mass of about @xmath15 approaches that of globular clusters ( brandner _ et al .
_ 2008 ) . of the nine clusters in her list , six have ages of 3 myr or less and can be considered as part of the embedded , young cluster population . to estimate the numbers of massive , embedded clusters that may inhabit the milky way , ascenso ( 2008 ) scaled the 2.5 kpc ecms to correspond to that expected for the entire area of the galactic disc ( @xmath16 kpc ) .
figure 3 shows the scaled ecms , along with the six massive clusters from ascenso s list .
the predicted distribution of clusters is also plotted , calculated for a mass spectrum index @xmath17 and normalized to the scaled observed distribution at cluster masses @xmath18 ( ascenso 2008 ) . a significant number of young clusters are predicted for cluster masses in excess of @xmath19 , but only a few such clusters have been identified .    understanding the nature and origin of high - mass ( @xmath20 ) clusters is of considerable interest because observations indicate that such large clusters are a major component of the vigorous star - formation activity in merging and starburst galaxies ( e.g. , bastian _
_ 2005 ; see also accompanying reviews by de grijs 2010 ; larsen 2010 ) . are there examples of such clusters being produced in the milky way at the present epoch or is the ecms truncated at high mass ?
the question of whether there exists a high - mass cutoff or truncation to the ecms is difficult to answer .
there are no known embedded clusters with reliably determined masses in excess of @xmath21 .
one would think that such clusters would be easily identified if they existed .
however , the predicted numbers are low enough that small - number statistics could be responsible , in part , for the observed lack of such systems in the galaxy today . the present - day masses of classical open clusters offer additional constraints on the possibility of a truncation of the upper end of the primordial ecms .
for example , consider the most massive open clusters in our galaxy with ages @xmath22 myr .
those systems all have masses @xmath23 ( e.g. , bruch & sanders 1983 ; mermilliod 2000 ) .
recently , astronomers identified three red - supergiant - rich clusters with ages between 10 and 20 myr and masses between 2 and @xmath24 ( figer _ et al .
_ 2006 ; davies _ et al .
_ 2007 ; clark _ et al .
if the cluster - formation rate were constant over time , figure 4 suggests that we might expect at least 100 objects with masses @xmath25 and with ages of 300 myr or less .
so far , only three such objects are known .
therefore , either the ecms is truncated at masses @xmath26 or such clusters did not survive emergence from molecular clouds and dissolved very rapidly into the field .
the theory of stellar structure and evolution predicts fairly precisely the mass - dependent evolution of stellar luminosities and temperatures with age .
this enables the age of a stellar population to be determined by placing the stars on a hertzsprung
russell ( hr ) diagram , which plots the distribution of stellar luminosities versus temperatures .
the locations of stars are then compared to the predictions of stellar - evolution theory .
this has proved to be a powerful tool for age dating mature clusters whose more massive stars have finished their main hydrogen - burning stage and experience significant luminosity and temperature evolution .
however , age dating young , embedded clusters has proved to be considerably more difficult .
the reason for this is twofold .
first , the majority of stars in an embedded cluster are pre - main - sequence ( pms ) stars , that is , stars that have not yet evolved to the point where nuclear burning commences in their interiors . for these stars , the theoretical trajectories on the hr diagram can be highly uncertain , especially for stars with ages of 1 myr or less .
second , it is difficult to measure accurate luminosities and temperatures for pms stars since they are obscured by interstellar dust , exhibit infrared excess emission and are variable . nonetheless , a number of embedded clusters have been age dated in this fashion and typical median ages of stellar members in clusters range between 1 and 3 myr , with the age spread in a given cluster of comparable magnitude to its age ( lada & lada 2003 )
. these results are consistent with the early work of leisawitz _ et al . _
( 1989 ) , who found that clusters older than about 5 myr exhibit no evidence of interstellar molecular gas within or near their boundaries .
thus , the embedded - cluster phase of evolution lasts somewhere between 3 and 5 myr .
[ it is interesting to note here that recent analysis using only main - sequence cluster members on the hr diagram by naylor ( 2009 ) yields age estimates 1.5 to 2 times older than derived from the pms stars . ]        assuming typical ages of embedded clusters of between 1 and 2 myr , lada & lada ( 2003 ) estimated their birthrate within 2 kpc from the sun at @xmath27 to 4 myr@xmath28 kpc@xmath29 . upon emergence from their dusty wombs ,
embedded clusters become optically visible open clusters . yet
, the embedded - cluster birthrate is an order of magnitude higher than that ( 0.45 myr@xmath28 kpc@xmath29 ) , estimated for classical open clusters within 2 kpc of the sun ( battinelli & capuzzo
dolcetta 1991 ) .
the large difference in the estimated birthrates is key in understanding the formation and early evolution of galactic clusters .
the distribution of ages for a combined sample of embedded and open clusters within 2 kpc is shown in figure 4 .
embedded clusters fall into the first bin and classical open clusters have ages that span a range from roughly @xmath30 to @xmath31 yr .
the distribution is surprisingly flat .
also plotted is the predicted age distribution for a constant rate of cluster formation , assuming all clusters can live to ages of at least a few times @xmath31 yr .
there is a large and increasing discrepancy between the expected and observed numbers .
this figure suggests a high infant - mortality rate for embedded clusters .
most of these clusters do not survive emergence from their parent molecular cloud .
roughly 90% dissolve in less than 10 myr .
less than about 4% of the initial cluster population reach an age of 100 myr .
lada & lada ( 2003 ) suggest that only the most massive embedded clusters ( @xmath32 ) survive to ages of 100 myr or more .
it has been understood for some time that passing massive interstellar clouds could tidally disrupt open clusters with mass densities @xmath33 pc@xmath2 within 200 myr ( spitzer 1958 ; gieles _ et al . _ 2006 ) .
these tidal interactions with giant molecular clouds ( gmcs ) likely account for the flat age distribution of clusters with ages @xmath34 myr .
however , the rapid disintegration of most clusters in the first 10 myr of life requires another mechanism . to understand this process
, we need to first understand the conditions under which embedded clusters form .
the intimate physical association with interstellar gas and dust is the defining characteristic of embedded clusters .
such objects are immersed in dusty , gaseous material to varying degrees .
the most deeply embedded clusters are believed to be the youngest and least evolved objects of this class .
they are found within massive , dense structures or cores inside gmcs .
gmcs are , along with globular clusters , the most massive objects in the milky way .
they are complex structures , often defined by long filamentary features .
such clouds are primarily composed of molecular hydrogen gas and characterized by a wide range of gas densities ( @xmath35 cm@xmath2 ) . star formation has long been known to take place in the dense ( @xmath37 ) component of molecular clouds ( lada 1992 ) and this component only accounts for between 1 and 10% of the total mass of gmcs observed in the vicinity of the sun .
there are also indications that the star - formation rate in local clouds is closely related to their total content of dense gas ( lada _ et al .
the latter is found in the form of filaments and also in more compact cores .
individual stars can form along dense filaments , usually in compact , dense cores within the filamentary structure .
clusters of stars typically form in the most massive and dense cores within a gmc
. often , these massive cores are compact structures located at the ends or at an intersection of filaments , as illustrated in figure 5 ( myers 2009 ) .
the dense cores that form clusters range in mass from @xmath38 to 2000 and have sizes on the order of 12 pc in diameter ( e.g. , higuchi et al . 2009 ) .
the spatial extents of embedded clusters are comparable to the sizes of the dense cores within which they reside .
this suggests that the stars may be in approximate virial equilibrium with the gas in the core .
the cluster members may thus not be colder dynamically than the surrounding gas . on the other hand , there are some indications that the stars within clusters may have started out dynamically colder than the gas .
starless dense molecular cores , the precursors to individual new stars , are characterized by subvirial velocities in a number of clouds ( e.g. , walsh _ et al .
_ 2004 ; peretto _ et al .
_ 2006 ; lada _ et al .
_ 2008 ) and protostellar objects have also been found to be spatially more confined than more evolved young stellar objects ( e.g. , teixeira _ et al .
_ 2006 ; muench _ et al .
such observations suggest that , at birth , new stars in clusters are dynamically colder than the bulk of the gas surrounding them .
such stars are then expected to violently relax and collapse to smaller configurations . in doing so
, dynamical interactions can produce segregation of the most massive stars , similar to observations ( allison _ et al .
comparison of the stellar content of embedded clusters with the mass of their surrounding cores yields star - formation efficiencies of typically between 10 and 30% .
if the cluster stars were indeed initially dynamical colder than the bulk molecular gas , these measured ` final ' efficiencies will be higher than the initial birth efficiencies ( lada _ et al .
_ 1984 ; goodwin & bastian 2006 ) . the observed low to modest final star - formation efficiencies are key to understanding the early dynamical evolution and infant mortality of such objects .
m with the _
spitzer space telescope _ ( rathborne _
this cloud contains in excess of 1000 within a radius of a parsec and yet exhibits no indication of active star formation .
it is possibly a future site of cluster formation . ]
-1.0 in    the most promising sites to investigate the formation of clusters may be found in a class of objects known as infrared dark clouds ( irdcs ) ( e.g. , rathborne _ et al .
they are identified as high - extinction dark clouds when viewed against bright , diffuse mid - infrared ( 10@xmath39 m ) galactic background emission . to be regions of high extinction in
the mid - infrared requires them to have enormous column densities ( @xmath40 cm@xmath29 ) and correspondingly high volume densities ( @xmath42 cm@xmath2 ) , just the type of conditions that are ripe for star and cluster formation .
the more massive of these clouds tend to be concentrated in the so - called 5 kpc ring ( simon _ et al .
_ 2006 ) , a large - scale galactic structure rich in massive molecular clouds and active star formation .
millimetre - wave studies show many irdcs to contain cores of sufficient mass ( 1001000 ) to form modest to rich clusters . yet , many of these objects show no evidence of significant star or cluster formation .
figure 6 shows one of these irdcs imaged with the _ spitzer space telescope_. this particular cloud contains more than @xmath13 , but exhibits no evidence for star formation .
its compact size and large mass suggest that it is gravitationally bound and very likely a site of future cluster formation .
( 2009 ) estimate that a 300 cluster could ultimately emerge from this cloud .
other massive irdcs have exhibited evidence for the formation of massive stars within their boundaries . since
massive stars are rarely formed in isolation , these clouds may also be the sites of incipient cluster formation ( rathborne _ et al . _ 2006 ) .
because these clouds are identified from mid - infrared observations , which can only be obtained from space , their existence has only been known for about a decade .
since the massive irdcs tend to be relatively distant , they are challenging to observe and , as a consequence , little is known about their detailed physical properties .
unfortunately , too little is known about the detailed conditions in the clusterless irdcs to provide any strong constraints on the problem of cluster formation .
however , future observational capabilities provided by such facilities as the _ atacama large millimeter array _ and the _ james webb space telescope _ hold the promise for obtaining a more detailed description of the physical conditions in these objects and for deeper insights into the cluster - formation problem .
knowledge of the internal structure of an embedded cluster can provide additional clues relating to its origin and early dynamical history .
for the youngest clusters , the spatial distribution of member stars may be closely related to the primordial structure of the cold , dense molecular gas from which the clusters formed . in their review of embedded - cluster properties , lada & lada ( 2003 ) identified two structural types of clusters , those that exhibit extended , irregular surface - density distributions , often with multiple peaks , and those that appear to be relatively more compact and have well - defined centrally condensed surface - density distributions with one dominant peak .
more recently , ferreira & lada ( 2009 , personal communication ) used ground - based infrared imaging observations to systematically investigate the structure of approximately 40 previously catalogued embedded clusters .
they found that about 60% of these objects were of the centrally condensed variety .
the remainder exhibit relatively flat surface - density profiles with more irregular boundaries and more elongated shapes than the centrally condensed clusters .
observations with the _ spitzer space telescope _ at infrared wavelengths of 324 @xmath6 m provide an improved census of cluster membership and thus enable better measurements of cluster shapes and sizes .
these observations have shown that most embedded clusters are elongated , with a typical aspect ratio of about 2 , very similar to the aspect ratios of the filamentary molecular clouds in which the clusters are born ( allen _ et al .
_ 2007 ; gutermuth _ et al .
these observations imply that the structures of most embedded clusters probably reflect those of the original molecular clouds from which they formed .
they also indicate that understanding the origin of such clusters is intimately tied to understanding the origin and evolution of the dense molecular material that produces them .
a potentially significant clue concerning the nature of the physical process of cluster formation is found in analysis of the distributions of spatial separations of the youngest members of some embedded clusters .
figure 7 shows the distributions of the nearest - neighbour separations of protostars in an embedded cluster in the ngc 2264 star - forming complex ( teixeira _ et al .
this distribution is strongly peaked at an angular separation of 20 arcsec and clearly different from that expected for a random distribution of separations .
the characteristic separation corresponds to the jeans length for the density and temperature of the cloud core containing the cluster and suggests that jeans fragmentation may be the physical process that governs the formation of the cluster from the massive cloud core .
_ spitzer space telescope _ studies of other clusters have found such evidence for jeans fragmentation to be a more common property of embedded clusters than previously thought ( gutermuth _ et al .
because stars form in dense gas , it is the initial distribution of the dense material in a cloud that controls the degree to which a cloud forms an extended population of distributed stars and/or compact clusters .
that the mass spectrum of embedded clusters resembles the shape of the mass spectrum of massive cores in gmcs appears to support this surmise ( lada & lada 2003 ) .
once dense gas is present , star formation likely proceeds quickly .
the timescale is somewhat uncertain but the fact that numerous examples of massive , clusterless cores in irdcs exist probably indicates a time for the onset of cluster formation on the order of a million years . on the basis of the ages of stars in embedded clusters , the time needed to construct the full cluster is likely 24 myr from the onset of significant star formation .
the physical mechanism that causes a massive core to fragment and form a cluster of stars is not clear . that the jeans scale appears to describe
the separations of the young stellar objects in embedded clusters suggests that a jeans - like fragmentation process is responsible for the formation of a stellar cluster from a massive , dense core .
support for this idea also derives from measurements of the initial mass function ( imf ) of stars within individual clusters .
the stellar imf has a characteristic mass of about 0.3 , which appears invariant among clusters and even stars in the field .
this suggests a mass scale for star formation that is consistent with thermal jeans fragmentation ( larson 1992 ) .
the jeans instability describes the situation where gravity overcomes the internal thermal pressure of a gas cloud and leads to both its fragmentation and the collapse of its individual fragments .
however , it is not clear how this mechanism would operate in a very massive core where internal turbulence dominates thermal motions .
moreover , since the critical conditions of density and temperature which determine the value of the jeans length and mass likely vary among regions , one would expect the characteristic stellar mass of clusters to likewise vary , but it does not .
recent observations by the _
chandra _ and _ spitzer space telescopes _ have revived the idea that external triggering could represent another mode of cluster formation in the galaxy . in the cepheus b ( getman _
_ 2009 ) and w5 ( koenig _ et al .
_ 2009 ) star - forming complexes , spatiotemporal patterns of star formation have been found in the distribution of recently formed stars .
such patterns of sequential star formation have long been regarded as evidence for triggered star formation ( e.g. , elmegreen & lada 1977 ) .
triggering produces enhanced compression of existing molecular clouds , raising their mean densities and , in doing so , increasing the fraction of cloud mass above the critical densities necessary for star formation .
one way to compress a cloud is through increasing its external pressure .
it is possible that conditions in merging and interacting galaxies , and in the nuclear regions of starburst galaxies , produce globally increased external pressures that significantly increase the fraction of dense material in a cloud and lead to the formation of massive dense cores and clusters .
lack of such conditions ( i.e. , high external pressures ) may explain why very massive clusters are not presently forming in the milky way .
although it is beyond the scope of this particular review to discuss the theory of cluster formation , it is useful to briefly mention the results of recent calculations since some insight into the physical process of cluster formation can be provided by numerical simulations .
a more detailed discussion of the theory of cluster formation can be found in the accompanying review by clarke ( 2010 ) .
simulations of cluster formation follow the collapse and fragmentation of massive , initially turbulent clouds from some specified initial gaseous state down to the formation of stellar objects .
the results of these calculations yield impressive details concerning the clusters they form but are very sensitive to the input physics .
in particular , calculations starting with roughly similar initial conditions ( i.e. , @xmath43 turbulent clouds ) can have drastically different outcomes , depending on whether feedback from newly formed stars is included .
models without feedback ( e.g. , bate & bonnell 2005 ) produce highly filamentary cloud configurations reminiscent of observed clouds .
these models produce extensive ( efficient ) fragmentation and form rich clusters of primarily low - mass stars whose characteristic mass is directly related to the jeans mass and thus varies with initial conditions ( i.e. , density ) .
models that include feedback ( e.g. , krumholtz _ et al .
_ 2007 ; bate 2009 ) produce more massive cloud cores , smaller clusters , more massive stars and an invariant characteristic stellar mass , all in better agreement with observations of embedded clusters .
most importantly perhaps , the models which include feedback from the newly formed stars are also characterized by low star - formation efficiencies , again similar to observations .
once embedded clusters are formed , the process of their subsequent dynamical evolution and emergence from molecular clouds is relatively well understood .
analytic and numerical calculations have shown that the key factor in determining the evolution of an embedded cluster is the ( final ) star - formation efficiency characterizing its formation ( e.g. , lada _ et al . _ 1984 ; geyer & burkert 2001 ; kroupa & boily 2002 ) .
observations indicate that the star - formation efficiencies for embedded clusters range between 10 and 30% ( lada & lada 2003 ; higuchi _ et al .
_ 2009 ) with some indications that those with the highest efficiencies are the most evolved objects ( higuchi _ et al .
numerical simulations of cluster formation suggest that these low efficiencies are caused by the radiative feedback of the recently formed stars on their surrounding dense gas . in any case , the gravitational glue that binds the stars and gas together in an embedded cluster system is largely provided by the gas .
we expect that , as more and more stars form in the cluster , the feedback from the stars due to radiation and energetic jets and outflows , generated in the early stellar - accretion process , becomes destructive and begin to disperse the dust and gas in the core .
massive stars are particularly destructive and a single o star can completely disrupt an entire cloud in a very short period of time . because of the low efficiency of conversion of gaseous to stellar mass in these systems , they can easily be completely disrupted if the gas is removed on timescales comparable to the dynamical time of the system .
thus , low efficiency coupled with relatively rapid gas removal can explain why 90% of embedded clusters do not survive their emergence from molecular clouds as bound systems ( lada & lada 2003 ) .
paradoxically , the clusters that survive emergence from molecular clouds as bound systems are likely the most massive systems at the upper end of the embedded - cluster mass function .
this is despite the fact that these systems are the ones most likely to produce destructive o stars .
simulations indicate that , as clusters are disrupted by gas - removal processes , this disruption is rarely 100% complete and bound remnants containing a small fraction of the original stellar mass are often left behind ( lada _ et al . _ 1984 ; kroupa & boily 2002 ) .
the size of the bound remnant depends on the initial size of the cluster and on how rapidly gas removal occurs .
lada & lada ( 2003 ) suggest that the field population of open clusters in the disc of the milky way was likely provided by embedded clusters with masses in excess of 500 .
this may also explain the relative lack of truly massive ( i.e. , @xmath20 ) clusters in the milky way with ages in excess of 10 myr .
the few that are formed in the milky way do not emerge completely unscathed .
only smaller remnants of these giants survive as bound systems beyond 10 myr .
more than two centuries after herschel ( 1785 ) first considered the question , we now understand that clusters are formed in cold , dark molecular clouds .
stellar clusters begin their lives deeply buried in dense molecular gas and dust as embedded infrared protoclusters . over the past two decades ,
advances in detector and telescope technology have led to steady advances in our knowledge of embedded protoclusters .
these extremely young stellar systems appear to account for a significant fraction of all star formation presently taking place in the milky way and , thus , they are tracers of the current epoch of star formation in the galaxy . although considerable progress has been made towards understanding the basic physical properties of embedded clusters , a physical theory of cluster formation eludes us . in the milky way
, the primary mode of cluster formation at the present epoch appears to be one in which relatively compact , massive dense cores in gmcs undergo a process of jeans - like fragmentation that transforms cold , dense interstellar material into stellar form .
triggering of clouds by shocks or other mechanisms that lead to an increase in external pressure may represent an additional mode of galactic cluster formation .
although its relative contribution to present - day cluster formation in the milky way is unclear , this latter mechanism may be significant in interacting galaxies and nuclear starbursts , and perhaps even in the early history of the milky way itself .
it is now clear that to develop a predictive theory of cluster formation will require both a better understanding of the process of star formation and a comprehensive understanding of the physical mechanisms that organize the dense material of a gmc into massive , compact cores which , for a typical gmc in the milky way , occupy only a small fraction ( @xmath44% ) of its volume and account for only a small fraction ( 110% ) of its total mass .
i thank joana ascenso , ken janes , lisa townsley , eric feigelson , elizabeth lada and bruno ferreira for informative discussions , joana ascenso for preparing figure 3 and joo alves for providing figure 5 .
i am grateful to simon goodwin and richard de grijs for constructive comments that improved this review .
allen , l. , megeath , s. t. , gutermuth , r. , myers , p. c. , wolk , s. , adams , f. c. , muzerolle , j. , young , e. & pipher , j. l. 2007 the structure and evolution of young clusters , in _ protostars & planets v _ ( eds b. reipurth , d. jewitt & k. keil ) , pp .
361376 , tucson : university of arizona press .
bastian , n. , trancho , g. , konstantopoulos , i. s. & miller , b. 2009 _ gemini _ spectroscopic survey of young star clusters in merging / interacting galaxies : iii .
the antennae .
_ astrophys .
j. _ , * 701 * , 607619 .
clark , j. s. , negueruela , i. , davies , b. , larionov , v. m. , ritchie , b. w. , figer , d. f. , messineo , m. , crowther , p. a. & arkharov , a. a. 2009 a third red supergiant rich cluster in the scutum  crux arm . _
_ , * 498 * , 109114 .
feigelson , e. , townsley , l. , gdel , m. & stassun , k. 2007 x - ray properties of young stars and stellar clusters . in _ protostars & planets v _
( eds . b. reipurth , d. jewitt & k. keil ) , pp . 313328 , tucson : university of arizona press .
figer , d. f. , mackenty , j. w. , massimo , r. , smith , k. , najarro , f. , kudritzki , r. p. & herrero , a. 2006 discovery of an extraordinarily massive cluster of supergiants .
_ astrophys .
j. _ , * 643 * , 11661179 .
getman , k. v. , feigelson , e. d. , luhman , k. l. , sicilia
aguilar , a. , wang , j. & garmire , g. p. 2009
protoplanetary disk evolution around the triggered star - forming region cepheus b. _ astrophys .
j. _ , * 699 * , 14541472 .
gieles , m. , portegies zwart , s. f. , baumgardt , h. , athanassoula , e. , lamers , h. j. g. l. m. , sipior , m. & leenaarts , j. 2006 star cluster disruption by giant molecular clouds .
_ mon . not .
r. astron .
_ , * 371 * , 793804 .
gdel , m. , padgett , d. l. & dougados , c. 2007 the taurus cloud : multiwavelength surveys with _
xmm  newton _ ,
spitzer space telescope _ and the _ canada  france
hawaii telescope_. in _ protostars & planets v _
b. reipurth , d. jewitt & k. keil ) , pp . 329344 , tucson : university of arizona press .
gutermuth , r. a. , megeath , s. t. , myers , p. c. , allen , l. e. , pipher , j. l. & fazio , g. g. 2009 a _ spitzer _ survey of young stellar clusters within one kiloparsec of the sun : cluster core extraction and basic structural analysis . _ astrophys . j. suppl
_ , * 184 * , 1883 .
higuchi , a. e. , kurono , y. , saito , m. & kawabe , r. 2009 a mapping survey of dense clumps associated with embedded clusters : evolutionary stages of cluster - forming clumps .
_ astrophys .
j. _ , in press ( arxiv:0909.1132 ) .
hernndez , j. , hartmann , l. , calvet , n. , jeffries , r. d. , gutermuth , r. , muzerolle , j. & stauffer , j. 2008 a _
view of protoplanetary disks in the @xmath45 velorum cluster .
_ astrophys .
, * 686 * , 11951208 .
koenig , x. p. , allen , l. e. , gutermuth , r. a. , hora , j. l. , brunt , c. m. & muzerolle , j. 2009 clustered and triggered star formation in w5 : observations with _
spitzer_. _ astrophys .
j. _ , * 688 * , 11421158 .
mermilliod , j .- c . 2000 massive clusters in the milky way and magellanic clouds . in _ massive stellar clusters _
a. lanon & c. boily ) , astron .
* 211 * , 4348 , san francisco : astron .
pac .
piskunov , a. e. , kharchenko , n. v. , schilbach , e. , roser , s. , sholtz , r. d. & zinnecker , h. 2008 the initial luminosity and mass functions of the galactic open clusters .
_ , * 487 * , 557566 .
peretto , n. , andr , p. & belloche , a. 2006 probing the formation of intermediate - to - high mass stars in protoclusters . a detailed millimeter study of the ngc 2264 clumps .
_ , * 445 * , 979998 .
porras , a. , christopher , m. , allen , l. , di francesco , j. , megeath , s. t. & myers , p. c. 2003 a catalog of young stellar groups and clusters within 1 kiloparsec of the sun .
j. _ , * 126 * , 19161924 .
sicilia  aguilar , a. , hartmann , l. , calvet , n. , megeath , s. t. , muzerolle , j. , allen , l. , dalessio , p. , mern , b. , stauffer , j. , young , e. & lada , c. 2006 disk evolution in cep ob2 : results from the _ spitzer space telescope_. _ astrophys . j. _ , * 638 * , 897919 .
teixeira , p. s. , lada , c. j. , young , e. t. , marengo , m. , muench , a. , muzerolle , j. , siegler , n. , rieke , g. , hartmann , l. , megeath , s. t. & fazio , g. 2006 identifying primordial substructure in ngc 2264 .
_ astrophys . j. _ , * 636 * , l45l48 .    telleschi , a. , gdel , m. , briggs , k.  r. , audard , m. & palla , f. 2007 x - ray emission from t tauri stars and the role of accretion : inferences from the _ xmm  newton _ extended survey of the taurus molecular cloud .
_ , * 468 * , 425442 .
winston , e. , megeath , s. t. , wolk , s. j. , muzerolle , j. , gutermuth , r. , hora , j. l. , allen , l. e. , spitzbart , b. , myers , p. & fazio , g. g. 2007 a combined _ spitzer _ and _ chandra _ survey of young stellar objects in the serpens cloud core . _ astrophys .
j. _ , * 669 * , 493518 .
white , r. j. , greene , t. p. , doppmann , g. , covey , k. r. & hillenbrand , l. a. 2007 stellar properties of embedded protostars , in _ protostars & planets v _
b. reipurth , d. jewitt & k. keil ) , pp . 117130 , tucson : university of arizona press .


eclipsing binaries have historically contributed a wealth to stellar astrophysics .
they have been used to determine distances , compute fundamental stellar parameters , and test stellar evolution models .
the _ kepler _ mission @xcite and its unprecedented precise photometry of @xmath0 stars , has allowed us to create a catalog of 2605 eclipsing binaries ( hereafter ebs ) in the _ kepler _ field @xcite .
this catalog is rich in interesting objects for individual study and also presents a large sample of ebs for statistical analysis . in studying this sample , we can attempt to determine the occurrence rate of ebs , circumbinary planets , and multiple star systems .    some theories for short - period binary star formation call for the presence of a third - body .
in these scenarios , the close binary was not created _ in situ _ , but rather at a larger separation as a part of a multiple star system @xcite .
tidal friction and kozai cycles between the inner - binary and a companion can cause the inner - orbit to shrink over time @xcite , and result in a hierarchical multiple system @xcite .
the spectroscopy and imaging studies by @xcite and @xcite have found 40% of binaries with periods less than 10 days , and 96% with periods less than 3 days , have a wide tertiary companion .
the general interpretation of these findings is that the tightest binaries likely became hardened over time through interactions with the tertiary companion , and the system evolves toward an increasingly hierarchical configuration .
indeed , the slowpokes study of ultra - wide binaries in the sloan digital sky survey @xcite found that the widest visual pairs with physical separations of 0.011 pc , in fact contain a tight binary @xmath180% of the time @xcite , again confirming the general picture that tight binaries are nearly always accompanied by wide tertiaries and that the tightest binaries are accompanied by the widest tertiaries .
discovery and study of these multiple systems gives new insight into the physics of ebs .
statistically , we can compare observed rates of multiple systems to theoretical models for short - period binary formation .
we can also model each system individually to study the disruptive dynamical effects seen in some cases .
@xcite determines ephemerides for the entire _ kepler _ eclipsing binary catalog .
if there are no external effects , a linear ephemeris will correctly predict all eclipse times of an eb . by measuring the exact time of each eclipse for a particular binary and comparing it to the calculated time from the linear ephemeris
, we can create an etv curve ( ` eclipse timing variations ' ; sometimes also referred to as an o - c diagram ) .
any trend in these timing residuals may be the result of one or more physical effects occurring in the system .    using transit timings and eclipse timings to find exoplanets is a well - known method @xcite .
@xcite , @xcite and @xcite used transit timings to detect and study multiple planetary systems , while _ kepler _ 16 @xcite , 34 , and 35 @xcite were validated , in part , through their eclipse timing variations .
the processes that can induce etvs , which are the focus of this paper , include the following :    * light time travel effect ( ltte ) : a third - body perturbing the center of mass of the binary system creates a light - time delay along the line of sight which can cause eclipses to appear earlier or later than expected . *
non - hierarchical third - body : the presence of a third - body actually changes the period of the binary over time . * mass transfer : mass transfer between the components in the binary changes the period . * gravitational quadrupole coupling ( applegate effect ) : spin - orbit transfer of angular momentum in a close binary due to one of the stars being active produces period changes up to @xmath2 times the binary period @xcite . *
apsidal motion : the rotation of the line of apsides causes a change in the time between primary and secondary eclipses even though the period remains unchanged ( requires an eccentric orbit ) .
* spurious signals : due to spots and other effects that distort the eb light curve .
@xcite previously published a list of 39 candidate third - body _ kepler _ systems using eclipse times and @xcite published a preliminary study on timing variations in 41 _ kepler _ eclipsing binaries .
@xcite will provide eclipse times for detached binaries , and this paper provides eclipse times for close binaries .
together , these two papers will comprehensively cover all 2605 binaries in the catalog .
s essentially uninterrupted observing over a long time baseline presents the opportunity to precisely time the eclipses and detect any underlying signals due to third bodies , apsidal motion , dynamical interaction , etc . due to the large number of ebs in the entire catalog , it is necessary to create an automated method for timing eclipses across the catalog .
short period and overcontact systems present a particular challenge due to spot activity and data convolution , due to a relatively long integration time .    in this paper
we discuss our method for automating eclipse timings for close _ kepler _ ebs in section 2 .
eclipse timings are reported for 1279 binaries in the catalog in section 3 . in section 4 ,
light time travel effect models for the 236 that are flagged as potential third - body candidates are also provided .
we discuss our findings in section 5 in the context of binary formation and evolution theory , and summarize our conclusions as well as information for accessing the products of our comprehensive eclipsing timing measurements in section 6 .
@xcite will update the _ kepler _ eclipsing binary catalog , raising the count of ebs from 2165 to 2605 .
the database is kept up - to - date with future data and revisions at ` http://keplerebs.villanova.edu ` . as changes and updates are made to the catalog , etvs
are being recomputed and updated automatically and made available in real - time through the online catalog .
@xcite will provide eclipse times for binaries with flat out - of - eclipse regions , covering most of the detached binaries with periods greater than 1 day .
there we locally detrend each eclipse and use a piecewise hermite spline template to determine the time of mid - eclipse .
this technique performs well on the set of detached systems but is not optimal for overcontact systems , systems with strong reflection effects or tidal distortion , or short - period binaries with only a few points in each eclipse due to _ kepler _
s 30 minute cadence .
for this reason , we divide the catalog based on the morphology parameter as described in @xcite .
this parameter is a value between 0 and 1 which describes the `` detachedness '' of an eclipsing binary , with 0 being completely detached and 1 being overcontact or ellipsoidal .
@xcite report timings for binaries with a morphology parameter less than 0.5 .
our method addresses and determines eclipse times for the remainder of the _ kepler _ eclipsing binary catalog .
the distribution of the catalog between these two methods is shown in fig .
[ etvcat ] , with 1279 binaries in the sample for this paper .
we detrend and phase `` sap '' ( simple aperture photometry ) _ kepler _ data through q16 as described in @xcite .
the upper envelope of the raw data is fit with a chain of legendre polynomials using a sigma - clipping technique and manually setting the breaks between sections and orders of the polynomials .
the data are then divided by this fit , resulting in a flat baseline .
these detrended data are then phased on the linear ephemeris as reported in @xcite , and used as input into the etv code , described below .
we fit a polynomial chain to the phased light curve data as described in @xcite .
this analytic function is a chain of four polynomials that is continuous , but not necessarily differentiable , at knots which were optimized to find the best overall solution .
this function does not represent a physical model , but rather analytically describes the mean phased shape of the binary light curve , an example of which can be seen in fig .
[ ecl_bounds ] .
we then take this analytical representation and , using a combination of heuristic and bisection approaches , determine the horizontal shift required to minimize the @xmath3 ( cost function ) for each individual eclipse as shown in fig .
[ bisection ] . in order to minimize the effect due to spots or imperfect detrending ,
a vertical shift is first determined using linear least squares for each eclipse and is applied before computing cost functions for horizontal shifts .
the cost function is initially sampled at 20 evenly - spaced phase shifts between -0.05 and 0.05 phase .
the minimum of this sampling is then used as the center of the bisection algorithm to quickly find the local minimum of the cost function .
the resulting @xmath4 values are unusually large because the errors on the _ kepler _ data are only formal and do not include any absolute calibration effect @xcite .
therefore , for each eclipse , we normalize the entire cost function such that the minimum cost is set to @xmath5 , where @xmath6 is the number of data points used for that eclipse and @xmath7 is the degrees of freedom , which we take to be 1 .
this reduced cost function is then used to compute 1-sigma errors on each timing to correspond to the @xmath8 contour .    for the shortest binaries in the catalog ,
however , the long - cadence data result in significant phase - smearing and limits our method to a very minimal number of points per cycle to determine a fit .
if there were to be a third - body , the signal would likely be buried in the noise induced by these factors .
for this reason , we include as many data points as possible in each eclipse timing . each data point is considered to belong to an eclipse if its phase as determined by the initial linear ephemeris is within bounds .
we initially set these bounds to be the mid - point between polyfit knots in the out - of - eclipse region as shown in fig .
[ ecl_bounds ] . to improve results for particular objects being studied individually , changing these bounds to use the knots ( instead of the mid - points ) can sometimes lower the systematics in the signal . for
any given eclipse , if the region between these bounds is not fully sampled or does not have at least 3 data points , then timings are not computed for that eclipse .
eclipse timings are then compared to the values expected from the linear ephemeris as reported by @xcite to compute the residuals and test for the presence of an etv signal .      due to a typically small number of points per eclipse ,
our timings are sensitive to various imperfections in the data processing , affecting the measured eclipse time and potentially introducing noise and/or fictitious signals in the etv signal . instrumental or astrophysical pulsations on top of the binary signal can change the shape of a single eclipse which can mimic a timing variation .
the detrending process attempts to remove these additional signals , but is not perfect , struggles at removing signals that happen during eclipse , and can also introduce spurious signals . also , all polyfits in the current version of the catalog use chains of four second - order polynomials , which does not always result in the ideal fit and can leave slight phase - dependent residuals . for the purpose of pipeline processing , we limit ourselves to second order polynomials , but note that , for special cases and in - depth studies , higher precision timings can be obtained by increasing the order of the fit .
until all polyfits are updated to a higher order in the future , we will use the second - order fits and manually run etv signals with a higher order for any individual etv signals that warrant further study . in the cases
when a binary has a period that is near - commensurate to _ kepler _ s 29.44 minute cadence , the period and cadence may beat , which results in a separate spurious signal .
any combination of these effects can cause issues in determining true and precise eclipse times when dealing with only a few data points .    fig .
[ anti - phase ] demonstrates how the cost function for the phase shift is affected by the vertical discrepancy in the out - of - eclipse region , creating a fictitious signal in which the etvs of the primary and secondary eclipse are in anti - phase .
the left of fig .
[ anti - phase ] plots four different eclipses , showing that over time the data in the out - of - eclipse region can be higher on either the right or the left .
when measuring timings for the primary and secondary separately , the cost function will artificially be minimized by `` pulling '' the analytic function towards the region with lower flux . since this will affect the primary and secondary in the opposite direction , we can mitigate for this effect by also running the fit over the entire phase .
this effectively averages out the anti - phase effect in the primary and secondary eclipses , projecting the real etv signal of the entire system .
[ anti - phase - beforeafter ] shows two cases where the anti - phase signal was removed , clearly showing whether there is a presence of any underlying etv signal .
these signals that show a `` random walk '' nature are discussed by @xcite .
unfortunately , since there is no rigorous way to discriminate between true and fictitious anti - phase signals , this process would also hide a physical etv signal such as apsidal motion .
since we are dealing with short - period binaries , most of these orbits will be quite circular so we do not expect to be able to detect any systems with apsidal motion anyway .
if a binary has short - cadence data available , they are usually limited to a short time baseline .
since we are generally looking for long period trends in the etv signal , we measure timings from the full long cadence dataset .
we also run timings on any available short cadence data in case some signal can be detected .    a few short period binaries in the catalog that appear to be overcontact or ellipsoidal variables
could actually be detached systems whose light curves are convolved by _ kepler _
s 30 minute long cadence exposure .
since these systems are less prone to timing noise due to spots or mass transfer than true overcontact systems , the phase smearing and limited data points per eclipse are the main issues preventing us from recovering any etv signal .
for this reason , short cadence data were requested and obtained via director s discretionary time for 31 short - period detached ebs in the catalog without previous short - cadence observations in the hope of detecting third - body etv signatures which were not visible in the long - cadence data .
unfortunately , none seem to exhibit any significant etvs .
eclipse times are computed for all available short - cadence data , but due to the longer baseline , long - cadence timings are used for detection and fitting of potential third - body orbits .
eclipse timing variations on the individual eclipses and the entire phase have been run for all objects with a morphology parameter greater ( less detached ) than 0.5 in the latest _ kepler _ eclipsing binary catalog .
our method requires at least 3 data points per timing , which allows us to get primary and secondary eclipse timings individually for long cadence data of binaries with periods as short as 3 hours and full phase timings for binaries as short as 1.5 hours .
we are able to determine timings for any binary with short cadence data , but since short cadence data availability is sparse and generally not for the whole length of the mission , short cadence etvs are determined separately .
plots and data for detrended light curves and eclipse times for the entire sample are available as a part of the _ kepler _ eclipsing binary online catalog at ` http://keplerebs.villanova.edu ` .
an excerpt of the eclipse times is shown in table [ tableetvs ] . with the third version of the catalog
being released shortly , the database will be updated as new data become available and ephemerides are further refined .
this etv code is incorporated into the pipeline : as objects are updated or added , their etvs are recomputed and updated in real - time .
as etvs are computed , the ephemerides in the catalog are refined by fitting a linear fit through the entire - phase timings and adjusting the values as necessary to get a `` flat '' trend . for any etv with a long - term sinusoidal trend ,
this could introduce systematics depending on the part of the sine curve observed and used to fit the linear trend .
in particular , for very long etv signals ( of the order of 1000 days and more ) , the measured orbital period of the binary will be anomalous because the variation can not be accounted for from available data .
all etv measurements were examined by eye for the presence of any interesting signal , discarding any that seem to be spurious based on their individual primary and secondary eclipse timings .
we do not expect to see evidence of apsidal motion in many of our targets due to their short periods and , consequently , circular orbits .
we also do not expect to be able to detect any signals due to gravitational quadrupole coupling .
this mechanism is able to create period changes with amplitudes on the order of @xmath2 times the period of the binary , meaning a maximum of 3.5 seconds for a binary with a period of 4 days , falling well within our noise limits .
etv signals that are sinusoidal in nature or show any sign of curvature are flagged and fit for both a third - body signal and a parabolic mass transfer model . for the cases where we only see a sign of curvature and not a full cycle
, we could either be seeing a section in a long period third body signal or mass transfer . to determine whether we consider the signal as a candidate third body or mass transfer
, we compare the two models using the bayesian information criterion @xcite : @xmath9 where @xmath10 are the data , @xmath11 the model , @xmath12 the number of data points , and @xmath13 is the number of parameters used in the fit . in the case of the eccentric ltte model @xmath14 , for the circular ltte model @xmath15 , and for the mass transfer model @xmath16 .
the better fit , as determined by the lower bic value , then determines whether we consider the signal as a candidate third body or mass transfer .
thirty - one etv signals were better fit by a parabola than an ltte orbit , and are possibly caused by mass transfer or the applegate effect instead of the presence of a third body @xcite . a selection of these signals are shown in fig .
[ etv_mt ] , all kics are listed in table [ tablemt ] , all of which are available on the online catalog .      236 binaries ( of the sample )
were flagged as candidate third bodies .
the results of the model fits are reported in table [ tablerapp ] with a selection plotted in fig .
[ etv_ltte ] . based on the fitted period
, we then divided these third body candidates into three sections .
the first group contains third body signals with periods less than 700 days , such that there are at least two full cycles of the signal present in data through q16 .
these systems have the highest confidence and are most likely due to the presence of a tertiary component .
the second group contains signals with periods between 700 and 1400 days , such that there is at least one full cycle present .
the last group contains signals with periods longer than 1400 days .
often these detections merely show some sign of curvature in the etv signal and so a full sinusoidal signal can not yet be confirmed .
for this reason the fits generally have large errors and many of these may not even be true triple systems , particularly the signals on the closest binaries which are more likely to be due to mass transfer .
@xcite presented an initial study of eclipse timings in 41 _ kepler _ binaries .
of their entire sample of 41 binaries , 40 are still in the _ kepler _ eb catalog ( kic 4678873 has since been removed from the catalog as a false positive ) , 32 fall under the scope of this paper ( have a morphology less detached than 0.5 ) , and 9 appear in our list of third - body signals .
they identified 14 out of their original 41 as candidate third - body systems , with others being identified as likely caused by starspots , pulsations , and apsidal motion . of their 14 candidate third - body systems ,
all 14 are still in the _ kepler _ eb catalog , 12 fall under the scope of this paper , and 9 appear in our list of third - body signals . those that appear in both lists
are noted in table [ tableresults ] .
binaries that they list as candidate third - body systems , but we do not , either show significant noise or would be very long period ltte orbits .
@xcite recently reported 39 triple - star _
binaries due to etv signatures . of these 39 , 21 fall under the scope of this paper , 19 of which also appear in our list of third - body signals , with the other two determined to be unlikely caused by a third - body due to their very short periods and notable spot activity .
the detections that overlap both of these studies are also noted in table [ tablerapp ] . in most cases ,
the model fits from both studies are consistent . in general ,
due to our treatment of the full _ kepler _ dataset now available , tertiary parameters should now be more precise and longer period third body signals are now more apparent .
any disagreement is likely due to a slightly differ inner - binary ephemeris or the addition of the physical delay in their models ( discussed further below ) .
@xcite presents analytic functions for the light time travel effect ( ltte ) component of the etv residual signal . using the same form as @xcite ,
the timings can be expressed by    @xmath17\ ] ]    where @xmath18 p_3^{2/3 } \end{aligned}\ ] ] and @xmath19 is a time offset , @xmath20 is the mass of the third body , @xmath21 is the mass of the entire system , and @xmath22 , @xmath23 , @xmath24 , @xmath25 , @xmath26 , and @xmath27 are the period , inclination , eccentricity , argument of periastron , eccentric anomaly , and mean anomaly of the third body orbit , respectively .
this expression was then used to fit all etvs flagged as potential third body signals .
the period was first estimated using lomb - scargle periodogram and used as input into a series of levenberg - marquardt fits , each using a different starting guess for eccentricity .
the fit with the lowest chi - squared was then kept and the errors estimated from the covariance matrix . if the final fit had an eccentricity consistent with 0 , then @xmath24 was set to 0 , @xmath25 to @xmath28 , and the fitting was redone with circular constraints to get appropriate error estimates on the remaining parameters .
this gives values and estimated errors for @xmath22 , @xmath24 , and @xmath29 ( table [ tableresults ] ) .
a sample of some of these etv signals and their respective fits can be seen in fig .
[ etv_ltte ] , with fits for all candidate third body signals available in the online version of the _ kepler _ eb catalog
. we can only provide estimate periods for the sample of etv signals with less than one full cycle in the data .
even these periods can be significantly biased based on the section of the cycle that is in the observed baseline and should be treated with reservation .
these very long period cases are provided separately at the end of table [ tableresults ]      @xcite included physical delays in their models of 39 _ kepler _ binaries with possible third - body etvs , sometimes contributing largely to the overall model .
this dynamical effect occurs when the presence of a third body changes the period of the inner binary .
[ phys_hist ] shows the distribution in their targets and the overlapping targets of the ratio of the amplitude of the physical delay compared to the total amplitude in the etv signal .
21 of their targets overlap with ours , but due to the short - period inner - binary , the physical delay rarely contributes significantly to these model . from their results , it seems that the ltte effect dominates over the physical delay for binaries with periods less than 3 days , which covers the vast majority of our targets .
for some of these binaries with ltte signals , tertiary eclipses have also been found that confirm the presence and third body period , and significantly constrain the inclination of the third body .
any binary which was identified to have a possible third body due to its etv signal and also has a detected tertiary eclipse is noted as such in table [ tableresults ] .
kic 2856960 , for instance , has an inner - binary period of 0.259 days with an etv signal resulting in a ltte fit with a period of @xmath30 days .
this period is consistent with the previously determined period for the tertiary events of 204.25 days ( fig .  [
2856960 ] ) .
this is also consistent with the ltte period of @xmath31 days reported by @xcite , and the tertiary event period of @xmath32 days by @xcite .    in the case of kic 2835289 ( fig .
[ 2835289 ] ) , we have only observed one potential tertiary event in q9 . without at least three consecutive events
, we can not rigorously confirm that the eclipse is a third - body as opposed to a blended eclipsing binary .
however , the eclipse seems to show the eclipse of both stars in the inner - binary and etv signal shows a possible long - term third body orbit suggesting a period of approximately 800 days .
if this proves to be a true third body , then _ kepler _ just missed an event before the beginning of the mission and may have observed another event in q17 , which has yet to be processed .
kic 6543674 also shows a single tertiary eclipse in q2 .
a second tertiary eclipse was missed during a break in the _ kepler _ data , but we were able to observe an additional tertiary event from the ground , giving a third body period of @xmath33 days @xcite .
in this case , we do not have a full orbit of the etv signal and the ltte model period is quite uncertain .
nine binaries that show third - body etv signals ( kic 3936357 , 4069063 , 5310387 , 6629588 , 7375612 , 8122124 , 8758716 , 10014830 , and 10855535 ) also show constant changes in their eclipse depths ( fig .  [ dv ] ) , which could either be caused by a change in inclination or apsidal motion perhaps induced by the third body .
we plan to follow these up later with full photodynamical models .
it is also possible that some of these etvs could be composed of multiple signals .
kic 5310387 , 6144827 , 8145477 , 11612091 , and 11825204 , for example , may have both an ltte and quadratic component or two ltte signals as is shown in the residuals in fig .
in general , the stronger signal is fitted and noted .
in this study we find a third body rate of in our sample of close binaries , nearly all of which have inner binary periods shorter than 3 days ( fig .  [ ptrip ] ) .
this is much lower than the third body rate of 96% found by the previous studies mentioned .
however , our identification of tertiary companions is certainly a lower limit for several reasons .
first , our ability to detect a third body is very sensitive to both inclination and mass of the third body , such that low - mass tertiaries and/or tertiaries whose orbital planes are highly inclined relative to the inner binary orbital plane do not present detectable ltte effects . of our total sample of 1279 binaries , 3 ( )
show an ltte orbit and visible tertiary eclipses .
111 ( ) have ltte orbits with periods shorter than the span of our photometric data but do not show tertiary eclipses , suggesting that the eclipses fell in a gap in the data or the orbits are not well enough aligned to show eclipses .
thus there is evidence from these examples that in a few percent of cases we are indeed missing true third bodies because of inclination non - alignments .
94 ( ) have ltte orbits with periods longer than the photometric baseline . in these cases
we do not have well constrained periods and our chances of detecting a tertiary eclipse are slim .
a second reason that our determination of the third - body occurrence is likely a lower limit is that the very close binaries that comprise our sample here generally present more noise in the etv signal , which could easily bury a weak ltte signal .
we have employed a method that minimizes false positives due to spurious etv signals , and thus necessarily have eliminated some potentially true ltte signals .
third , and perhaps most important , the limited timespan of the currently available _ kepler _ data ( @xmath11400 days ) significantly restricts us to detect third bodies with orbital periods comparable to or shorter than 1400 days . relative to the full span of tertiary separations found in previous works @xcite , with separations as large as @xmath11 pc , we are at present sampling only the relatively closest tertiary companions .
indeed , @xcite found among tight binaries that the rate of third bodies with orbital periods less than @xmath13 years ( comparable to our limit based on the duration of the available _ kepler _ data ) is 15% @xmath34 3% .
thus our finding of a third - body occurrence rate with a period less than 1400 days of is compatible with the expected rate , though it appears we are likely still missing a fraction of some systems for the reasons already mentioned .    the distribution of periods of potential third body orbits is also shown in fig .  [ ptrip ] .
we can clearly see a falloff in detection past the current length of the _ kepler _ mission of @xmath1 1400 days , as expected .
however , for third - body periods shorter than @xmath1 1400 days , for which our detectability is relatively good , the occurrence rate does appear to increase toward longer third - body periods , consistent with the period distribution of third bodies among tight binaries found by @xcite .
furthermore , we find that the triples on the widest orbits are found around the shortest period binaries , which is consistent with models that tighten the inner binary orbit through the presence , and gradual widening , of a companion .
we presented our technique for computing precise eclipse timings for 1279 close eclipsing binaries in the _ kepler _ eclipsing binary catalog .
these precise eclipse timings are complemented by the eclipse timings to be reported by @xcite for longer period , detached ebs . for the ebs
whose timings are reported here , our method has been developed specifically to deal with the challenge of constantly changing light levels arising from spots and other phenomena that distort the light curves and could cause spurious eclipse timing variation ( etv ) signals .
ebs with etv signals suggesting the possible presence of a third body have been identified and have been fit with a ltte orbit model in order to determine the likely parameters of the third bodies . in the current sample of 1279 close ebs , we have identified 236 that likely have tertiary companions .
the parameters of these fits are also available online and are updated as new data become available .    our measured occurrence rate of of close binaries with tertiary companions with periods up to @xmath35 days ( limited by the current timespan of _ kepler _ data ) , appears to be broadly consistent with the expectation that @xmath36% of close binaries will have tertiaries of such periods @xcite .
indeed , we already find in our data that the periods of third bodies rise among the tightest binaries , consistent with previous work that has found a very high rate of third bodies in very wide orbits around the tightest binaries , presumably the result of dynamical tightening of inner binaries through widening of the tertiary .
eclipse timings for all ebs are updated in real - time and are freely available as a community resource at ` http://keplerebs.villanova.edu ` .
we thank nathan de lee and phil cargile for helpful discussions and darin ragozzine , eric ford , and joshua pepper for their feedback .    this project is supported through the _ kepler _ participating scientist award nsr303065 .
kec and kgs acknowledge support from nasa adap grant nnx12ae22 g .
jao and wfw gratefully acknowledge support from the nsf via grant ast-1109928 , and also nasa via the _ kepler _ psp grant nnx12ad23 g .
applegate , j.  h.  1992 , , 385 , 621 armstrong , d. , pollacco , d. , watson , c.  a. , et al .
2012 , , 545 , l4 batalha , n.  m. , borucki , w.  j. , koch , d.  g. , et al .
2010 , , 713 , l109 bonnell , i.  a.  2001 , the formation of binary stars , 200 , 23 borkovits , t. , csizmadia , s. , forgcs - dajka , e. , & hegeds , t.  2011 , , 528 , a53 borucki , w.  j. , koch , d. , basri , g. , et al .
2010 , science , 327 , 977 dhital , s. , west , a.  a. , stassun , k.  g. , & bochanski , j.  j.  2010 , , 139 , 2566 doyle , l.  r. , carter , j.  a. , fabrycky , d.  c. , et al .
2011 , science , 333 , 1602 fabrycky , d.  c. , ford , e.  b. , steffen , j.  h. , et al .
2012 , , 750 , 114 fabrycky , d. , & tremaine , s.  2007 , , 669 , 1298 ford , e.  b. , ragozzine , d. , rowe , j.  f. , et al .
2012 , , 756 , 185 gies , d.  r. , williams , s.  j. , matson , r.  a. , et al .
2012 , , 143 , 137 hilditch , r.  w.  2001 , an introduction to close binary stars , by r.  w.  hilditch , pp .
isbn 0521241065 .
cambridge , uk : cambridge university press , march 2001 .
jenkins , j.  m. , caldwell , d.  a. , chandrasekaran , h. , et al .
2010 , , 713 , l120 law , n.  m. , dhital , s. , kraus , a. , stassun , k.  g. , & west , a.  a.  2010 , , 720 , 1727 lee , j.  w. , kim , s .- l . , lee , c .- u . , et al .
2013 , , 763 , 74 kirk et .
al  2013 ( in preparation ) matijevi , g. , pra , a. , orosz , j.  a. , et al .
2012 , , 143 , 123 orosz et al .  2013 ( in preparation ) pra , a. , batalha , n. , slawson , r.  w. , et al .
2011 , , 141 , 83 pra , a. , guinan , e.  f. , devinney , e.  j. , et al .
2008 , , 687 , 542 rappaport , s. , deck , k. , levine , a. , et al .
2013 , arxiv:1302.0563 reipurth , b. , & mikkola , s.  2012 , , 492 , 221 schwarz , g.  1978 .
the annals of statistics , 6 , 461 schwarz , r. , haghighipour , n. , eggl , s. , pilat - lohinger , e. , & funk , b.  2011 , , 414 , 2763 slawson , r.  w. , pra , a. , welsh , w.  f. , et al .
2011 , , 142 , 160 steffen , j.  h. , ford , e.  b. , rowe , j.  f. , et al .
2012 , , 756 , 186 thackeray - lacko , b. , hill , m. , orosz , j.  a. , et al .  2013 , american astronomical society meeting abstracts , 221 , # 142.09 tokovinin , a.  a.  1997 , astronomy letters , 23 , 727 tokovinin , a. , thomas , s. , sterzik , m. , & udry , s.  2006 , , 450 , 681 tran , k. , levine , a. , rappaport , s. , et al .
2013 , , 774 , 81 welsh , w.  f. , orosz , j.  a. , carter , j.  a. , et al .
2012 , , 481 , 475    rrrrr 1433410 & 56107.275 & -6.307 & 39.658 & primary + 1433410 & 56107.275 & -4.493 & 20.390 & entire + 1433410 & 56107.416 & -47.606 & 40.954 & secondary + 1433410 & 56107.558 & -10.541 & 18.835 & entire + 1433410 & 56107.558 & 21.341 & 30.499 & primary + 1433980 & 55740.148 & -443.232 & 143.078 & secondary + 1433980 & 55740.944 & -283.392 & 70.502 & primary + 1433980 & 55740.944 & -170.899 & 69.552 & entire + 1433980 & 55741.741 & -161.482 & 90.634 & secondary + 1433980 & 55742.537 & -248.659 & 61.776 & primary +    rrr [ tablemt ] 2305372 & 3104113 & 3765708 + 4074532 & 4851217 & 4853067 + 5020034 & 5471619 & 5770860 + 5792093 & 6044064 & 6044543 + 6066379 & 6213131 & 6314173 + 6464285 & 6677225 & 7696778 + 7938468 & 7938870 & 8758161 + 9087918 & 9402652 & 9840412 + 9934052 & 10030943 & 10292413 + 10736223 & 11097678 &
11144556 + 11924311 & & +    rrrrrrr 3228863 & 668.4 & @xmath37 & 0.08\{0.06,0.12 } & @xmath38 & 189\{187,194 } & @xmath39 + 4647652 & 753.5 & @xmath40 & 0.35\{0.10,0.44 } & @xmath41 & 228\{183,274 } & @xmath42 + 4909707 & 505.3 & @xmath43 & 0.54\{0.31,0.66 } & @xmath44 & 493\{378,627 } & @xmath45 + 5128972 & 447.8 & @xmath46 & 0.33\{0.25,0.41 } & @xmath47 & 259\{244,271 } & @xmath48 + 5264818 & 296.3 & @xmath49 & 0.37\{0.13,0.53 } & @xmath50 & 145\{107,196 } & @xmath51 + 5310387 & 214.2 & @xmath52 & 0.53\{0.34,0.61 } & @xmath53 & 31\ { 27 , 37 } & @xmath54 + 5376552 & 334.5 & @xmath55 & 0.40\{0.35,0.45 } & @xmath56 & 94\ { 91 , 98 } & @xmath57 + 6370665 & 285.9 & @xmath58 & 0.22\{0.07,0.33 } & @xmath59 & 67\ { 61 , 74 } & @xmath60 + 6531485 & 48.3 & & 0.44\{0.33,0.63 } & & 72\ { 31,109 } & + 7690843 & 74.3 & @xmath61 & 0.25\{0.08,0.42 } & @xmath62 & 71\ { 51 , 91 } & @xmath63 + 8043961 & 476.7 & @xmath64 & 0.25\{0.14,0.33 } & @xmath65 & 194\{179,213 } & @xmath66 + 8192840 & 803.9 & @xmath67 & 0.63\{0.52,0.70 } & @xmath68 & 208\{187,223 } & @xmath69 + 8386865 & 293 & @xmath70 & 0.38\{0.27,0.48 } & @xmath71 & 171\{156,210 } & @xmath72 + 8394040 & 394.8 & @xmath73 & 0.61\{0.50,0.67 } & @xmath74 & 369\{345,391 } & @xmath75 + 8904448 & 548.1 & @xmath76 & 0.59\{0.50,0.66 } & @xmath77 & 171\{158,192 } & @xmath78 + 9451096 & 106.7 & @xmath79 & 0.24\{0.10,0.36 } & @xmath80 & 90\ { 59,144 } & @xmath81 + 9722737 & 443.9 & @xmath82 & 0.22\{0.16,0.27 } & @xmath83 & 230\{225,236 } & @xmath84 + 9912977 & 753.7 & @xmath85 & 0.31\{0.16,0.39 } & @xmath86 & 105\ { 94,117 } & @xmath87 + 10226388 & 934.9 & @xmath88 & 0.32\{0.24,0.39 } & @xmath89 & 465\{434,493 } & @xmath90 + 10991989 & 554.2 & @xmath91 & 0.30\{0.21,0.37 } & @xmath92 & 256\{239,274 } & @xmath93 + 11042923 & 839 & @xmath94 & 0.17\{0.09,0.25 } & @xmath95 & 223\{213,230 } & @xmath96 +    rrrrrr 2856960 & 0.60 & 0.259 & @xmath97 & @xmath98 & @xmath99 + 3228863 & 0.65 & 0.731 & @xmath37 & @xmath38 & @xmath39 + 3245776 & 0.96 & 1.492 & @xmath100 & @xmath101 & @xmath102 + 3641446 & 0.95 & 2.100 & @xmath103 & @xmath104 & @xmath105 + 4037163 & 0.58 & 0.635 & @xmath106 & @xmath107 & @xmath108 + 4909707 & 0.72 & 2.302 & @xmath43 & @xmath44 & @xmath45 + 5128972 & 0.74 & 0.505 & @xmath46 & @xmath47 & @xmath48 + 5264818 & 0.91 & 1.905 & @xmath49 & @xmath50 & @xmath51 + 5310387 & 0.96 & 0.442 & @xmath52 & @xmath53 & @xmath54 + 5376552 & 0.82 & 0.504 & @xmath55 & @xmath56 & @xmath57 + 5459373 & 0.97 & 0.287 & @xmath109 & @xmath110 & @xmath111 + 5560831 & 0.60 & 0.868 & @xmath112 & @xmath113 & @xmath114 + 6302592 & 0.93 & 1.578 & @xmath115 & @xmath116 & @xmath117 + 6370665 & 0.94 & 0.932 & @xmath58 & @xmath59 & @xmath60 + 7362751 & 0.73 & 0.338 & @xmath118 & @xmath119 & @xmath120 + 7657914 & 0.72 & 0.475 & @xmath121 & @xmath122 & @xmath123 + 7685689 & 0.77 & 0.325 & @xmath124 & @xmath125 & @xmath126 + 7690843 & 0.64 & 0.786 & @xmath61 & @xmath62 & @xmath63 + 8043961 & 0.63 & 1.559 & @xmath64 & @xmath65 & @xmath66 + 8145477 & 0.88 & 0.566 & @xmath127 & @xmath128 & @xmath129 + 8190491 & 0.95 & 0.778 & @xmath130 & @xmath38 & @xmath131 + 8211618 & 0.73 & 0.337 & @xmath132 & @xmath133 & @xmath134 + 8330092 & 0.80 & 0.322 & @xmath135 & @xmath136 & @xmath137 + 8386865 & 0.99 & 1.258 & @xmath70 & @xmath71 & @xmath72 + 8394040 & 0.77 & 0.302 & @xmath73 & @xmath74 & @xmath75 + 8904448 & 0.74 & 0.866 & @xmath76 & @xmath77 & @xmath78 + 9075704 & 0.68 & 0.513 & @xmath138 & @xmath139 & @xmath140 + 9451096 & 0.54 & 1.250 & @xmath79 & @xmath80 & @xmath81 + 9706078 & 0.55 & 0.614 & @xmath141 & @xmath142 & @xmath143 + 9722737 & 0.78 & 0.419 & @xmath82 & @xmath83 & @xmath84 + 9994475 & 0.76 & 0.318 & @xmath144 & @xmath145 & @xmath146 + 10014830 & 0.61 & 3.031 & @xmath147 & @xmath148 & @xmath149 + 10855535 & 0.98 & 0.113 & @xmath150 & @xmath151 & @xmath152 + 10991989 & 0.54 & 0.974 & @xmath91 & @xmath92 & @xmath93 + 11247386 & 0.75 & 0.394 & @xmath153 & @xmath154 & @xmath155 + 2302092 & 0.89 & 0.295 & @xmath156 & @xmath157 & @xmath158 + 2450566 & 0.98 & 1.845 & @xmath159 & @xmath160 & @xmath161 + 2835289 & 0.94 & 0.858 & @xmath162 & @xmath163 & @xmath164 + 3839964 & 0.78 & 0.256 & @xmath165 & @xmath166 & @xmath167 + 4069063 & 0.56 & 0.504 & @xmath168 & @xmath169 & @xmath170 + 4138301 & 0.90 & 0.253 & @xmath171 & @xmath172 & @xmath173 + 4244929 & 0.91 & 0.341 & @xmath174 & @xmath175 & @xmath176 + 4451148 & 0.82 & 0.736 & @xmath177 & @xmath178 & @xmath179 + 4547308 & 0.88 & 0.577 & @xmath180 & @xmath38 & @xmath181 + 4647652 & 0.68 & 1.065 & @xmath40 & @xmath41 & @xmath42 + 4670267 & 0.60 & 2.006 & @xmath182 & @xmath183 & @xmath184 + 4681152 & 0.55 & 1.836 & @xmath185 & @xmath186 & @xmath187 + 4762887 & 0.95 & 0.737 & @xmath188 & @xmath65 & @xmath189 + 4859432 & 0.76 & 0.385 & @xmath190 & @xmath191 & @xmath192 + 4937217 & 0.82 & 0.429 & @xmath193 & @xmath194 & @xmath195 + 4945857 & 0.74 & 0.335 & @xmath196 & @xmath56 & @xmath197 + 5269407 & 0.53 & 0.959 & @xmath198 & @xmath38 & @xmath199 + 5478466 & 0.97 & 0.483 & @xmath200 & @xmath201 & @xmath202 + 5611561 & 0.74 & 0.259 & @xmath203 & @xmath65 & @xmath204 + 5790912 & 0.77 & 0.383 & @xmath205 & @xmath206 & @xmath207 + 5791886 & 0.76 & 0.325 & @xmath208 & @xmath209 & @xmath210 + 5975712 & 0.87 & 1.136 & @xmath211 & @xmath212 & @xmath213 + 6050116 & 0.77 & 0.240 & @xmath214 & @xmath215 & @xmath216 + 6118779 & 0.90 & 0.364 & @xmath217 & @xmath218 & @xmath219 + 6281103 & 0.98 & 0.363 & @xmath220 & @xmath221 & @xmath222 + 6469946 & 0.51 & 0.716 & @xmath223 & @xmath224 & @xmath225 + 6516874 & 0.60 & 0.916 & @xmath226 & @xmath227 & @xmath228 + 6543674 & 0.53 & 2.391 & @xmath229 & @xmath230 & @xmath231 + 6615041 & 0.75 & 0.340 & @xmath232 & @xmath233 & @xmath234 + 6629588 & 0.55 & 2.264 & @xmath235 & @xmath236 & @xmath237 + 6671698 & 0.73 & 0.472 & @xmath238 & @xmath239 & @xmath240 + 6766325 & 0.92 & 0.440 & @xmath241 & @xmath242 & @xmath243 + 7035139 & 0.79 & 0.310 & @xmath244 & @xmath245 & @xmath246 + 7119757 & 0.64 & 0.743 & @xmath247 & @xmath248 & @xmath249 + 7272739 & 0.75 & 0.281 & @xmath250 & @xmath251 & @xmath252 + 7385478 & 0.54 & 1.655 & @xmath253 & @xmath254 & @xmath255 + 7518816 & 0.65 & 0.467 & @xmath256 & @xmath257 & @xmath258 + 7877062 & 0.81 & 0.304 & @xmath259 & @xmath260 & @xmath81 + 8045121 & 0.99 & 0.263 & @xmath261 & @xmath262 & @xmath263 + 8192840 & 0.95 & 0.434 & @xmath67 & @xmath68 & @xmath69 + 8242493 & 0.73 & 0.283 & @xmath264 & @xmath38 & @xmath265 + 8563964 & 1.00 & 0.338 & @xmath266 & @xmath262 & @xmath267 + 8690104 & 0.77 & 0.409 & @xmath268 & @xmath269 & @xmath270 + 8739802 & 0.93 & 0.275 & @xmath271 & @xmath272 & @xmath273 + 8957887 & 0.76 & 0.347 & @xmath274 & @xmath275 & @xmath276 + 8982514 & 0.84 & 0.414 & @xmath277 & @xmath262 & @xmath278 + 9091810 & 0.69 & 0.480 & @xmath279 & @xmath38 & @xmath280 + 9101279 & 0.58 & 1.811 & @xmath281 & @xmath282 & @xmath283 + 9272276 & 0.78 & 0.281 & @xmath284 & @xmath285 & @xmath286 + 9283826 & 0.84 & 0.357 & @xmath287 & @xmath56 & @xmath288 + 9353234 & 0.86 & 1.487 & @xmath289 & @xmath290 & @xmath291 + 9412114 & 0.85 & 0.250 & @xmath292 & @xmath293 & @xmath294 + 9532219 & 0.74 & 0.198 & @xmath295 & @xmath110 & @xmath296 + 9592145 & 0.65 & 0.489 & @xmath297 & @xmath298 & @xmath299 + 9612468 & 1.00 & 0.133 & @xmath300 & @xmath301 & @xmath302 + 9665086 & 0.67 & 0.297 & @xmath303 & @xmath304 & @xmath305 + 9821923 & 0.95 & 0.350 & @xmath306 & @xmath307 & @xmath308 + 9838047 & 0.84 & 0.436 & @xmath309 & @xmath310 & @xmath311 + 9912977 & 0.59 & 1.888 & @xmath85 & @xmath86 & @xmath87 + 10226388 & 0.77 & 0.661 & @xmath88 & @xmath89 & @xmath90 + 10275197 & 0.78 & 0.391 & @xmath312 & @xmath313 & @xmath314 + 10322582 & 0.86 & 0.291 & @xmath315 & @xmath316 & @xmath317 + 10383620 & 0.64 & 0.735 & @xmath318 & @xmath319 & @xmath320 + 10388897 & 0.99 & 0.344 & @xmath321 & @xmath322 & @xmath323 + 10724533 & 0.75 & 0.745 & @xmath324 & @xmath325 & @xmath326 + 10727655 & 0.73 & 0.353 & @xmath327 & @xmath262 & @xmath328 + 10848807 & 0.74 & 0.346 & @xmath329 & @xmath330 & @xmath331 + 10905804 & 0.68 & 0.751 & @xmath332 & @xmath333 & @xmath334 + 10916675 & 0.87 & 0.419 & @xmath335 & @xmath336 & @xmath337 + 10934755 & 0.68 & 0.786 & @xmath338 & @xmath339 & @xmath340 + 11042923 & 0.76 & 0.390 & @xmath94 & @xmath95 & @xmath96 + 11246163 & 0.77 & 0.279 & @xmath341 & @xmath342 & @xmath343 + 11347875 & 0.86 & 3.455 & @xmath344 & @xmath345 & @xmath346 + 11604958 & 0.72 & 0.299 & @xmath347 & @xmath348 & @xmath349 + 11716688 & 0.94 & 0.301 & @xmath350 & @xmath351 & @xmath352 + 11825204 & 0.98 & 0.210 & @xmath353 & @xmath354 & @xmath355 + 12019674 & 0.76 & 0.355 & @xmath356 & @xmath357 & @xmath358 + 12055255 & 0.90 & 0.221 & @xmath359 & @xmath360 & @xmath361 + 12071741 & 0.94 & 0.314 & @xmath362 & @xmath363 & @xmath364 + 12458133 & 0.76 & 0.333 & @xmath365 & @xmath366 & @xmath367 + 3114667 & 0.52 & 0.889 & @xmath368 & & + 2983113 & 0.89 & 0.395 & @xmath369 & & + 4066203 & 0.93 & 0.363 & @xmath370 & & + 4074708 & 0.73 & 0.302 & @xmath33 & & + 4904304 & 0.89 & 0.389 & @xmath369 & & + 5282464 & 0.72 & 0.496 & @xmath33 & & + 5956776 & 0.61 & 0.569 & @xmath370 & & + 6153219 & 0.62 & 0.530 & @xmath370 & & + 6265720 & 0.92 & 0.312 & @xmath369 & & + 6794131 & 0.81 & 1.613 & @xmath369 & & + 7506164 & 0.88 & 0.558 & @xmath369 & & + 7590728 & 0.69 & 0.477 & @xmath371 & & + 9544350 & 0.92 & 2.260 & @xmath372 & & + 9776718 & 0.87 & 0.544 & @xmath373 & & + 9788457 & 0.60 & 0.963 & @xmath372 & & + 10095469 & 0.60 & 0.678 & @xmath373 & & + 10796477 & 0.74 & 0.485 & @xmath33 & & + 2159783 & 0.87 & 0.374 & @xmath374 & & + 3127873 & 0.91 & 0.672 & @xmath375 & & + 3221207 & 0.81 & 0.474 & @xmath376 & & + 3342425 & 0.93 & 0.393 & @xmath377 & & + 3766353 & 0.53 & 2.667 & @xmath378 & & + 3848042 & 0.99 & 0.411 & @xmath379 & & + 3935319 & 0.75 & 0.353 & @xmath35 & & + 3936357 & 0.76 & 0.369 & @xmath380 & & + 4077442 & 0.58 & 0.693 & @xmath381 & & + 4464999 & 0.77 & 0.434 & @xmath382 & & + 4563150 & 0.79 & 0.275 & @xmath383 & & + 4758368 & 0.57 & 3.750 & @xmath384 & & + 4945588 & 0.99 & 1.129 & @xmath384 & & + 5008287 & 0.94 & 0.292 & @xmath385 & & + 5015926 & 0.75 & 0.363 & @xmath35 & & + 5097446 & 0.60 & 1.288 & @xmath378 & & + 5296877 & 0.95 & 0.377 & @xmath374 & & + 5353374 & 0.78 & 0.393 & @xmath376 & & + 5389616 & 0.99 & 0.407 & @xmath386 & & + 5513861 & 0.57 & 1.510 & @xmath387 & & + 5770431 & 0.89 & 0.392 & @xmath388 & & + 5820209 & 0.81 & 0.656 & @xmath389 & & + 5951553 & 0.95 & 0.432 & @xmath385 & & + 6144827 & 0.79 & 0.235 & @xmath390 & & + 6187893 & 0.59 & 0.789 & @xmath391 & & + 6287172 & 0.95 & 0.204 & @xmath35 & & + 6370361 & 0.84 & 0.455 & @xmath377 & & + 7137798 & 0.51 & 2.254 & @xmath392 & & + 7269843 & 0.79 & 0.268 & @xmath377 & & + 7339345 & 0.74 & 0.260 & @xmath393 & & + 7367833 & 0.76 & 0.286 & @xmath394 & & + 7375612 & 0.97 & 0.160 & @xmath386 & & + 7680593 & 0.97 & 0.276 & @xmath384 & & + 7697065 & 0.69 & 0.273 & @xmath386 & & + 7709086 & 0.78 & 0.409 & @xmath383 & & + 7773380 & 0.94 & 0.308 & @xmath35 & & + 7871200 & 0.74 & 0.243 & @xmath395 & & + 7878402 & 0.73 & 0.374 & @xmath396 & & + 8016214 & 0.53 & 3.175 & @xmath386 & & + 8122124 & 1.00 & 0.249 & @xmath389 & & + 8189196 & 0.98 & 2.304 & @xmath397 & & + 8222945 & 0.99 & 0.451 & @xmath385 & & + 8231231 & 0.89 & 0.712 & @xmath389 & & + 8257903 & 0.76 & 0.515 & @xmath398 & & + 8265951 & 0.81 & 0.780 & @xmath383 & & + 8703528 & 0.74 & 0.400 & @xmath384 & & + 8715667 & 0.85 & 0.406 & @xmath377 & & + 8758716 & 1.00 & 0.107 & @xmath388 & & + 9026766 & 0.75 & 0.272 & @xmath399 & & + 9083523 & 0.65 & 0.918 & @xmath400 & & + 9097798 & 0.99 & 0.334 & @xmath401 & & + 9181877 & 0.74 & 0.321 & @xmath379 & & + 9347868 & 0.82 & 0.318 & @xmath402 & & + 9657096 & 0.94 & 2.138 & @xmath35 & & + 9724080 & 0.94 & 1.174 & @xmath382 & & + 9882280 & 0.75 & 0.289 & @xmath403 & & + 9956124 & 0.91 & 0.363 & @xmath404 & & + 10007533 & 0.88 & 0.648 & @xmath405 & & + 10135584 & 0.96 & 0.391 & @xmath375 & & + 10216186 & 0.64 & 0.606 & @xmath406 & & + 10228991 & 0.97 & 2.799 & @xmath389 & & + 10229723 & 0.83 & 0.629 & @xmath407 & & + 10481912 & 0.96 & 0.442 & @xmath408 & & + 10485137 & 0.71 & 0.445 & @xmath403 & & + 10557008 & 0.77 & 0.265 & @xmath384 & & + 10711938 & 0.95 & 0.358 & @xmath383 & & + 11091082 & 0.82 & 0.385 & @xmath35 & & + 11151970 & 0.87 & 0.312 & @xmath409 & & + 11305087 & 0.78 & 0.309 & @xmath410 & & + 11496078 & 0.87 & 0.300 & @xmath411 & & + 11509282 & 0.83 & 0.634 & @xmath389 & & + 11566174 & 0.76 & 0.277 & @xmath412 & & + 11717798 & 0.83 & 0.375 & @xmath413 & & + 11805235 & 0.79 & 0.395 & @xmath406 & & + 11910076 & 0.96 & 0.348 & @xmath414 & & + 12055421 & 0.96 & 0.386 & @xmath389 & & + 12267718 & 0.89 & 0.545 & @xmath415 & & + 12554536 & 0.63 & 0.684 & @xmath394 & & +