When solving problems of inhomogeneous gas dynamics in the slip regime, it is necessary to know the boundary conditions for the velocity, temperature, heat fluxes, etc., that is, the boundary conditions for the gas macroparameters. In particular, such problems arise in developing the theory of thermophoresis of moderately large aerosol particles [1]. The problem of monatomic and molecular (di- and polyatomic) gas slip along a boundary surface is considered in many publications (see, for example, [2-8]). The first-order effects include the isothermal and thermal gas slips characterized by the coefficients C m and K TS, respectively. In contrast to a monatomic gas, the molec ules of diatomic and polyatomic gases have internal degrees of freedom, which considerably complicates the kinetic equation [9]. The lack of reliable models for the intermolecular interaction potential predetermines the need to construct model kinetic equations [10]. In this study, for a diatomic g as whose molecules have rotational degrees of freedom, we propose a model kinetic equation obtained by developing the approach described in [6]. With the use of this model equation, the problem of diatomic gas slip along a plane surface is solved. As a result, for diatomic gases the coefficients C m and K TS, which depend on the thermophysical gas parameters and the intensity of inelastic collisions, are obtained. © Springer Science+Business Media, Inc. 2004.