 The ability of various interatomic potentials to reproduce the properties of silicene, that is, 2D single-layer silicon, polymorphs was examined. Structural and mechanical properties of flat, low-buckled, trigonal dumbbell, honeycomb dumbbell, and large honeycomb dumbbell silicene phases were obtained using density functional theory and molecular statics calculations with Tersoff, NAM, Stillinger Weber, EDIP, Rieck's FF, COM, and machine learning-based interatomic potentials. A quantitative systematic comparison and a discussion of the results obtained are reported. This article was authored by Marcin Mazziars.