
We systematically study the highresolution and polarized Raman spectra of multilayer (ML) MoTe2. The layerbreathing (LB) and shear (C) modes are observed in the ultralowfrequency region, which are used to quantitatively evaluate the interlayer coupling in ML MoTe2 based on the linear chain model, in which only the nearest interlayer coupling is considered. The Raman spectra on three different substrates verify the negligible substrate effect on the phonon frequencies of ML MoTe2. Ten excitation energies are used to measure the highfrequency modes of Nlayer MoTe2 (NL MoTe2; N is an integer). Under the resonant excitation condition, we observe Ndependent Davydov components in ML MoTe2 , originating from the Ramanactive A'1(A21g) modes at ~172 cm1. More than two Davydov components are observed in NL MoTe2 for N larger than 4 by Raman spectroscopy. The Ndependent Davydov components are further investigated based on the symmetry analysis. A van der Waals model only considering the nearest interlayer coupling has been proposed to well understand the Davydov splitting of highfrequency A'1(A21g) modes. The different resonant profiles for the two Davydov components in 3L MoTe2 indicate that proper excitation energy of ~1.82.2 eV must be chosen to observe the Davydov splitting in ML MoTe2 . Our work presents a simple way to identify layer number of ultrathin MoTe2 flakes by the corresponding number and peak position of Davydov components. Our work also provides a direct evidence from Raman spectroscopy of how the nearest van der Waals interactions significantly affect the frequency of the highfrequency intralayer phonon modes in multilayer MoTe2 and expands the understanding on the lattice vibrations and interlayer coupling of transition metal dichalcogenides and other twodimensional materials.