This document aims at reviewing the different types of clustering algorithms and substructures detection techniques in order to study the spatial and kinematic clustering of stars and detect the gas components in molecular clouds. It is the deliverable: Report on Optimal Substructure Techniques for Stellar, Gas and Combined Samples, for the EU H2020 (COMPET-5-2015 - Space) project (A Gaia and Herschel Study of the Density Distribution and Evolution of Young Massive Star Clusters), Grant Agreement Number: 687528, with abbreviated code name StarFormMapper (SFM) project. The document is organized in the following sections: 1. General Introduction 2. Clustering of Discrete Distributions 3. Clustering of Continuous Distributions 4. Clustering in Astrophysics 5. StarFormMapper 6. Summary and Conclusions

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.

Multiplicity and clustering of young pre-main sequence stars appear as critical clues to constrain the star formation process. Taurus is the archetypical example of the most quiescent star forming regions that may still retain primeval signatures of star formation. This work identifies local overdense stellar structures at the 99.8\% confidence level above random expectation using the DBSCAN algorithm, and setting its free parameters based on the one-point correlation function and the k-nearest neighbor statistics. Nearly half of the entire stellar population in Taurus is found to be concentrated in 20 dense, tiny and prolate regions called NESTs (for Nested Elementary STructures). They are regularly spaced ($\approx 2$ pc) and mainly oriented along the gas filaments axes. Each NEST contains between 4 and 23 stars. Inside NESTs, the surface density of stars may be as high as 2500 pc$^{-2}$. Nearly half (11) of these NESTs contain about 75\% of the class 0/I objects. The balance between Class I, II, and, III fraction within the NESTs suggests that they may be ordered as an evolutionary temporal scheme. We have inferred that only 20\% of stars in Taurus do not belong to any kind of stellar groups. The mass of NESTs covers the same range as that of dense molecular cores. The distribution in size is bimodal peaking at 12.5 and 50 kAU and the distribution of the number of YSOs in NESTs as a function of size exhibits two regimes. We propose that the NESTs in their two size regimes represent the spatial imprints of stellar distribution at birth as they may have emerged within few millions years from their natal cloud either from a single core or from a chain of cores. We have identified them as the preferred sites of star formation in Taurus. These NESTs are the regions of highest stellar density and intermediate spatial scale structures between ultra-wide pairs and loose groups.

Using new proper motion data from recently published catalogs, we revisit the membership of previously identified members of the Upper Scorpius association. We confirmed 750 of them as cluster members based on the convergent point method, compute their kinematic parallaxes and combined them with Gaia parallaxes to investigate the 3D structure and geometry of the association using a robust covariance method. We find a mean distance of $146\pm 3\pm 6$~pc and show that the morphology of the association defined by the brightest (and most massive) stars yields a prolate ellipsoid with dimensions of $74\times38\times32$~pc$^{3}$, while the faintest cluster members define a more elongated structure with dimensions of $98\times24\times18$~pc$^{3}$. We suggest that the different properties of both populations is an imprint of the star formation history in this region.

We collected the multiplicity data of stars in Taurus to build an up-to-date stellar/multiplicity catalog. After a general study of nearest-neighbor statistics on spatial random distribution, we introduce the one-point correlation $\Psi$ function to complement the pair correlation function and define the spatial regimes departing from randomness in Taurus. We then perform a set of statistical studies to characterize the binary regime that prevails in Taurus. The $\Psi$ function in Taurus has a scale-free trend with a similar exponent as the correlation function at small scale. It extends almost 3 decades up to $\sim 60$ kAU showing a potential extended wide binary regime. This was hidden in the correlation function due to the clustering pattern blending. Distinguishing two stellar populations, single stars versus multiple systems (separation $ \leq 1$ kAU), within Class II/III stars observed at high angular resolution, we highlight a major spatial neighborhood difference between the two populations using nearest-neighbor statistics. The multiple systems are three times more likely to have a distant companion within 10 kAU when compared to single stars. We show that this is due to the presence of most probable physical ultra-wide pairs. These UWPs are biased towards high multiplicity and higher-stellar-mass components at shorter separations. The multiplicity fraction per ultra-wide pair with separation less than 10 kAU may be as high as 83.5 $\pm$ 19.6\%. We suggest that these young pre-main sequence UWPs may be pristine imprints of their spatial configuration at birth resulting from a cascade fragmentation scenario of the natal molecular core. They could be the older counterparts, at least for those separated by less than 10 kAU, to the $\le$ 0.5 Myr prestellar cores/Class 0 multiple objects observed at radio/mm wavelengths.

We provide a new view on the Cygnus-X north complex by accessing for the first time the low mass content of young stellar populations in the region. CFHT/WIRCam camera was used to perform a deep near-IR survey of this complex, sampling stellar masses down to ~0.1 M$_\odot$. Several analysis tools, including a extinction treatment developed in this work, were employed to identify and uniformly characterise a dozen unstudied young star clusters in the area. Investigation of their mass distributions in low-mass domain revealed a relatively uniform log-normal IMF with a characteristic mass of 0.32$\pm$0.08 M$_\odot$ and mass dispersion of 0.40$\pm$0.06. In the high mass regime, their derived slopes showed that while the youngest clusters (age < 4 Myr) presented slightly shallower values with respect to the Salpeter's, our older clusters (4 Myr < age < 18 Myr) showed IMF compliant values and a slightly denser stellar population. Although possibly evidencing a deviation from an 'universal' IMF, these results also supports a scenario where these gas dominated young clusters gradually 'build up' their IMF by accreting low-mass stars formed in their vicinity during their first ~3 Myr, before the gas expulsion phase, emerging at the age of ~4 Myr with a fully fledged IMF. Finally, the derived distances to these clusters confirmed the existence of at least 3 different star forming regions throughout Cygnus-X north complex, at distances of 500-900 pc, 1.4-1.7 kpc and 3.0 kpc, and revealed evidence of a possible interaction between some of these stellar populations and the Cygnus-OB2 association.