A direct method and ab initio force constants were used to calculate phonon dispersion curves and phonon density in Al. The force constants were determined from the Hellmann-Feynman forces induced by the displacement of an atom in the 2 x 2 x 2 fcc crystallographic supercell. This size of the supercell gives exact phonon frequencies at Γ, X, L, W points of the Brillouin zone. The calculated phonon dispersion curves are in good agreement with the experimental data.PACS numbers: 63.20.-e, 71.15.Nc .The vibrational properties determine a wide range of macroscopic behaviour of solids, e.g. specific heat and the sound velocity. In addition, very low-frequency modes can be associated with phase transformations, while imaginary frequencies provide an indication that the calculated structure is not the most stable. Finally, the phonon spectrum enables a good approximation to free energies to be made via the quasiharmonic approximation. In view of these, derivation of phonons from αb initio calculations has become a very important topic [1]. Generally the αb initio calculations of phonon frequencies fall into two methods: the linear response and the direct approach. Works still appear in the literature, in which various combinations of materials/methods/approximations and codes are tested. Aluminium is often considered to be a representative free-electron like metal and in recent years several calculations of phonons spectra for this material have been published. For example, the different approaches based on the phenomenological and the pseudopotential force constants have been compared in Ref. [2]. The αb initio cohesive energies were used to evaluate the phonon dispersion relations with and without three-body interaction [3]. Calculations based on the linear response theory which