Based on the theory of surface lattice dynamics, the surface phonon spectrums along three symmetrical directions of $\bar \varGamma \bar L$ , $\bar L\bar M$ and $\bar \varGamma \bar M$ are simulated for the W(100) surface by using the modified analytic embedded atom method. The polarization vectors at different symmetrical points are also calculated. According to the criterion and marking method of surface mode, the surface modes along different symmetrical directions are drawn, the distribution range and mode coupling of surface modes are discussed as well. The vibration frequencies of surface modes calculated by us have been compared to available experimental datum and some theoretical values correspondingly. The results display that the present results are general agreement with the referenced experimental or theoretical results. Based on the calculated polarization vector, the surface vibration states are constructed for the atomic layers in the neighboring surface. And the polarization and local features of the surface modes along different symmetrical directions are analyzed. The results show that there are some coupling phenomena between surface mode dispersion, such as avoid crossing and independence crossing. The avoid crossing is found between the surface-mode branch S ₁and the surface-mode branch S ₂near ${\bar \zeta _y} = 0.32$ along $\bar L\bar M$ direction. In the region, going from $\bar L$ to $\bar M$ , S ₁changes from ypolarization to zpolarization, and S ₂changes from zpolarization to ypolarization. The independence crossings exist between surface-mode branch S ₁and surface-mode branch S ₂at ${\bar \zeta _x} = 0.5$ along $\bar \varGamma \bar L$ direction, and surface-mode branch S ₂and surface-mode branch S ₃at ${\bar \zeta _x} = 0.5$ along $\bar L\bar M$ direction, respectively. Before and after the crossings, the polarization and local features of the surface modes have not changed. Inspection of the polarization vectors, the coupling phenomena are iconically demonstrated.