Work function and surface energy per unit area were calculated in the framework of density functional theory (DFT) with Linearized A ug mented Plane Wave Plus Local Orbital method in full potential for a clean symmetric slab of silicon containing two (top and bottom) surfaces. The surfaces were theoretically modeled using supercell technique by stacking a variety of silicon layers along (111)...