The influence of diffusion and defects of crystal surface on the crystal vibration mode are an important and basic subject in surface physics research. The frequency of lattice vibration corresponds to the energy band of the system. Since the vibrations of the atoms in the crystal lattice are not isolated from each other, and the crystal lattice is periodic, thereby forming a lattice wave in the crystal. The lattice wave represents that all the atoms in the crystal vibrate at an identical frequency, which is often called a vibration mode. The lattice chain model has been studied as the vibrating mode of phonon and the energy-band in solid state physics. The vibrating modes of the lattice chain model have been analyzed with the Newton equation and the Born-von-Karman boundary condition in the literaure. In general, it is difficult to solve this problem due to the complex nonlinear characteristic of the interactions between the matter particles and the environment. Noting the complicacy in directly diagonalizing quantum Hamiltonian operator of a long chain, we introduce the invariant eigenoperator method (IEO) for deriving the energy gap of a given crystal lattice without solving its eigenstates in the Heisenberg picture. The Heisenberg equation is as important as the Schrödinger equation. However, it has been seldom used for directly deriving the energy-gap in previous studies. Following the Heisenberg's original idea that most observable physical quantity in quantum mechanics is energy spectrum, Hong-yi Fan, one of the authors of the present paper, developed the IEO method. This method provides a natural result of combining both the Schrödinger operator and the Heisenberg equation. Using the IEO method, we study the vibration modes of crystal lattice, which are affected by absorbing an atom with mass m0, which is different from the mass of atom in the crystal. Moreover, the attractive potential constantβ0 of the lattice surface differs from the inner constantβ. With the help of invariant eigen-operator method, we deduce the vibration mode ω=√(2β(1-cosh α))/ħm, where α=ln-(mβ0+m0(-2β+β0)+√β0√-4 mm0β+(m+m0)2β0)/2m0β. Our numerical results show that vibration mode ω depends not only on the absorption potential and the mass of the absorbed atom, but also on the mass of the lattice atom and the inner potential. In general, by discussing the vibration modes via some numerical solutions or approximate methods, we show the relations between the system vibration modes with different parameters which describe the environment influences. These results can deepen our understanding of quantum Brownian motion and demonstrate the applicability of the IEO method.