many years of effort have been spent on trying to understand the structure and evolution of the super star cluster ngc2070 , at the centre of the 30doradus nebula in the large magellanic cloud ( otherwise known as the tarantula nebula ) since its discovery as such by lacaille in 1751 . in this case
, the prize is clearly worth the effort as this object may hold the key to answers to many fundamental questions regarding the physics of the formation of stars , stellar clusters and galaxies and their role in the cosmos ( e.g. brandl 2005 ) .
thanks also to this effort , we now know much about the brightest stars more massive than @xmath5m@xmath6that can be more easily detected and resolved . in spite of the rather large distances ( @xmath7kpc ) involved ( hunter et al .
1995a ) , these objects can be studied even in the very compact core of ngc2070 , referred to as r136 , which was first resolved unambiguously by the hubble space telescope ( hst ) using the faint object camera ( weigelt et al .
1991 ; de marchi et al .
but we still know very little about the other part of this population consisting of the fainter , low mass main sequence ( ms ) and pre - ms ( pms ) stars and brown dwarfs . these might be present in even larger numbers if the tapered power law imf proposed recently by de marchi , paresce & portegies zwaart ( 2010 ) for young and very old clusters continues to the hydrogen burning limit and beyond .
their number and spatial distribution are critical for understanding the low mass star formation process in a cluster of very massive stars and their role in determining the ultimate fate of the cluster itself .
the current observational situation regarding low mass stars in ngc2070 is contradictory .
although there is good evidence for a significant population of these stars , the precise nature of their spatial and mass distributions is not well understood .
sirianni et al ( 2000 ) , using @xmath8- and @xmath9-band hst observations claim that the mass function ( mf ) flattens below @xmath10m@xmath6 , while zinnecker et al .
( 2002 ) and more recently andersen et al .
( 2010 ) , using hst observations in the h band , find no such feature in their data down to their completeness limit around 1m@xmath6 .
so far , little can be said of what happens below this limit , mainly because these stars in ngc2070 and its compact core r136 are very difficult to detect and characterize accurately .
first , because of their inherent faintness especially with respect to the high background due to the numerous much brighter stars in the field .
second , because they scatter widely across the colour  magnitude diagram ( cmd ) making accurate subtraction of field stars and older cluster pms stars and subdwarfs near the ms unreliable without useful separation criteria .
finally , the patchy extinction across the face of the nebula ( e.g. indebetouw et al . 2009 ) needs to be properly accounted for in order to accurately determine their true brightness .    for these reasons , use of the wide field camera 3 ( wfc3 ) , the new panchromatic camera on hst , is , in principle , ideal for the task of establishing the main characteristics of the low mass population in ngc2070 . with its high spatial resolution over a wide field encompassing most of the central part of the cluster and its unprecedented sensitivity in the near uv and ir , the wfc3 allows us to tackle for the first time most of the problems summarised above with some hope of resolving them . in this paper , we present the first results of a deep survey of ngc2070 obtained in 2009 october as part of the early release science programme with the wfc3 camera .
this project has as its ultimate objective that of determining as accurately as possible the properties of the low - mass initial mf in this cluster , which could be a prototype of the much larger and distant unresolved super star clusters in starburst galaxies . in order for this type of study to become possible ,
we must first carefully identify and distinguish from one another stars of different ages , since several generations of objects are known to populate the 30dor region ( e.g. walborn & blades 1997 ) .
so far , studies of the star formation history of 30dor have been limited to the most massive stars , typically above 20m@xmath6 , and hence to the past 10myr or so , because these works rely on spectroscopy of individual objects that can not easily be obtained for less massive stars ( e.g. selman et al . 1999 ) .
clearly , access to a wider mass range is crucial , so the goal of this paper is to understand how star formation has proceeded for much smaller objects , i.e. stars with masses as low as @xmath11m@xmath6 .
therefore , in this work we concentrate on a systematic search for low - mass stars in the pms phase , which we can identify in a reliable way and which allow us to study how star formation has proceeded in this area over the past @xmath3myr .
only then will it be possible to address the properties of the mf , which will be the subject of future work .
the structure of this paper is as follows .
we describe the observations in section2 , the reduction of the data in section3 and the photometric analysis in section4 .
the correction for extinction and the resultant corrected cmds are presented in section5 .
the hertzsprung
russell ( hr ) diagram and its analysis with recent pms models are the topic of section6 , where we also discuss the preliminary physical implications of our results on the star formation history in these regions .
a summary is provided in section7 .
the observations discussed in this paper were obtained in 2009 october ( 20  27 ) with the wfc3 on board the hst .
we refer the reader to dressel et al .
( 2010 ) for information about the wfc3 instrument and its performances .
six broad - band and one narrow - band filters were used in these observations , as shown in table[tab1 ] , where the number of exposures and the total exposure times in each filter are also given .
note that an extensive dithering pattern was employed , combining both long and short exposures , in order to improve the sampling of the telescope s point spread function ( psf ) , to minimise the effects of detector blemishes , to cover the gap between the two detectors in the field of view , and to recover the magnitudes of the brightest objects that would otherwise saturate in long exposures .
lcr[bh ] f336w & 24 & @xmath12 + f438w & 16 & @xmath13 + f555w & 20 & @xmath14 + f656n & 8 & @xmath15 + f814w & 20 & @xmath16 + f110w & 9 & @xmath17 + f160w & 12 & @xmath18 [ tab1 ]    for illustration purposes , figure[fig0 ] presents a four - colour composite image of the observed field .
the colours result from assigning different hues to each monochromatic image , as indicated in the figure caption .
the bright compact star cluster in the upper right - hand side of the frame is r136 , the core of the massive ngc2070 association .
a cursory glance at figure[fig0 ] already shows that the centre of the cluster is dominated by bright blue stars , although a handful of bright red objects are seen to the east of it , suggesting the presence of an older generation of stars .
this is consistent with the conclusions of selman et al .
( 1999 ) , who found at least three generations of massive stars , with the youngest and bluest members more centrally concentrated .
we will return in more detail to the spatial distribution of these objects in section4 .
throughout this paper we will concentrate on the observations in the optical wavelength range , and in particular on the f555w , f656n and f814w bands ( hereafter referred to as @xmath8 , @xmath19 and @xmath9 , respectively ) , whereas the entire wfc3 dataset including the observations at ir wavelengths will be covered in more detail in a forthcoming paper ( beccari et al .
, in preparation ) .
as regards the data reduction , the long and short exposures were treated in a different way .
the deep datasets were combined using the iraf _ multidrizzle _ task , with the appropriate calibration files and geometric distortion coefficients for the uvis channel of the wfc3 . given the large number of deep frames available in the uvis filters and the careful dithering strategy adopted for these observations , we were able to derive one single `` drizzled '' image per band , free of cosmic rays and blemishes and with a rather sharp and well sampled psf , ideal for the detection of faint stars .
producing a deep , high - quality combined image per filter is particularly useful in a field like 30dor where stars are embedded in extensive nebulosity associated with the gas that makes the photometry and the determination of the background more difficult on individual images .
star detection and identification was done by running the automated iraf _ daofind _
routine separately on the @xmath8- and @xmath9-band drizzled frames , with the detection threshold set at @xmath20 above the local background level . in order for
objects to be classified as real stars , we required that they be detected in both filters .
the presence of numerous gas filaments and irregularities in the background can cause spurious detections , particularly in the h@xmath1 band , where considerable emission is present ( see figure[fig0 ] .
for this reason , we visually inspected the detected stars and excluded from the bona - fide list all the detections not clearly associated with point - like sources .
more specifically , we applied an unsharp - masking filter ( e.g. malin 1977 ) , with a smoothing radius of 15 pixel meant to subtract the low - frequency signal .
unsharp masking makes structures and filaments in the nebular gas easier to see than in the direct image . in this way we compared the structures and filaments with the positions of the detected stars , in order to avoid spurious detections ( see beccari et al .
2010 for an application of this method to the case of the massive galactic cluster ngc3603 observed with the same camera ) .
once the master list was defined in this way , we proceeded to measuring the magnitudes of each star with the iraf _ phot _ aperture photometry task .
since the central regions are rather crowded , we opted for a photometric aperture of 3 pixel radius and a background annulus enclosing the area between 4 and 7 pixel radius around the stars . following de marchi et al .
( 1993 ) , we took as background value the mode of the counts within the annulus , in such a way to minimise the contamination due to neighbouring objects .
this procedure was not applied directly on the combined drizzled frames , but rather on a subset of those .
more precisely , having at our disposal 16 deep exposures both in the @xmath8 and @xmath9 bands , we used multidrizzle to build for each filter four drizzled images , each comprising four individual exposures .
we then ran the iraf phot task on each of them , using as input coordinates the master list derived before . in this way
, we could derive four independent measurements of the magnitude for each object in both the @xmath8 and @xmath9 band and the standard deviation of the four measurements immediately provides the photometric uncertainty .
we followed the same procedure in all other visible bands , including h@xmath1 for which the four combined images comprise each two individual exposures .
the approach followed for the short exposures ( ranging from a @xmath21s to 30s in duration ) is different , owing to the limited number of exposures available in each filter . in that case
we performed the photometry on each individual flat - fielded image using psf fitting , having first built a proper psf model on each frame using isolated and well sampled stars , which was then fitted to each object in the images by means of the standard daophot / allstar routine .
the reason why we did not use psf - fitting photometry also for the deep exposures is that image combination with multidrizzle does not necessarily preserve the uniformity of the psf .
the two photometric catalogues obtained in this way for the long and short exposures were independently calibrated following the specific recipes developed by kalirai et al .
( 2009 ) for the uvis channel of the wfc3 .
briefly , all magnitudes were first rescaled to an aperture of @xmath22 and then transformed into the vegamag system using the synthetic zero points given by kalirai et al . ( 2009 ) .
the accuracy of the photometry in both catalogues is very good . the typical photometric uncertainty over the entire sample ranges from @xmath23mag at @xmath24 or @xmath25 to @xmath26mag at @xmath27 or @xmath28 ( note that in the following sections we will limit our study to objects with the smallest uncertainties , as indicated ) .
finally , the two photometric catalogues were combined and merged so as to avoid duplications .
we used stars in the magnitude range @xmath29 , where errors are small , to verify the integrity of the photometry by comparing the magnitudes derived for the same star in the long and short exposures .
we found the match to be excellent , with any differences in magnitude being comparable to the photometric uncertainties in the short exposures ( @xmath30mag ) .
the final combined catalogue contains 22291 stars with well defined fluxes in the @xmath8 and @xmath9 bands . of these
, 18142 objects are also detected in h@xmath1 .
the cmd of all objects with combined photometric uncertainty in the @xmath8 and @xmath9 bands not exceeding @xmath31mag is shown in figure[fig1 ] .
the cmd shown in figure[fig1 ] is the deepest so far obtained for this cluster . at @xmath32
it reaches about 3mag deeper than the cmds previously obtained with the wfpc2 camera on board the _ hst _ in these bands ( sirianni et al .
2000 ; hunter et al .
note that this is not only due to the improved sensitivity of the wfc3 camera , but also to the considerably shorter exposure times in the i band and the smaller field of view of the wfpc2 camera , which limited the observations to the most crowded regions of the cluster .
the striking difference between the cmd of figure[fig1 ] and that of sirianni et al .
( 2000 ) is the prominent lower ms occupying the region fainter than @xmath33 and bluer than @xmath32 and revealing the presence in this field of a considerably older stellar population than the @xmath34myr old stars associated with r136 ( e.g. de marchi et al . 1993 ; hunter et al .
1995a ; massey & hunter 1998 ) .
we will discuss the age distribution of the stars in this field in sections5 and 6 .
the other interesting feature seen in figure[fig1 ] is the large number of objects in the range @xmath35 that appear brighter than the lower ms .
their location is consistent with that of the pms population studied by sirianni et al .
( 2000 ; see also andersen et al .
2009 ) , but their number is much larger than those detected before .
we discuss the properties of these stars in section5 , but we first need to address an important issue that affects the cmd , namely interstellar extinction .
the width of the upper ms seen in figure[fig1 ] , at @xmath36 , is much larger than the photometric uncertainty at those magnitudes ( indicated by the horizontal bars on the left - hand side of the figure ) , revealing that there is a wide spread of extinction in this field .
actually , there is general consensus that differential reddening is present in these regions . from the analysis of high - resolution narrow - band images , hunter et al .
( 1995b ) concluded that the nebulosity in this area is not uniform even on scales of order @xmath37 .
however , in a subsequent study of the stellar populations in this region hunter et al .
( 1996 ) opted for a uniform reddening correction across the whole field , since they concluded that uncertainties ( e.g. red leak ) in the response of the wfpc2 filters did not justify a star - by - star correction .
also sirianni et al .
( 2000 ) used the same @xmath38 value for all stars , although this may have lead to larger uncertainties in their mass determination of individual pms stars , as selman et al .
( 1999 ) and more recently andersen et al .
( 2009 ) have pointed out . in order to overcome this problem , in section4.1
we will use a robust method based on the determination of the individual reddening corrections for a large group of bright stars ( romaniello 1998 ) , which are then used to derive a correction also for fainter neighbouring objects .
the young population in this field , responsible for the upper ms in figure[fig1 ] , has an age of @xmath34myr ( e.g. de marchi et al .
1993 ; hunter et al .
1995a ) or perhaps younger ( @xmath10myr ) if wolf rayet stars are still burning hydrogen ( massey & hunter 1998 ; de koter , heap & hubeny 1998 ; walborn et al .
one can therefore determine the reddening towards each of those stars by measuring their colour and magnitude displacements with respect to a theoretical isochrone for the appropriate age , metallicity and distance .
as for the metallicity , throughout this paper we will assume @xmath39 or about one third @xmath40 , since this is a typical value for the lmc ( e.g. hill , andrievsky & spite 1995 ; geha et al . 1998 ) .
as regards the distance to 30dor , we will use @xmath41kpc since this is the geometrical distance to the neighbouring sn1987 object ( panagia et al .
1991 , later updated in panagia 1999).the resulting distance modulus of @xmath42 is fully consistent with that adopted by ( walborn & blades 1997 ) .
we show in figure[fig2 ] as a solid line ( in blue in the online version ) an isochrone for an age of 3myr and a metallicity of @xmath39 , from marigo et al .
( 2008 ) , and the adopted distance modulus of @xmath42 .
having assumed a younger age , i.e. 1 or 2myr , would have only affected the shape of the isochrone in the top @xmath43mag range , which we will not consider in this analysis .
the isochrone is also already reddened by the amount corresponding to the intervening absorption along the line of sight due to the milky way , which fitzpatrick & savage ( 1984 ) quantified in @xmath44 , in turn corresponding to @xmath45 .
the other curves on the left - hand side of the figure ( in blue in the online version ) represent the same isochrone reddened , for purpose of illustration , by an additional amount of extinction , with the long - dashed , short - dashed and dot - dashed lines corresponding respectively to @xmath46 , @xmath47 and @xmath48 . as for the extinction law , for the portion in addition to the @xmath44 due to the milky way we have used the law determined spectroscopically by scuderi et al .
( 1996 ) for star2 in the field of sn1987a , which we take as representative of the 30dor region as well .
the arrow in the figure shows the corresponding reddening vector for @xmath49 . the long - dashed line , corresponding to @xmath46 , forms a natural left envelope to the upper ms at @xmath50 and leaves very few objects in the range @xmath51 at those magnitudes .
this is expected because of the considerable amount of foreground extinction known to be present in the lmc towards 30dor : although the @xmath52 value that we find is somewhat larger than the @xmath53 reported by fitzpatrick ( 1985 ) , the latter value has been determined for the brightest stars at the very centre of 30dor that generally lie in evacuated cavities of hii gas of their own making and should therefore have a lower extinction .
note that , at magnitudes fainter than @xmath33 , where the photometric uncertainty is still rather small ( see bars in figure[fig1 ] ) , there is a large number of objects between the @xmath46 and @xmath44 solid curves .
in fact , the latter offers a much better fit to the lower ms , suggesting that these stars belong to an older foreground population . at least some of them could be members of the much older lmc field , which is known to have ages in excess of @xmath43gyr ( e.g. panagia et al .
2000 ; harris & zaritsky 2004 ) .
interestingly , there are also a few objects that appear bluer than the @xmath44 curve
. although they are in general consistent with our photometric uncertainties , some of them could actually be very young stars with a circumstellar disc seen at high inclination ( @xmath54 ) .
objects of this type would appear somewhat bluer than their photospheric colour due to light scattering on the circumstellar disc , as well as several magnitudes fainter than their photospheric brightness due to extinction along the line of sight caused by an almost edge - on disc .
however , according to the models of robitaille et al .
( 2006 ) for the spectral energy distribution of young stars seen at various viewing angles , objects of this type can only account for a few percent of the total young population , so the vast majority of stars near the @xmath44 solid curve must belong to the field in the lmc foreground .    the presence of a foreground population of older objects becomes even more obvious when one looks at the other set of curves on the right - hand side of the figure ( shown in orange in the online version ) , corresponding to isochrones for an age of 2gyr and the same metallicity , distance , and reddening values used for the 3myr isochrones , also from the models of marigo et al .
the 2gyr isochrones are shown up to a mass of @xmath55m@xmath6 , corresponding to the approximate location of the red giant clump ( rc ) , a characteristic feature of the cmd of low - mass , metal - rich stars experiencing their he - burning phase .
it represents the counterpart at high metallicity of the horizontal branch seen in cmds of metal- poor globular clusters .
the rc usually appears as a very tight , well defined feature in the cmd of stellar populations that could in principle be used as a distance indicator ( e.g. stanek & garnavich 1998 ; paczynski & stanek 1998 ; cole 1998 ; girardi et al .
1998 ; girardi & salaris 2001 ) .
when there is patchy absorption , like in our case , and the amount of extinction varies from place to place , the rc becomes in practice an elongated ellipse running parallel to the direction of the reddening vector .
this is precisely the nature of the apparent sequence of points running diagonally between @xmath56 and @xmath57 .
this feature can also be used to derive the extinction law ( de marchi , panagia & romaniello 2007 ; panagia , de marchi & romaniello 2008 ) .
the implications of this finding for the properties of the absorbing material in this and other lmc fields will be addressed in detail elsewhere ( de marchi & panagia 2011 , in preparation ) , but it is immediately clear from figure[fig2 ] that the young and old stars present in this field have rather different extinction ranges . while stars in the rc span the entire range @xmath58 , there are in practice no objects in the upper ms with @xmath59 . therefore
, the reddening correction that we can derive using the objects in the upper ms will only be valid for those stars and , in a statistical sense , for low - mass stars of similar age if they have a similar spatial distribution , but it will not be appropriate for older objects , for which we will have to resort to an average correction .
the procedure that we followed for reddening correction is explained in detail in the following subsection .
we have used all the upper ms objects inside the box shown in figure[fig2 ] as reference stars for our determination of the reddening distribution across the field , with a total of about 800 stars selected in this way . assuming that they are all ms stars , we have measured their displacement from the solid isochrone along the reddening vector , whose direction is given by the arrow and the tilted sides of the box (
as mentioned above , we adopted the extinction law of scuderi et al .
we have explicitly excluded the brightest objects where an age difference of 1 or 2myr could have an effect , albeit small , on the shape of the isochrone .
the reddening distribution ranges from @xmath60 or @xmath61 to @xmath62 or @xmath63 , corresponding respectively to the 17% and 83% distribution limits , with a median value of @xmath64 or @xmath49 .
the difference between the 17% and 83% limits indicates a @xmath65 spread of @xmath66mag in @xmath67 , indicating that reddening varies considerably in this field .
the reddening values towards each of the selected ms stars were then used to derive a reddening correction for the objects in their vicinities .
this procedure implicitly assumes that stars projected within a relatively small region of the field have similar extinction values , even though no information is available on the depth of their distribution nor on the uniformity of the absorbing material . in order to verify the validity of this assumption , we first applied the reddening correction to the upper ms stars themselves .
each star was corrected using the average of the extinction values of the nearest upper ms stars , excluding the star itself ( no distance weighing was used , since we have no information on the depth of the distribution ) . after experimenting with a different number of neighbours ( namely 1 , 5 , 10 , and 20 )
, we concluded that using the 5 closest neighbours produces the tightest upper ms .
this choice corresponds to using reference stars within a radius of less than @xmath68 , with a typical effective radius of @xmath69 .
after reddening correction , the @xmath70 spread in the upper ms drops to @xmath71mag and the upper ms itself appears considerably tighter in the cmd than in the raw data ( see figure[fig6 ] ) .
this result convinced us that the method , albeit statistical in nature , can be applied to other objects in the field , provided that they have an age and spatial distribution consistent with that of the upper ms stars .
the cmd and isochrones in figure[fig2 ] help us to identify stars of similar age .
this is the case of ms stars roughly brighter than @xmath72 and bluer than @xmath73 and of young pms objects roughly redder and fainter than the dot - dashed 2gyr isochrone ( these objects are shown as dark dots in the figure ) .
however , the cmd alone does not allow us to assign an age to objects falling in the region between the solid and dot - dashed 2gyr isochrones ( shown as light dots ) , since they could be pms stars of various ages , low - mass ms stars , or red giants .
therefore , in order to identify bona - fide young objects in this region we have to resort to a more reliable youth indicator , namely a strong h@xmath1 excess emission .
this feature is characteristic of low - mass star formation and is attributed to the accretion process , whereby gravitational energy , released by infalling matter , ionises and excites the surrounding gas ( e.g. knigl 1991 ; shu et al .
1994 ) .
a detailed account on how to reliably detect these objects using multi - colour photometry is given in de marchi et al .
( 2010 ) and de marchi et al .
the method , already employed also by beccari et al .
( 2010 ) to search for objects with h@xmath1 excess in ngc3603 , will be used in section5 for an accurate selection of bona - fide pms stars after reddening correction .
briefly , the method uses the median @xmath74 dereddened colour of stars with small ( @xmath75mag ) photometric uncertainties in the three bands ( @xmath76 and @xmath19 ) to define a reference template that is used to identify objects with excess h@xmath1 emission .
however , as de marchi et al .
( 2010 ) have shown , even before reddening correction the colour  colour diagram @xmath77 , @xmath78 provides a robust identification of stars with h@xmath1 excess , since in these bands the reddening vector runs almost parallel to the median photospheric colours of normal stars .
therefore , following that procedure , we looked initially for stars whose @xmath77 colour exceeds by at least three times the photometric uncertainty the @xmath77 colour of normal stars ( i.e. those with no excess ) with the same @xmath78 colour .
a more stringent selection will be applied after reddening correction ( see section5 ) .
the pms candidates selected in this way , a total of 2374 , are shown as thick dots in figure[fig3 ] , where the dashed line represents the median photospheric @xmath77 colour as a function of @xmath78 . for comparison
, we also show as a thin dotted line the corresponding colours derived from the model atmospheres of bessell , castelli & pletz ( 1998 ) for stars with effective temperatures in the range 3500 k @xmath79 40000k , surface gravity @xmath80 and a metallicity index @xmath81=-0.5 $ ] , as appropriate for the lmc ( dufour 1984 ) .
the theoretical colours were calculated for the specific filters of the wfc3 camera and were reddened by @xmath82 , following the extinction law of scuderi et al .
( 1996 ) mentioned above for the purpose of comparison with the data .
note that , due to the still preliminary calibration of the wfc3 photometric system in the h@xmath1 band , we had to apply a small correction to the @xmath77 colour of @xmath83 mag to force the @xmath77 colour of vega to be 0 at @xmath84 , as is appropriate for the vegamag photometric system .
the objects with h@xmath1 excess emission are also shown as thick symbols in the cmd of figure[fig4 ] .
besides the 3myr and 2gyr isochrones , we show there also the isochrones for pms stars of age 5myr ( upper thin line ) and 10myr ( lower thin line ) from the models of deglinnocenti et al .
( 2008 ) for the assumed lmc metallicity ( @xmath39 ) and having taken @xmath46 for the reddening .
these isochrones allow us to distinguish in an approximate way between younger and older pms stars ( a more accurate age determination requires the reddening correction and will be possible only at the end of this process ) and to study their spatial distribution . using the 5myr pms isochrone as a reference ,
we have marked as thick dots and crosses respectively the objects above and below that line , and hence approximately younger and older than 5myr .    in figure[fig5 ] we show the radial distribution of these low - mass stars with h@xmath1 excess , drawn from the nominal centre of r136 .
the short - dashed line corresponds to stars younger than @xmath85myr , the long - dashed line to stars older than 10myr , while the distribution of upper ms stars is shown by the solid line , rescaled by a factor of @xmath86 to make the comparison easier .
younger and older pms stars have markedly different distributions : older stars have a rather uniform distribution and are almost absent inside a radius of @xmath87 , while younger objects show a density gradient that follows remarkably well that of upper ms stars down to @xmath37 , where photometric incompleteness starts to dominate as indicated by the rapidly increasing photometric uncertainties in these regions .
the choice of which approximate age ( 5myr or 10myr ) to adopt for separating `` younger '' and `` older '' stars is somewhat arbitrary , but the dotted line in figure[fig5 ] , corresponding to stars with age in the range @xmath88myr , suggests that objects with h@xmath1 excess older than @xmath89myr have a radial distribution remarkably different from that of upper ms stars .
this is perhaps not surprising , considering that 5myr is about twice the age of the massive members of r136 ( @xmath90myr ; de marchi et al .
1993 ; hunter et al . 1995a ; massey & hunter 1998 ; de koter et al . 1998
; walborn et al .
1999 ) .
we will discuss these differences in the radial distribution in more detail in section6 .
for now , figure[fig5 ] confirms that it is appropriate to use the reddening values derived from upper ms stars also for similarly young low - mass objects , but it suggests that this type of correction would not be suitable for older stars outside of the inner @xmath91 . as figure[fig5 ] implies , these objects are uniformly distributed over the field and are thus consistent with being in the foreground of r136 . in order to obtain at least an estimate of their reddening , we can use the stars with h@xmath1 excess that appear bluer than the isochrone for @xmath46 , i.e. the long dashed line in figures[fig2 ] and [ fig4 ] .    as mentioned earlier
, @xmath46 is the lower limit to the reddening towards the upper ms stars of r136 , so objects bluer than the long - dashed line in those figures are most likely in the foreground .
a total of 144 stars satisfy this condition and , assuming that they are in the late stages of pms evolution , i.e. close to the ms , we derive for them a median reddening value of @xmath82 .
if these objects were in a considerably earlier stage of their pms evolution , this procedure would give us an upper limit to their reddening .
on the other hand , the fact that they are uniformly distributed over the field indicates that they can not be as young as the other much more concentrated objects with h@xmath1 excess ( see figure[fig5 ] ) .
we therefore take @xmath82 as a representative value for the reddening towards all these stars and in general for all objects with h@xmath1 excess shown as crosses in figure[fig4 ] .
we can now correct for extinction the magnitudes of all objects shown in figure[fig4 ] in the following way .
all stars located within @xmath92 of the centre of r136 are corrected individually using the average extinction value of the nearest five upper ms stars . as regards objects outside @xmath92 , those with an age lower than @xmath85myr ( i.e. thick and thin dots in figure[fig4 ] ) are corrected in the same way , whereas older stars with excess ( crosses in figure[fig4 ] , red in the online version ) are corrected using the average @xmath82 value .
no correction is attempted for other objects , since we can not determine their ages on the basis of the cmd alone , and they are not considered further in this work .
note that when the magnitude of an object is corrected for extinction using its nearest neighbours , it is assumed that the placement of the object with respect to the absorbing material is similar to that of the reference stars in its vicinities .
this process is accompanied by some statistical fluctuations , since as we have seen there is an intrinsic spread of @xmath71mag in the @xmath67 distribution of our upper ms reference stars .
however , the uncertainties that this assumption can introduce in the final magnitudes of the stars are relatively small . the typical standard deviation of the @xmath67 values of the five nearest neighbours is @xmath93mag and , more generally , over the entire field these statistical uncertainties are not expected to exceed the reddening spread measured for the upper ms reference stars before reddening correction , i.e. @xmath66mag in @xmath67 or @xmath94mag in @xmath38 .
the cmd corrected for extinction in this way is shown in figure[fig6 ] .
it must be understood that it does not cover all objects that were observed , but only those for which we could determine a reliable correction for reddening .
the symbol types are the same as used in figure[fig4 ] , the solid line represents the theoretical isochrone from marigo et al .
( 2008 ) for a ms age of 3myr , a metallicity of @xmath39 and a distance modulus of @xmath42 , whereas the dashed line corresponds to the isochrone of deglinnocenti et al .
( 2008 ) for a pms age of 5myr and the same metallicity and distance modulus .
the pms isochrone of deglinnocenti et al .
( 2008 ) has been translated from the theoretical plane ( @xmath95 , @xmath96 ) to the plane of the observations using the model atmospheres of bessell et al .
( 1998 ) folded through the specific wfc3 filters employed here .
the pms isochrone meets the ms at @xmath97 and it is reassuring that the two curves are in excellent agreement at brighter magnitudes , since they are based on independent models and different input physics .
the cmd corrected for extinction in figure[fig6 ] reveals that at dereddened @xmath98 magnitudes fainter than @xmath99 there is a conspicuous number of objects that are above ( i.e. brighter than ) the ms and , as such , they are consistent with a population of pms stars .
these objects extend down to @xmath100 , or more than three magnitudes deeper than the limit reached by sirianni et al .
( 2000 ) and andersen et al .
( 2009 ) , who detected and studied the brightest members of this population .    as already mentioned in section4 , having corrected the photometry for extinction allows us to better clarify and confirm the pms nature of these stars by looking at the specific signatures of the mass accretion process that is supposed to accompany this evolutionary phase .
the accretion process is believed to be the cause of the strong excess emission observed in these objects ( e.g. knigl 1991 ; shu et al .
for example , the presence of a strong h@xmath1 emission line ( with equivalent width in excess of @xmath101 ) in young stellar objects is normally interpreted as the signature of ongoing accretion ( e.g. feigelson & montmerle 1999 ; white & basri 2003 ) .
note that , although it is customary to use negative equivalent widths for emission lines , in this work we will use @xmath102 to refer to the equivalent width of the h@xmath1 emission line , defined as a positive quantity .
in section4 we used this method for a preliminary coarse identification of stars with h@xmath1 excess emission , but here we employ very conservative criteria to refine the selection of bona - fide pms stars . in particular , we now limit the search only to stars with a combined photometric uncertainty in the @xmath76 and @xmath19 bands of less than @xmath103mag and fainter than @xmath104 ( corresponding to @xmath105l@xmath106 ) and we look for objects with dereddened @xmath107 colour excess at least four times larger than the photometric uncertainty in those bands ( equivalent to a @xmath108 detection limit ) .    in section4 we used as a reference template the median photospheric colours of stars with small photometric uncertainty , the majority of which are older ms stars and have no h@xmath1 excess .
this approach can not be followed here since many of the objects for which we could compute a reddening correction have h@xmath1 excess . as a reference template
we use instead the theoretical photospheric colours derived from model atmospheres . as figure[fig3 ] shows , the models of bessell et al .
( 1998 ) for surface gravity @xmath109 and metallicity @xmath81 = -0.5 $ ] are in excellent agreement with the observed colours .
this was also shown to be the case for stars in the field of sn1987a ( de marchi et al .
2010 ) and in ngc346 ( de marchi et al .
2011 ) .    besides selecting stars with excess emission above the @xmath108 level
, we also imposed constraints on the equivalent width of the h@xmath1 emission , @xmath102 . following the procedure outlined in de marchi et al .
( 2010 ) to derive @xmath102 from the observed colours , we further restricted the selection to objects with @xmath110 in order to avoid contamination by stars with significant chromospheric activity
. furthermore , since at temperatures @xmath111k or colours @xmath112 the sample could be contaminated by be stars that are evolving off the ms , for those objects we set a more stringent condition , namely @xmath113 . in a survey of about 100 stars of type
be in the galaxy ( cot & waters 1987 ) , only one star is found with @xmath113  and the largest majority have values in the range from @xmath114  to @xmath115 .
a total of 1159 objects satisfy these conservative detection conditions : we will take them hereafter as bona - fide pms stars .
obviously , less stringent limits would result in a larger number of pms stars , albeit with lower significance , but it must be understood that , even without changing the selection criteria , the actual number of pms stars in this field must be considerably larger . this is because our detection and classification relies on the presence of h@xmath1 excess emission at the time of the observations , but very young populations are known to show large variations in their h@xmath1 emission over hours or days ( e.g. fernandez et al .
1995 ; smith et al . 1999 ; alencar et al .
2001 ) , with only about one third of them at any given time being active h@xmath1 emitters with @xmath116 ( panagia et al .
2000 ; de marchi et al . 2010
; de marchi et al . 2011 ) .
therefore , the 1159 stars identified as bona - fide pms objects are those that were undergoing active mass accretion and therefore showed h@xmath1 excess emission at the time of these observations . on the other hand , since the duty cycles of any two stars are completely independent of one another , the sample of bona - fide pms stars identified in this way is representative of the entire pms population in this field .
hereafter , unless otherwise indicated , when talking about pms stars we will refer exclusively to the 1159 bona - fide pms objects identified through their h@xmath1 excess emission .
we show in figure[fig7 ] their locations in the cmd ( thick dots ) .
the majority occupy the region above the ms where young pms objects are expected , but there are many located at or near the ms , suggesting an older age .
a similar situation has already been noticed in other star forming regions , namely the sn1987a field in the lmc ( de marchi et al .
2010 ) , ngc3603 in the milky way ( beccari et al .
2010 ) , and ngc346 in the small magellanic cloud ( de marchi et al .
2011 ) .
the skeptical reader might be concerned that the h@xmath1 excess emission of many of the objects near the ms could in reality be due to nebular contamination of the gas present in the field . however , as mentioned in section3 , we specifically excluded from our photometry all objects for which the photometric aperture would be contaminated by gas filaments . if the filament completely fills both the core aperture used to measure the flux of the star and the surrounding background annulus , it will not alter the photometry of the star , since it will be equivalent to a uniform background .
obviously , if the background annulus does not properly sample the contamination , for instance because the filament only covers part of it , the emission present in the gas may be erroneously associated to the object ( for an illustration of this effect see beccari et al .
precisely for this reason , as discussed in section3 , we have removed from the initial star lists all objects in the vicinity of a filament .
therefore , we are confident that the h@xmath1 excess emission of our bona - fide pms stars is an intrinsic feature of those objects and we must interpret their distribution across the cmd as a sign of an age spread .
more quantitative information on the age of these objects can be obtained from the comparison with pms isochrones , shown in figure[fig8 ] , where the effective temperatures ( @xmath117 ) and luminosities ( @xmath118 ) of these stars are arranged in the hr diagram .
the value of @xmath117 was derived for each star directly from the dereddened @xmath78 colour , since this is an excellent index for temperature determinations in the range @xmath119k ( e.g. von braun et al .
1998 ; bessell et al .
1998 ) . for the conversion we made use of the models of bessell et al .
( 1998 ) already discussed in section4 .
as for the bolometric luminosity , it was obtained from the dereddened @xmath8 magnitude and the same models , having assumed a distance modulus of @xmath42 .
the solid lines in the figure corresponds to the pms evolutionary track from the pisa group ( deglinnocenti et al .
2008 ; tognelli et al .
2011 ) for metallicity @xmath39 .
the tracks correspond to stellar masses from @xmath120m@xmath6to @xmath121m@xmath6and reveal the presence of bona - fide pms stars over the entire mass range .
the isochrones ( dashed lines ) correspond , from right to left , to ages ranging from @xmath122myr through to 32myr with a constant age step of a factor of 2 , selected in such a way that our photometric uncertainties are slightly smaller than the typical separation between the isochrones .
pms objects are seen across the entire age range .
the comparison of the hr diagram with pms evolutionary tracks and isochrones allows us to quantify the age spread of the pms stars that was already evident in the cmd .
individual stellar masses and ages were derived via interpolation over a finer grid than the one shown in the figure .
we followed the interpolation procedure developed by romaniello ( 1998 ) , which does not make assumptions on the properties of the population , such as the functional form of the imf . on the basis of the measurement errors
, this procedure provides the probability distribution for each individual star to have a given value of the mass and age ( the method is conceptually identical to the one presented recently by da rio et al . 2010 ) .    an important aspect to consider
is the uncertainty associated with the age determination of our pms stars from their location in the hr diagram . besides photometric errors in the @xmath8 and @xmath9-band data ( which are always less than @xmath31mag )
, the major sources of uncertainty are the effects of reddening and the comparison with stellar models . as regards the latter , de marchi et al .
( 2010 ) show that , excluding possible inaccuracies in the models input physics and interpolation errors , the largest source of uncertainty comes from the use of models that do not properly describe the stellar population under study , e.g. because of the wrong metallicity , and from differences between models of various authors .
they conclude that systematic uncertainties of order 20% are to be expected for the mass and possibly higher for the age .    as for reddening
, dust extinction systematically displaces stars in the hr diagram towards lower luminosities and effective temperatures . therefore ,
if the photometry of a pms star is not corrected for reddening , its age will typically be underestimated .
a study currently being conducted ( de marchi & panagia , in prep . ) shows that for reddening @xmath123 there is a rather linear relationship of the type @xmath124 , where @xmath125 and @xmath126 are the ages before and after reddening correction .
for example , as discussed in section4 , the spread of the upper ms in figure[fig1 ] corresponds to a reddening spread of @xmath127mag in @xmath67 .
therefore , if we had used the average reddening value of @xmath49 to correct the magnitudes of all young stars , instead of using the values of the nearest neighbours , we would have introduced an age uncertainty of a factor of @xmath55 .
similarly , if reddening correction for older stars with excess ( i.e. those indicated by crosses in figure[fig4 ] ) had been done using the average @xmath38 value of the nearest upper ms neighbours instead of @xmath82 , the derived ages would be typically twice as old . as a further example , the double arrow in figure[fig8 ] indicates the amount and direction of the displacement in the hr diagram caused by the typical uncertainty of @xmath128mag in @xmath67 associated with our reddening correction ( see section4.1 ) . in the specific case shown ,
the corresponding uncertainty on the age is a factor of @xmath55 while the uncertainty on the mass is less than 15% .
this analysis suggests that the typical uncertainty on the ages that we derive is less than a factor of 2 .
for this reason , the age distribution of pms stars is shown in figure[fig9 ] using as bin size a constant logarithmic step of a factor of 2 , which does not suffer from the relative age uncertainties .
the solid line gives the number of stars inside each age bin as a function of time , whereas the dashed line provides an apparent value of the star formation rate in units of stars per myr ( i.e. the number of stars in each bin divided by the size of the bin ) .
the dashed line seems to suggest a rather uniform star formation rate over the past 4myr , of order @xmath129m@xmath6myr@xmath130 ( the mass of our pms stars younger than 4myr ranges from @xmath131m@xmath6to @xmath132m@xmath6 , respectively the 17 and 83 percentiles , with an average of 2m@xmath6 ) , and a progressive drop at earlier times .
obviously , this is necessarily a lower limit to the star formation rate , since our measurements only account for pms stars in the range @xmath133m@xmath6that had h@xmath1 excess emission at the time of the observations .
furthermore , since the fraction of pms stars in an active state of accretion is expected to decrease with age ( see de marchi , panagia & sabbi 2011 ) , the apparent drop in the star formation beyond 2myr can at least in part be explained .    keeping in mind that we can only provide a lower limit to the star formation rate , near the peak of the distribution , at @xmath43myr
, this limit corresponds to @xmath134m@xmath6myr@xmath130 ( the average mass of those objects is @xmath135m@xmath6 ) , while at @xmath136myr it drops by over an order of magnitude to @xmath101m@xmath6myr@xmath130 ( average mass @xmath137m@xmath6 ) .
a more detailed study of the star formation rate of 30 dor will be the topic of a future paper , but the histograms in figure[fig9 ] already show quite convincingly that star formation has been present in this region for an extended period of time , exceeding 30myr .      from a spectroscopic study of massive stars ( @xmath138m@xmath6 ) within a radius of @xmath139 of r136 , selman et al .
( 1999 ) concluded that recent star formation in this region shows at least three distinct bursts , respectively @xmath140 and @xmath55myr ago , with increasing strength .
they also found that the younger generation is considerably more concentrated than the others and is responsible for most of the star - formation within 6pc or @xmath141 of the cluster centre .
it is therefore useful to compare this picture with the spatial distribution of our pms stars .
although we do not have the age resolution needed to establish unambiguously whether the formation has been continuous during the past 30myr or whether repeated bursts have occurred , we can look at the current age distribution for clues as to the relationships between stars of different ages .
the age distribution of these pms stars is characterised by a median of @xmath142myr , with about @xmath143 of the stars being younger than @xmath2myr and an equal number older than @xmath4myr .
hereafter we will refer to these two groups , comprising 415 objects each , as `` younger '' and `` older '' pms stars , respectively .
note that , with a typical age uncertainty not exceeding a factor of 2 , these two groups are well separated in age , since the former has a median age of @xmath10myr and the latter @xmath144myr .
interestingly , also their projected spatial density distributions are very different .
they are shown in figure[fig10 ] by contour lines of different colour , with logarithmic scaling , plotted over an image of the field in the h@xmath1 band .
the contour plots have been derived after smoothing the distribution with a gaussian beam with size @xmath145 or @xmath43pc , as indicated by the circle in the lower left corner of the figure .
the lowest contour level corresponds to a local density of pms stars twice as high as the average pms stars density over the entire field .
the step between levels corresponds to a factor of @xmath146 .
there is a pronounced difference in the spatial distribution of younger and older pms stars .
the former are concentrated near the centre of the r136 cluster , where they form a density plateau , while older pms stars occupy preferentially the region to the east of r136 and are remarkably absent from the cluster centre .
since these stars have a rather uniform distribution across the field , we should expect to find them also near the location of the young r136 cluster , contrary to what we observe .
furthermore , with a median age of @xmath144myr , these stars are considerably older than the @xmath10myr old members of r136 , so we should not expect any spatial relationship ( either a correlation or anti - correlation ) between the distributions of the two types of objects , unless something else related to the younger population prevents us from detecting older pms stars near the cluster centre .
photometric incompleteness due to the young massive stars could seem an obvious culprit , since the younger pms stars are on average brighter than the older objects .
however , our analysis shows that the difference in the distribution of younger and older pms stars remains the same as shown in figure[fig10 ] even when we restrict the choice of objects to the same magnitude range ( e.g. @xmath147 ) . therefore
, the apparent paucity of older pms stars near the centre is not caused by selection effects in their detection .
a more likely explanation is the erosion of circumstellar discs via photo - evaporation caused by the far - ultraviolet ( fuv ) radiation of the young massive members of r136 .
if the discs of pms stars are disrupted , the accretion process dwindles and these objects would no longer be detectable through their h@xmath1 excess emission . to test this hypothesis ,
we have looked at how older pms stars of various h@xmath1 luminosity @xmath148 are distributed in the field . in figure[fig11 ] we show , as an example , all pms stars with ages in the range @xmath149myr ( 115 objects in total ) using different symbols according to their h@xmath1 luminosity .
dots correspond to stars with @xmath150l@xmath151 , squares to @xmath152l@xmath151 , and triangles to intermediate values .
these values have been selected in such a way to have an equal number of objects in each group . in the background
we show as a reference a negative image of the field in the h@xmath1 band .
the distribution of the three classes of objects is not random : stars with higher @xmath148 are typically closer to the centre of r136 and located in areas where the background has a higher emission , whereas stars with lower @xmath148 are further away and preferentially in regions where the background emission is low .
although systematic differences in the mass of the three classes of objects could produce the observed distribution of @xmath148 , this can be excluded since the median mass of the three groups of stars is the same at @xmath153m@xmath6 .
we can also exclude that this is an artefact of our photometry since the background is subtracted locally from each star ( see section3 ) and therefore the diffuse h@xmath1 nebular emission does not affect the intrinsic h@xmath1 brightness of individual objects .
however , a fraction of the h@xmath1 luminosity of the stars closest to the centre of r136 could in principle be due to fluorescence of their circumstellar discs caused by the very strong ionising radiation of the massive stars .
it is , therefore , important to quantify this effect .
according to massey & hunter ( 1998 ) , adopting the temperature calibration of vacca , garmany & shull ( 1996 ) , the total bolometric luminosity of the @xmath3 brightest and most massive stars in r136 is of order @xmath154l@xmath106 and most of them are of spectral type o4 or earlier .
based on this information and following panagia ( 1973 ) , we have determined at which distance from the centre of r136 the effects of the ionising radiation impinging face - on on a circumstellar disc of radius 100au would contribute significantly to the h@xmath1 luminosity of its star .
the median @xmath148 value of the pms stars shown in figure[fig11 ] is @xmath155l@xmath106 , so we have considered three cases of contamination , at approximately the 10% , 50% and 100% level . the dotted , dot - dashed and solid circles shown in figure[fig11 ] correspond respectively to the distances at which the contribution to the h@xmath1 luminosity due to fluorescence is @xmath156 and @xmath157l@xmath106 ( note that these distances are necessarily upper limits since they are projected on the plane of the sky and we have no information on the distribution of these objects along the line of sight ) .
it is immediately clear that the effects of fluorescence can only be important for a handful of stars , leaving the largest majority of old pms objects in this field unaffected .
the origin of the apparent correlation between their positions and h@xmath1 luminosities has , therefore , to be found elsewhere .
a natural explanation would be that stars in regions of low emission , where the density of the gas is lower , such as those to the east of r136 , are much less shielded from the radiation of massive stars than objects in denser regions , such as those to the north and west of r136 .
therefore , their discs have been exposed to photo - evaporation for a longer time and their h@xmath1 luminosity is today further reduced .    in support of this hypothesis ,
the photometry of all stars enclosed in the low - background region resembling a `` christmas tree '' and centered at ( @xmath158 ) , which is better seen in figure[fig0 ] , is consistent with a very low extinction value , close to the @xmath44 indicated by fitzpatrick & savage ( 1984 ) as the contribution of the milky way to the intervening absorption along the line of sight .
furthermore , this region is characterised by very low emission at ir ( e.g. indebetouw et al . 2009 ) and radio wavelengths ( e.g. johansson et al . 1998 ; rubio et al .
1998 ) , confirming that it is devoid of both dust and gas .
the region also corresponds with a high - velocity structure seen by chu & kennicutt ( 1994 ) in their detailed h@xmath1 echelle study of the kinematics of 30dor .
these authors find a bright redshifted knot @xmath159 east of r136 , with a velocity offset of @xmath160 km s@xmath130 and very low emission .
the knot is part of a fast expanding shell to the east of the cluster that correlates well with a region of diffuse x - ray emission seen in images obtained with the _ einstein _ satellite .
these authors conclude that both the high - velocity features and the x - ray emission are related to high - velocity shocks from fast stellar winds and supernovae .
this picture is supported by the later study of 30dor by townsley et al .
( 2006 ) with the _ chandra x - ray observatory_. their observations reveal rather strong diffuse emission in the area of the `` christmas tree '' ( their region 6 , coincident with the east fast shell of chu & kennicutt 1994 ) with an absorption - corrected x - ray surface brightness of @xmath161erg s@xmath130 pc@xmath162 in the @xmath163kev band .
a spectral fit to this emission shows a soft thermal plasma with @xmath164kev ( townsley et al .
therefore , in such an evacuated region ,    the discs of pms stars must be much more exposed to photo - evaporation .
conversely , as discussed in section4 , extinction towards the young massive stars at the centre of r136 is both higher and highly variable , indicating a much denser environment where gas and dust can intercept a significant fraction of the uv radiation from the massive central stars , thereby reducing the efficiency of the photo - evaporation process .
indeed , most of the pms stars with @xmath165l@xmath106 in figure[fig11 ] are in regions of significant @xmath166co@xmath167 emission , suggesting a high density of molecular gas , and the column density in hi in front of r136 is known to be very high ( @xmath168 cm@xmath162 ; de boer , koornneef & savage 1980 ) .    thus , although this picture is still necessarily qualitative , it provides a plausible explanation of the observed distribution of h@xmath1 luminosities of older pms stars across the field via the effects of photo - evaporation .
these effects have been discussed by de marchi et al .
( 2010 ) in the field of sn1987a . in that environment , located @xmath170 or 300 pc sw of 30dor ,
the pms stars around the hottest massive young stars are both less numerous and fainter in h@xmath1 emission , suggesting that their circumstellar discs have been eroded more efficiently by enhanced photo - evaporation .    unlike the relatively simple case of the sn1987a field where little gas is present , the interstellar medium in 30 dor is known to be much denser and highly fragmented , with considerable density variations across the field ( e.g. poglitsch et al .
1995 ) .    from the h@xmath1 echelle study mentioned above ,
chu & kennicutt ( 1994 ) conclude that the kinematics of 30dor is dominated by a large number of expanding structures , over scales of 2  20pc and with expansion velocities of 100  300 km s@xmath130 and kinetic energies of @xmath171erg .
these structures are probably associated with supernova remnants embedded in supershells produced by stellar winds and supernovae explosions .
these processes can release a large amount of kinetic energy into the surrounding environment , and as such they could in principle contribute to the dissipation of circumstellar discs , adding to the effects of photo - evaporation . to better understand this , it is useful to calculate the amount of kinetic energy possibly deposited on such a disc by a nearby supernova explosion .
we will assume hereafter a total energy release of @xmath172foe @xmath173erg per supernova explosion , a circumstellar radius @xmath174au around a pms star located at a distance of @xmath175pc from the supernova .
the maximum fraction of energy deposited on such a disc hit face - on by the blast can be written as @xmath176 and amounts to @xmath177erg for the specific parameters assumed here . the motion that this energy release can impart onto the material of the disc depends on the total mass and size of the disc and hence on its density , which in turn could depend e.g. on the radial location . in the simplified case of a homogeneous disc of mass @xmath178 that we express as a fraction @xmath179 of the total stellar mass @xmath180 , i.e.
@xmath181 , the velocity @xmath182 imparted by the ejecta to the disc material can be expressed as :    @xmath183^{-1 } \ ,      \bigg [ \frac{r}{100 \ , { \rm au } } \bigg ] \ ,      \left [ \phi\ , \frac{m}{{\rm m}_\odot } \right]^{-1/2 } \ , { \rm km~s}^{-1}.\ ] ]    de marchi et al .
( 2011 ) have shown that stars in the magellanic clouds accrete a substantial fraction of their total mass during the pms phase , particularly at low masses , so for an order - of - magnitude calculation it is not unreasonable to assume @xmath184@xmath185 .
we can see in this way that for a star with @xmath186m@xmath6at a distance of 10pc , the supernova ejecta hitting the disc face - on would impart a motion to its material with a velocity of @xmath43 km s@xmath130 .
this value is about a factor of 5 smaller than the escape velocity at the edge of the disc for the typical pms stars that we detect in the `` christmas tree '' ( median mass @xmath187m@xmath6 ) . for a disc seen under a smaller solid angle
, the effect would be proportionally smaller .
on the other hand , less massive discs around less massive stars located closer to the supernova would obviously experience a much stronger effect .
for instance for a star of mass @xmath188m@xmath6at a distance of 1pc from the supernova , the impulse received by the disc material would be sufficient to exceed the escape velocity at @xmath189au .
this could significanly affect or possibly even stop the accretion process , eventually limiting the final mass of the star .
on the other hand , with an average surface density of @xmath190pms stars per pc@xmath191 , only a very small fraction of the discs around objects in this field would be severely affected by any one supernova explosion .    in summary ,
it is clear that the interplay between shocks and photo - evaporation effects due to the massive members of r136 can have important consequences on disc erosion , potentially affecting the formation of low - mass stars in these regions .
however , without more specific kinematical information for individual objects , it remains    difficult to quantify the expected effects on the circumstellar discs of pms stars using photometry alone .
for a better understanding of this complex region and a more quantitative characterisation of the differences between the regions preferentially occupied by older pms stars with high and low @xmath148 , we may have to await high spectral and spatial - resolution spectroscopic observations to obtain a map of the amount and kinematics of the gas in the interstellar medium along the lines of sight towards these objects .
finally , we briefly address the implications of our findings for the current understanding of whether and how different generations of stars in 30 dor are causally connected .
walborn & blades ( 1987 ) and walborn ( 1991 ) reported three early o stars embedded in dense nebular knots to the north and west of r136 and suggested that these objects could belong to a younger generation , triggered by the energy released by r136 . using ground - based spectral classification of a large number of ob stars in the 30 dor nebula , walborn & blades ( 1997 ) studied the ages and spatial distribution of these objects , discovering an additional eight early o stars embedded in dense nebular knots . from a morphological analysis of these regions
, they concluded that 30 dor has a two - stage starburst structure in which the energetic activity of the @xmath90myr old r136 has triggered the formation of a younger ( @xmath192myr ) generation of stars at a typical projected distance of @xmath193pc , inside the dense molecular clouds that surround r136 to the northeast and west .
later ir observations with nicmos ( walborn et al . 1999 ; brandner et al . 2001 ) resolved many of these knots into many new ir sources , including multiple systems , clusters , and nebular structures , indicating active star formation .    from mainly morphological considerations in this and other galactic star forming regions , walborn et al .
( 1999 ; see also walborn 2002 ) suggested that this mode of two - stage starburst , with generations separated by @xmath10myr , is ubiquitous .
they concluded that in 30 dor the first generation of stars , still present , is now triggering the birth of the second , peripheral generation of very massive stars , which in a few million years will be the only one left , after the most massive stars in the core of r136 have disappeared .
these authors predict 30 dor to become a giant shell h ii region like n11 and ngc 604 , with its most massive stars around the periphery .
our discovery of a considerable spread in the ages of pms stars in this region , with the younger generation confined to the cluster s centre and older objects uniformly distributed across the field ( see figure[fig10 ] ) , raises some questions as to the validity of this picture at least as regards the formation of low - mass stars triggered by older ( massive ) objects in their vicinities .
for instance , if we compare the projected distribution of pms stars with ages of less than 2myr with that of objects in the age - range @xmath194myr we see that they have the same centre , with the former slightly more concentrated , rather than distributed in the periphery of the latter as suggested in walborn s ( 2002 ) triggered massive star formation scenario .
as already shown in figure[fig10 ] , the stars around the periphery of r136 are systematically older and not younger than the members of r136 itself
. the same conclusion had been reached by selman et al .
( 1999 ) from a spectroscopic study of the massive stars ( @xmath195m@xmath6 ) in these regions .
our findings suggest that if causal effects are present between different generations of stars , their interpretation can not rely only on the relative ages and projected distributions of the objects , but must consider proper kinematical and possibly dynamical information . while it is undeniable that we are seeing a case of sequential star formation , in the sense that there are separate generations of stars over time , no clear signs of causal effects are seen .
therefore , although the concept of triggered star formation may well be valid , its relevance in the case of 30 dor has yet to be demonstrated .
we have studied the properties of the stellar populations in the central regions of the 30dor cluster in the lmc using observations obtained with the wfc3 camera on board the hst in the @xmath8 , @xmath9 and @xmath19 bands .
the main results of this work can be summarised as follows .    1 .
our photometry delivers the deepest cmd so far obtained for this cluster , reaching down to @xmath196 or about 3 mag deeper than previous hst observations , and reveals for the first time a prominent lower ms fainter than @xmath33 .
we attribute this feature to a considerably older stellar population than the @xmath34myr old stars associated with r136 .
the cmd shows also a large number of objects in the range @xmath35 that appear brighter than the lower ms and that are consistent with , but extend much deeper than , the population of pms stars detected and studied by sirianni et al .
( 2000 ) and andersen et al .
the observations also reveal the presence of considerable differential extinction across the field .
we quantify the total extinction towards massive ms stars younger than @xmath34myr to be in the range @xmath197 ( respectively the 17% and 83% levels ) with a median value of @xmath49 .
we show that these stars can be used to derive a statistical reddening correction also for stars in their vicinities , if they have a similar age and spatial distribution .
however , this is not the case for older stars , which show a broader reddening distribution starting already at @xmath44 ( corresponding to the contribution of the milky way along the line of sight towards 30dor ) .
a proper correction for extinction , therefore , requires at least an approximate knowledge of the age of the stars , which is not possible for all objects in the cmd .
3 .   in order to distinguish relatively young objects from older stars
, we looked for stars with h@xmath1 excess emission , since this feature identifies candidate pms stars , and we studied their distribution both across the field and in the cmd .
younger pms stars follow the spatial distribution of massive stars and this allows us to correct them for reddening individually , using the extinction information derived for the massive objects . on the other hand ,
older pms objects are uniformly distributed over the field and appear to be systematically less reddened than younger pms stars .
we derive for them an average extinction correction of @xmath82 .
4 .   reddening correction allows us to proceed to a more accurate identification of bona - fide pm stars .
we searched for stars with @xmath77 colour exceeding that of normal stars at the @xmath108 level or more and with @xmath198   ( or @xmath199  for stars with @xmath200k ) , in order to avoid contamination by stars with chromospheric activity or by the winds of be stars .
a total of 1159 objects satisfy these stringent conditions and we take them as bona - fide pms stars .
we derived the age of pms stars through the comparison with theoretical pms isochrones for metallicity @xmath39 , as appropriate for the lmc .
about @xmath201 of these objects are younger than 4myr and an equal number are older than 12myr , indicating that star formation in this field has proceeded over a long time , although our age resolution does not allow us to discriminate between an extended episode and a series of short and frequent bursts .
as regards the past 4myr , we obtain a lower limit to the star formation rate of @xmath129m@xmath6myr@xmath130 for objects in the range @xmath202m@xmath6 .
we find a pronounced difference in the spatial distribution of younger ( @xmath203myr ) and older ( @xmath204myr ) pms stars .
the former are concentrated near the centre of r136 , whereas older pms objects are mostly located to the east of it and are remarkably absent near the cluster centre .
this effect , which is not due to photometric incompleteness , is in general consistent with the photo - evaporation of the circumstellar discs of older stars caused by the uv radiation of the massive members of r136 .
however , not all older pms stars are affected in the same way , since some are in regions where the high density of the molecular gas can reduce the efficiency of the photo - evaporation process .
our results show that several generations of stars are present in this field , since there are objects with ages ranging from less than 1 to over 30myr , although no clear signs of causal effects or of triggered star formation are seen .
it appears , however , that the presence of the massive r136 cluster is sculpting the region and affecting low - mass neighbouring objects through the effects of photo - evaporation .
we thank an anonymous referee , whose suggestions have helped us to improve the presentation of this work .
this paper is based on early release science observations made by the wfc3 scientific oversight committee .
we are grateful to the director of the space telescope science institute for awarding director s discretionary time for this programme .
fp is grateful to the space science department of esa for their hospitality via the visitor programme .
np acknowledges partial support by hst
nasa grants go11547.06a and go11653.12a , and stsci  ddrf grant d0001.82435 .
zinnecker , h. , andersen , m. , brandl , b. , brandner , w. , hunter , d. , larson , r. , mccaughrean , m. , meylan , g. , moneti , a. 2002 , in `` extragalactic star clusters , '' iau symp .
207 , eds .
d. geisler , e. grebel , d. minniti ( san francisco : asp ) , 531