Classification of Magnetohydrodynamic Simulations Using Wavelet Scattering Transforms
Andrew K. Saydjari, Stephen K. N. Portillo, Zachary Slepian, Sule Kahraman, Blakesley Burkhart, Douglas P. Finkbeiner
Abstract
Abstract The complex interplay of magnetohydrodynamics, gravity, and supersonic turbulence in the interstellar medium (ISM) introduces a non-Gaussian structure that can complicate a comparison between theory and observation. In this paper, we show that the wavelet scattering transform (WST), in combination with linear discriminant analysis (LDA), is sensitive to non-Gaussian structure in 2D ISM dust maps. WST-LDA classifies magnetohydrodynamic (MHD) turbulence simulations with up to a 97% true positive rate in our testbed of 8 simulations with varying sonic and Alfvénic Mach numbers. We present a side-by-side comparison with two other methods for non-Gaussian characterization, the reduced wavelet scattering transform (RWST) and the three-point correlation function (3PCF). We also demonstrate the 3D-WST-LDA, and apply it to the classification of density fields in position–position–velocity (PPV) space, where density correlations can be studied using velocity coherence as a proxy. WST-LDA is robust to common observational artifacts, such as striping and missing data, while also being sensitive enough to extract the net magnetic field direction for sub-Alfvénic turbulent density fields. We include a brief analysis of the effect of point-spread functions and image pixelization on 2D-WST-LDA applied to density fields, which informs the future goal of applying WST-LDA to 2D or 3D all-sky dust maps to extract hydrodynamic parameters of interest.