Quantifying kinematic substructure in star-forming regions with statistical tests of spatial autocorrelation

We investigate whether spatial-kinematic substructure in young star-forming regions can be quantified using Moran's I statistic. Its presence in young star clusters would provide an indication that the system formed from initially substructured conditions, as expected by the hierarchical model of star cluster formation, even if the cluster were spatially smooth and centrally concentrated. Its absence, on the other hand, would be evidence that star clusters form monolithically. The Moran's I statistic is applied to N-body simulations of star clusters with different primordial spatial-velocity structures, and its evolution over time is studied. It is found that this statistic can be used to reliably quantify spatial-kinematic substructure, and can be used to provide evidence as to whether the spatial-kinematic structure of regions with ages ≲6 Myr is best reproduced by the hierarchical or monolithic models of star formation. Moran's I statistic is also able to conclusively say whether the data are not consistent with initial conditions that lack kinematic substructure, such as the monolithic model, in regions with ages up to, and potentially beyond, 10 Myrs. This can therefore provide a kinematic signature of the star cluster formation process that is observable for many Myr after any initial spatial structure has been erased.