The thermal behavior of an ammonia-covered Si͑100͒ surface is investigated by infrared spectroscopy and density functional methods. Upon adsorption at room temperature, ͑Si͒NH 2 and Si-H species are formed on the surface. Comparison of the vibrational studies with density functional calculations suggests that the ͑Si͒NH 2 structures are preferentially located on the same side along the silicon dimer row on a ͑2 ϫ 1͒ reconstructed Si͑100͒ surface, although a mixture of different long-range configurations is likely formed. Decomposition of these ͑Si͒NH 2 species is observed to start at temperatures as low as 500 K. Theoretical predictions of the vibrational modes indicate that at this point, the spectrum is composed of a combination of ͑Si͒ 2 NH and ͑Si͒ 3 N vibrational signatures, which result from insertion of N into Si-Si bonds. Our computational study of the formation of ͑Si͒ 2 NH structures indicates that subsurface insertion is more feasible if the strain imposed during the insertion in a Si dimer is attenuated by a ͑Si͒ 2 NH structure already inserted in the neighboring dimer along the same silicon dimer row. This cooperative reaction lowers the energetic requirements for subsurface insertion, providing a theoretical explanation for the mechanism of thermal decomposition of NH 3 on Si͑100͒ and for other systems where subsurface migration is observed experimentally.