We have investigated the dynamics of the relaxed GaAs(llO) surface using an ab initio linearresponse approach in the framework of density-functional theory. The relaxation geometry was found by minimizing the total energy with the help of the Hellmann-Feynman forces. In terms of the electronic ground-state properties we have calculated the full phonon dispersion of GaAs(llO) along high symmetry lines of the surface Brillouin zone by means of a self-consistent first-order perturbation scheme without any adjustable parameters. Our results are in excellent agreement with all available experimental data.PACS numbers: 68.35.Ja, 71.10.+X, 81.60.CpThe last years have witnessed a substantial progress in studies of the vibrations in crystal surfaces. Inelastic He atom scattering and electron energy loss spectroscopy have reached such a stage of sophistication that surface phonon dispersion can be measured with high precision for a great variety of systems. The need for microscopic models which describe correctly not only the structural but also the dynamical properties of crystal surfaces has intensified theoretical investigations on this topic. A short summary of what has been done so far is given in [1], While model calculations can be used for a proper description of the force constants in the surfaces of ionic crystals, a self-consistent treatment of the electrons is necessary in metals and semiconductors. Up to now computationally demanding ab initio calculations of surface phonons have been done only for some metals [2-4] and elemental semiconductors [5]. In order to extend the self-consistent description to the surface dynamics of binary semiconductors, we have applied the densityfunctional linear-response approach proposed in [6,7] to GaAs(llO). The method used is different from all former surface calculations of this type. By treating the electronic response with the help of the one-particle Green's function, the advantages of frozen phonon calculations [2] and of the dielectric function approach [3][4][5] can be combined in one formalism. The method has been applied very successfully to the bulk dynamics of a large number of elemental and binary semiconductors [6] and to other materials [8]. Especially in the case of GaAs all of the details of the bulk phonon spectrum (not only in high symmetry directions) have been reproduced. Therefore we have a reliable basis for our surface phonon calculation.Among the surfaces of binary semiconductors, GaAs(110) has been studied the most, both experimentally [9-15] and theoretically [16][17][18][19][20][21][22][23][24]. It is now well established that the surface relaxation is mainly characterized by a bond-length-conserving rotation of the surface chains by a tilt angle of about 30° [9-11,16-19], with As shifted above the ideal (110) plane and Ga shifted towards the bulk. In contrast to the detailed examinations of the structural and electronical properties, up to now only a few investigations have been made on the surface phonons of GaAs(llO).On the experimental side inel...