MNRAS 000,1–23 (2020) Preprint 12 February 2021 Compiled using MNRAS LATEX style file v3.0 High contrast and resolution near infrared photometry of the core of R136¢ Zeinab Khorrami,1y Maud Langlois,2 Paul C. Clark,1 Farrokh Vakili,3,9 Anne S. M. Buckner,4 Marta Gonzalez,5 Paul Crowther,6 Richard Wünsch,7 Jan Palouš,7 Stuart Lumsden,8 Estelle Moraux5 1School of Physics and Astronomy, Cardiff University, The Parade, Cardiff CF24 3AA, UK 2 Universite de Lyon, Universite Lyon 1, CNRS, CRAL UMR5574, Saint-Genis Laval, France 3 Universite Cote d’Azur, OCA, CNRS, Lagrange, France 4 School of Physics and Astronomy, University of Exeter, Stocker Road, Exeter, EX4 4QL, UK 5 Universite Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France 6 Department of Physics and Astronomy, Hounsfield Road, University of Sheffield, Sheffield, S3 7RH, UK 7 Astronomical Institute of the Czech Academy of Sciences, BočníII 1401/1a, 141 00 Praha 4, Czech Republic 8 School of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK 9 Department of Physics, Shahid Beheshti University, G.C., Tehran, Iran Accepted XXX. Received YYY; in original form ZZZ ABSTRACT We present the sharpest and deepest near infrared photometric analysis of the core of R136, a newly formed massive star cluster at the centre of the 30 Doradus star forming region in the Large Magellanic Cloud. We used the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO Very Large Telescope and operated in its IRDIS imaging mode, for the second time with longer exposure time in the H- and K filters. Our aim was to (i) increase the number of resolved sources in the core of R136, and (ii) to compare with the first epoch to classify the properties of the detected common sources between the two epochs. Within the field of view (FOV) of 10.8” × 12.1” (2.7pc × 3.0pc), we detected 1499 sources in both H and K filters, for which 76% of these sources have visual companions closer than 0.2". The larger number of detected sources, enabled us to better sample the mass function (MF). The MF slopes are estimated at ages of 1, 1.5 and 2 Myr, at different radii, and for different mass ranges. The MF slopes for the mass range of 10-300 M are about 0.3 dex steeper than the mass range of 3-300 M , for the whole FOV and different radii. Comparing the JHK colours of 790 sources common in between the two epochs, 67% of detected sources in the outer region (A¡ 3”) are not consistent with evolutionary models at 1 − 2Myr and with extinctions similar to the average cluster value, suggesting an origin from ongoing star formation within 30 Doradus, unrelated to R136. Key words: open clusters and associations: individual: R136 – stars: luminosity function, mass function – stars: massive – instrumentation: adaptive optics arXiv:2102.05972v1 [astro-ph.SR] 11 Feb 2021 1 INTRODUCTION standing of the true nature of R136 has significantly improved with increasing telescope resolution. The combination of photom- R136 is a very massive young star cluster that lies at the centre of etry (Hunter et al. 1995; Brandl et al. 1996; Andersen et al. 2009; the Tarantula nebula in the Large Magellanic Cloud (LMC). Host- De Marchi et al. 2011; Cignoni et al. 2015), ultraviolet spectrom- ing the most massive stars known in the Local Universe (Crowther etry (STIS/MAMA, Crowther et al.(2016)), visible (HST/FOS, et al. 2010; Crowther et al. 2016), R136 provides a unique op- Massey & Hunter(1998)) and near infrared (NIR) observations portunity to observationally study the formation of massive stars (VLT/SINFONI, Schnurr et al.(2009)) has resulted in more con- and clusters in the earliest stages of their evolution. Our under- straints on the R136 stellar population and its most luminous stars. ¢ Based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programme ID 0102.D0271 In 2015, R136 was observed for the first time in the NIR y E-mail: KhorramiZ@cardiff.ac.uk by the second generation Spectro-Polarimetric High-contrast Exo- © 2020 The Authors 2 Z. Khorrami et al. 1 planet Research (SPHERE, Beuzit et al. 2019) instrument of the Table 1. Exposure time log of VLT/SPHERE observations on R136. Very Large Telescope (VLT). Thanks to SPHERE’s high contrast and extreme adaptive optics (XAO), the sharpest images from the Obs. date Filter Single/Total _24=[nm] Δ_ [nm] central region of R136 were recorded in J and K. For the first time, Exposure[s] more than one thousand sources were detected in these bandpasses 2015-09-22 BB-J 4.0/1200 1245 235 within the small field of view (FOV) of IRDIS (10.9” × 12.3”) cov- 2015-09-22 BB-K 4.0/1200 2182 294 ering almost 2.7 × 3.1 pc of R136 core (Khorrami et al. 2017). 2018-10-10 BB-K 4.0/4352 2182 294 The SPHERE data on R136 core were used to partially resolve and 2018-11-14 BB-H 2.0/4352 1626 286 study the initial mass function (IMF) covering a mass range of 3 to 300 M , at ages of 1 and 1.5 Myr. The density in the core of R136 (A < 1.4 pc) was estimated and extrapolated in 3D for larger Table 2. Observing condition for R136 data in 2015 and 2018 radii (up to 6pc). Even at high angular resolution in the NIR the Data SR Seeing["] Airmass N?ℎ>C stars in the core are still unresolved due to crowding and central concentration of bright sources. Using evolutionary models of stars 2015-J 0.40 ± 0.05 0.63 ± 0.1 1.54-1.59 1110 more massive than 40 solar masses, Bestenlehner et al.(2020) finds 2015-K 0.75 ± 0.03 0.63 ± 0.1 1.61-1.67 1059 that the initial mass function (IMF) of massive stars in R136 is top- 2018-H 0.71 ± 0.05 0.58 ± 0.05 1.70-1.55 1658 heavy with a power law exponent W = 2.0 ± 0.3. They also estimate 2018-K 0.83 ± 0.03 0.69 ± 0.10 1.52-1.45 2528 the age of R136 between 1 and 2 Myr with a median age of around 1.6 Myr. and the filters information. Figure1 shows the final reduced images This paper presents the second epoch observation of R136 in the H (left) and K (right). These images are the deepest and in the NIR by VLT/SPHERE. These observations, obtained with sharpest images taken from the core of R136 so far. The Wolf-Rayet longer exposures in different NIR filters, aimed first to increase star R136a1 was used for guiding the AO loop of SPHERE confirm- the number of resolved sources in the core of R136, and second to ing the high level of performances even for faint guide stars which enable comparison with the first epoch to classify the property of the surpass both NACO and MAD performances. Figure2 compares detected common sources between these two epochs. We divided the highest resolution available images from the very central part of this paper as follows. In Section2 we explain in detail the data and R136 taken by HST (WFPC2 by Hunter et al. 1995) in left (top in the observing conditions. We then present the photometric analysis V, bottom in H), VLT/MAD (Campbell et al. 2010) in the middle in Section3 and the method we use to correct for completeness (top in K, bottom in H), and IRDIS in right (top in K, bottom in H). in Section4. In Section5 we describe the details of extinction These images shows the effect of better contrast, pixel-sampling, estimation. On the basis of the results we investigate the stellar mass sensitivity, and AO correction on resolving the stellar sources in the functions and the density in Section6. In Section7, we compare the most crowded part of R136. photometric analysis between the two epochs and between the two The range of airmass during these observations was 1.52 to sets of filters. We conclude by summarising our results in Sec.8. 1.45 in K, and 1.70 to 1.55 in the H. A log of observations is provided in Table2. For comparison we include the observing log for our first epoch data-sets taken in 2015. 2 OBSERVATIONS We used the SPHERE pipeline package to correct for dark (current), flat (fielding), (spatial) distortion, bad pixels and thermal We obtained the data via ESO Time Observation run (0102.D-0271) background due to instrument and sky (in K). Furthermore and to image R136 using the “classical imaging" mode of IRDIS Lan- in order to reach the highest sensitivity and the largest number of glois et al.(2014) in H- and K-bandpass filters. For our purpose detectable sources, additional corrections were carried out on the we used the same spectral band split into the two IRDIS chan- images. Based on a Gaussian fit using selected stars, we estimated nels to correct for residual detector artefacts such as hot pixels and and corrected the subpixel image drifts before combining the indi- uncorrelated detector noise among other instrumental effects. We vidual images. This allowed us to correct for residual tip-tilt errors maximized the image sharpness of the total exposure by discard- with an accuracy of a few mas. The final step is performed to com- ing the single frames with poorer Strehl ratio and by correcting bine the parallel images on the left and right side of the detector, "a posteriori" the residual tip-tilt image motion on each short ex- after the anamorphism correction.