Simulated analogues II: a new methodology for non-parametric matching of models to observations

        ABSTRACT
        Star formation is a multiscale problem, and only global simulations that account for the connection from the molecular cloud-scale gas flow to the accreting protostar can reflect the observed complexity of protostellar systems. Star-forming regions are characterized by supersonic turbulence, and as a result, it is not possible to simultaneously design models that account for the larger environment and in detail reproduce observed stellar systems. Instead, the stellar inventories can be matched statistically, and the best matches found that approximate specific observations. Observationally, a combination of single-dish telescopes and interferometers are now able to resolve the nearest protostellar objects on all scales from the protostellar core to the inner $10\, \mathrm{au}$. We present a new non-parametric methodology which uses high-resolution simulations and post-processing methods to match simulations and observations using deep learning. Our goal is to perform a down-selection from large data sets of synthetic images to a ranked list of best-matching candidates with respect to the observation. This is particularly useful for binary and multiple stellar systems that form in turbulent environments. The objective is to accelerate the rate at which we can do such comparisons, remove biases from hand-picking matches, and contribute to identifying the underlying physical processes that drive the creation and evolution of observed protostellar systems.