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 present a survey of variable stars detected in K2 Campaign 13 within the massive intermediate age ($\sim1$ Gyr) open cluster NGC 1817. We identify a complete sample of 44 red clump stars in the cluster, and have measured asteroseismic quantities ($\nu_{\rm max}$ and/or $\Delta \nu$) for 29 of them. Five stars showed suppressed dipole modes, and the occurrence rates indicate that mode suppression is unaffected by evolution through core helium burning. A subset of the giants in NGC 1817 (and in the similarly aged cluster NGC 6811) have $\nu_{\rm max}$ and $\Delta \nu$ values at or near the maximum observed for core helium burning stars, indicating they have core masses near the minimum for fully non-degenerate helium ignition. Further asteroseismic study of these stars can constrain the minimum helium core mass in red clump stars and the physics that determines this limit. Two giant stars show photometric variations on timescales similar to previously measured spectroscopic orbits. Thirteen systems in the field show eclipses, but only five are probable cluster members. We identify 32 $\delta$ Sct pulsators, 27 $\gamma$ Dor candidates, and 7 hybrids that are probable cluster members, with most new detections. We used the ensemble properties of the $\delta$ Sct stars to identify stars with possible radial pulsation modes. Among the oddities we have uncovered are: an eccentric orbit for a short-period binary containing a $\delta$ Sct pulsating star; a rare subgiant within the Hertzsprung gap showing $\delta$ Sct pulsations; and two hot $\gamma$ Dor pulsating star candidates.

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.

We identify a likely member of the NGC 752 open cluster (based on Gaia DR2 proper motions and parallax) having photometry consistent with being a white dwarf and a relatively low-mass main sequence star. If confirmed, this would be the first white dwarf known in this nearby cluster, and the brightest.

We use deep multi-epoch near-IR images of the VISTA Variables in the V\'ia L\'actea (VVV) Survey to search for RR Lyrae stars towards the Southern Galactic plane. Here we report the discovery of a group of RR Lyrae stars close together in VVV tile d025. Inspection of the VVV images and PSF photometry reveals that most of these stars are likely to belong to a globular cluster, that matches the position of the previously known star cluster FSR\,1716. The stellar density map of the field yields a $>100$ sigma detection for this candidate globular cluster, that is centered at equatorial coordinates $RA_{J2000}=$16:10:30.0, $DEC_{J2000}=-$53:44:56; and galactic coordinates $l=$329.77812, $b=-$1.59227. The color-magnitude diagram of this object reveals a well populated red giant branch, with a prominent red clump at $K_s=13.35 \pm 0.05$, and $J-K_s=1.30 \pm 0.05$. We present the cluster RR Lyrae positions, magnitudes, colors, periods and amplitudes. The presence of RR Lyrae indicates an old globular cluster, with age $>10$ Gyr. We classify this object as an Oosterhoff type I globular cluster, based on the mean period of its RR Lyrae type ab, $<P>=0.540$ days, and argue that this is a relatively metal-poor cluster with $[Fe/H] = -1.5 \pm 0.4$ dex. The mean extinction and reddening for this cluster are $A_{K_s}=0.38 \pm 0.02$, and $E(J-K_s)=0.72 \pm 0.02$ mag, respectively, as measured from the RR Lyrae colors and the near-IR color-magnitude diagram. We also measure the cluster distance using the RR Lyrae type ab stars. The cluster mean distance modulus is $(m-M)_0 = 14.38 \pm 0.03$ mag, implying a distance $D = 7.5 \pm 0.2$ kpc, and a Galactocentric distance $R_G=4.3$ kpc.

Discerning the nature of open cluster candidates is essential for both individual and statistical analyses of cluster properties. Here we establish the nature of thirteen cluster candidates from the FSR cluster list using photometry from the 2MASS and deeper, higher resolution UKIDSS-GPS and VISTA-VVV surveys. These clusters were selected because they were flagged in our previous studies as expected to contain a large proportion of pre-main sequence members or are at unusually small/large Galactocentric distances. We employ a decontamination procedure of JHK photometry to identify cluster members. Cluster properties are homogeneously determined and we conduct a cross comparative study of our results with the literature (where available). Seven of the here studied clusters were confirmed to contain PMS stars, one of which is a newly confirmed cluster. Our study of FSR1716 is the deepest to date and is in notable disagreement with previous studies, finding that it has a distance of about 7.3kpc and age of 10-12Gyr. As such, we argue that this cluster is a potential globular cluster candidate.

We present the extended source catalogue for the UKIRT Widefield Infrared Survey for H2 (UWISH2). The survey is unbiased along the inner Galactic Plane from l \approx 357deg to l \approx 65deg and |b| < 1.5deg and covers 209 square degrees. A further 42.0 and 35.5 square degrees of high dust column density regions have been targeted in Cygnus and Auriga. We have identified 33200 individual extended H2 features. They have been classified to be associated with about 700 groups of jets and outflows, 284 individual (candidate) Planetary Nebulae, 30 Supernova Remnants and about 1300 Photo-Dissociation Regions. We find a clear decline of star formation activity (traced by H2 emission from jets and photo-dissociation regions) with increasing distance from the Galactic Centre. About 60% of the detected candidate Planetary Nebulae have no known counterpart and 25% of all Supernova Remnants have detectable H2 emission associated with them.

Until now it has been impossible to observationally measure how star cluster scale height evolves beyond 1Gyr as only small samples have been available. Here we establish a novel method to determine the scale height of a cluster sample using modelled distributions and Kolmogorov-Smirnov tests. This allows us to determine the scale height with a 25% accuracy for samples of 38 clusters or more. We apply our method to investigate the temporal evolution of cluster scale height, using homogeneously selected sub-samples of Kharchenko et al. (MWSC), Dias et al. (DAML02), WEBDA, and Froebrich et al. (FSR). We identify a linear relationship between scale height and log(age/yr) of clusters, considerably different from field stars. The scale height increases from about 40pc at 1Myr to 75pc at 1Gyr, most likely due to internal evolution and external scattering events. After 1Gyr, there is a marked change of the behaviour, with the scale height linearly increasing with log(age/yr) to about 550pc at 3.5Gyr. The most likely interpretation is that the surviving clusters are only observable because they have been scattered away from the mid-plane in their past. A detailed understanding of this observational evidence can only be achieved with numerical simulations of the evolution of cluster samples in the Galactic Disk. Furthermore, we find a weak trend of an age-independent increase in scale height with galactocentric distance. There are no significant temporal or spatial variations of the cluster distribution zero point. We determine the Sun's vertical displacement from the Galactic Plane as $Z_\odot=18.5\pm1.2$pc.

Determining star cluster distances is essential to analyse their properties and distribution in the Galaxy. In particular it is desirable to have a reliable, purely photometric distance estimation method for large samples of newly discovered cluster candidates e.g. from 2MASS, UKIDSS-GPS and VISTA-VVV. Here, we establish an automatic method to estimate distances and reddening from NIR photometry alone, without the use of isochrone fitting. We employ a decontamination procedure of JHK photometry to determine the density of stars foreground to clusters and a galactic model to estimate distances. We then calibrate the method using clusters with known properties. This allows us to establish distance estimates with better than 40% accuracy. We apply our method to determine the extinction and distance values to 378 known open clusters and 397 cluster candidates from the list of Froebrich, Scholz and Raftery (2003). We find that the sample is biased towards clusters of a distance of approximately 3kpc, with typical distances between 2 and 6kpc. Using the cluster distances and extinction values, we investigate how the average extinction per kiloparsec distance changes as a function of Galactic longitude. We find a systematic dependence that can be approximated by A_H(l)[mag/kpc]=0.10+0.001*|l-180deg|/deg for regions more than 60deg from the Galactic Centre.