Abstract: The effect of vacancy defects on the structural properties and the thermal stability of free standing silicene – a buckled structure of hexagonally arranged silicon atoms – is studied using reactive molecular dynamics simulations. Pristine silicene is found to be stable up to 1500 K, above which the system transits to a three-dimensional amorphous configuration. Vacancy defects result in local structural changes in the system and considerably reduce the thermal stability of silicene: depending on the size of the vacancy defect, the critical temperature decreases by more than 30%. However, the system is still found to be stable well above room temperature within our simulation time of 500 ps. We found that the, stability of silicene can be increased by saturating the dangling bonds at the defect edges by foreign atoms (e.g., hydrogen).