A semiclassical wave packet model that permits simultaneous study of surface diffraction, Debye–Waller effects, and inelastic energy transfer processes at the gas–surface interface is presented. The incident atomic beam is represented by a quantum mechanical wave packet whose momentum-space probability density corresponds to that employed in the molecular-beam experiments being simulated. A semiclassical forced-oscillator-type approximation couples the surface motion to that of the wave packet through a time-varying potential obtained from a solution of Hamilton’s equations of motion for the lattice. Explicit integration procedures are used to evolve the wave packet through this time-varying field. The average final kinetic energy of the beam, the angular distributions and diffraction patterns, Debye–Waller broadening effects, and the surface phonon spectrum are obtained from the scattered wave packet and its Fourier transform. The model has been applied to the case of an atomic hydrogen beam having a square distribution of incident momenta impinging upon a simple two-dimensional lattice. The results yield diffraction patterns that correlate with the known grating of the lattice. The average final kinetic energy of the beam is found to vary linearly with incident energy and with surface temperature in accord with the results of recent molecular beam experiments, and a surface phonon spectrum that exhibits one-, two-, and three-phonon processes is obtained.