Herein, we performed first principle calculation and classical molecular dynamics simulation to study structural optimization, band structure, and mechanical properties of differently stacked multilayer silicene. Several local energy minima have been identified as metastable conformation with different stacking mode and layer number. Bandstructure of low buckled AA bilayer silicene optimized with SCAN+rvv10 presents semiconducting behavior with a bandgap of 0.4419ev. Young's modulus of multilayer silicene shows low dependency on layer number or stacking mode. Whereas, fracture stress and strain is sensitive to the number of layers, specific stacking mode, and chirality. Furthermore, bending modulus of multilayer silicene (e.g., 0.44ev for monolayer silicene) is even lower than that of graphene, which may attribute to the flexibility of bond angle. Figure. 1 Top and side view of three 3 × 3 basic stacking mode: (a) AA; (b) AA'; (c) AB. Atoms with larger radius correspond to bottom layer. Different colours denote atoms in different buckling directions.