Formation of negative ions as a result of the interactions of a heavy ion beam with a surface under coverage of an alkali element is calculated by developing a Monte Carlo simulation code. Discussions are made on the relation among the main parameters of the physical process. The effects of the incident ion energy and cesium coverage on velocity distribution of the sputtered atoms, and consequently the negative ionization efficiency, are calculated. The simulation results show good agreement with the available experimental results.