“…Modeling such devices is challenging due to their high plasma density, reactive environment, and low pressure [1,3]. The modular code, hybrid plasma equipment model (HPEM) [4], has been successfully applied to the numerical investigations of IMPVD [1,5,6,7]. However, as the development of plasma devices progresses towards the lower pressure regime (e.g., ≤ 5 mTorr), the breakdown of the continuum assumption limits the application of the fluid model in HPEM and suggests the incorporation of the kinetic method, such as the direct simulation Monte Carlo (DSMC), in plasma simulations.…”