In rarefied gas dynamics scattering kernels deserve special attention since they contain all the essential information about the effects of physical and chemical properties of the gas-solid surface interface on the gas scattering process. However, to study the impact of the gas-surface interactions on the large-scale behavior of fluid flows, these scattering kernels need to be integrated in larger-scale models like Direct Simulation Monte Carlo (DSMC). In this work, the Gaussian mixture (GM) model, an unsupervised machine learning approach, is utilized to establish a scattering kernel for monoatomic (Ar) and diatomic (H\textsubscript{2}) gases directly from Molecular Dynamics (MD) simulations data. The GM scattering kernel is coupled to a pure DSMC solver to study isothermal and non-isothermal rarefied gas flows in a system with two parallel walls. To fully examine the coupling mechanism between the GM scattering kernel and the DSMC approach, a one-to-one correspondence between MD and DSMC particles is considered here. Benchmarked by MD results, the performance of the GM-DSMC is assessed against the Cercignani-Lampis-Lord (CLL) kernel incorporated into DSMC simulation (CLL-DSMC). The comparison of various physical and stochastic parameters shows the better performance of the GM-DSMC approach. Especially for the diatomic system, the GM-DSMC outperforms the CLL-DSMC approach. The fundamental superiority of the GM-DSMC approach confirms its potential as a multi-scale simulation approach for accurately measuring flow field properties in systems with highly nonequilibrium conditions.