Plane-wave density functional calculations are used to investigate the pressure dependence of the geometry and Gamma-point phonons of FeS2 pyrite up to 150 GPa. The linear response method is employed to calculate the vibrational properties. Raman-active modes are in excellent agreement with the experimental data available up to 50 GPa,(1) and we predict the evolution with pressure of the IR-active modes for which no high-pressure spectroscopic data have been reported so far. Over the wide pressure range investigated here, all vibrational frequencies depend nonlinearly on pressure; their pressure dependence is quantified by determining the full set of mode Gruneisen parameters and their pressure derivatives.