Adaptive phase field modeling is used in order to model the formation mechanism of a silicon faceted interface in three dimensions. We investigate the faceting condition for equilibrium shapes and dynamic situations. In this study, we propose a new anisotropic function of surface energy for the phase-field simulations in three-dimension, and negative stiffness is further considered. The morphological evolutions are presented and compare well with experimental findings. The growth mechanism is further discussed.