Radiative heating and cooling of solid particles are studied numerically by considering the vibrational motions of solid particles and their relationships to the radiative absorption and emission with the molecular dynamics method. The radiative absorption and emission are assumed as a square-shaped function of the vibrational frequency and the Lennard-Jones potential is used for estimating the forces between particles. The change in the solid temperature is proportional to the difference between the products of the vibrational probability and the energy intensity of the absorption state and that of the emission state.