Aims. To demonstrate that `INDICATE' is a powerful spatial analysis tool which when combined with kinematic data from Gaia DR2 can be used to robustly probe star formation history. Methods. We compared the dynamic & spatial distributions of young stellar objects (YSOs) at various evolutionary stages in NGC 2264 using Gaia DR2 proper motion data and INDICATE. Results. The dynamic & spatial behaviours of YSOs at different evolutionary stages are distinct. Dynamically, Class II YSOs predominately have non-random trajectories that are consistent with known substructures, whereas Class III YSOs have random trajectories with no clear expansion or contraction patterns. Spatially, there is a correlation between the evolutionary stage and source concentration: 69.4% of Class 0/I, 27.9% of Class II, and 7.7% of Class III objects are found to be clustered. The proportion of YSOs clustered with objects of the same class also follows this trend. Class 0/I objects are both found to be more tightly clustered with the general populous/objects of the same class than Class IIs and IIIs by a factor of 1.2/4.1 and 1.9/6.6, respectively. An exception to these findings is within 0.05deg of S Mon where Class III objects mimic the behaviours of Class II sources across the wider cluster region. Our results suggest (i) current YSOs distributions are a result of dynamical evolution, (ii) prolonged star formation has been occurring sequentially, and (iii) stellar feedback from S Mon is causing YSOs to appear as more evolved sources. Conclusions. Designed to provide a quantitative measure of clustering behaviours, INDICATE is a powerful tool with which to perform rigorous spatial analyses. Our findings are consistent with what is known about NGC 2264, effectively demonstrating that when combined with kinematic data from Gaia DR2 INDICATE can be used to robustly study the star formation history of a cluster.

We introduce our new code MYOSOTIS (Make Your Own Synthetic ObservaTIonS) which is designed to produce synthetic observations from simulated clusters. The code can synthesise observations from both ground- and spaced-based observatories, for a range of different filters, observational conditions and angular/spectral resolution. In this paper, we highlight some of the features of MYOSOTIS, creating synthetic observations from young massive star clusters. Our model clusters are simulated using nbody6 code and have different total masses, half-mass radii, and binary fractions. The synthetic observations are made at the age of 2 Myr with Solar metallicity and under different extinction conditions. For each cluster, we create synthetic images of the Hubble Space Telescope (HST) in the visible (WFPC2/F555W) as well as Very Large Telescopes (VLT) in the nearIR (SPHERE/IRDIS/Ks). We show how MYOSOTIS can be used to look at mass function (MF) determinations. For this aim we re-estimate stellar masses using a photometric analysis on the synthetic images. The synthetic MF slopes are compared to their actual values. Our photometric analysis demonstrate that depending on the adopted filter, extinction, angular resolution and pixel sampling of the instruments, the power-law index of the underlying MFs can be shallower than the observed ones by at least 0.25 dex which is in agreement with the observed discrepancies reported in the literature, specially for young star clusters.

Aims. To present the new novel statistical clustering tool 'INDICATE' which assesses and quantifies the degree of spatial clustering of each object in a dataset, discuss its applications as a tracer of morphological stellar features in star forming regions, and to look for these features in the Carina Nebula (NGC 3372). Results. We successfully recover known stellar structure of the Carina Nebula, including the 5 young star clusters in this region. Four sub-clusters contain no, or very few, stars with a degree of association above random which suggests they may be fluctuations in the field rather than real clusters. In addition we find: (1) Stars in the NW and SE regions have significantly different clustering tendencies, which is reflective of differences in the apparent star formation activity in these regions. Further study is required to ascertain the physical origin of the difference; (2) The different clustering properties between these two regions are even more pronounced for OB stars; (3) There are no signatures of classical mass segregation present in the SE region - massive stars here are not spatially concentrated together above random; (4) Stellar concentrations are more frequent around massive stars than typical for the general population, particularly in the Tr14 cluster; (5) There is a relation between the concentration of OB stars and the concentration of (lower mass) stars around OB stars in the centrally concentrated Tr14 and Tr15, but no such relation exists in Tr16. We conclude this is due to the highly sub-structured nature of Tr16. Conclusions. INDICATE is a powerful new tool employing a novel approach to quantify the clustering tendencies of individual objects in a dataset within a user-defined parameter space. As such it can be used in a wide array of data analysis applications.

This paper presents the sharpest near-IR images of the massive cluster R136 to date, based on the extreme adaptive optics of the SPHERE focal instrument implemented on the ESO/VLT and operated in its IRDIS imaging mode. Stacking-up a few hundreds of short exposures in J and Ks spectral bands over a FoV of 10.9"x12.3" centered on the R136a1 stellar component, enabled us to carry a refined photometric analysis of the core of R136. We detected 1110 and 1059 sources in J and Ks images respectively with 818 common sources. Thanks to better angular resolution and dynamic range, we found that more than 62.6% (16.5%) of the stars, detected both in J and Ks data, have neighbours closer than 0.2" (0.1"). Among resolved and detected sources R136a1 and R136c have visual companions and R136a3 is resolved as two stars separated by 59mas. The new set of detected sources were used to re-assess the age and extinction of R136 based on 54 spectroscopically stars that have been recently studied with HST slit-spectroscopy. Over 90% of these 54 sources identified visual companions (closer than 0.2"). We found the most probable age and extinction for these sources within the photometric and spectroscopic error-bars. Additionally, using PARSEC evolutionary isochrones and tracks, we estimated the stellar mass range for each detected source (common in J and K data) and plotted the generalized histogram of mass (MF with error-bars). Using SPHERE data, we have gone one step further and partially resolved and studied the IMF covering mass range of (3-300) Msun at the age of 1 and 1.5 Myr. The density in the core of R136 is estimated and extrapolated in 3D and larger radii (up to 6pc). We show that the stars in the core are still unresolved due to crowding, and the results we obtained are upper limits. Higher angular resolution is mandatory to overcome these difficulties.

MATISSE represents a great opportunity to image the environment around massive and evolved stars. This will allow one to put constraints on the circumstellar structure, on the mass ejection of dust and its reorganization , and on the dust-nature and formation processes. MATISSE measurements will often be pivotal for the understanding of large multiwavelength datasets on the same targets collected through many high-angular resolution facilities at ESO like sub-millimeter interferometry (ALMA), near-infrared adaptive optics (NACO, SPHERE), interferometry (PIONIER, GRAVITY), spectroscopy (CRIRES), and mid-infrared imaging (VISIR). Among main sequence and evolved stars, several cases of interest have been identified that we describe in this paper.

Continuous progress in data quality from HST, recent multiwavelength high resolution spectroscopy and high contrast imaging from ground adaptive optics on large telescopes need modeling of R136 to understand its nature and evolutionary stage. To produce the best synthesized multiwavelength images of R136 we need to simulate the effect of dynamical and stellar evolution, mass segregation and binary stars fraction on the survival of young massive clusters with the initial parameters of R136 in the LMC, being set to the present knowledge of this famous cluster. We produced a series of 32 young massive clusters using the NBODY6 code. Each cluster was tracked with adequate temporal samples to follow the evolution of R136 during its early stages. To compare the NBODY6 simulations with observational data, we created the synthetic images from the output of the code. We used the TLUSTY and KURUCZ model atmospheres to produce the fluxes in HST/ WFPC2 filters. GENEVA isochrones were used to track the evolution of stars. Then, we derived the observable parameters from synthetic scenes at the spatial resolution of HST/WFPC2 in the F814W filter (790.48nm). Surface brightness profile of the cluster, half-light radius, mass function and neighbor radius were used to select the best representation of R136. We compared the simulations of R136 to its HST imagery by creating synthetic scenes at the same resolution, pixel scale and FOV of the HST. We applied the same photometric analysis of the images as of the real ones. Having extracted the stellar sources, we estimated the mass-function, the surface brightness profile, the half-light radius and the neighbor radius across R136. The interpretation of these criteria point to the fact that an initially non-segregated cluster scenario is more representative of R136. This result pleads for the formation of massive stars by accretion instead of collision.

We present new near-infrared photometric measurements of the core of the young massive cluster NGC 3603 obtained with extreme adaptive optics. The data were obtained with the SPHERE instrument mounted on ESO Very Large Telescope, and cover three fields in the core of this cluster. We applied a correction for the effect of extinction to our data obtained in the J and K broadband filters and estimated the mass of detected sources inside the field of view of SPHERE/IRDIS, which is 13.5"x13.5". We derived the mass function (MF) slope for each spectral band and field. The MF slope in the core is unusual compared to previous results based on Hubble space telescope (HST) and very large telescope (VLT) observations. The average slope in the core is estimated as -1.06^{+0.26}_{-0.26} for the main sequence stars with 3.5 Msun < M < 120 Msun.Thanks to the SPHERE extreme adaptive optics, 814 low-mass stars were detected to estimate the MF slope for the pre-main sequence stars with 0.6 Msun< M < 3.5 Msun , Gamma = -0.54^{+0.11}_{-0.11} in the K-band images in two fields in the core of the cluster. For the first time, we derive the mass function of the very core of the NGC 3603 young cluster for masses in the range 0.6 - 120 Msun. Previous studies were either limited by crowding, lack of dynamic range, or a combination of both.

In this thesis we study several aspects of dynamical evolution of stellar clusters. The results of more than 200 simulations of single-mass star clusters with different initial total mass, half-mass radius and galactocentric distance, are reported. Recent studies of star clusters show a linear relation between a star cluster's dissolution time and its two-body relaxation time in logarithmic scale. We found that the single-mass star clusters do not show such a linear relation. We present new modified initial parameters to obtain a linear relation for single-mass star clusters. Also the evolution of multi-mass clusters and their lifetime, in the presence of the Galaxy is investigated. We simulate about 90 multi-mass star clusters with the Nbody6 code. These clusters have different initial total mass, half-mass radius and galactocentric distance. Finally we investigate the evolution of the stellar mass function and show that the slopes of the mass functions decrease with time. In addition we study the effect of galactocentric distance of star clusters on the evolution of the mass function.